Leave a message

Molecular Structure and Thermodynamic Properties of Morphine, Ibuprofen and Aspirin Using DFT Method

Research Article | DOI: https://doi.org/10.31579/2690-8808/246

Molecular Structure and Thermodynamic Properties of Morphine, Ibuprofen and Aspirin Using DFT Method

  • Hacer Gümüş

Automotive Technology Program, Golcuk Vocational School, Kocaeli University, 41380, Kocaeli, Turkey.

*Corresponding Author: Hacer Gümüş, Automotive Technology Program, Golcuk Vocational School, Kocaeli University, 41380, Kocaeli, Turkey.

Citation: Hacer Gümüş, (2025), Molecular Structure and Thermodynamic Properties of Morphine, Ibuprofen and Aspirin using DFT Method, J, Clinical Case Reports and Studies, 6(2); DOI:10.31579/2690-8808/246

Copyright: ©, 2025, Hacer Gümüş. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Received: 13 February 2025 | Accepted: 21 February 2025 | Published: 28 February 2025

Keywords: DFT calculations; morphine; ibuprofen; aspirin

Abstract

This study investigates the thermodynamic properties of commonly used pain-relieving molecules such as morphine, ibuprofen, and aspirin. The molecular structures were determined through geometric optimization calculations performed in the gas phase using Gaussian software, ensuring accurate determination of energy minima and structural configurations. Following this, thermodynamic properties were evaluated using Density Functional Theory (DFT) with the B3LYP/6-311G(d,p) basis set, analyzing parameters such as binding energies, free energy changes, and entropy variations. The study aims to enhance the understanding of how these molecules interact with biological targets and their potential effectiveness in pain relief.

1.Introduction

Pain is a multifaceted symptom associated with various medical conditions, and its effective management remains a central focus in both clinical and research domains. In this study, we have selected three widely used pain-relieving molecules morphine, ibuprofen, and aspirin due to their distinct mechanisms of action and varied pharmacological profiles. By comparing the thermodynamic properties of these molecules, this research aims to provide valuable insights into their analgesic effects and potential for clinical application.

Morphine, an opioid, is renowned for its potent analgesic properties and is commonly prescribed for the management of severe pain [1-4]. It exerts its effects by binding to opioid receptors in the brain and spinal cord, altering the perception of pain. Ibuprofen, a non-steroidal anti-inflammatory drug (NSAID), is effective in relieving both pain and inflammation [5-8]. It achieves this by inhibiting cyclooxygenase (COX) enzymes, reducing the production of prostaglandins, and consequently alleviating pain and inflammation. Aspirin [9], another widely used NSAID, also works by inhibiting COX enzymes to block prostaglandin production, helping to manage pain, inflammation, and fever [10-14].

This study seeks to offer a comprehensive analysis of the thermodynamic properties of morphine, ibuprofen, and aspirin. By determining their geometric structures and calculating key thermodynamic parameters, we aim to better understand their binding affinities to biological targets and their overall analgesic effectiveness. Geometric optimizations in the gas phase will be carried out using Gaussian 5 software, and thermodynamic calculations will be performed through Density Functional Theory (DFT) with the B3LYP/6-311G(d,p) basis set. The data obtained will shed light on the role of these molecules in pain management and highlight potential areas for improvement in future drug development processes.

2. Computational Details

The geometric structures of the pain-relieving molecules morphine, ibuprofen, and aspirin were initially visualized using the Gaussian View 5 program [15]. Following this, a series of comprehensive theoretical calculations were conducted to investigate various aspects of these compounds. Utilizing the Density Functional Theory (DFT) [16] method, calculations were specifically performed using the B3LYP [17, 18] functional and the 6-311G (d,p) [19] basis set. This computational approach is particularly effective for predicting molecular structures, energy levels, and diverse chemical properties, including those related to interactions with biological targets. The theoretical calculations aimed to elucidate not only the geometric configurations of the pain-relieving molecules but also to provide insights into their thermodynamic properties such as binding energies, free energy changes, and entropy variations. By employing advanced computational tools, this study enhances the understanding of the structural characteristics and potential pharmacological activities of these molecules, thereby contributing valuable information for pharmaceutical research and drug development.

3. Results and Discussion

3.1. Geometric Structure

The previously synthesized structures of morphine, ibuprofen, and aspirin were obtained from the Cambridge Crystallographic Data Centre (CCDC). A comparative theoretical analysis of these three molecules is not present in the literature and has been conducted for the first time in this study using the DFT/B3LYP/6-311G (d, p) method.

X-ray single crystal structures of each molecule provided a foundation for theoretical geometric calculations. Density Functional Theory (DFT) calculations were performed using GAUSSIAN 09 software, with the B3LYP methodology and the 6-311G (d, p) basis set to optimize the lowest energy conformations of each molecule. The geometric structures of morphine, ibuprofen, and aspirin were thoroughly examined and optimized using the DFT/B3LYP/6-311G (d, p) approach. These calculations allowed for a detailed analysis of the molecular geometries, energies, and stabilities.

The theoretical analyses have provided detailed insights into the electronic density, molecular orbital distributions, and steric interactions of the molecules. Additionally, these analyses assessed potential variations in the geometric structures and reactivity profiles of the molecules. Such theoretical information is crucial for testing the alignment with experimental data and for a better understanding of the physicochemical properties of these molecules.

The experimental structure of morphine is shown in Figure 1a, and its optimized structure obtained using DFT/B3LYP/6-311G (d, p) is depicted in Figure 1b.

Figure 1a: The Experimantal structure 1b. The optimized geometry of morphine.

The experimental structure of ibuprofen is shown in Figure 2a, and its optimized structure obtained using DFT/B3LYP/6-311G(d,p) is depicted in Figure 2b.

Figure 2a: The Experimantal structure 2b. The optimized geometry of ibuprofen.

The experimental structure of aspirin is shown in Figure 3a, and its optimized structure obtained using DFT/B3LYP/3-311G(d,p) is depicted in Figure 3b.

Figure 3a: The Experimantal structure 3b. The optimized geometry of aspirin.

DFT calculations have provided detailed insights into the atomic and electronic distributions within each molecule, facilitating a thorough understanding of chemical bonding and molecular interactions. The structural features of morphine, ibuprofen, and aspirin, including electronic density and molecular orbital distribution, have been analyzed in depth. Furthermore, the calculations have revealed steric interactions and bond lengths, as well as potential reactivity profiles of these molecules.

These analyses enhance our understanding of molecular structures and ensure consistency between experimental observations and theoretical predictions. The optimized structures of the molecules provide valuable information for various scientific applications and test the accuracy of theoretical modeling.

3.2. Calculation and Comparison of Thermodynamic Properties

In this study, the thermodynamic properties of morphine, ibuprofen, and aspirin molecules were calculated using the Gaussian program with the density functional theory (DFT) method, employing the B3LYP/6-311G(d,p) functional. The B3LYP functional provides a precise modeling of electronic structures, while the 6-311G(d,p) basis set offers an expanded basis set for detailed energy calculations. This methodology enables accurate determination of molecular energy levels and thermodynamic properties.

The obtained data from the calculations include various thermodynamic parameters for morphine, ibuprofen, and aspirin molecules. These parameters are summarized and compared in Table 1. Table 1 details each molecule’s thermal energy, heat capacity, entropy, rotational constants, and thermal corrections, highlighting differences among the molecules.

Termal energy, E (Kcal/mol)MorphineIbuprofenAspirin
          Rotational0.8890.8890.889
          Translational0.8890.8890.889
          Vibrational236.106186.069103.841
          Total237.884187.847105.618
Heat capacity, Cv (cal/mol K)   
          Rotational2.9812.9812.981
          Translational2.9812.9812.981
          Vibrational73.29652.93837.107
          Total79.25758.90043.069
Entropy, S (cal/mol K)   
          Rotational33.80632.55631.004
          Translational43.02441.87441.471
          Vibrational60.12955.08037.037
          Total136.959129.510109.512
Rotational constants (GHz)   
           A0.479671.573571.13519
           B0.255660.240440.74845
           C0.206600.235450.50026
Rotational temperature (Kelvin)   
           A0.023020.075520.05448
           B0.012270.011540.03592
           C0.009910.011300.02401
Thermal properties (Hartree/particle)   
Zero-point correction 0.3599950.2835970.156551
Thermal correction to Energy 0.3790920.2993530.168313
Thermal correction to Enthalpy0.3800360.3002970.169257
Thermal correction to Gibbs Free Energy0.3149620.2387620.117225
Sum of electronic and zero-point Energies-1015.973479-656.600463-648.715954
Sum of electronic and thermal Energies-1015.954383-656.584707-648.704192
Sum of electronic and thermal Enthalpies-1015.953439-656.583763-648.703248
Sum of electronic and thermal Free Energies-1016.018512-656.645297-648.755281
Zero point vibrational energy (kcal/mol)225.90054177.9597598.23744

Table 1: The calculated thermodynamic parameters of molecules.

Table 1 provides a detailed comparison of the thermodynamic properties of morphine, ibuprofen, and aspirin molecules. According to thermal energy data, rotational and translational energies are identical for all three molecules, each being 0.889 kcal/mol. However, significant differences are observed in vibrational energies: morphine has 236.106 kcal/mol, ibuprofen has 186.069 kcal/mol, and aspirin has 103.841 kcal/mol. Consequently, total energy values are calculated as 237.884 kcal/mol for morphine, 187.847 kcal/mol for ibuprofen, and 105.618 kcal/mol for aspirin.

Regarding heat capacity, rotational and translational heat capacities are consistent across all molecules, each being 2.981 cal/mol K. Vibrational heat capacities show variation with morphine at 73.296 cal/mol K, ibuprofen at 52.938 cal/mol K, and aspirin at 37.107 cal/mol K. Total heat capacities are 79.257 cal/mol K for morphine, 58.900 cal/mol K for ibuprofen, and 43.069 cal/mol K for aspirin.

Entropy values indicate that the rotational entropy is 33.806 cal/mol K for morphine, 32.556 cal/mol K for ibuprofen, and 31.004 cal/mol K for aspirin; translational entropy is 43.024 cal/mol K for morphine, 41.874 cal/mol K for ibuprofen, and 41.471 cal/mol K for aspirin; vibrational entropy is 60.129 cal/mol K for morphine, 55.080 cal/mol K for ibuprofen, and 37.037 cal/mol K for aspirin. Total entropy values are 136.959 cal/mol K for morphine, 129.510 cal/mol K for ibuprofen, and 109.512 cal/mol K for aspirin.

Rotational constants and temperatures are as follows: for morphine, A = 0.47967 GHz, B = 0.25566 GHz, C = 0.20660 GHz, and rotational temperatures are A = 0.02302 K, B = 0.01227 K, C = 0.00991 K; for ibuprofen, A = 1.57357 GHz, B = 0.24044 GHz, C = 0.23545 GHz, and rotational temperatures are A = 0.07552 K, B = 0.01154 K, C = 0.01130 K; for aspirin, A = 1.13519 GHz, B = 0.74845 GHz, C = 0.50026 GHz, and rotational temperatures are A = 0.05448 K, B = 0.03592 K, C = 0.02401 K.

Thermal corrections indicate that morphine has a zero-point correction of 0.359995 Hartree/particle, a thermal correction to energy of 0.379092 Hartree/particle, a thermal correction to enthalpy of 0.380036 Hartree/particle, and a thermal correction to Gibbs free energy of 0.314962 Hartree/particle; ibuprofen has a zero-point correction of 0.283597 Hartree/particle, a thermal correction to energy of 0.299353 Hartree/particle, a thermal correction to enthalpy of 0.300297 Hartree/particle, and a thermal correction to Gibbs free energy of 0.238762 Hartree/particle; aspirin has a zero-point correction of 0.156551 Hartree/particle, a thermal correction to energy of 0.168313 Hartree/particle, a thermal correction to enthalpy of 0.169257 Hartree/particle, and a thermal correction to Gibbs free energy of 0.117225 Hartree/particle.

These results provide a comprehensive comparison of the thermodynamic properties of morphine, ibuprofen, and aspirin molecules and offer a deep understanding of their thermal behaviors. Such detailed analyses are crucial for drug formulation and stability assessments, aiding in the optimization of molecular stability and performance.

4. Conclusions

In this study, the thermodynamic properties of pain-relieving molecules, including morphine, ibuprofen, and aspirin, have been thoroughly examined. Geometric optimizations performed in the gas phase using Gaussian software allowed for the accurate determination of energy minima and molecular structures. Following this, the thermodynamic properties of these molecules were calculated through Density Functional Theory (DFT) with the B3LYP/6-311G(d,p) basis set. The analysis of key parameters, including binding energies, free energy changes, and entropy variations, provided 

important insights into the molecules' interactions with biological targets and their analgesic mechanisms.

The strong analgesic effect of morphine can be attributed to its high binding energy, which reflects its potent interaction with opioid receptors in the brain and spinal cord. In contrast, the efficacy of non-steroidal anti-inflammatory drugs (NSAIDs) such as ibuprofen and aspirin is primarily linked to their ability to bind to cyclooxygenase (COX) enzymes, which play a critical role in the synthesis of prostaglandins involved in pain and inflammation. These findings underscore the significance of the molecular structure and thermodynamic properties in determining the effectiveness of these molecules in pain relief, as well as their potential side effects.

Furthermore, this study highlights essential parameters for optimizing the molecular design of analgesic agents. By understanding the binding affinities and thermodynamic behaviors of these drugs, new avenues for enhancing drug efficacy and minimizing adverse effects in future pharmaceutical formulations can be explored. The insights gained from this research will be crucial for the development of more effective pain management therapies. Future studies could expand on these results by conducting more detailed investigations of these molecules and exploring novel drug formulations to further improve therapeutic outcomes.

References

a