Diffusional interaction behavior of NSAIDs in lipid bilayer membrane using molecular dynamics (MD) simulation: Aspirin and Ibuprofen

Insight into the effect of ibuprofen on the permeability of the membrane: a molecular dynamic simulation study

Studying interactions between drugs and cell membranes is of great interest to designing novel drugs, optimizing drug delivery, and discerning drug mechanism action. In this study, we investigated the physical properties of the bilayer membrane model of POPC upon interaction with ibuprofen (IBU) using molecular dynamics simulations. The area per lipid (APL) was calculated to describe the effect of ibuprofen on the packing properties of the lipid bilayer. The APL was 0.58 nm2 and 0.63 nm2 for the membrane in low and high IBU respectively, and 0.57 nm2 for the membrane without IBU. Our finding showed that the mean square deviation (MSD) increased with increased ibuprofen content. In addition, the order parameter for the hydrocarbon chain of lipids increased with increased ibuprofen content. There was an increment in the transfer free energy after the head group region while it was maximum in the hydrophobic core for hydrogen peroxide (H2O2) (∼6.2 kcal.mol-1) and H2O (∼3.4 kcal.mol-1) which then decreased to respective values of (∼4.6 kcal.mol-1), and (∼2.3 kcal.mol-1) at the center of the bilayer in the presence of IBU. It seems that in the presence of ibuprofen, the free energy profile of the permeability of water and H2O2 significantly decreased. These findings show that ibuprofen significantly influences the physical properties of the bilayer by decreasing the packing and intermolecular interaction in the hydrocarbon chain region and increasing the water permeability of the bilayer. These results may provide insights into the local cytotoxic side effects of ibuprofen and its underlying molecular mechanisms.Communicated by Ramaswamy H. Sarma.

An in silico study of protein-protein interactions and design of novel peptides for TrkA in ameloblastoma

Ameloblastoma is a benign odontogenic jawbone tumor. The binding of Nerve growth factor (NGF) to receptor tyrosine kinase A (TrkA) promotes cell survival, proliferation, and differentiation via PI3K/AKT and Ras/MAPK signaling. Although the exact cause of ameloblastoma remains unknown, elevated levels of NGF and TrkA expression in ameloblastoma are associated with aggressive tumor behavior and poor patient outcomes. It is previously demonstrated that His 4, Arg 9, and Glu 11 residues of NGF made crucial interactions with the TrkA subunit. The main aim of our present study to develop potential therapeutic strategies by identifying promising peptide candidates. The objectives include starting with a detailed in silico analysis to identify a crucial peptide sequence of NGF that is bound by TrkA, creating a library of novel peptides from the identified peptide sequence through a single-point mutation on interacting residues (His 4, Arg 9, and Glu 11), and selecting the top peptides based on docking score, interactions analysis, and desirable pose analysis. The study ultimately designed a hybrid peptide candidate through the simultaneous and continuous mutation of the top residues, resulting in a peptide that exhibited a more specific interaction with TrkA, blocking the binding site and preventing the interaction between NGF and TrkA.Communicated by Ramaswamy H. Sarma.

Molecular dynamics simulation of phosphatidylcholine membrane in low ionic strengths of sodium chloride

The one-microsecond molecular dynamics simulations of a membrane-protein complex investigate the influence of the aqueous sodium chloride solutions on the structure and dynamics of a palmitoyl-oleoyl-phosphatidylcholine bilayer membrane. The simulations were performed on five different concentrations (40, 150, 200, 300, and 400 mM) in addition to a salt-free system by using the charmm36 force field for all atoms. Four biophysical parameters, (membrane thicknesses of annular and bulk lipids, and the area per lipid of both leaflets), were computed separately. Nevertheless, the area per lipid was expressed by using the Voronoi algorithm. All time-independent analyses were carried out for the last 400 ns trajectories. Different concentrations revealed dissimilar membrane dynamics before equilibration. The biophysical properties of the membrane (thickness, area-per-lipid, and order parameter) have non-significant changes with increasing ionic strength, however, the 150 mM system had exceptional behavior. Sodium cations were dynamically penetrating the membrane forming weak coordinate bonds with single or multiple lipids. Nevertheless, the binding constant was unaffected by the cation concentration. The electrostatic and Van der Waals energies of lipid-lipid interactions were influenced by the ionic strength. On the other hand, the Fast Fourier Transform was performed to figure out the dynamics at the membrane-protein interface. The nonbonding energies of membrane-protein interactions and order parameters explained the differences in the synchronization pattern. All results were consensus with experimental and theoretical works.Communicated by Ramaswamy H. Sarma.

Interactions between β-cyclodextrin as a carrier for anti-cancer drug delivery: a molecular dynamics simulation study

A series of molecular dynamics simulations were performed on 5-fluorouracil (5-Fu), Alendronate (Ald), and Temozolomide (TMZ) anticancer drugs in the presence and absence of β-cyclodextrin (βCD) as a carrier. Thermodynamic investigations showed that the van der Waals interaction energy was dominant in loading all drugs inside the βCD cavity. The sum of the interaction energies illustrated that the highest affinity was related to Ald (-136.5 kJ/mol), which in turn was due to the presence of bulky and charged atoms of phosphorus and oxygen, although TMZ (-115.92 kJ/mol) showed a very high affinity as well. At the same time, the hydrogen bond analysis also represented that Ald had the most hydrogen bond (1.97) with the highest half-life (3.13 ps) with βCD. Investigation of the root mean fluctuation (RMSF) indicated that all the drugs had a relatively rigid structure and maintain this rigidity during loading in the βCD cavity, and in the meantime, Ald was slightly more flexible than 5-Fu and TMZ. The area of ​the primary hydroxyl rim decreased in all drug-containing systems, which in turn was caused by the attractive interaction of drugs with oxygens in the primary hydroxyl rim. Especially for those drugs that were able to penetrate to the end of the primary hydroxyl rim of the βCD, that means TMZ and 5-Fu. Meanwhile, due to the lack of Ald penetration to the end of the primary hydroxyl rim, the area change in the Ald-containing system was less than in the two others.Communicated by Ramaswamy H. Sarma.

Molecular Dynamics Study on Structure, Vibrational Properties, and Transport Coefficients of Liquid Alumina

The structure, vibrational density of states, and transport coefficients of liquid alumina were studied using molecular dynamics simulations. At the temperature of 2500 K, 3000 K, 3500 K, and 4000 K, systems with three different densities were constructed, respectively, including the configurations with densities of 2.81 g/cm³ and 3.17 g/cm³, and the relaxed ones with nearly zero pressure at each temperature. With the changes in temperature or density, the transformations on the structural, vibrational and transport properties were discussed. The Born-Mayer-Huggins type of atomic interactions was used, with newly optimized parameters. The analysis of the interatomic correlations indicated that the short-range order of liquid alumina was mainly constructed by AlO₄ tetrahedra, also a certain number of AlO₃ and AlO₅ was present. Meanwhile, the structural transitions on the elemental units occurred as either the temperature or density increased. Two primary frequency bands were observed in each vibrational density of states spectrum, with the higher frequency bands produced by the O atom vibrations, and the lower frequency ones generated by the Al atom vibrations. Self-diffusion coefficients were estimated using the linear behavior of the mean-squared displacement at long time, while by using the Green-Kubo relation during equilibrium molecular dynamics simulations, thermal conductivities and viscosities were calculated. Significantly, the viscosity at 2500 K with a density of 2.81 g/cm³ was equal to 25.23 mPa s, which was very close to the experimental finding.

Molecular Insights into Factors Affecting the Generation Behaviors, Dynamic Properties, and Interfacial Structures of Methane Gas Bubbles

Molecular dynamics simulations were performed to study the effects of temperatures, pressures, and methane mole fractions on the generation behaviors, dynamic properties, and interfacial structures of methane gas bubbles. Methane gas bubbling can be promoted by high temperatures and high mole fractions of methane, which come from the generation of larger methane clusters in solution. Bubbles were found to be highly dynamic, with more methane molecules exchanging between bubbles and the surrounding solution at high pressures and in systems with high mole fractions of methane. The interfacial structures between bubbles and the surrounding solution were rough at a molecular level, and the roughness of the outermost methane and water molecules was high at high temperatures, low pressures, and in systems with high methane mole fractions. The dissolution of methane molecules depended on the interactions between the outermost methane and water molecules, which would become stronger with decreasing temperatures, increasing pressures, and decreasing methane mole fractions. The results obtained can help in understanding both the generation behaviors of bubbles when gas hydrates decompose and the re-nucleation behaviors of gas hydrates in the presence of bubbles.

Research progress on the potential novel analgesic BU08028

Pain is a common symptom accompanying several clinical conditions and causes serious distress to patients. Addressing pain management is an important aspect of disease treatment, including cancer therapy. Opioid analgesics used to manage pain in human and veterinary medicine have been associated with substance dependence and other adverse effects, thereby limiting their application. Thus, the development of opioid analgesics with good safety profiles with minimal adverse effects and no addictive effects, is presently the focus of pain research. As a new potential analgesic, (2S)-2-[(5R,6R,7R,14S)-N-cyclopropylmethyl-4,5-epoxy-6,14-ethano-3-hydroxy-6-methoxymorphinan-7-yl]-3,3-dimethylpentan-2-ol (BU08028) has fewer adverse effects than other analgesics and is expected to be a safer alternative. In this review, we discuss the development of the opioid analog BU08028 and summarize its analgesic effects and biological characteristics, including efficiency, safety, and tolerance. Furthermore, we elaborate on studies showing the bifunctional effect of BU08028, which targets both mu opioid peptide and nociceptin-orphanin FQ peptide receptors, as well as the unique advantages of using BU08028 over single-target opioid agonists. Previous studies have suggested that BU08028 can not only weaken the reward and abuse effects of opioids and other drugs, but also enhance the anti-nociceptive effect of the mu opioid peptide receptors, making it a potent analgesic. Besides, we describe studies suggesting that BU08028 inhibits the effects of alcohol, making it a candidate drug for the management of alcohol addiction. Our review suggests that BU08028 is a potential novel medicine for managing pain and addiction.

Mechanistic Understanding from Molecular Dynamics in Pharmaceutical Research 2: Lipid Membrane in Drug Design

We review the use of molecular dynamics (MD) simulation as a drug design tool in the context of the role that the lipid membrane can play in drug action, i.e., the interaction between candidate drug molecules and lipid membranes. In the standard "lock and key" paradigm, only the interaction between the drug and a specific active site of a specific protein is considered; the environment in which the drug acts is, from a biophysical perspective, far more complex than this. The possible mechanisms though which a drug can be designed to tinker with physiological processes are significantly broader than merely fitting to a single active site of a single protein. In this paper, we focus on the role of the lipid membrane, arguably the most important element outside the proteins themselves, as a case study. We discuss work that has been carried out, using MD simulation, concerning the transfection of drugs through membranes that act as biological barriers in the path of the drugs, the behavior of drug molecules within membranes, how their collective behavior can affect the structure and properties of the membrane and, finally, the role lipid membranes, to which the vast majority of drug target proteins are associated, can play in mediating the interaction between drug and target protein. This review paper is the second in a two-part series covering MD simulation as a tool in pharmaceutical research; both are designed as pedagogical review papers aimed at both pharmaceutical scientists interested in exploring how the tool of MD simulation can be applied to their research and computational scientists interested in exploring the possibility of a pharmaceutical context for their research.

Thermodynamics and Intermolecular Interactions during the Insertion of Anionic Naproxen into Model Cell Membranes

The insertion process of Naproxen into model dimyristoylphosphatidylcholine (DMPC) membranes is studied by resorting to state-of-the-art classical and quantum mechanical atomistic computational approaches. Molecular dynamics simulations indicate that anionic Naproxen finds an equilibrium position right at the polar/nonpolar interphase when the process takes place in aqueous environments. With respect to the reference aqueous phase, the insertion process faces a small energy barrier of ≈5 kJ mol^-1 and yields a net stabilization of also ≈5 kJ mol^-1. Entropy changes along the insertion path, mainly due to a growing number of realizable microstates because of structural reorganization, are the main factors driving the insertion. An attractive fluxional wall of noncovalent interactions is characterized by all-quantum descriptors of chemical bonding (natural bond orbitals, quantum theory of atoms in molecules, noncovalent interaction, density differences, and natural charges). This attractive wall originates in the accumulation of tiny transfers of electron densities to the interstitial region between the fragments from a multitude of individual intermolecular contacts stabilizing the tertiary drug/water/membrane system.

Electric-field-induced electroporation and permeation of reactive oxygen species across a skin membrane

Electroporation processes affect the permeability of cell membranes, which can be utilized for the delivery of plasma species in cancer therapy. By means of computational dynamics, many aspects of membrane electroporation have been unveiled at the atomic level for lipid membranes. Herein, a molecular dynamics simulation study was performed on native and oxidized membrane systems with transversal electric fields. The simulation result shows that the applied electric field mainly affects the membrane properties so that electroporation takes place and these pores are lined by hydrophilic headgroups of the lipid components. The calculated hydrophobic thickness, lateral diffusion and pair correlation revealed the role of 5α-CH in creation of water-pore in an oxidized membrane. Additionally, the permeability of reactive oxygen species was examined through these electroporated systems. The permeability study suggested that water pores in the membrane facilitate the penetration of these species across the membrane to the interior of the cell. These findings may have significance in experimental applications in vivo as once the reactive oxygen species reaches the interior of the cell, they may cause oxidative stress and induce apoptosis.Communicated by Ramaswamy H. Sarma.