Molecular dynamics study of the mechanical properties of drug loaded model systems: A comparison of a polymersome with a bilayer.

In this work we implement a new methodology to study structural and mechanical properties of systems having spherical and planar symmetries throughout Molecular Dynamics simulations. This methodology is applied here to a drug delivery system based in polymersomes, as an example. The chosen model drug was the local anesthetic prilocaine due to previous parameterization within the used coarse grain scheme. In our approach, mass density profiles (MDPs) are used to obtain key structural parameters of the systems, and pressure profiles are used to estimate the curvature elastic parameters. The calculation of pressure profiles and radial MPDs required the development of specific methods, which were implemented in an in-house built version of the GROMACS 2018 code. The methodology presented in this work is applied to characterize poly(ethylene oxide)-poly(butadiene) polymersomes and bilayers loaded with the model drug prilocaine. Our results show that structural properties of the polymersome membrane could be obtained from bilayer simulations, with significantly lower computational cost compared to whole polymersome simulations, but the bilayer simulations are insufficient to get insights on their mechanical aspects, since the elastic parameters are canceled out for the complete bilayer (as consequence of the symmetry). The simulations of entire polymersomes, although more complex, offer a complementary approach to get insights on the mechanical behavior of the systems.

Mechanical properties of drug loaded diblock copolymer bilayers: A molecular dynamics study.

In this work, we present results of coarse-grained simulations to study the encapsulation of prilocaine (PLC), both neutral and protonated, on copolymer bilayers through molecular dynamics simulations. Using a previously validated membrane model, we have simulated loaded bilayers at different drug concentrations and at low (protonated PLC) and high (neutral PLC) pH levels. We have characterized key structural parameters of the loaded bilayers in order to understand the effects of encapsulation of PLC on the bilayer structure and mechanical properties. Neutral PLC was encapsulated in the hydrophobic region leading to a thickness increase, while the protonated species partitioned between the water phase and the poly(ethylene oxide)-poly(butadiene) (PBD) interface, relaxing the PBD region and leading to a decrease in the thickness. The tangential pressures of the studied systems were calculated, and their components were decomposed in order to gain insights on their compensation. In all cases, it is observed that the loading of the membrane does not significantly decrease the stability of the bilayer, indicating that the system could be used for drug delivery.

Diblock copolymer bilayers as model for polymersomes: A coarse grain approach.

This paper presents a new model for polymersomes developed using a poly(ethylene oxide)-poly(butadiene) diblock copolymer bilayer. The model is based on a coarse-grained approach using the MARTINI force field. Since no MARTINI parameters exist for poly(butadiene), we have refined these parameters using quantum mechanical calculations and molecular dynamics simulations. The model has been validated using extensive molecular dynamics simulations in systems with several hundred polymer units and reaching up to 6 µs. These simulations show that the copolymer coarse grain model self-assemble into bilayers and that NPT and NPNγT ensemble runs reproduce key structural and mechanical experimental properties for different copolymer length chains with a similar hydrophilic weight fraction.

Study of polymorphism in imatinib mesylate: a quantum chemical approach using electronic and vibrational spectra.

Imatinib mesylate, 4-(4-methyl-piperazin-1-ylmethyl)-N-u[4-methyl-3-(4-pyridin-3-yl)pyrimidine-2-ylamino)phenyl]benzamide methanesulfonate is a therapeutic drug that is approved for the treatment of chronic myelogeneous leukemia (CML) and gastrointestinal stromal tumors (GIST). It is known that imatinib mesylate exists in two polymorphic forms α and ß. However, ß-form is more stable than the α-form. In this work, we present a detailed vibrational spectroscopic investigation of ß-form by using FT-IR and FT-Raman spectra. These data are supported by quantum mechanical calculations using DFT employing 6-311G(d,p) basis set, which allow us to characterize completely the vibrational spectra of this compound. The FT-IR spectrum of α-form has also been discussed. The importance of hydrogen-bond formation in the molecular packing arrangements of both forms has been examined with the vibrational shifts observed due to polymorphic changes. The red shift of the NH stretching bands in the infrared spectrum from the computed wavenumber indicates the weakening of the NH bond. The UV-vis spectroscopic studies along with the HOMO-LUMO analysis of both polymorphs (α and ß) were performed and their chemical activity has been discussed. The TD-DFT method was used to calculate the electronic absorption spectra in the gas phase as well as in the solvent environment using IEF-PCM model and 6-31G basis set. Finally, the results obtained complements to the experimental findings.

Conformational polymorphism on imatinib mesylate: grinding effects.

Crystal structures of polymorphs α and ß of imatinib mesylate were obtained. Thermal behavior and grinding effects were studied by X-ray powder diffraction and differential scanning calorimetry techniques. Molecules in forms α and ß exhibit significant conformational differences due to dissimilar intramolecular interactions, which stabilize their molecular conformations. In spite of that, both crystal structures present a dimer-chain arrangement. Dimers are mainly determined by hydrogen bonding interactions and some weak π-π interactions. Connections between dimers are provided by mesylate ions to determine chains of dimers. Neighboring chains are linked by very weak interactions: C-H···π interactions in form α and π-π interactions in form ß. At room temperature, thermal disorder was observed in the mesylate ion in form α, which could be removed at low temperatures (-123°C). Form ß was found to be the more stable form at room temperature. Both polymorphs exhibit a tendency to generate amorphous material by grinding, which can be converted to a crystalline phase by either temperature or aging. When amorphous crystallization is kinetically studied at room temperature, form ß is obtained after a week. Conversely, when the crystallization is activated by temperature, the final obtained crystal form depends on the starting material, proving the importance of seeding.

Towards crystal structure prediction of complex organic compounds--a report on the fifth blind test.

Following on from the success of the previous crystal structure prediction blind tests (CSP1999, CSP2001, CSP2004 and CSP2007), a fifth such collaborative project (CSP2010) was organized at the Cambridge Crystallographic Data Centre. A range of methodologies was used by the participating groups in order to evaluate the ability of the current computational methods to predict the crystal structures of the six organic molecules chosen as targets for this blind test. The first four targets, two rigid molecules, one semi-flexible molecule and a 1:1 salt, matched the criteria for the targets from CSP2007, while the last two targets belonged to two new challenging categories - a larger, much more flexible molecule and a hydrate with more than one polymorph. Each group submitted three predictions for each target it attempted. There was at least one successful prediction for each target, and two groups were able to successfully predict the structure of the large flexible molecule as their first place submission. The results show that while not as many groups successfully predicted the structures of the three smallest molecules as in CSP2007, there is now evidence that methodologies such as dispersion-corrected density functional theory (DFT-D) are able to reliably do so. The results also highlight the many challenges posed by more complex systems and show that there are still issues to be overcome.