Grigorios Megariotis

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 proteases 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.

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 AT₁ 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.

Equation of State Based Slip Spring Model for Entangled Polymer Dynamics

A mesoscopic, mixed particle- and field-based Brownian dynamics methodology for the simulation of entangled polymer melts has been developed. Polymeric beads consist of several Kuhn segments, and their motion is dictated by the Helmholtz energy of the sample, which is a sum of the entropic elasticity of chain strands between beads, slip springs, and nonbonded interactions. Following earlier works in the field [Phys. Rev. Lett. 2012, 109, 148302], the entanglement effect is introduced by the slip springs, which are springs connecting either nonsuccessive beads on the same chain or beads on different polymer chains. The terminal positions of slip springs are altered during the simulation through a kinetic Monte Carlo hopping scheme, with rate-controlled creation/destruction processes for the slip springs at chain ends. The rate constants are consistent with the free energy function employed and satisfy microscopic reversibility at equilibrium. The free energy of nonbonded interactions is derived from an app...

Comparative study of interactions of aliskiren and AT1 receptor antagonists with lipid bilayers.

The renin-angiotensin-aldosterone system (RAAS) plays a key role in the regulation of blood pressure. Renin is the rate limiting enzyme of the RAAS and aliskiren is a highly potent and selective inhibitor of the human renin. Renin is known to be active both in the circulating blood stream as well as locally, when bound to the (pro)-renin receptor ((P)RR). In this study we have investigated a possible mechanism of action of aliskiren, in which its accumulation in the plasma membrane is considered as an essential step for effective inhibition. Aliskirens interactions with model membranes (cholesterol rich and poor) have been investigated by applying different complementary techniques: differential scanning calorimetry (DSC), Raman spectroscopy, magic angle spinning (MAS) nuclear magnetic resonance (NMR) spectroscopy and small- and wide-angle X-ray scattering (SAXS and WAXS). In addition, in silico molecular dynamics (MD) calculations were applied for further confirmation of the experimental data. Aliskirens thermal effects on the pre- and main transition of dipalmitoyl-phosphatidylcholine (DPPC) membranes as well as its topographical position in the bilayer show striking similarities to those of angiotensin II type 1 receptor (AT1R) antagonists. Moreover, at higher cholesterol concentrations aliskiren gets expelled from the membrane just as it has been recently demonstrated for the angiotensin receptor blocker (ARB) losartan. Thus, we propose that both the AT1R and the (P)RR-bound renin active sites can be efficiently blocked by membrane-bound ARBs and aliskiren when cholesterol rich membrane rafts/caveolae are formed in the vicinity of the receptors.

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.

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.

Multiscale simulations of hexa-peri-hexabenzocoronene and hexa-n-dodecyl-hexa-peri-hexabenzocoronene

This study concerns atomistic and coarse-grained molecular dynamics simulations of two disk-shaped molecules in the crystalline phase. The coarse-grained models were developed by applying the Iterative Boltzmann Inversion method which is a systematic coarse-graining method. In particular, a set of radial distribution functions and intramolecular distributions are reproduced at the coarse-grained level. Before applying coarse-graining, a reliable atomistic model was developed to reproduce the main experimental properties of these molecules. The crystalline phases are analyzed in terms of the Saupe ordering tensor.

Chapter 8:Theoretical Studies of Interactions in Nanomaterials and Biological Systems

In this article, we review some of our recent work and related articles by other authors, which deal with various aspects of nanomaterials and biological complexes and in particular, the effects of interactions on selected properties.

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.