Interactions of a Novel Anthracycline with Oligonucleotide DNA and Cyclodextrins in an Aqueous Environment.

Berubicin, a chemotherapy medication belonging to the class of anthracyclines, is simulated in double-stranded DNA sequences and cyclodextrins in an aqueous environment via full-atom molecular dynamics simulations on the time scale of microseconds. The drug is studied in both the neutral and protonated states so as to better comprehend the role of its charge in the formed complexes. The noncovalent berubicin-DNA and berubicin-cyclodextrin complexes are investigated in detail, paying special attention to their thermodynamic description by employing the double decoupling method, the solvent balance method, the weighted solvent accessible surface model, and the linear interaction energy method. A novel approach for extracting the desolvation thermodynamics of the binding process is also presented. Both the binding and desolvation Gibbs energies are decomposed into entropic and enthalpic contributions so as to elucidate the nature of complexation and its driving forces. Selected structural and geometrical properties of all the complexes, which are all stable, are analyzed. Both cyclodextrins under consideration are widely utilized for drug delivery purposes, and a comparative investigation between their bound states with berubicin is carried out.

Molecular simulations of fluoxetine in hydrated lipid bilayers, as well as in aqueous solutions containing ß-cyclodextrin.

Fluoxetine, which is a well-known antidepressant drug, is studied in hydrated cholesterol-free and cholesterol-containing lipid bilayers through unbiased and biased atomistic molecular dynamics simulations. The latter are conducted for the calculation of the potential of mean force (PMF) of fluoxetine along an axis perpendicular to the two leaflets of the bilayer. The PMF indicates that the drug prefers to reside inside the lipid phase and allows us to calculate important thermodynamic properties, such as the Gibbs energy difference of partitioning from the water to the lipid phase and the Gibbs energy barrier for hopping events between the two leaflets of the bilayer. The results from the biased simulations are in accord with the mass density profiles calculated from the unbiased simulations. Moreover, we estimate the effect of fluoxetine mole fraction on the order parameters of the lipid alkyl chains and on the area per lipid. It is also found that fluoxetine forms a hydrogen bond network with lipids and water molecules penetrating into the lipid phase. In addition, fluoxoetine is studied in detail in aqueous solutions containing ß-cyclodextrin. It is observed from unbiased molecular dynamics simulations that the two aforementioned molecules form a noncovalent complex spontaneously and the calculated binding free energy is in agreement with the literature.

In silico study of levodopa in hydrated lipid bilayers at the atomistic level.

This article presents atomistic molecular dynamics and umbrella sampling simulations of levodopa at various concentrations in hydrated cholesterol-free 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) and cholesterol-containing 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) bilayers. Levodopa is the standard medication for Parkinson's disease and is marketed under various trade names; in the context of this article, the levodopa molecule is mostly studied in its zwitterionic form but some results concerning the neutral levodopa are presented as well for comparison purposes. The motivation is to study in detail how levodopa behaves in different hydrated lipid membranes, primarily from the thermodynamic point of view, and reveal aspects of mechanism of its permeation through them. Dependencies of properties on the levodopa concentration are also investigated. Special attention is paid to the calculation of mass density profiles, order parameters and self-diffusion coefficients. Levodopa zwitterions, which form a hydrogen bond network with water and phospholipid molecules, are found to be preferentially located at the water/lipid interface, as well as in the aqueous phase surrounding the cholesterol-free and cholesterol-containing bilayers. This is concluded from the potentials of mean force calculated by umbrella sampling simulations as levodopa is transferred from the lipid to the aqueous phase along an axis perpendicular to the two leaflets of the membranes.

Slip Spring-Based Mesoscopic Simulations of Polymer Networks: Methodology and the Corresponding Computational Code.

In previous work by the authors, a new methodology was developed for Brownian dynamics/kinetic Monte Carlo (BD/kMC) simulations of polymer melts. In this study, this methodology is extended for dynamical simulations of crosslinked polymer networks in a coarse-grained representation, wherein chains are modeled as sequences of beads, each bead encompassing a few Kuhn segments. In addition, the C++ code embodying these simulations, entitled Engine for Mesoscopic Simulations for Polymer Networks (EMSIPON) is described in detail. A crosslinked network of cis-1,4-polyisoprene is chosen as a test system. From the thermodynamic point of view, the system is fully described by a Helmholtz energy consisting of three explicit contributions: entropic springs, slip springs and non-bonded interactions. Entanglements between subchains in the network are represented by slip springs. The ends of the slip springs undergo thermally activated hops between adjacent beads along the chain backbones, which are tracked by kinetic Monte Carlo simulation. In addition, creation/destruction processes are included for the slip springs at dangling subchain ends. The Helmholtz energy of non-bonded interactions is derived from the Sanchez⁻Lacombe equation of state. The isothermal compressibility of the polymer network is predicted from equilibrium density fluctuations in very good agreement with the underlying equation of state and with experiment. Moreover, the methodology and the corresponding C++ code are applied to simulate elongational deformations of polymer rubbers. The shear stress relaxation modulus is predicted from equilibrium simulations of several microseconds of physical time in the undeformed state, as well as from stress-strain curves of the crosslinked polymer networks under deformation.

Exploring the interactions of irbesartan and irbesartan-2-hydroxypropyl-ß-cyclodextrin complex with model membranes.

The interactions of irbesartan (IRB) and irbesartan-2-hydroxypropyl-ß-cyclodextrin (HP-ß-CD) complex with dipalmitoyl phosphatidylcholine (DPPC) bilayers have been explored utilizing an array of biophysical techniques ranging from differential scanning calorimetry (DSC), small angle X-ray scattering (SAXS), ESI mass spectrometry (ESI-MS) and solid state nuclear magnetic resonance (ssNMR). Molecular dynamics (MD) calculations have been also conducted to complement the experimental results. Irbesartan was found to be embedded in the lipid membrane core and to affect the phase transition properties of the DPPC bilayers. SAXS studies revealed that irbesartan alone does not display perfect solvation since some coexisting irbesartan crystallites are present. In its complexed form IRB gets fully solvated in the membranes showing that encapsulation of IRB in HP-ß-CD may have beneficial effects in the ADME properties of this drug. MD experiments revealed the topological and orientational integration of irbesartan into the phospholipid bilayer being placed at about 1nm from the membrane centre.

Rational design, efficient syntheses and biological evaluation of N,N'-symmetrically bis-substituted butylimidazole analogs as a new class of potent Angiotensin II receptor blockers.

A series of symmetrically bis-substituted imidazole analogs bearing at the N-1 and N-3 two biphenyl moieties ortho substituted either with tetrazole or carboxylate functional groups was designed based on docking studies and utilizing for the first time an extra hydrophobic binding cleft of AT1 receptor. The synthesized analogs were evaluated for their in vitro antagonistic activities (pA2 values) and binding affinities (-logIC50 values) to the Angiotensin II AT1 receptor. Among them, the potassium (-logIC50 = 9.04) and the sodium (-logIC50 = 8.54) salts of 4-butyl-N,N'-bis{[2'-(2H-tetrazol-5-yl)biphenyl-4-yl]methyl}imidazolium bromide (12a and 12b, respectively) as well as its free acid 11 (-logIC50 = 9.46) and the 4-butyl-2-hydroxymethyl-N,N'-bis{[2'-(2H-tetrazol-5-yl)biphenyl-4-yl]methyl}imidazolium bromide (14) (-logIC50 = 8.37, pA2 = 8.58) showed high binding affinity to the AT1 receptor and high antagonistic activity (potency). The potency was similar or even superior to that of Losartan (-logIC50 = 8.25, pA2 = 8.25). On the contrary, 2-butyl-N,N'-bis{[2'-[2H-tetrazol-5-yl)]biphenyl-4-yl]methyl}imidazolium bromide (27) (-logIC50 = 5.77) and 2-butyl-4-chloro-5-hydroxymethyl-N,N'-bis{[2'-[2H-tetrazol-5-yl)]biphenyl-4-yl]methyl}imidazolium bromide (30) (-logIC50 = 6.38) displayed very low binding affinity indicating that the orientation of the n-butyl group is of primary importance. Docking studies of the representative highly active 12b clearly showed that this molecule has an extra hydrophobic binding feature compared to prototype drug Losartan and it fits to the extra hydrophobic cavity. These results may contribute to the discovery and development of a new class of biologically active molecules through bis-alkylation of the imidazole ring by a convenient and cost effective synthetic strategy.

Conformational analysis of two novel cytotoxic C2-substituted pyrrolo[2,3-f]quinolines in aqueous media, organic solvents, membrane bilayers and at the putative active site.

We have performed: (i) conformational analysis of two novel cytotoxic C2-substituted pyrrolo[2,3-f]quinolines 5e and 5g in deuterated dimethylsulfoxide (DMSO-d(6)) utilizing NOE results from NMR spectroscopy; (ii) molecular dynamics (MD) calculations in water, DMSO and dimyristoyl phosphatidylcholine bilayers and (iii) molecular docking and MD calculations on DNA nucleotide sequences. The obtained results for the two similar in structure molecules showed differences in: (i) their conformational properties in silico and in media that reasonably simulate the biological environment; (ii) the way they are incorporated into the lipid bilayers and therefore their diffusion ability and (iii) molecular docking capacity as it is depicted from their different binding scores.

Comparative study of the AT1 receptor prodrug antagonist candesartan cilexetil with other sartans on the interactions with membrane bilayers.

Drug-membrane interactions of the candesartan cilexetil (TCV-116) have been studied on molecular basis by applying various complementary biophysical techniques namely differential scanning calorimetry (DSC), Raman spectroscopy, small and wide angle X-ray scattering (SAXS and WAXS), solution ¹H and ¹³C nuclear magnetic resonance (NMR) and solid state ¹³C and ³¹P (NMR) spectroscopies. In addition, ³¹P cross polarization (CP) NMR broadline fitting methodology in combination with ab initio computations has been applied. Finally molecular dynamics (MD) was applied to find the low energy conformation and position of candesartan cilexetil in the bilayers. Thus, the experimental results complemented with in silico MD results provided information on the localization, orientation, and dynamic properties of TCV-116 in the lipidic environment. The effects of this prodrug have been compared with other AT1 receptor antagonists hitherto studied. The prodrug TCV-116 as other sartans has been found to be accommodated in the polar/apolar interface of the bilayer. In particular, it anchors in the mesophase region of the lipid bilayers with the tetrazole group oriented toward the polar headgroup spanning from water interface toward the mesophase and upper segment of the hydrophobic region. In spite of their localization identity, their thermal and dynamic effects are distinct pointing out that each sartan has its own fingerprint of action in the membrane bilayer, which is determined by the parameters derived from the above mentioned biophysical techniques.

Dual inhibitors for aspartic proteases HIV-1 PR and renin: advancements in AIDS-hypertension-diabetes linkage via molecular dynamics, inhibition assays, and binding free energy calculations.

Human immunodeficiency virus type 1 protease (HIV-1 PR) and renin are primary targets toward AIDS and hypertension therapies, respectively. Molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) free-energy calculations and inhibition assays for canagliflozin, an antidiabetic agent verified its effective binding to both proteins (ΔG(pred) = -9.1 kcal mol(-1) for canagliflozin-renin; K(i,exp)= 628 nM for canagliflozin-HIV-1 PR). Moreover, drugs aliskiren (a renin inhibitor) and darunavir (an HIV-1 PR inhibitor) showed high affinity for HIV-1 PR (K(i,exp)= 76.5 nM) and renin (K(i,pred)= 261 nM), respectively. Importantly, a high correlation was observed between experimental and predicted binding energies (r(2) = 0.92). This study suggests that canagliflozin, aliskiren, and darunavir may induce profound effects toward dual HIV-1 PR and renin inhibition. Since patients on highly active antiretroviral therapy (HAART) have a high risk of developing hypertension and diabetes, aliskiren-based or canagliflozin-based drug design against HIV-1 PR may eliminate these side-effects and also facilitate AIDS therapy.

Computational studies of darunavir into HIV-1 protease and DMPC bilayer: necessary conditions for effective binding and the role of the flaps.

Human immunodeficiency virus type 1 protease (HIV-1 PR) is one of the main targets toward AIDS therapy. We have selected the potent drug darunavir and a weak inhibitor (fullerene analog) as HIV-1 PR substrates to compare protease's conformational features upon binding. Molecular dynamics (MD), molecular mechanics Poisson-Boltzmann surface area (MM-PBSA), and quantum-mechanical (QM) calculations indicated the importance of the stability of HIV-1 PR flaps toward effective binding: a weak inhibitor may induce flexibility to the flaps, which convert between closed and semiopen states. A water molecule in the darunavir-HIV-1 PR complex bridged the two flap tips of the protease through hydrogen bonding (HB) interactions in a stable structure, a feature that was not observed for the fullerene-HIV-1 PR complex. Additionally, despite that van der Waals interactions and nonpolar contribution to solvation favored permanent fullerene entrapment into the cavity, these interactions alone were not sufficient for effective binding; enhanced electrostatic interactions as observed in the darunavir-complex were the crucial component of the binding energy. An alternative pathway to the usual way of a ligand to access the cavity was also observed for both compounds. Each ligand entered the binding cavity through an opening between the one flap of the protease and a neighboring loop. This suggested that access to the cavity is not necessarily regulated by flap opening. Darunavir exerts its biological action inside the cell, after crossing the membrane barrier. Thus, we also initiated a study on the interactions between darunavir and the DMPC bilayer to reveal that the drug was accommodated inside the bilayer in conformations that resembled its structure into HIV-1 PR, being stabilized via HBs with the lipids and water molecules.