Molecular dynamics (MD) simulations and biased MD simulation were carried out for the neutral form of Paracetamol inserted in fully hydrated dipalmitoylphosphatidylcholine (DPPC) and dimyristoylphosphatidylcholine (DMPC) lipid bilayers. For comparison, fully hydrated DMPC and DPPC lipid bilayers were also simulated separately without Paracetamol. The simulation time for each system was 50 ns. At two concentrations of Paracetamol, various properties of the lipid bilayer such as area per lipid, order parameter, diffusion coefficient, radial distribution function, electrostatic potential, mass density and hydrogen bonds have been calculated. Also, the convergence in time of the free energy profile of the Paracetamol along a DPPC bilayer normal was calculated by umbrella sampling method. From the obtained results, it can be concluded that neutral form of Paracetamol shows a generally similar behaviour in DPPC and DMPC lipid bilayers. It was shown that the addition of Paracetamol causes a decrease in tail order parameter of both DPPC and DMPC lipid bilayers and the tail of Paracetamol adopts an inward orientation in the lipid bilayers. Also from the free energy profile, the high penetration barrier in the bilayer centre was determined.

Naproxen and relafen, as nonsteroidal antiinflammatory drugs, were simulated in neutral and charged forms and their effects on a lipid bilayer membrane were investigated by molecular dynamics simulation using Groningen machine for chemical simulations software (GROMACS). Simulation of 10 systems was performed, which included different dosages of the drug molecules, naproxen and Relafen, in charged and neutral forms, and a mixture of naproxen and Relafen in neutral forms. The effects of the mixture and the individual drugs' dosages on membrane properties, such as electrostatic potential, order parameter, diffusion coefficients, and hydrogen bond formation, were analyzed. Hydration of the drugs in the membrane system was investigated using radial distribution function analysis. Using fully hydrated dimyristoylphosphatidylcholine (DMPC) as a reference system, 128 lipid molecules and water molecules were simulated exclusively, and the same simulation technique was performed on 10 other systems, including drug mixtures and a DMPC membrane. Angular distributions of lipid chains of the membrane were calculated, and the effects of the drug insertion and chain orientation in the membrane were evaluated.

No abstract

The interaction of PEGylated anti-hypertensive drugs, amlodipine, atenolol and lisinopril with lipid bilayer membrane dimyristoylphosphatidylcholine (DMPC) has been studied in nine different simulation systems consisting of 128 lipid molecules and appropriate number of water molecules by molecular dynamics method and by utilizing GROMACS software. The influences of PEGylation on the mentioned drugs and the differences in application of two types of spacer molecules on the performance of drugs and DMPC membrane have been evaluated and mass density of the components in the simulation box, mean square displacement (MSD), electrostatic potential, hydrogen bonding, radial distribution function (RDF), area per lipid, order parameter, and angle distribution of the component molecules including drug, DMPC and PEG has been investigated. Furthermore, umbrella sampling analysis indicated that, PEGylation of the drugs made amlodipine to behave more hydrophilic, whereas in case of lisinopril and atenolol, PEGylation made these drugs to behave more hydrophobic. In almost all of the simulated systems, PEGylation increased the diffusion coefficient of the drugs.

The anti-hypertensive drugs amlodipine, atenolol and lisinopril, in ordinary and PEGylated forms, with different combined-ratios, were studied by molecular dynamics simulations using GROMACS software. Twenty simulation systems were designed to evaluate the interactions of drug mixtures with a dimyristoylphosphatidylcholine (DMPC) lipid bilayer membrane, in the presence of water molecules. In the course of simulations, various properties of the systems were investigated, including drug location, diffusion and mass distribution in the membrane; drug orientation; the lipid chain disorder as a result of drug penetration into the DMPC membrane; the number of hydrogen bonds; and drug surface area. According to the results obtained, combined drugs penetrate deeper into the DMPC lipid bilayer membrane, and the lipid chains remain ordered. Also, the combined PEGylated drugs, at a combination ratio of 1:1:1, enhance drug penetration into the DMPC membrane, reduce drug agglomeration, orient the drug in a proper angle for easy penetration into the membrane, and decrease undesirable lipotoxicity due to distorted membrane self-assembly and thickness. Graphical abstract ᅟ.

Articaine, as a local anesthetic drug has been simulated in neutral and charged forms, and its interaction with the dimyristoylphosphatidylcholine (DMPC) lipid bilayer membrane is investigated by molecular dynamics simulation using GROMACS software. In order to obtain the optimum location of the drug molecules, as they penetrate into the membrane, umbrella sampling is applied and the free energy is calculated. The effect of protein binding to DMPC membrane on the process of drug diffusion through the membrane is considered. Five simulation systems are designed and by applying the potential of mean force, the molecular dynamics simulation on the system is performed. In light of the obtained results, the electrostatic potential, variation of lipid bilayer's order parameter and the diffusion coefficient of drug are discussed.

The interactions of the drugs amlodipine and paroxetine, which are prescribed respectively for treatment of hypertension and depression, with the metabolizing enzyme cytochrome CYP2B4 as the drug target, have been studied by molecular dynamics (MD) simulation. Poly ethylene glycol was used to control the drugs' interactions with each other and with the target CYP2B4. Thirteen simulation systems were carefully designed, and the results obtained from MD simulations indicated that amlodipine in the PEGylated form prescribed with paroxetine in the nonPEGylated form promotes higher cytochrome stability and causes fewer fluctuations as the drugs approach the target CYP2B4 and interact with it. The simulation results led us to hypothesize that the combination of the drugs with a specific drug ratio, as proposed in this work, manifests more effective diffusivity and less instability while metabolizing with enzyme CYP2B4. Also, the active residues in the CYP2B4 enzyme that interact with the drugs were determined by MD simulation, which were consistent with the reported experimental results. Graphical Abstract Efficient drug-enzyme interactions, as a result of PEGylation.

In this work, the effects of the anti-hypertensive drug amlodipine in native and PEGylated forms on the malfunctioning of negatively charged lipid bilayer cell membranes constructed from DMPS or DMPS + DMPC were studied by molecular dynamics simulation. The obtained results indicate that amlodipine alone aggregates and as a result its diffusion into the membrane is retarded. In addition, due to their large size aggregates of the drug can damage the cell, rupturing the cell membrane. It is shown that PEGylation of amlodipine prevents this aggregation and facilitates its diffusion into the lipid membrane. The interaction of the drug with negatively charged membranes in the presence of an aqueous solution of NaCl, as the medium, is investigated and its effects on the membrane are considered by evaluating the structural properties of the membrane such as area per lipid, thickness, lipid chain order and electrostatic potential difference between bulk solution and lipid bilayer surface. The effect of these parameters on the diffusion of the drug into the cell is critically examined and discussed.