1
|
Buczek A, Rzepiela K, Kupka T, Broda MA, Kar T. Uracil-water interaction revisited - in search of single H-bonded secondary minima. Phys Chem Chem Phys 2024; 26:5169-5182. [PMID: 38261334 DOI: 10.1039/d3cp04057g] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2024]
Abstract
Monohydrated uracil (UW) complexes are stabilized by both O⋯HO and NH⋯O hydrogen bonds (H-bonds), simultaneously participating in forming three stable cyclic structures. The role and contribution of these individual H-bonds (O⋯HO and NH⋯O) to the stability of the three UW complexes are still not understood, because of the technical problems in obtaining their optimized structures by standard geometry optimization. The present study explores a non-standard approach to identify three single H-bonded local minima structures without imaginary frequency using DFT (M06-2X, B3LYP and B3LYP-D3), MP2 and CCSD(T) theories and Dunning's correlation-consistent aug-cc-pVTZ basis set, in both vacuum and aqueous media (CPCM method). The results reveal that these new structures are very shallow secondary minima between two deep wells or next to a deep well of primary minima (double H-bonded structures) in the potential energy surface. The H-bond energy of these single H-bonded complexes is found to be less sensitive to a wide range (about 15-20 degrees) of O⋯HO and NH⋯O angles, and the linearity is preferred in the stable three single H-bonded structures. The technical method used to locate such a shallow minimum is described in detail and may be useful for identifying local minima in other cases where consecutive multiple H-bonded structures are global minima. Energy decomposition (using symmetry adapted perturbation theory (SAPT)) of interaction energy, electron redistribution, and relevant vibrational modes are discussed.
Collapse
Affiliation(s)
- Aneta Buczek
- Faculty of Chemistry, University of Opole, Oleska 48, 45-052 Opole, Poland.
| | - Kacper Rzepiela
- Faculty of Chemistry, University of Opole, Oleska 48, 45-052 Opole, Poland.
| | - Teobald Kupka
- Faculty of Chemistry, University of Opole, Oleska 48, 45-052 Opole, Poland.
| | - Małgorzata A Broda
- Faculty of Chemistry, University of Opole, Oleska 48, 45-052 Opole, Poland.
| | - Tapas Kar
- Department of Chemistry and Biochemistry, Utah State University, Logan, UT 84322, USA.
| |
Collapse
|
2
|
Rzepiela K, Kaminský J, Buczek A, Broda MA, Kupka T. Electron Correlation or Basis Set Quality: How to Obtain Converged and Accurate NMR Shieldings for the Third-Row Elements? Molecules 2022; 27:molecules27238230. [PMID: 36500321 PMCID: PMC9737175 DOI: 10.3390/molecules27238230] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 11/02/2022] [Revised: 11/17/2022] [Accepted: 11/22/2022] [Indexed: 11/29/2022]
Abstract
The quality of theoretical NMR shieldings calculated at the quantum-chemical level depends on various theoretical aspects, of which the basis set type and size are among the most important factors. Nevertheless, not much information is available on the basis set effect on theoretical shieldings of the NMR-active nuclei of the third row. Here, we report on the importance of proper basis set selection to obtain accurate and reliable NMR shielding parameters for nuclei from the third row of the periodic table. All calculations were performed on a set of eleven compounds containing the elements Na, Mg, Al, Si, P, S, or Cl. NMR shielding tensors were calculated using the SCF-HF, DFT-B3LYP, and CCSD(T) methods, combined with the Dunning valence aug-cc-pVXZ, core-valence aug-cc-pCVXZ, Jensen polarized-convergent aug-pcSseg-n and Karlsruhe x2c-Def2 basis set families. We also estimated the complete basis set limit (CBS) values of the NMR parameters. Widely scattered nuclear shieldings were observed for the Dunning polarized-valence basis set, which provides irregular convergence. We show that the use of Dunning core-valence or Jensen basis sets effectively reduces the scatter of theoretical NMR results and leads to their exponential-like convergence to CBS. We also assessed the effect of vibrational, temperature, and relativistic corrections on the predicted shieldings. For systems with single bonds, all corrections are relatively small, amounting to less than 4% of the CCSD(T)/CBS value. Vibrational and temperature corrections were less reliable for H3PO and HSiCH due to the high anharmonicity of the molecules. An abnormally high relativistic correction was observed for phosphorus in PN, reaching ~20% of the CCSD(T)/CBS value, while the correction was less than 7% for other tested molecules.
Collapse
Affiliation(s)
- Kacper Rzepiela
- Faculty of Chemistry, University of Opole, 48 Oleska Street, 46-052 Opole, Poland
| | - Jakub Kaminský
- Institute of Organic Chemistry and Biochemistry of the CAS, Flemingovo nám. 2, 166 10 Prague, Czech Republic
- Correspondence: (J.K.); (T.K.)
| | - Aneta Buczek
- Faculty of Chemistry, University of Opole, 48 Oleska Street, 46-052 Opole, Poland
| | - Małgorzata A. Broda
- Faculty of Chemistry, University of Opole, 48 Oleska Street, 46-052 Opole, Poland
| | - Teobald Kupka
- Faculty of Chemistry, University of Opole, 48 Oleska Street, 46-052 Opole, Poland
- Correspondence: (J.K.); (T.K.)
| |
Collapse
|
3
|
Najgebauer P, Staś M, Wrzalik R, Broda MA, Wieczorek PP, Andrushchenko V, Kupka T. Muscimol hydration and vibrational spectroscopy – The impact of explicit and implicit water. J Mol Liq 2022. [DOI: 10.1016/j.molliq.2022.119870] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
|
4
|
Kupka T, Gajda T, Ochędzan-Siodłak W, Buczek A, Broda MA. On the impact of side methyl groups on the structure and vibrational properties of β-carotenoids. The case of butadiene and isoprene. Food Chem 2022; 369:130880. [PMID: 34438344 DOI: 10.1016/j.foodchem.2021.130880] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/21/2021] [Revised: 08/02/2021] [Accepted: 08/14/2021] [Indexed: 11/04/2022]
Abstract
Theoretical consideration about the impact of methyl groups on the structure and vibrational properties of β-carotenoids, using medium size molecules of trans-butadiene and trans-isoprene, are reported. Density functional theory (DFT) calculations with correlation-consistent and polarization-consistent basis sets were applied to trans-1,3-butadiene and trans-isoprene as the smallest building bricks of β-carotenoids. Their structure and harmonic vibrations were estimated in the complete basis set limit (CBS) using the non-linear least square fit. Optimized geometries and harmonic frequencies, obtained with B3LYP and BLYP density functionals and large basis sets, were favorably reproduced by a significantly faster approach, using a recently modified STO(1M)-3G Slater-type basis set. Selected density functionals with STO(1M)-3G and 6-311++G** basis sets were also successful in predicting β-carotene structures and harmonic vibrations. This work demonstrates the potential applicability of the proposed level of theory for larger molecules, including β-carotenoids, present in numerous natural food sources. The proposed scheme of molecular modeling, applied to biologically active compounds in food, could provide a deeper insight into their function in vivo, which is directly related to their structure and spectroscopic properties. It could also support the experimental qualitative analysis, based on peak assignment of β-carotenoids in various food sources.
Collapse
Affiliation(s)
- Teobald Kupka
- University of Opole, Faculty of Chemistry, 48, Oleska Street, 45-052 Opole, Poland.
| | - Tomasz Gajda
- University of Opole, Faculty of Chemistry, 48, Oleska Street, 45-052 Opole, Poland
| | | | - Aneta Buczek
- University of Opole, Faculty of Chemistry, 48, Oleska Street, 45-052 Opole, Poland
| | - Małgorzata A Broda
- University of Opole, Faculty of Chemistry, 48, Oleska Street, 45-052 Opole, Poland
| |
Collapse
|
5
|
Buczek A, Kupka T, Broda MA, Maślanka S, Pentak D. Liposomes as nonspecific nanocarriers for 5-Fluorouracil in the presence of cyclodextrins. J Mol Liq 2021. [DOI: 10.1016/j.molliq.2021.117623] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
|
6
|
Staś M, Broda MA, Siodłak D. Thiazole-amino acids: influence of thiazole ring on conformational properties of amino acid residues. Amino Acids 2021; 53:673-686. [PMID: 33837859 PMCID: PMC8128816 DOI: 10.1007/s00726-021-02974-0] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/05/2020] [Accepted: 03/29/2021] [Indexed: 12/29/2022]
Abstract
Post-translational modified thiazole-amino acid (Xaa-Tzl) residues have been found in macrocyclic peptides (e.g., thiopeptides and cyanobactins), which mostly inhibit protein synthesis in Gram + bacteria. Conformational study of the series of model compounds containing this structural motif with alanine, dehydroalanine, dehydrobutyrine and dehydrophenylalanine were performed using DFT method in various environments. The solid-state crystal structure conformations of thiazole-amino acid residues retrieved from the Cambridge Structural Database were also analysed. The studied structural units tend to adopt the unique semi-extended β2 conformation; which is stabilised mainly by N-H⋯NTzl hydrogen bond, and for dehydroamino acids also by π-electron conjugation. The conformational preferences of amino acids with a thiazole ring were compared with oxazole analogues and the role of the sulfur atom in stabilising the conformations of studied peptides was discussed.
Collapse
Affiliation(s)
- Monika Staś
- Faculty of Chemistry, University of Opole, 45-052, Opole, Poland.
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Science, Flemingovo Náměstí 2, 166 10, Praha 6, Czech Republic.
| | | | - Dawid Siodłak
- Faculty of Chemistry, University of Opole, 45-052, Opole, Poland.
| |
Collapse
|
7
|
Nackiewicz J, Kołodziej Ł, Poliwoda A, Broda MA. Oxidation of diclofenac in the presence of iron(II) octacarboxyphthalocyanine. Chemosphere 2021; 265:129145. [PMID: 33302208 DOI: 10.1016/j.chemosphere.2020.129145] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/04/2020] [Revised: 10/29/2020] [Accepted: 11/26/2020] [Indexed: 06/12/2023]
Abstract
This paper presents the results of the research on the influence of catalytic activity of iron(II) octacarboxyphthalocyanines (FePcOC) on the transformation of diclofenac (DCF) which is the most popular anti-inflammatory analgesic. Diclofenac poses a serious threat to the natural environment. The paper demonstrates that diclofenac, in the presence a monomeric form of iron octacarboxyphthalocyanine and hydroxyl radicals (HO•) (from H2O2), undergoes a transformation into diclofenac-2,5-iminoquinone (DCF-2,5-IQ), causing distinct changes in the UV-Vis absorption spectrum. In the presence of iron octacarboxyphthalocyanine and H2O2, the previously colourless diclofenac solution becomes intense orange. As a result, a new band at approx. 450 nm appears in the absorption spectrum. HPLC analysis has shown that the concentration of diclofenac decreases with time. TD-DFT calculations using the CAM-B3LYP/6-31+G (d, p) method have been conducted to confirm experimental data concerning the formation of a new band at λmax = 450 nm.
Collapse
Affiliation(s)
- Joanna Nackiewicz
- Faculty of Chemistry, University of Opole, Oleska 48, Opole, 45-052, Poland.
| | - Łukasz Kołodziej
- Faculty of Chemistry, University of Opole, Oleska 48, Opole, 45-052, Poland
| | - Anna Poliwoda
- Faculty of Chemistry, University of Opole, Oleska 48, Opole, 45-052, Poland
| | - Małgorzata A Broda
- Faculty of Chemistry, University of Opole, Oleska 48, Opole, 45-052, Poland
| |
Collapse
|
8
|
Rzepiela K, Buczek A, Kupka T, Broda MA. Factors Governing the Chemical Stability and NMR Parameters of Uracil Tautomers and Its 5-Halogen Derivatives. Molecules 2020; 25:E3931. [PMID: 32872098 PMCID: PMC7504704 DOI: 10.3390/molecules25173931] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/23/2020] [Revised: 08/21/2020] [Accepted: 08/26/2020] [Indexed: 01/10/2023] Open
Abstract
We report on the density functional theory (DFT) modelling of structural, energetic and NMR parameters of uracil and its derivatives (5-halogenouracil (5XU), X = F, Cl, Br and I) in vacuum and in water using the polarizable continuum model (PCM) and the solvent model density (SMD) approach. On the basis of the obtained results, we conclude that the intramolecular electrostatic interactions are the main factors governing the stability of the six tautomeric forms of uracil and 5XU. Two indices of aromaticity, the harmonic oscillator model of aromaticity (HOMA), satisfying the geometric criterion, and the nuclear independent chemical shift (NICS), were applied to evaluate the aromaticity of uracil and its derivatives in the gas phase and water. The values of these parameters showed that the most stable tautomer is the least aromatic. A good performance of newly designed xOPBE density functional in combination with both large aug-cc-pVQZ and small STO(1M)-3G basis sets for predicting chemical shifts of uracil and 5-fluorouracil in vacuum and water was observed. As a practical alternative for calculating the chemical shifts of challenging heterocyclic compounds, we also propose B3LYP calculations with small STO(1M)-3G basis set. The indirect spin-spin coupling constants predicted by B3LYP/aug-cc-pVQZ(mixed) method reproduce the experimental data for uracil and 5-fluorouracil well.
Collapse
Affiliation(s)
| | | | - Teobald Kupka
- Department of Physical Chemistry and Molecular Modeling, Faculty of Chemistry, University of Opole, 48 Oleska Street, 45−052 Opole, Poland; (K.R.); (A.B.)
| | - Małgorzata A. Broda
- Department of Physical Chemistry and Molecular Modeling, Faculty of Chemistry, University of Opole, 48 Oleska Street, 45−052 Opole, Poland; (K.R.); (A.B.)
| |
Collapse
|
9
|
Kupka T, Broda MA, Wieczorek PP. What is the form of muscimol from fly agaric mushroom (Amanita muscaria) in water? An insight from NMR experiment supported by molecular modeling. Magn Reson Chem 2020; 58:584-593. [PMID: 31912552 DOI: 10.1002/mrc.4990] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/22/2019] [Revised: 12/30/2019] [Accepted: 01/02/2020] [Indexed: 06/10/2023]
Abstract
The biologically active alkaloid muscimol is present in fly agaric mushroom (Amanita muscaria), and its structure and action is related to human neurotransmitter γ -aminobutyric acid (GABA). The current study reports on determination of muscimol form present in water solution using multinuclear 1 H and 13 C nuclear magnetic resonance (NMR) experiments supported by density functional theory molecular modeling. The structures of three forms of free muscimol molecule both in the gas phase and in the presence of water solvent, modeled by polarized continuous model, and nuclear magnetic isotropic shieldings, the corresponding chemical shifts, and indirect spin-spin coupling constants were calculated. Several J-couplings observed in proton and carbon NMR spectra, not available before, are reported. The obtained experimental spectra, supported by theoretical calculations, favor the zwitterion form of muscimol in water. This structure differs from NH isomer, previously determined in dimethyl sulfoxide (DMSO) solution. In addition, positions of signals C3 and C5 are reversed in both solvents.
Collapse
Affiliation(s)
- Teobald Kupka
- Faculty of Chemistry, University of Opole, Opole, Poland
| | | | | |
Collapse
|
10
|
Dyrda G, Broda MA, Hnatejko Z, Pędziński T, Słota R. Adducts of free-base meso-tetraarylporphyrins with trihaloacetic acids: Structure and photostability. J Photochem Photobiol A Chem 2020. [DOI: 10.1016/j.jphotochem.2020.112445] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
11
|
Kupka T, Mnich A, Broda MA. Performance of revised STO(1M)-3G basis set for prediction of 5-fluorocytosine chemical shifts. Magn Reson Chem 2019; 57:489-498. [PMID: 31013546 DOI: 10.1002/mrc.4879] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/14/2019] [Revised: 04/14/2019] [Accepted: 04/15/2019] [Indexed: 06/09/2023]
Abstract
Nuclear shieldings and chemical shifts of 5-fluorocytosine (5FC) were predicted in the gas phase and DMSO solution modeled by polarizable continuum model using B3LYP density functional and revised STO(1M)-3G basis set. For comparison, eight arbitrary selected basis sets including STO-3G and medium-size Pople-type and larger dedicated Jensen-type ones were applied. The former basis sets were significantly smaller, but the calculated structural parameters, harmonic vibrational frequencies, were very accurate and close to those obtained with larger, polarization-consistent ones. The predicted 13 C and 1 H chemical shieldings of 5FC and cytosine, selected as parent molecule, were acceptable (root mean square for 13 C chemical shifts in DMSO of about 5 ppm and less) though less accurate than those calculated with large basis sets, dedicated for prediction of nuclear magnetic resonance parameters.
Collapse
Affiliation(s)
- Teobald Kupka
- Faculty of Chemistry, University of Opole, 48, Oleska Street, 46-052, Opole, Poland
| | - Adrianna Mnich
- Faculty of Chemistry, University of Opole, 48, Oleska Street, 46-052, Opole, Poland
| | - Małgorzata A Broda
- Faculty of Chemistry, University of Opole, 48, Oleska Street, 46-052, Opole, Poland
| |
Collapse
|
12
|
Kupka T, Gajda Ł, Stobiński L, Kołodziej Ł, Mnich A, Buczek A, Broda MA. Local aromaticity mapping in the vicinity of planar and nonplanar molecules. Magn Reson Chem 2019; 57:359-372. [PMID: 31034627 DOI: 10.1002/mrc.4880] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/01/2019] [Accepted: 04/23/2019] [Indexed: 06/09/2023]
Abstract
We report on nucleus-independent magnetic shielding (NICS) scans over the centers of six- and five-membered rings in selected metal phthalocyanines (MPc) and fullerene C60 for more accurate characterization of local aromaticity in these compounds. Detailed tests were conducted on model aromatic molecules including benzene, pyrrole, indole, isoindole, and carbazole and subsequently applied to H2 Pc, ZnPc, Al(OH)Pc, and CuPc. Similar behavior of three selected magnetic probes, Bq, 3 He, and 7 Li+ , approaching perpendicularly the ring centers, was observed. For better visualization of shielding zone over the centers of aromatic rings, we introduced a simple mathematical procedure: the first and second derivatives of scan curves with respect to magnetic probe position enabled their additional examination. It allowed an easier localization of curve minimum and discrimination between areas in space varying by the magnetic field magnitude and to illustrate local aromaticity of two different kinds of rings in MPc with better resolution. Our results supported earlier reports on very low aromaticity indexes of pyrrole ring incorporated into MPc and significant aromaticity of the central macrocycle. No direct dependence between harmonic oscillator model of aromaticity and NICS was observed. Instead, a correlation between position of scan curve minimum and its magnitude were observed. In addition, the NICS values and 3 He chemical shifts in the middle of neutral C60 and C606- anion agreed well with the reported experimental NMR values for He@C60 and He@C606- .
Collapse
Affiliation(s)
- Teobald Kupka
- Faculty of Chemistry, University of Opole, Opole, Poland
| | - Łukasz Gajda
- Faculty of Chemistry, University of Opole, Opole, Poland
| | - Leszek Stobiński
- Faculty of Process and Chemical Engineering, Warsaw University of Technology, 1, Waryńskiego Street, 00-645, Warsaw, Poland
- nanoMAT Ltd., Waryńskiego 1 (FL. 4, RM. 401), 00-645, Warsaw, Poland
| | | | - Adrianna Mnich
- Faculty of Chemistry, University of Opole, Opole, Poland
| | - Aneta Buczek
- Faculty of Chemistry, University of Opole, Opole, Poland
| | | |
Collapse
|
13
|
Kupka T, Buczek A, Broda MA, Gajda Ł, Ignasiak M. Convergence of nuclear magnetic shieldings and one-bond 1 J( 11 B 1 H) indirect spin-spin coupling constants in small boron molecules. Magn Reson Chem 2018; 56:338-351. [PMID: 29361201 DOI: 10.1002/mrc.4714] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/24/2017] [Revised: 01/10/2018] [Accepted: 01/13/2018] [Indexed: 06/07/2023]
Abstract
Self-consistent field Hartree-Fock, density functional theory, and coupled-cluster calculations of the nuclear magnetic shielding constants of BH and BH3 molecules have been conducted to characterize the convergence of individual results obtained with correlation-consistent and polarization-consistent basis sets. The individual 11 B and 1 H NMR parameters were estimated in the complete basis set limit and compared with benchmark results. Only the KT3 density functional accurately reproduced 11 B shielding in BH molecule.
Collapse
Affiliation(s)
- Teobald Kupka
- Faculty of Chemistry, University of Opole, 48, Oleska Street, 46-052, Opole, Poland
| | - Aneta Buczek
- Faculty of Chemistry, University of Opole, 48, Oleska Street, 46-052, Opole, Poland
| | - Małgorzata A Broda
- Faculty of Chemistry, University of Opole, 48, Oleska Street, 46-052, Opole, Poland
| | - Łukasz Gajda
- Faculty of Chemistry, University of Opole, 48, Oleska Street, 46-052, Opole, Poland
| | - Monika Ignasiak
- Faculty of Chemistry, University of Opole, 48, Oleska Street, 46-052, Opole, Poland
| |
Collapse
|
14
|
Gajda Ł, Kupka T, Broda MA, Leszczyńska M, Ejsmont K. Method and basis set dependence of the NICS indexes of aromaticity for benzene. Magn Reson Chem 2018; 56:265-275. [PMID: 29211311 DOI: 10.1002/mrc.4690] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/24/2017] [Revised: 11/21/2017] [Accepted: 11/24/2017] [Indexed: 05/24/2023]
Abstract
The role of theory level in prediction of benzene magnetic indexes of aromaticity is analysed and compared with calculated nuclear magnetic shieldings of 3 He used as NMR probe. Three closely related nucleus-independent chemical shift (NICS) based indexes were calculated for benzene at SCF-HF, MP2, and DFT levels of theory and the impact of basis set on these quantities was studied. The changes of benzene NICS(0), NICS(1), and NICS(1)zz parameters calculated using SCF-HF, MP2 and several density functionals were within 1 to 3 ppm. Similar deviations between magnetic indexes of aromaticity were observed for values calculated with selected basis sets. Only very small effect of polar solvent on benzene aromaticity was predicted. The 3 He nuclear magnetic isotropic shielding (σ) and its zz-components (σzz ) of helium atom approaching the centre of benzene ring from above produced similar curves versus benzene-He distance to NICS parameters calculated for similarly moving Bq ghost atom. We also propose an experimental verification of NICS calculations by designing the 3 He NMR measurement for benzene saturated with helium gas or in low temperature matrices.
Collapse
Affiliation(s)
- Łukasz Gajda
- Faculty of Chemistry, University of Opole, 48, Oleska Street, 45-052, Opole, Poland
| | - Teobald Kupka
- Faculty of Chemistry, University of Opole, 48, Oleska Street, 45-052, Opole, Poland
| | - Małgorzata A Broda
- Faculty of Chemistry, University of Opole, 48, Oleska Street, 45-052, Opole, Poland
| | | | - Krzysztof Ejsmont
- Faculty of Chemistry, University of Opole, 48, Oleska Street, 45-052, Opole, Poland
| |
Collapse
|
15
|
Abstract
The conformational propensities of the Aib residue on the example of two model peptides Ac-Aib-NHMe (1) and Ac-Aib-NMe2 (2), were studied by B3LYP and M06-2X functionals, in the gas phase and in the polar solvents. To verify the reliability of selected functionals, we also performed MP2 calculations for the tested molecules in vacuum. Polarizable continuum models (PCM and SMD) were used to estimate the solvent effect. Ramachandran maps were calculated to find all energy minima. Noncovalent intramolecular interactions due to hydrogen-bonds and dipole attractions between carbonyl groups are responsible for the relative stabilities of the conformers. In order to verify the theoretical results, the available conformations of similar X-ray structures from the Cambridge Crystallographic Data Center (CCDC) were analyzed. The results of the calculations show that both derivatives with the Aib residue in the gas phase prefer structures stabilized by intramolecular N–H⋯O hydrogen bonds, i.e., C5 and C7 conformations, while polar solvent promotes helical conformation with φ, ψ values equal to +/−60°, +/−40°. In addition, in the case of molecule 2, the helical conformation is the only one available in the polar environment. This result is fully consistent with the X-ray data. Effect of solvent on the Ramachandran maps of the model peptides with Aib residue ![]()
Collapse
Affiliation(s)
- Roksana Wałęsa
- Faculty of Chemistry, University of Opole, 48, Oleska St., 45-052, Opole, Poland
| | - Małgorzata A Broda
- Faculty of Chemistry, University of Opole, 48, Oleska St., 45-052, Opole, Poland.
| |
Collapse
|
16
|
Kusakiewicz-Dawid A, Porada M, Ochędzan-Siodłak W, Broda MA, Bujak M, Siodłak D. Pyrazole amino acids: hydrogen bonding directed conformations of 3-amino-1H-pyrazole-5-carboxylic acid residue. J Pept Sci 2017; 23:716-726. [PMID: 28608410 DOI: 10.1002/psc.3018] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/12/2017] [Revised: 05/24/2017] [Accepted: 05/25/2017] [Indexed: 11/06/2022]
Abstract
A series of model compounds containing 3-amino-1H-pyrazole-5-carboxylic acid residue with N-terminal amide/urethane and C-terminal amide/hydrazide/ester groups were investigated by using NMR, Fourier transform infrared, and single-crystal X-ray diffraction methods, additionally supported by theoretical calculations. The studies demonstrate that the most preferred is the extended conformation with torsion angles ϕ and ψ close to ±180°. The studied 1H-pyrazole with N-terminal amide/urethane and C-terminal amide/hydrazide groups solely adopts this energetically favored conformation confirming rigidity of that structural motif. However, when the C-terminal ester group is present, the second conformation with torsion angles ϕ and ψ close to ±180° and 0°, respectively, is accessible. The conformational equilibrium is observed in NMR and Fourier transform infrared studies in solution in polar environment as well as in the crystal structures of other related compounds. The observed conformational preferences are clearly related to the presence of intramolecular interactions formed within the studied residue. Copyright © 2017 European Peptide Society and John Wiley & Sons, Ltd.
Collapse
Affiliation(s)
| | - Monika Porada
- Faculty of Chemistry, University of Opole, Oleska 48, 45-052, Opole, Poland
| | | | - Małgorzata A Broda
- Faculty of Chemistry, University of Opole, Oleska 48, 45-052, Opole, Poland
| | - Maciej Bujak
- Faculty of Chemistry, University of Opole, Oleska 48, 45-052, Opole, Poland
| | - Dawid Siodłak
- Faculty of Chemistry, University of Opole, Oleska 48, 45-052, Opole, Poland
| |
Collapse
|
17
|
Staś M, Bujak M, Broda MA, Siodłak D. Conformational preferences and synthesis of isomers Z and E of oxazole-dehydrophenylalanine. Biopolymers 2017; 106:283-94. [PMID: 27061820 DOI: 10.1002/bip.22852] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/10/2016] [Revised: 03/30/2016] [Accepted: 04/05/2016] [Indexed: 12/17/2022]
Abstract
Dehydrophenylalanine, ΔPhe, is the most commonly studied α,β-dehydroamino acid. In nature, further modifications of the α,β-dehydroamino acids were found, for example, replacement of the C-terminal amide group by oxazole ring. The conformational properties of oxazole-dehydrophenylalanine residue (ΔPhe-Ozl), both isomers Z and E, were investigated. To determine all possible conformations, theoretical calculations were performed using Ac-(Z/E)-ΔPhe-Ozl(4-Me) model compounds at M06-2X/6-31++G(d,p) level of theory. Ac-(Z/E)-ΔPhe-Ozl-4-COOEt compounds were synthesized and the conformational preferences of each isomer, Z and E, were investigated using FTIR and NMR-NOE in solutions of increasing polarity (CHCl3 , DMSO-d6). The solid-state low-temperature structures of Ac-(Z)-ΔPhe-Ozl-4-COOEt and its intermediate analog Ac-(Z)-ΔPhe-Ozn(4-OH)-4-COOEt were also determined. In a weakly polar environment, the ΔPhe-Ozl residue has a tendency to adopt the conformation β2 with the calculated φ and ψ angles of -127° and 0° for the isomer Z and -170° and 26° for the isomer E. The increase of environment polarity favors the helical conformation α and the beta-turn like conformation β, but the conformation β2 seems to be still accessible. The (E)-ΔPhe-Ozl residue can be obtained from the isomer Z in photoisomerization reaction. However, hydroxyl-oxazoline-dehydrophenylalanine ΔPhe-Ozn(4-OH) decomposes in such conditions. Alternatively, (E)-ΔPhe-NH2 can be applied as a substrate in the Hantzsch reaction. © 2016 Wiley Periodicals, Inc. Biopolymers (Pept Sci) 106: 283-294, 2016.
Collapse
Affiliation(s)
- Monika Staś
- Faculty of Chemistry, University of Opole, Oleska 48, Opole, 45-052, Poland
| | - Maciej Bujak
- Faculty of Chemistry, University of Opole, Oleska 48, Opole, 45-052, Poland
| | - Małgorzata A Broda
- Faculty of Chemistry, University of Opole, Oleska 48, Opole, 45-052, Poland
| | - Dawid Siodłak
- Faculty of Chemistry, University of Opole, Oleska 48, Opole, 45-052, Poland
| |
Collapse
|
18
|
Buczek A, Makowski M, Jewgiński M, Latajka R, Kupka T, Broda MA. Toward engineering efficient peptidomimetics. Screening conformational landscape of two modified dehydroaminoacids. Biopolymers 2016; 101:28-40. [PMID: 23606332 DOI: 10.1002/bip.22264] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/15/2013] [Accepted: 04/12/2013] [Indexed: 11/11/2022]
Abstract
Effective peptidomimetics should posses structural rigidity and appropriate interaction pattern leading to potential spatial and electronic matching to the target receptor site. Rational design of such small bioactive molecules could push chemical synthesis and molecular modeling toward faster progress in medicinal chemistry. Conformational properties of N-t-butoxycarbonyl-glycine-(E/Z)-dehydrophenylalanine N',N'-dimethylamides (Boc-Gly-(E/Z)-ΔPhe-NMe2 ) in chloroform were studied by NMR and IR spectroscopy. The experimental findings were supported by extensive calculations at DFT(B3LYP, M06-2X) and MP2 levels of theory and the β-turn tendency for both isomers of the studied dipeptide were determined in vacuum and in solution. The theoretical data and experimental IR results were used as an additional information for the NMR-based determination of the detailed solution conformations of the peptides. The obtained results reveal that N-methylation of C-terminal amide group changes dramatically the conformational properties of studied dehydropeptides. Theoretical conformational analysis reveals that the tendency to adopt β-turn conformations is much weaker for the N-methylated Z isomer (Boc-Gly-(Z)-ΔPhe-NMe2 ), both in vacuum and in polar environment. On the contrary, N-methylated E isomer (Boc-Gly-(E)-ΔPhe-NMe2 ) can easily adopt β-turn conformation, but the backbone torsion angles (φ1, ψ1, φ2, ψ2) are off the limits for common β-turn types.
Collapse
Affiliation(s)
- Aneta Buczek
- Department of Physical Chemistry and Molecular Modeling, Faculty of Chemistry, University of Opole, Oleska 48, 45-052, Opole, Poland
| | | | | | | | | | | |
Collapse
|
19
|
Kupka T, Buczek A, Broda MA, Stachów M, Tarnowski P. DFT studies on the structural and vibrational properties of polyenes. J Mol Model 2016; 22:101. [PMID: 27048200 PMCID: PMC4821862 DOI: 10.1007/s00894-016-2969-1] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/22/2015] [Accepted: 03/16/2016] [Indexed: 11/29/2022]
Abstract
Detailed density functional theory (DFT) calculations on the structure and harmonic frequencies of model all-trans and all-cis polyenes were undertaken. For the first time, we report on the convergence of selected B3LYP/6-311++G** and BLYP/6-311++G** calculated structural parameters resulting from a systematic increase in polyene size (chains containing 2 to 14 C = C units). The limiting values of the structural parameters for very long chains were estimated using simple three-parameter empirical formulae. BLYP/6-311++G** calculated ν(C = C) and ν(C–C) frequencies for all-trans and all-cis polyenes containing up to 14 carbon–carbon double bonds were used to estimate these values for very long chains. Correction of raw, unscaled vibrational data was performed by comparing theoretical and experimental wavenumbers for polyenes chains containing 3 to 12 conjugated C = C units with both ends substituted by tert-butyl groups. The corrected ν(C = C) and ν(C–C) wavenumbers for all-trans molecules were used to estimate the presence of 9 – 12 C = C units in all-trans polyene pigment in red coral. Detailed density functional theory (DFT) calculations on the structure and harmonic frequencies of model all-trans and all-cis polyenes were undertaken. For the first time, we report on the convergence of selected B3LYP/6-311++G** and BLYP/6-311++G** calculated structural parameters resulting from a systematic increase in polyene size (chains containing 2 to 14 C=C units). The limiting values of the structural parameters for very long chains were estimated using simple three-parameter empirical formulae. ![]()
Collapse
Affiliation(s)
- Teobald Kupka
- Faculty of Chemistry, University of Opole, 48 Oleska Street, 45-052, Opole, Poland.
| | - Aneta Buczek
- Faculty of Chemistry, University of Opole, 48 Oleska Street, 45-052, Opole, Poland.
| | - Małgorzata A Broda
- Faculty of Chemistry, University of Opole, 48 Oleska Street, 45-052, Opole, Poland
| | - Michał Stachów
- Faculty of Chemistry, University of Opole, 48 Oleska Street, 45-052, Opole, Poland
| | - Przemysław Tarnowski
- Faculty of Chemistry, University of Opole, 48 Oleska Street, 45-052, Opole, Poland
| |
Collapse
|
20
|
Kliber M, Broda MA, Nackiewicz J. Interactions of zinc octacarboxyphthalocyanine with selected amino acids and with albumin. Spectrochim Acta A Mol Biomol Spectrosc 2016; 155:54-60. [PMID: 26580510 DOI: 10.1016/j.saa.2015.11.003] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/24/2015] [Revised: 09/26/2015] [Accepted: 11/02/2015] [Indexed: 06/05/2023]
Abstract
Effect of selected amino acids (glycine, l-histidine, l-cysteine, l-serine, l-tryptophan) and albumin on the spectroscopic properties and photostability of zinc octacarboxyphthalocyanine (ZnPcOC) was explored in the phosphate buffer at a pH of 7.0. The photodegradation of ZnPcOC alone and in the presence of amino acids or albumin has been investigated in aqueous phase using UV-366nm and daylight irradiation. Kinetic analysis showed that the interaction with amino acids or albumin enhances the photostability of ZnPcOC. To answer the question of how zinc phthalocyanine interacts with amino acids extensive DFT calculations were performed. Analysis of the optimized geometry features of ZnPcOC: amino acids complexes in the gas phase and in water environment as well as the BSSE corrected interaction energies indicates that the more likely is the formation of equatorial complexes in which H-bonds are formed between the COOH groups of the phthalocyanine and carboxyl or amino groups of amino acids. UV-Vis spectra calculated by employing time dependent density functional theory (TD-DFT) are also consistent with this conclusion.
Collapse
Affiliation(s)
- Marta Kliber
- Department of Physical Chemistry and Molecular Modeling, Faculty of Chemistry, University of Opole, Oleska 48, Opole 45-052, Poland
| | - Małgorzata A Broda
- Department of Physical Chemistry and Molecular Modeling, Faculty of Chemistry, University of Opole, Oleska 48, Opole 45-052, Poland
| | - Joanna Nackiewicz
- Department of Physical Chemistry and Molecular Modeling, Faculty of Chemistry, University of Opole, Oleska 48, Opole 45-052, Poland.
| |
Collapse
|
21
|
Buczek A, Kupka T, Broda MA, Żyła A. Predicting the structure and vibrational frequencies of ethylene using harmonic and anharmonic approaches at the Kohn-Sham complete basis set limit. J Mol Model 2016; 22:42. [PMID: 26800989 PMCID: PMC4723629 DOI: 10.1007/s00894-015-2902-z] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/24/2015] [Accepted: 12/22/2015] [Indexed: 11/25/2022]
Abstract
In this work, regular convergence patterns of the structural, harmonic, and VPT2-calculated anharmonic vibrational parameters of ethylene towards the Kohn–Sham complete basis set (KS CBS) limit are demonstrated for the first time. The performance of the VPT2 scheme implemented using density functional theory (DFT-BLYP and DFT-B3LYP) in combination with two Pople basis sets (6-311++G** and 6-311++G(3df,2pd)), the polarization-consistent basis sets pc-n, aug-pc-n, and pcseg-n (n = 0, 1, 2, 3, 4), and the correlation-consistent basis sets cc-pVXZ and aug-cc-pVXZ (X = D, T, Q, 5, 6) was tested. The BLYP-calculated harmonic frequencies were found to be markedly closer than the B3LYP-calculated harmonic frequencies to the experimentally derived values, while the calculated anharmonic frequencies consistently underestimated the observed wavenumbers. The different basis set families gave very similar estimated values for the CBS parameters. The anharmonic frequencies calculated with B3LYP/aug-pc-3 were consistently significantly higher than those obtained with the pc-3 basis set; applying the aug-pcseg-n basis set family alleviated this problem. Utilization of B3LYP/aug-pcseg-n basis sets instead of B3LYP/aug-cc-pVXZ, which is computationally less expensive, is suggested for medium-sized molecules. Harmonic BLYP/pc-2 calculations produced fairly accurate ethylene frequencies. In this study, the performance of the VPT2 scheme implemented using density functional theory (DFT-BLYP and DFT-B3LYP) in combination with the polarization-consistent basis sets pc-n, aug-pc-n, and pcseg-n (n = 0, 1, 2, 3, 4), and the correlation-consistent basis sets cc-pVXZ and aug-cc-pVXZ (X = D, T, Q, 5, and 6) was tested. For the first time, we demonstrated regular convergence patterns of the structural, harmonic, and VPT2-calculated anharmonic vibrational parameters of ethylene towards the Kohn–Sham complete basis set (KS CBS) limit ![]()
Collapse
Affiliation(s)
- Aneta Buczek
- Faculty of Chemistry, University of Opole, 48, Oleska Street, 45-052, Opole, Poland.
| | - Teobald Kupka
- Faculty of Chemistry, University of Opole, 48, Oleska Street, 45-052, Opole, Poland.
| | - Małgorzata A Broda
- Faculty of Chemistry, University of Opole, 48, Oleska Street, 45-052, Opole, Poland
| | - Adriana Żyła
- Faculty of Chemistry, University of Opole, 48, Oleska Street, 45-052, Opole, Poland.,Department of Biosystematics, University of Opole, Oleska 22, 45-052, Opole, Poland.,Institute of Physics, Adam Mickiewicz University in Poznan, 65, Umultowska Street, 61-614, Poznań, Poland
| |
Collapse
|
22
|
Dyrda G, Słota R, Broda MA, Mele G. meso-Aryl-substituted free-base porphyrins: formation, structure and photostability of diprotonated species. Res Chem Intermed 2015. [DOI: 10.1007/s11164-015-2245-5] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
23
|
Wałęsa R, Man D, Engel G, Siodłak D, Kupka T, Ptak T, Broda MA. The Impact of Model Peptides on Structural and Dynamic Properties of Egg Yolk Lecithin Liposomes - Experimental and DFT Studies. Chem Biodivers 2015; 12:1007-24. [DOI: 10.1002/cbdv.201400179] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/09/2014] [Indexed: 12/21/2022]
|
24
|
Jaremko Ł, Jaremko M, Buczek A, Broda MA, Kupka T, Jackowski K. 1H and 13C shielding measurements in comparison with DFT calculations performed for two 2-(acetyloamino)-N,N-dimethyl-3-phenylacrylamide isomers. Chem Phys Lett 2015. [DOI: 10.1016/j.cplett.2015.03.021] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
|
25
|
Wałęsa R, Ptak T, Siodłak D, Kupka T, Broda MA. Experimental and theoretical NMR studies of interaction between phenylalanine derivative and egg yolk lecithin. Magn Reson Chem 2014; 52:298-305. [PMID: 24639342 DOI: 10.1002/mrc.4064] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/19/2013] [Revised: 01/18/2014] [Accepted: 02/28/2014] [Indexed: 06/03/2023]
Abstract
The interaction of phenylalanine diamide (Ac-Phe-NHMe) with egg yolk lecithin (EYL) in chloroform was studied by (1)H and (13)C NMR. Six complexes EYL-Ac-Phe-NHMe, stabilized by N-H···O or/and C-H···O hydrogen bonds, were optimized at M06-2X/6-31G(d,p) level. The assignment of EYL and Ac-Phe-NHMe NMR signals was supported using GIAO (gauge including atomic orbital) NMR calculations at VSXC and B3LYP level of theory combined with STO-3Gmag basis set. Results of our study indicate that the interaction of peptides with lecithin occurs mainly in the polar 'head' of the lecithin. Additionally, the most probable lecithin site of H-bond interaction with Ac-Phe-NHMe is the negatively charged oxygen in phosphate group that acts as proton acceptor.
Collapse
Affiliation(s)
- Roksana Wałęsa
- Faculty of Chemistry, University of Opole, Oleska 48 Str., 45-052, Opole, Poland
| | | | | | | | | |
Collapse
|
26
|
Abstract
Dehydroamino acids with the methylated N-terminal peptide group occur in natural small cyclic peptides. The structural analysis was used to investigate the cis-trans isomerization of the N-terminal tertiary amide group of diamides: Ac-(Z)-Δ(Me)Abu-NHMe (1), Ac-(Z)-Δ(Me)Phe-NHMe (2), Ac-(E)-Δ(Me)Phe-NHMe (3), Ac-Δ(Me)Ala-NHMe (4), and Ac-(Me)Ala-NHMe (5). The compounds were analyzed in the solid state by an X-ray crystallography (1-3), and in the solution by FTIR (MeCN and CHCl(3) ) and NMR (DMSO-d6 and CDCl(3) ) methods (1-5). In the solid state, the studied compounds adopt the cis configuration of N-terminal amide. In solution, this configuration also prevails for the dehydroamino acids 1-4, in contrast to the saturated analog 5. The results indicate that N-methyldehydroamino acids present a promising tool to induce the cis configuration of the amide bond.
Collapse
Affiliation(s)
- Dawid Siodłak
- Faculty of Chemistry, University of Opole, Oleska 48, Opole 45-052, Poland.
| | | | | | | | | |
Collapse
|
27
|
Siodłak D, Staś M, Broda MA, Bujak M, Lis T. Conformational properties of oxazole-amino acids: effect of the intramolecular N-H···N hydrogen bond. J Phys Chem B 2014; 118:2340-50. [PMID: 24528177 DOI: 10.1021/jp4121673] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
Oxazole ring occurs in numerous natural peptides, but conformational properties of the amino acid residue containing the oxazole ring in place of the C-terminal amide bond are poorly recognized. A series of model compounds constituted by the oxazole-amino acids occurring in nature, that is, oxazole-alanine (L-Ala-Ozl), oxazole-dehydroalanine (ΔAla-Ozl), and oxazole-dehydrobutyrine ((Z)-ΔAbu-Ozl), was investigated using theoretical calculations supported by FTIR and NMR spectra and single-crystal X-ray diffraction. It was found that the main feature of the studied oxazole-amino acids is the stable conformation β2 with the torsion angles φ and ψ of -150°, -10° for L-Ala-Ozl, -180°, 0° for ΔAla-Ozl, and -120°, 0° for (Z)-ΔAbu-Ozl, respectively. The conformation β2 is stabilized by the intramolecular N-H···N hydrogen bond and predominates in the low polar environment. In the case of the oxazole-dehydroamino acids, the π-electron conjugation that is spread on the oxazole ring and C(α)═C(β) double bond is an additional stabilizing interaction. The tendency to adopt the conformation β2 clearly decreases with increasing the polarity of environment, but still the oxazole-dehydroamino acids are considered to be more rigid and resistant to conformational changes.
Collapse
Affiliation(s)
- Dawid Siodłak
- Faculty of Chemistry, University of Opole , Oleska 48, 45-052 Opole, Poland
| | | | | | | | | |
Collapse
|
28
|
Wałęsa R, Broda MA. Solvent effects on the conformational preferences of model peptoids. MP2 study. J Pept Sci 2014; 20:203-11. [PMID: 24520015 DOI: 10.1002/psc.2601] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/23/2013] [Revised: 11/12/2013] [Accepted: 11/22/2013] [Indexed: 01/04/2023]
Abstract
The influence of aqueous environment on the main-chain conformation (ω0 , ϕ, and ψ dihedral angles) of two model peptoids: N-acetyl-N-methylglycine N'-methylamide (Ac-N(Me)-Gly-NHMe) (1) and N-acetyl-N-methylglycine N',N'-dimethylamide (Ac-N(Me)-Gly-NMe₂) (2) was investigated by MP2/6-311++G(d,p) method. The Ramachandran maps of both studied molecules with cis and trans configuration of the N-terminal amide bond in the gas phase and in water environment were obtained and all energy minima localized. The polarizable continuum model was applied to estimate the solvation effect on conformation. Energy minima of the Ac-N(Me)-Gly-NHMe and Ac-N(Me)-Gly-NMe₂ have been analyzed in terms of the possible hydrogen bonds and C = O dipole attraction. To validate the theoretical results obtained, conformations of the similar structures gathered in the Cambridge Crystallographic Data Centre were analyzed. Obtained results indicate that aqueous environment in model peptoids 1 and 2 favors the conformation F (ϕ and ψ = -70º, 180º), and additionally significantly increases the percentage of structures with cis configuration of N-terminal amide bond in studied compounds.
Collapse
|
29
|
|
30
|
Buczek A, Wałęsa R, Broda MA. β-turn tendency in N-methylated peptides with dehydrophenylalanine residue: DFT study. Biopolymers 2012; 97:518-28. [DOI: 10.1002/bip.22034] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/05/2010] [Revised: 02/02/2011] [Accepted: 02/04/2011] [Indexed: 01/10/2023]
|
31
|
Słota R, Broda MA, Dyrda G, Ejsmont K, Mele G. Structural and molecular characterization of meso-substituted zinc porphyrins: a DFT supported study. Molecules 2011; 16:9957-71. [PMID: 22134401 PMCID: PMC6264692 DOI: 10.3390/molecules16129957] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/02/2011] [Revised: 11/22/2011] [Accepted: 11/22/2011] [Indexed: 11/16/2022] Open
Abstract
Structural parameters of a range of over 100 meso-substituted zinc porphyrins were reviewed and compared to show how far the nature of the functional group may affect the interatomic distances and bond angles within the porphyrin core. It was proved that even despite evident deformations of the molecular structure, involving twisting of the porphyrin's central plane, the coupled π-bonding system remains flexible and stable. DFT calculations were applied to a number of selected porphyrins representative for the reviewed compounds to emphasize the relevance of theoretical methods in structural investigations of complex macrocyclic molecular systems. Experimental and DFT-simulated IR spectral data were reported and analyzed in context of the individual molecular features introduced by the meso substituents into the porphyrin moiety base. Raw experimental spectral data, including 1H- and 13C-NMR, UV-Vis, FTIR, XRD, and other relevant physicochemical details have been provided for a specially chosen reference zinc porphyrin functionalized by tert-butylphenyl groups.
Collapse
Affiliation(s)
- Rudolf Słota
- Faculty of Chemistry, Opole University, ul. Oleska 48, 45-052 Opole, Poland.
| | | | | | | | | |
Collapse
|
32
|
Siodłak D, Grondys J, Broda MA. The conformational properties of α,β-dehydroamino acids with a C-terminal ester group. J Pept Sci 2011; 17:690-9. [PMID: 21805538 DOI: 10.1002/psc.1390] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/05/2011] [Revised: 05/31/2011] [Accepted: 06/03/2011] [Indexed: 11/06/2022]
Abstract
α,β-Dehydroamino acid esters occur in nature. To investigate their conformational properties, a systematic theoretical analysis was performed on the model molecules Ac-ΔXaa-OMe [ΔXaa = ΔAla, (E)-ΔAbu, (Z)-ΔAbu, ΔVal] at the B3LYP/6-311+ + G(d,p) level in the gas phase as well as in chloroform and water solutions with the self-consistent reaction field-polarisable continuum model method. The Fourier transform IR spectra in CCl(4) and CHCl(3) have been analysed as well as the analogous solid state conformations drawn from The Cambridge Structural Database. The ΔAla residue has a considerable tendency to adopt planar conformations C5 (ϕ, ψ ≈ - 180°, 180°) and β2 (ϕ, ψ ≈ - 180°, 0°), regardless of the environment. The ΔVal residue prefers the conformation β2 (ϕ, ψ ≈ - 120°, 0°) in a low polar environment, but the conformations α (ϕ, ψ ≈ - 55°, 35°) and β (ϕ, ψ ≈ - 55°, 145°) when the polarity increases. The ΔAbu residues reveal intermediate properties, but their conformational dispositions depend on configuration of the side chain of residue: (E)-ΔAbu is similar to ΔAla, whereas (Z)-ΔAbu to ΔVal. Results indicate that the low-energy conformation β2 is the characteristic feature of dehydroamino acid esters. The studied molecules constitute conformational patterns for dehydroamino acid esters with various side chain substituents in either or both Z and E positions.
Collapse
Affiliation(s)
- Dawid Siodłak
- Faculty of Chemistry, Opole University, Oleska 48, 45-052 Opole, Poland.
| | | | | |
Collapse
|
33
|
Buczek AM, Ptak T, Kupka T, Broda MA. Experimental and theoretical NMR and IR studies of the side-chain orientation effects on the backbone conformation of dehydrophenylalanine residue. Magn Reson Chem 2011; 49:343-349. [PMID: 21506144 DOI: 10.1002/mrc.2753] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/28/2011] [Revised: 02/09/2011] [Accepted: 02/16/2011] [Indexed: 05/30/2023]
Abstract
Conformation of N-acetyl-(E)-dehydrophenylalanine N', N'-dimethylamide (Ac-(E)-ΔPhe-NMe(2)) in solution, a member of (E)-α, β-dehydroamino acids, was studied by NMR and infrared spectroscopy and the results were compared with those obtained for (Z) isomer. To support the spectroscopic interpretation, the Φ, Ψ potential energy surfaces were calculated at the MP2/6-31 + G(d,p) level of theory in chloroform solution modeled by the self-consistent reaction field-polarizable continuum model method. All minima were fully optimized by the MP2 method and their relative stabilities were analyzed in terms of π-conjugation, internal H-bonds and dipole interactions between carbonyl groups. The obtained NMR spectral features were compared with theoretical nuclear magnetic shieldings, calculated using Gauge Independent Atomic Orbitals (GIAO) approach and rescaled to theoretical chemical shifts using benzene as reference. The calculated indirect nuclear spin-spin coupling constants were compared with available experimental parameters.
Collapse
Affiliation(s)
- Aneta M Buczek
- Faculty of Chemistry, University of Opole, 48, Oleska Street, 45-052 Opole, Poland
| | | | | | | |
Collapse
|
34
|
Buczek A, Siodłak D, Bujak M, Broda MA. Effects of side-chain orientation on the backbone conformation of the dehydrophenylalanine residue. Theoretical and X-ray study. J Phys Chem B 2011; 115:4295-306. [PMID: 21443240 DOI: 10.1021/jp200949t] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Two E isomers of α,β-dehydro-phenylalanine, Ac-(E)-ΔPhe-NHMe (1a) and Ac-(E)-ΔPhe-NMe(2) (2a), have been synthesized and their low temperature structures determined by single-crystal X-ray diffraction. A systematic theoretical analysis was performed on these molecules and their Z isomers (1b and 2b). The ϕ,ψ potential energy surfaces were calculated at the MP2/6-31+G(d,p) and B3LYP/6-31+G(d,p) levels in the gas phase and at the B3LYP/6-31+G(d,p) level in the chloroform and water solutions with the SCRF-PCM method. All minima were fully optimized by the MP2 and DFT methods, and their relative stabilities were analyzed in terms of π-conjugation, internal H-bonds, and dipole-dipole interactions between carbonyl groups. The results indicate that all the studied compounds can adopt the conformation H (ϕ, ψ ≈ ±40°, ∓120°) which is atypical for standard amino acids residues. A different arrangement of the side chain in the E and Z isomers causes them to have different conformational preferences. In the presence of a polar solvent both Z isomers of ΔPhe (1b and 2b) are found to adopt the 3(10)-helical conformation (left- and right-handed are equally likely). On the other hand, this conformation is not accessible or highly energetic for E isomers of ΔPhe (1a and 2a). Those isomers have an intrinsic inclination to have an extended conformation. The conformational space of the Z isomers is much more restricted than that of the E derivative both in the gas phase and in solution. In the gas phase the E isomers of ΔPhe have lower energies than the Z ones, but in the aqueous solution the energy order is reversed.
Collapse
Affiliation(s)
- Aneta Buczek
- Faculty of Chemistry, University of Opole, Poland
| | | | | | | |
Collapse
|
35
|
Abstract
Liposomes characterized by membranes featuring diverse fluidity (liquid-crystalline and/or gel phase), prepared from egg yolk lecithin (EYL) and dipalmitoylphosphatidylcholine (DPPC), were doped with selected metalloporphyrins and the time-related structural and dynamic changes within the lipid double layer were investigated. Porphyrin complexes of Mg(II), Mn(III), Fe(III), Co(II), Ni(II), Cu(II), Zn(II), and the metal-free base were embedded into the particular liposome systems and tested for 350 h at 24°C using the electron spin resonance (ESR) spin probe technique. 5-DOXYL, 12-DOXYL, and 16-DOXYL stearic acid methyl ester spin labels were applied to explore the interior of the lipid bilayer. Only the 16-DOXYL spin probe detected evident structural changes inside the lipid system due to porphyrin intercalation. Fluidity of the lipid system and the type of the porphyrin complex introduced significantly affected the intermolecular interactions, which in certain cases may result in self-assembly of metalloporphyrin molecules within the liposome membrane, reflected in the presence of new lines in the relevant ESR spectra. The most pronounced time-related effects were demonstrated by the EYL liposomes (liquid-crystalline phase) when doped with Mg and Co porphyrins, whereas practically no spectral changes were revealed for the metal-free base and both the Ni and Zn dopants. ESR spectra of the porphyrin-doped gel phase of DPPC liposomes did not show any extra lines; however, they indicated the formation of a more rigid lipid medium. Electronic configuration of the porphyrin's metal center appeared crucial to the degree of molecular reorganization within the phospholipid bilayer system.
Collapse
Affiliation(s)
- Dariusz Man
- Institute of Physics, Opole University, Oleska 48, 45-052 Opole, Poland
| | - Rudolf Słota
- Faculty of Chemistry, Opole University, Oleska 48, 45-052 Opole, Poland
| | | | - Giuseppe Mele
- Dipartimento di Ingegneria dell’Innovazione, Università del Salento, Via Arnesano, 73100 Lecce, Italy
| | - Jun Li
- Key Laboratory of Synthetic and Natural Functional Molecule Chemistry of Ministry of Education, The College of Chemistry and Materials Science, Northwest University, Xi’an, 710069 China
| |
Collapse
|
36
|
Buczek A, Kupka T, Broda MA. Extrapolation of water and formaldehyde harmonic and anharmonic frequencies to the B3LYP/CBS limit using polarization consistent basis sets. J Mol Model 2010; 17:2029-40. [PMID: 21153906 PMCID: PMC3143340 DOI: 10.1007/s00894-010-0913-3] [Citation(s) in RCA: 25] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/16/2010] [Accepted: 11/21/2010] [Indexed: 11/12/2022]
Abstract
The harmonic and anharmonic frequencies of fundamental vibrations in formaldehyde and water were successfully estimated using the B3LYP Kohn-Sham limit. The results obtained with polarization- and correlation-consistent basis sets were fitted with a two-parameter formula. Anharmonic corrections were obtained by a second order perturbation treatment (PT2). We compared the performance of the PT2 scheme on the two title molecules using SCF, MP2 and DFT (BLYP, B3LYP, PBE and B3PW91 functionals) methods combined with polarization consistent pc-n (n = 0, 1, 2, 3, 4) basis sets, Dunning’s basis sets (aug)-cc-pVXZ where X = D, T, Q, 5, 6 and Pople’s basis sets up to 6-311++G(3df,2pd). The influence of SCF convergence level and density grid size on the root mean square of harmonic and anharmonic frequency deviations from experimental values was tested. The wavenumber of formaldehyde CH2 anharmonic asymmetric stretching mode is very sensitive to grid size for large basis sets; this effect is not observed for harmonic modes. BLYP-calculated anharmonic frequencies consistently underestimate observed wavenumbers. On the basis of formaldehyde anharmonic frequencies, we show that increasing the Pople basis set size does not always lead to improved agreement between anharmonic frequencies and experimental values. Sensitivity of water B3LYP calculated harmonic and anharmonic vs(OH) frequencies on selected Pople and polarization consistent basis sets size. The results for pc-n basis sets were fitted with two parameter formula and the CBS(2,3,4) estimated ![]()
Collapse
Affiliation(s)
- Aneta Buczek
- Faculty of Chemistry, University of Opole, 48 Oleska Street, 45-052 Opole, Poland
| | | | | |
Collapse
|
37
|
Boniewska-Bernacka E, Man D, Słota R, Broda MA. Effect of tin and lead chlorotriphenyl analogues on selected living cells. J Biochem Mol Toxicol 2010; 25:231-7. [DOI: 10.1002/jbt.20380] [Citation(s) in RCA: 13] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/12/2010] [Revised: 09/22/2010] [Accepted: 10/19/2010] [Indexed: 11/05/2022]
|
38
|
Siodłak D, Grondys J, Lis T, Bujak M, Broda MA, Rzeszotarska B. The conformational properties of dehydrobutyrine and dehydrovaline: theoretical and solid-state conformational studies. J Pept Sci 2010; 16:496-505. [PMID: 20645424 DOI: 10.1002/psc.1267] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Abstract
Dehydrobutyrine is the most naturally occurring dehydroamino acid. It is also the simplest dehydroamino acid having the geometrical isomers E/Z. To investigate its conformational properties, a theoretical analysis was performed on N-acetyl-alpha,beta-dehydrobutyrine N'-methylamides, Ac-(E)-DeltaAbu-NHMe and Ac-(Z)-DeltaAbu-NHMe, as well as the dehydrovaline derivative Ac-DeltaVal-NHMe. The phi, psi potential energy surfaces and the localised conformers were calculated at the B3LYP/6-311 + + G(d,p) level of theory both in vacuo and with inclusion of the solvent (chloroform, water) effect (SCRF method). The X-ray crystal structures of Ac-(Z)-DeltaAbu-NHMe and Ac-DeltaVal-NHMe were determined at 85 and 100 K, respectively. The solid-state conformational preferences for the studied residues have been analysed and compared with the other related structures. Despite the limitations imposed by the C(alpha) = C(beta) double bond on the topography of the side chains, the main chains of the studied dehydroamino acids are more flexible than in standard alanine. The studied dehydroamino acids differ in their conformational preferences, which depend on the polarity of the environment. This might be a reason why the nature quite precisely differentiates between DeltaVal and each of the DeltaAbu isomers, and why, particularly so with the latter, they are used as a conformational tool to influence the biological action of usually small, cyclic dehydropeptides.
Collapse
Affiliation(s)
- Dawid Siodłak
- Faculty of Chemistry, University of Opole, Oleska 48, 45-052 Opole, Poland.
| | | | | | | | | | | |
Collapse
|
39
|
Broda MA, Buczek A, Siodłak D, Rzeszotarska B. The effect of beta-methylation on the conformation of alpha, beta-dehydrophenylalanine: a DFT study. J Pept Sci 2009; 15:465-73. [PMID: 19408247 DOI: 10.1002/psc.1137] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
Abstract
Dehydroamino acids are non-coded amino acids that offer unique conformational properties. Dehydrophenylalanine (DeltaPhe) is most commonly used to modify bioactive peptides to constrain the topography of the phenyl ring in the side chain, which commonly serves as a pharmacophore. The Ramachandran maps (in the gas phase and in CHCl(3) mimicking environments) of DeltaPhe analogues with methyl groups at the beta position of the side chain as well as at the C-terminal amide were calculated using the B3LYP/6-31 + G** method. Unexpectedly, beta-methylation alone results in an increase of conformational freedom of the affected DeltaPhe residue. However, further modification by introducing an additional methyl group at C-terminal methyl amide results in a steric crowding that fixes the torsion angle psi of all conformers to the value 123 degrees , regardless of the Z or E position of the phenyl ring. The number of conformers is reduced and the accessible conformational space of the residues is very limited. In particular, (Z)-Delta(betaMe)Phe with the tertiary C-terminal amide can be classified as the amino acid derivative that has a single conformational state as it seems to adopt only the beta conformation.
Collapse
|
40
|
Broda MA, Rzeszotarska B. Association of model peptides and dehydropeptides: N-acetyl-L-butyrine and (Z)-dehydrobutyrine N',N'-dimethylamides. Acta Biochim Pol 2006. [DOI: 10.18388/abp.2006_3381] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [What about the content of this article? (0)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
Abstract
These comparative studies on the aggregation behaviour of Ac-(Z)-DeltaAbu-NMe(2) and Ac-L-Abu-NMe(2) in carbon tetrachloride were performed by the analysis of their FTIR spectra and by theoretical calculations. The percentage of the monomeric form (alpha) decreased as concentration increased and this occurred to a higher degree for the (Z)-DeltaAbu derivative than for its saturated analogue. The dimerization constant K(D), calculated on the basis of the intensity of the monomer and associate bands in the nu(s)(N-H) vibration region, is by three orders of magnitude larger for Ac-(Z)-DeltaAbu-NMe(2) than for Ac-L-Abu-NMe(2). The obtained dimer geometries of the dehydro- compound were calculated by the B3LYP/6-31+G** method.
Collapse
|
41
|
Broda MA, Rzeszotarska B. Association of model peptides and dehydropeptides: N-acetyl-L-butyrine and (Z)-dehydrobutyrine N',N'-dimethylamides. Acta Biochim Pol 2006; 53:221-6. [PMID: 16528413] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/06/2005] [Revised: 02/16/2006] [Accepted: 03/02/2006] [Indexed: 05/07/2023]
Abstract
These comparative studies on the aggregation behaviour of Ac-(Z)-DeltaAbu-NMe(2) and Ac-L-Abu-NMe(2) in carbon tetrachloride were performed by the analysis of their FTIR spectra and by theoretical calculations. The percentage of the monomeric form (alpha) decreased as concentration increased and this occurred to a higher degree for the (Z)-DeltaAbu derivative than for its saturated analogue. The dimerization constant K(D), calculated on the basis of the intensity of the monomer and associate bands in the nu(s)(N-H) vibration region, is by three orders of magnitude larger for Ac-(Z)-DeltaAbu-NMe(2) than for Ac-L-Abu-NMe(2). The obtained dimer geometries of the dehydro- compound were calculated by the B3LYP/6-31+G** method.
Collapse
|
42
|
Broda MA, Siodłak D, Rzeszotarska B. Conformational investigation of alpha,beta-dehydropeptides. XV: N-acetyl-alpha,beta-dehydroamino acid N 'N '-dimethylamides: conformational properties from infrared and theoretical studies. J Pept Sci 2005; 11:546-55. [PMID: 15782429 DOI: 10.1002/psc.655] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Abstract
The FTIR spectra were analysed in the region of the nu(s)(N-H), AI(C=O) and nu(s)(Calpha=Cbeta) bands for a series of Ac-DeltaXaa-NMe2, where DeltaXaa = DeltaAla, (Z)-DeltaAbu, (Z)-DeltaLeu, (Z)-DeltaPhe and DeltaVal, to determine a predominant solution conformation of these alpha,beta-dehydropeptide-related molecules. Measurements were taken in CCl4, DCM and MeCN solutions. In the same way, spectra of saturated analogues Ac-Xaa-NMe2, where Xaa = Ala, Abu, Leu, Phe and Val, were investigated. To help interpret the spectroscopic results, conformational maps were calculated by the B3LYP/6-31+G** method. Also, the relative energies of all conformers of the dehydro compounds in vacuo as well as in the studied solvents in addition to the theoretical IR frequencies of these conformers were calculated. For comparison, molecules of two saturated analogues, Ac-L-Ala-NMe2 and Ac-L-Phe-NMe2, were calculated in a similar way. Both unsaturated and saturated compounds, which have an aliphatic side chain, occur in CCl4 and DCM mainly as a mixture of extended conformers with the C5 H-bond and open conformers. As solvent polarity increases, participation of the open conformers also increases, and in MeCN, the model amides are almost exclusively in the open form, except Ac-DeltaAla-NMe2, which shows a small amount of the H-bonded conformer. Ac-DeltaAla-NMe2 and Ac-DeltaAbu-NMe2 have stronger C5 hydrogen bonds than those of their saturated counterparts. As the calculations indicate, the open conformation of the unsaturated amides is conformer H/F with phi, psi -44 +/- 5 degrees, 127 +/- 4 degrees. This is the second lowest in energy conformer in vacuo and in CCl4 and the lowest one in more polar solvents. The open conformation of Ac-L-Ala-NMe2 constitutes conformer C with phi, psi -101.5 degrees, 112.7 degrees. For Ac-DeltaAla-NMe2 and Ac-DeltaAbu-NMe2, FTIR also reveals the presence of a third conformer. Calculations indicate that is the semiextended conformer D with the N1-H1...N2 hydrogen bond/contact. In all solvents, Ac-L-Phe-NMe2 and Ac-(Z)-DeltaPhe-NMe2 show only the extended E and the open H/F, respectively. In both there is an amide/pi(Ph) interaction.
Collapse
|
43
|
Broda MA, Siodłak D, Rzeszotarska B. Conformational investigation of alpha,beta-dehydropeptides. N-acetyl-(E)-dehydrophenylalanine N'-methylamide: conformational properties from infrared and theoretical studies, part XIV. J Pept Sci 2005; 11:235-44. [PMID: 15635637 DOI: 10.1002/psc.601] [Citation(s) in RCA: 13] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Abstract
N-Acetyl-(E)-dehydrophenylalanine N'-methylamide [Ac-(E)-DeltaPhe-NHMe], one of a few representative (E)-alpha,beta-dehydroamino acids, was studied by FTIR in dichloromethane and acetonitrile. To support spectroscopic interpretations and to gain some deeper insight into the Ac-(E)-DeltaPhe-NHMe molecule, the Ramachandran potential energy surface was calculated by the B3LYP/6-31G*//HF/3-21G method and the conformers localized were fully optimized at the B3LYP/6-31 + G** level. The spectra and calculations were compared with those of the related molecules Ac-DeltaAla-NHMe and Ac-(Z)-DeltaPhe-NHMe. The title compound assumes two conformational states in equilibrium in dichloromethane solution with a predominance of the extended conformer E. The Ac-(E)-DeltaPhe-NHMe spectrum is like that of Ac-DeltaAla-NHMe, particularly in the region of bands AI and AII, and unlike that of Ac-(Z)-DeltaPhe-NHMe. The positions of bands AI and II together with the nu(s)(N1--H1) band proves that the conformers E of both DeltaAla and (E)-DeltaPhe compounds are stabilized by the quite strong C5 hydrogen bonds N1--H1...O2. The same conclusion is drawn from the Ramachandran diagrams. The conformers E of both compounds are placed in the global minima and the gaps in energy order between them and the second conformer are large. The conformers E of DeltaAla and (E)-DeltaPhe, apart from the N1--H1...O2 hydrogen bond, show the Cbeta--H...O1 interaction, and Ac-(E)-DeltaPhe-NHMe displays the NH/pi interaction with the N2--H2 projecting in the first carbon atom of the phenyl ring. The C5 hydrogen bond is stronger in (E)-DeltaPhe than that in the DeltaAla compound. This is in agreement with interactions found in the calculated structures and can be explained by the influence of the phenyl ring in position (E). In acetonitrile, the molecule of Ac-(E)-DeltaPhe-NHMe loses its C5 hydrogen bond and becomes unfolded, whereas that of Ac-DeltaAla-NHMe does not vary practically. Adopting conformation E in a non-polar solvent seems to be a general feature of the (E)-DeltaXaa residues.
Collapse
|
44
|
Broda MA, Rospenk M, Siodłak D, Rzeszotarska B. Association of model peptides and dehydropeptides: N-acetyl-l-alanine- and dehydroalanine N′,N′-dimethylamides. J Mol Struct 2005. [DOI: 10.1016/j.molstruc.2004.12.045] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
|
45
|
Abstract
Ab initio/DFT analysis of the conformational properties of free Ac-Ala-NMe(2) (N-acetyl-L-alanine-N',N'-dimethylamide) in terms of the N-H.O, N-H.N, C-H.O hydrogen bonds and C(delta+) = O(delta-) dipole attractions was performed. The Ala residue combined with the C-terminal tertiary amide prefers an extended conformation and that characteristic of the (i + 1)th position of the betaVIb turn. These can be easily remodelled into a structure compatible with the (i + 1)th position of the betaII/betaVIa turn. The residue has also the potential to adopt the conformation accommodated at both central positions of the betaIII/betaIII' turn or the (i + 1)th position of the betaI/beta'I turn.
Collapse
|
46
|
Siodłak D, Rzeszotarska B, Broda MA, Kozioł AE, Kołodziejczyk E. Conformational investigation of alpha,beta-dehydropeptides. XIII. Conformational properties of N-acetyl-alpha,beta-dehydrovaline N',N'-dimethylamide. Acta Biochim Pol 2004. [DOI: 10.18388/abp.2004_3606] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
Abstract
The crystal structure of Ac-DeltaVal-NMe(2) (DeltaVal = alpha,beta-dehydrovaline) was determined by X-ray crystallography. The found angles phi = -60 degrees and psi = 125 degrees correspond exactly to the respective values of the (i + 1)th residue in idealised beta-turn II/VIa. Ab initio/DFT studies revealed that the molecule adopts the angle psi restricted only to about |130 degrees | and very readily attains the angle phi = about -50 degrees. This is in line with its solid-state conformation. Taken together, these data suggest that the DeltaVal residue combined with a C-terminal tertiary amide is a good candidate at the (i + 1)th position in a type II/VIa beta-turn.
Collapse
|
47
|
Siodłak D, Rzeszotarska B, Broda MA. Conformational properties of N-acetyl-L-alanine N',N'-dimethylamide. Acta Biochim Pol 2004; 51:137-143. [PMID: 15094834] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/31/2003] [Revised: 02/12/2004] [Accepted: 02/18/2004] [Indexed: 05/24/2023]
Abstract
Ab initio/DFT analysis of the conformational properties of free Ac-Ala-NMe(2) (N-acetyl-L-alanine-N',N'-dimethylamide) in terms of the N-H.O, N-H.N, C-H.O hydrogen bonds and C(delta+) = O(delta-) dipole attractions was performed. The Ala residue combined with the C-terminal tertiary amide prefers an extended conformation and that characteristic of the (i + 1)th position of the betaVIb turn. These can be easily remodelled into a structure compatible with the (i + 1)th position of the betaII/betaVIa turn. The residue has also the potential to adopt the conformation accommodated at both central positions of the betaIII/betaIII' turn or the (i + 1)th position of the betaI/beta'I turn.
Collapse
Affiliation(s)
- Dawid Siodłak
- Institute of Chemistry, University of Opole, Opole, Poland
| | | | | |
Collapse
|
48
|
Siodłak D, Rzeszotarska B, Broda MA, Kozioł AE, Kołodziejczyk E. Conformational investigation of alpha,beta-dehydropeptides. XIII. Conformational properties of N-acetyl-alpha,beta-dehydrovaline N',N'-dimethylamide. Acta Biochim Pol 2004; 51:145-152. [PMID: 15094835] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/31/2003] [Revised: 02/17/2004] [Accepted: 03/08/2004] [Indexed: 05/24/2023]
Abstract
The crystal structure of Ac-DeltaVal-NMe(2) (DeltaVal = alpha,beta-dehydrovaline) was determined by X-ray crystallography. The found angles phi = -60 degrees and psi = 125 degrees correspond exactly to the respective values of the (i + 1)th residue in idealised beta-turn II/VIa. Ab initio/DFT studies revealed that the molecule adopts the angle psi restricted only to about |130 degrees | and very readily attains the angle phi = about -50 degrees. This is in line with its solid-state conformation. Taken together, these data suggest that the DeltaVal residue combined with a C-terminal tertiary amide is a good candidate at the (i + 1)th position in a type II/VIa beta-turn.
Collapse
Affiliation(s)
- Dawid Siodłak
- Institute of Chemistry, University of Opole, Opole, Poland
| | | | | | | | | |
Collapse
|
49
|
Slodłak D, Broda MA, Rzeszotarska B, Dybała I, Kozioł AE. Conformational investigation of alpha, beta-dehydropeptides. XI. Molecular and crystal structure of Ac-(Z)-deltaPhe-NMe2 as compared to those of related molecules. J Pept Sci 2003; 9:64-74. [PMID: 12587884 DOI: 10.1002/psc.433] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
A series of three homologous dimethyldiamides Ac-(Z)-deltaPhe-NMe2, Ac-L-Phe-NMe2 and Ac-DL-Phe-NMe2 have been synthesized and their structures determined from single-crystal X-ray diffraction data. To learn more about the conformational preferences of the compounds studied, the fully relaxed phi, psi conformational energy maps on the free molecules of Ac-deltaAla-NMe2 and Ac-(Z)-deltaPhe-NMe2 were obtained with the HF/3-21G method and the calculated minima re-optimized with the DFT/B3LYP/6-31G** method. The crystal state results have been compared with the literature data. The studied dimethyldiamide Ac-deltaXaa-NMe2 combines the double bond in positions alpha, beta and the C-terminal tertiary amide within one molecule. As the representative probe with deltaXaa = deltaAla, (Z)-deltaLeu and (Z)-deltaPhe shows, in the solid state they adopt the conservative conformation with phi, psi approximately -45 degrees, approximately 130 degrees and with a non-planar tertiary amide bond, whatever the packing forces are. This conformation is located on the Ramachandran map in region H/F, which is of high-energy for common amino acids, but not so readily accessible to them. The free molecule calculations on Ac-deltaAla-NMe2 and Ac-(Z)-deltaPhe-NMe2 reveal that, in spite of dissimilar overall conformational profiles of these molecules, this structure is one of their low-energy conformers and for Ac-(Z)-deltaPhe-NMe2 it constitutes the global minimum. So, the theoretical results corroborate those experimental results proving that this structure is robust enough to avoid conformational distortion due to packing forces. In contrast to Ac-deltaXaa-NMe2, the saturated Ac-L/DL-Xaa-NMe2 shows the constancy of the associative patterns but do not prefer any molecular structure in the solid state.
Collapse
Affiliation(s)
- Dawid Slodłak
- Institute of Chemistry, University of Opole, 45-052 Opole, Poland
| | | | | | | | | |
Collapse
|
50
|
Siodłak D, Broda MA, Rzeszotarska B, Kozioł AE, Dybała I. Conformational properties of N',N'-dimethylamides of N-acetyldehydroalanine and N-acetyl-(Z)-dehydrophenylalanine. Acta Biochim Pol 2002; 48:1179-83. [PMID: 11995989] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/24/2023]
Abstract
Conformational preferences of Ac-deltaAla-NMe2 and Ac-(Z)-deltaPhe-NMe2 were studied and compared with those of their monomethyl counterparts as well as with those of their saturated analogues. X-Ray data and energy calculations revealed a highly conservative conformation of the dehydro dimethylamides, which is located in a high-energy region of the Ramachandran map.
Collapse
Affiliation(s)
- D Siodłak
- Institute of Chemistry, University of Opole, Poland
| | | | | | | | | |
Collapse
|