Insights into Transfer of Supramolecular Doxorubicin/Congo Red Aggregates through Phospholipid Membranes.

Doxorubicin (DOX) is a commonly used chemotherapeutic drug, from the anthracycline class, which is genotoxic to neoplastic cells via a DNA intercalation mechanism. It is effective and universal; however, it also causes numerous side effects. The most serious of them are cardiotoxicity and a decrease in the number of myeloid cells. For this reason, targeted DOX delivery systems are desirable, since they would allow lowering the drug dose and therefore limiting systemic side effects. Recently, synthetic dyes, in particular Congo red (CR), have been proposed as possible DOX carriers. CR is a planar molecule, built of a central biphenyl moiety and two substituted naphthalene rings, connected with diazo bonds. In water, it forms elongated ribbon-shaped supramolecular structures, which are able to selectively interact with immune complexes. In our previous studies, we have shown that CR aggregates can intercalate DOX molecules. In this way, they preclude DOX precipitation in water solutions and increase its uptake by MCF7 breast cancer cells. In the present work, we further explore the interactions between DOX, CR, and their aggregates (CR/DOX) with phospholipid membranes. In addition to neutral molecules, the protonated doxorubicin form, DXP, is also studied. Molecular dynamics simulations are employed to study the transfer of CR, DOX, DXP, and their aggregates through POPC bilayers. Interactions of CR, DOX, and CR/DOX with model monolayers are studied with Langmuir trough measurements. This study shows that CR may support the transfer of doxorubicin molecules into the bilayer. Both electrostatic and van der Waals interactions with lipids are important in this respect. The former promote the initial stages of the insertion process, the latter keep guest molecules inside the bilayer.

Chemometrics as a valuable tool for evaluating interactions between antiretroviral drugs and food.

AIMS: Clinically significant interactions with food occur for more than half of antiretroviral drugs. Different physiochemical properties deriving from the chemical structures of antiretroviral drugs may contribute to the variable food effect. Chemometric methods allow analysing a large number of interrelated variables concomitantly and visualizing correlations between them. We used a chemometric approach to determine the types of correlations among different features of antiretroviral drugs and food that may influence interactions. METHODS: Thirty-three antiretroviral drugs were analysed: ten nucleoside reverse transcriptase inhibitors, six non-nucleoside reverse transcriptase inhibitors, five integrase strand transfer inhibitors, ten protease inhibitors, one fusion inhibitor and one HIV maturation inhibitor. Input data for the analysis were collected from already published clinical studies, chemical records and calculations. We constructed a hierarchical partial least squares (PLS) model with three response parameters: postprandial change of time to reach maximum drug concentration (ΔTmax ), albumin binding (%) and logarithm of partition coefficient (logP). Predictor parameters were the first two principal components of principal component analysis (PCA) models for six groups of molecular descriptors. RESULTS: PCA models explained 64.4% to 83.4% of the variance of the original parameters (average: 76.9%), whereas the PLS model had four significant components and explained 86.2% and 71.4% of the variance in the sets of predictor and response parameters, respectively. We observed 58 significant correlations between ΔTmax , albumin binding (%), logP and constitutional, topological, hydrogen bonding and charge-based molecular descriptors. CONCLUSIONS: Chemometrics is a useful and valuable tool for analysing interactions between antiretroviral drugs and food.

Nucleoside Analog Reverse-Transcriptase Inhibitors in Membrane Environment: Molecular Dynamics Simulations.

The behavior of four drugs from the family of nucleoside analog reverse-transcriptase inhibitors (zalcitabine, stavudine, didanosine, and apricitabine) in a membrane environment was traced using molecular dynamics simulations. The simulation models included bilayers and monolayers composed of POPC and POPG phospholipids. It was demonstrated that the drugs have a higher affinity towards POPG membranes than POPC membranes due to attractive long-range electrostatic interactions. The results obtained for monolayers were consistent with those obtained for bilayers. The drugs accumulated in the phospholipid polar headgroup region. Two adsorption modes were distinguished. They differed in the degree of penetration of the hydrophilic headgroup region. Hydrogen bonds between drug molecules and phospholipid heads were responsible for adsorption. It was shown that apricitabine penetrated the hydrophilic part of the POPC and POPG membranes more effectively than the other drugs. Van der Waals interactions between S atoms and lipids were responsible for this.

The Molecular Bases of the Interaction between a Saponin from the Roots of Gypsophila paniculata L. and Model Lipid Membranes.

In view of the possible medical applications of saponins, the molecular structure of a GOTCAB saponin from the roots of Gypsophila paniculata L. was determined by NMR. The biological activity of saponins may depend on the interaction with cell membranes. To obtain more insight in the mechanism of membrane-related saponin function, an experimental and theoretical study was conducted. Ternary lipid systems composed of sphingomyelin, 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine, and cholesterol were used as models of mammalian cell membranes. The membrane-saponin interaction was studied experimentally by monitoring surface pressure in the monomolecular films formed at the air-aqueous subphase interface. The behavior of GOTCAB saponin in a water box and model monolayer systems was characterized by molecular dynamics simulations. The results obtained showed that, in the systems used, cholesterol had a decisive effect on the interaction between GOTCAB and phosphocholine or sphingomyelin as well as on its location within the lipid film.

Congo Red as a Supramolecular Carrier System for Doxorubicin: An Approach to Understanding the Mechanism of Action.

The uptake and distribution of doxorubicin in the MCF7 line of breast-cancer cells were monitored by Raman measurements. It was demonstrated that bioavailability of doxorubicin can be significantly enhanced by applying Congo red. To understand the mechanism of doxorubicin delivery by Congo red supramolecular carriers, additional monolayer measurements and molecular dynamics simulations on model membranes were undertaken. Acting as molecular scissors, Congo red particles cut doxorubicin aggregates and incorporated them into small-sized Congo red clusters. The mixed doxorubicin/Congo red clusters were adsorbed to the hydrophilic part of the model membrane. Such behavior promoted transfer through the membrane.

Charge distributions for molecular dynamics simulations from self-consistent polarization method.

Partial atomic charges are important force field parameters. They are usually computed by applying quantum-chemical calculations and the assumed population scheme. In this study polarization consistent scheme of deriving a charge distribution inside solute molecule is proposed. The environment effect is explicitly taken into account by distributing solvent molecules around the solute target. The performed analysis includes a few computational schemes (HF, MP2, B3LYP, and M026X), basis sets (cc-pvnz, n = 2, 3, …, 6), and electrostatically derived charge distributions (KS, CHELP, CHELPG, and HLY). It is demonstrated that the environment effect is very important and cannot be disregarded. The second solvation shell should be included to achieve the charge convergence. Huge corrections to charge distribution are due to induction and dispersion. The B3LYP/cc-pvqz level of theory is recommended for deriving the charges within self-consistent polarization scheme.

Elongation method with intermediate mechanical and electrostatic embedding for geometry optimizations of polymers.

The elongation method with intermediate mechanical and electrostatic embedding (ELG-IMEE) is proposed. The electrostatic embedding uses atomic charges generated by a charge sensitivity analysis (CSA) method and parameterized for three different population analyses, namely, the Merz-Singh-Kollman scheme, the charge model 5, and the atomic polar tensor. The obtained CSA models were tested on two model systems. Test calculations show that the electrostatic embedding provides several times of decrease in the difference of energies of testing and reference calculations in comparison with the conventional elongation approach (ELG). The mechanical embedding is implemented in a combination of the conventional elongation method and the ONIOM approach. Moreover, it was demonstrated that the geometry optimization with the ELG-IMEE reduces the errors in the optimized structures by about one order in root-mean-square deviation, when compared to ELG.

Mixed DPPC/POPC Monolayers: All-atom Molecular Dynamics Simulations and Langmuir Monolayer Experiments.

To elucidate the consequences of the saturated-unsaturated nature of lipid surface films, monolayers formed by an equimolar mixture of 1-palmitoyl-2-oleyl-sn-glycero-3-phosphocholine (POPC) and 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) lipids are investigated in a wide range of surface pressures. As such mixtures share some features with naturally-occurring surfactants, for example the lung surfactant, the systems are studied at the temperature relevant for human body. All-atom molecular dynamics simulations and Langmuir trough experiments are employed. The binary lipid mixture is compared with the corresponding one-component systems. Atomistic-level alterations of monolayer molecular properties upon lateral compression are scrutinized. These involve elevation of lateral ordering of lipid chains, modulation of chain and headgroup orientation, and reduction of lipid hydration. The presence of the unsaturated POPC in the DPPC/POPC mixture reduces the liquid expanded-liquid condensed coexistence region and moderates the phase transition. Simulations predict that nanoscale lipid de-mixing occurs with small transient DPPC clusters emerging due to local fluctuations of the lateral lipid arrangement. A vertical sorting of lipids induced by lateral compression is also observed, with DPPC transferred toward the water phase. Both the conformational lipid alterations due to monolayer compression as well as the existence of lateral dynamic inhomogeneities of the lipid film are potentially pertain to dynamic and non-homogeneous lipid interfacial systems.

Modeling Lung Surfactant Interactions with Benzo[a]pyrene.

By reducing the surface tension of the air-water interface in alveoli, lung surfactant (LS) is crucial for proper functioning of the lungs. It also forms the first barrier against inhaled pathogens. In this study we inspect the interactions of LS models with a dangerous air pollutant, benzo[a]pyrene (BaP). Dipalmitoylphosphatidylcholine (DPPC), 1-palmitoyl-2-oleoylphosphatidylcholine, and their 1:1 mixture are used as LS models. Pressure-area isotherms are employed to study macroscopic properties of the monolayers. We find that addition of BaP has a condensing effect, manifested by lowering the values of surface pressure and shifting the isotherms to smaller areas. Atomistic details of this process are examined by means of molecular dynamics simulations. We show that initially BaP molecules are accumulated in the monolayers. Upon compression, they are forced to the headgroups region and eventually expelled to the subphase. BaP presence results in reduction of monolayer hydration in the hydrophilic region. In the hydrophobic region it induces increased chain ordering, reduction of monolayer fluidity, and advances transition to the liquid condensed phase in the DPPC system.

The selective interactions of cationic tetra-p-guanidinoethylcalix[4]arene with lipid membranes: theoretical and experimental model studies.

Behavior of cationic tetra-p-guanidinoethylcalix[4]arene (CX1) and its building block, p-guanidinoethylphenol (mCX1) in model monolayer lipid membranes was investigated using all atom molecular dynamics simulations and surface pressure measurements. Members of two classes of lipids were taken into account: zwitterionic 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) and anionic 1,2-dimyristoyl-sn-glycero-3-phospho-l-serine sodium salt (DMPS) as models of eukaryotic and bacterial cell membranes, respectively. It was demonstrated that CX1 and mCX1 accumulate near the negatively charged DMPS monolayers. The adsorption to neutral monolayers was negligible. In contrast to mCX1, CX1 penetrated into the hydrophobic part of the monolayer. The latter effect, which is possible due to a flip-flop inversion of the CX1 orientation in the lipid layer compared to the aqueous phase, may be responsible for its antibacterial activity.

Intermediate electrostatic field for the generalized elongation method.

An intermediate electrostatic field is introduced to improve the accuracy of fragment-based quantum-chemical computational methods by including long-range polarizations of biomolecules. The point charge distribution of the intermediate field is generated by a charge sensitivity analysis that is parameterized for five different population analyses, namely, atoms-in-molecules, Hirshfeld, Mulliken, natural orbital, and Voronoi population analysis. Two model systems are chosen to demonstrate the performance of the generalized elongation method (ELG) combined with the intermediate electrostatic field. The calculations are performed for the STO-3G, 6-31G, and 6-31G(d) basis sets and compared with reference Hartree-Fock calculations. It is shown that the error in the total energy is reduced by one order of magnitude, independently of the population analyses used. This demonstrates the importance of long-range polarization in electronic-structure calculations by fragmentation techniques.

Atomic shell structure from the Single-Exponential Decay Detector.

The density of atomic systems is analysed via the Single-Exponential Decay Detector (SEDD). SEDD is a scalar field designed to explore mathematical, rather than physical, properties of electron density. Nevertheless, it has been shown that SEDD can serve as a descriptor of bonding patterns in molecules as well as an indicator of atomic shells [P. de Silva, J. Korchowiec, and T. A. Wesolowski, ChemPhysChem 13, 3462 (2012)]. In this work, a more detailed analysis of atomic shells is done for atoms in the Li-Xe series. Shell populations based on SEDD agree with the Aufbau principle even better than those obtained from the Electron Localization Function, which is a popular indicator of electron localization. A link between SEDD and the local wave vector is given, which provides a physical interpretation of SEDD.

Bond detectors for molecular dynamics simulations, Part I: Hydrogen bonds.

Charge sensitivity analysis in AMBER force-field resolution has been used in quest for detectors of hydrogen bonds (HBs). The process of HB formation was investigated on ab initio classical trajectories (B3LYP/6-31G*) of different nucleobase pairs. Several charge sensitivities, namely: electronegativity, hardness, Fukui function (FF), and polarization matrix, were analyzed. The global and constrained equilibria were considered. It was demonstrated that FF indices and polarization matrix elements are good detectors of HB formation.

Extracting information about chemical bonding from molecular electron densities via single exponential decay detector (SEDD).

The recently introduced molecular descriptor (Single Exponential Decay Detector - SEDD) [P. de Silva, J. Korchowiec, T. A. Wesolowski, ChemPhysChem 2012, 13, 3462] is used to visualize bonding patterns in molecules. In each point of space SEDD is simply related to the electron density: SEDD(r) = In[1/rho2(del(del rho/rho)2)2]. Either experimental or computed densities rho(r) can be used to evaluate SEDD. Here, maps of SEDD are obtained from theoretical densities and reveal such features as core electrons, chemical bonds, lone pairs and delocalization in aromatic systems. It is shown that SEDD provides fingerprints of aromaticity, which can be separated into geometric and electronic effects.

Glycolipid-cholesterol monolayers: towards a better understanding of the interaction between the membrane components.

In this work, the interaction between a synthetic analog of archaeal lipids and cholesterol was studied using Langmuir technique. The lipid, ß-Mal(3)O(C(16+4))(2), contained phytanyl chains attached via two ether bonds to the sn-2 carbon of the glycerol backbone. The preliminary studies showed that monolayers formed with the pure lipid have a liquid-like character; here, a hypothesis that admixing cholesterol to ß-Mal(3)O(C(16+4))(2) could confer a higher rigidity on the films was tested. To check this proposal, two-dimensional miscibility of cholesterol and ß-Mal(3)O(C(16+4))(2) in monomolecular films was studied using surface pressure and surface potential measurements, as well as Brewster angle microscopy and polarization-modulation infrared reflection absorption spectroscopy. The stability of the monomolecular films was evaluated based on thermodynamics of mixing of cholesterol and ß-Mal(3)O(C(16+4))(2). Atomic level information concerning the orientation of molecules and the degree of hydration of polar headgroups was obtained from molecular dynamics simulations.

Revealing the bonding pattern from the molecular electron density using single exponential decay detector: an orbital-free alternative to the electron localization function.

We introduce a new tool (single exponential decay detector: SEDD) to extract information about bonding and localization in atoms, molecules, or molecular assemblies. The practical evaluation of SEDD does not require any explicit information about the orbitals. The only quantity needed is the electron density (calculated or experimental) and its derivatives up to the second order.

Fast orbital localization scheme in molecular fragments resolution.

The method for localizing orbitals on a set of predefined molecular fragments is introduced. Regional localized molecular orbitals (RLMO) are obtained through block diagonalization of the one-electron density matrix and further refinement of the resulting eigenvectors. The algorithm is fast and reliable, as is illustrated by a few examples. Potential applications range from conceptual insight into a chemical bonding to reduced scaling computational techniques. RLMOs are particularly well suited for fragmentation computational methods and for exploiting the locality of electronic correlation in post-HF methods.

Application of explicitly localized molecular orbitals to electronic structure calculations.

We have recently generalized the method for localizing orbitals on a set of predefined molecular fragments [Phys. Chem. Chem. Phys. 2012, 14, 546]. The regional localized molecular orbitals (RLMO) are well suited for exploiting the locality of electronic correlation at post-Hartree-Fock level of theory. In this paper, the adequacy of RLMO representation is tested in the second-order local Møller-Plesset (LMP2) perturbation theory. Two model systems, namely, n-pentadecane and trans-retinal, are considered. Adequacy of RLMO/LMP2 method is discussed in conjunction with 'exact' MP2 and Pipek-Mezey LMP2 calculations. It is demonstrated that RLMO/MP2 method reduces correlation space and reproduces more than 99% of the correlation energy.

The lung surfactant activity probed with molecular dynamics simulations.

The surface of pulmonary alveolar subphase is covered with a mixture of lipids and proteins. This lung surfactant plays a crucial role in lung functioning. It shows a complex phase behavior which can be altered by the interaction with third molecules such as drugs or pollutants. For studying multicomponent biological systems, it is of interest to couple experimental approach with computational modelling yielding atomic-scale information. Simple two, three, or four-component model systems showed to be useful for getting more insight in the interaction between lipids, lipids and proteins or lipids and proteins with drugs and impurities. These systems were studied theoretically using molecular dynamic simulations and experimentally by means of the Langmuir technique. A better understanding of the structure and behavior of lung surfactants obtained from this research is relevant for developing new synthetic surfactants for efficient therapies, and may contribute to public health protection.

Energy partitioning scheme based on self-consistent method for subsystems: populational space approach.

The refinement of the previously proposed energy partitioning scheme, self-consistent charge and configuration method for subsystems (Korchowiec and Uchimaru, J Chem Phys 2000, 112, 1623), is proposed. Our new realization takes rigorously into account all the interactions between subsystems and guarantees proper symmetry of the intermediate wavefunctions. In addition, the scheme is supplemented with natural orbitals for chemical valence to trace the charge reorganization during polarization and charge transfer steps. The water dimer and ammonia borane are used to illustrate the proposed formalism.