Maciej Baginski

Comparative molecular dynamics simulations of amphotericin B–cholesterol/ergosterol membrane channels

Amphotericin B (AmB) is a very effective anti-fungal polyene macrolide antibiotic whose usage is limited by its toxicity. Lack of a complete understanding of AmBs molecular mechanism has impeded attempts to design less toxic AmB derivatives. The antibiotic is known to interact with sterols present in the cell membrane to form ion channels that disrupt membrane function. The slightly higher affinity of AmB toward ergosterol (dominant sterol in fungal cells) than cholesterol (mammalian sterol) is regarded as the most essential factor on which antifungal chemotherapy is based. To study these differences at the molecular level, two realistic model membrane channels containing molecules of AmB, sterol (cholesterol or ergosterol), phospholipid, and water were studied by molecular dynamics (MD) simulations. Comparative analysis of the simulation data revealed that the sterol type has noticeable effect on the properties of AmB membrane channels. In addition to having a larger size, the AmB channel in the ergosterol-containing membrane has a more pronounced pattern of intermolecular hydrogen bonds. The interaction between the antibiotic and ergosterol is more specific than between the antibiotic and cholesterol. These observed differences suggest that the channel in the ergosterol-containing membrane is more stable and, due to its larger size, would have a higher ion conductance. These observations are in agreement with experiments.

Amphotericin B and its new derivatives - mode of action.

Amphotericin B (AmB) is a well known antifungal and antiprotozoal antibiotic used in the clinic for several decades. Clinical applications of AmB, however, are limited by its nephrotoxicity and many other acute side effects which are not acceptable by patients when their life is not threaten. In order to improve the therapeutic index of this drug, lipid formulations have been introduced and many efforts have been made to obtain less toxic AmB derivatives by chemical modifications of the parent drug. This review presents concise knowledge about this fascinating compound and a critical review of the data published within last few years about the mechanism of action of this antibiotic. In particular, in the present work we discuss: i) structure and properties of AmB and its recently synthesized new derivatives; ii) antifungal and antileishmanial activity and toxicity of these compounds; and iii) mode of action of AmB and its derivatives at cellular and molecular levels, with particular attention paid to interactions of AmB and different components of cellular membranes.

Comparative conformational analysis of cholesterol and ergosterol by molecular mechanics

A comparative conformational analysis of cholesterol and ergosterol has been carried out using molecular mechanics methods. These studies are aimed at giving a better understanding of the molecular nature of the interaction of these sterols with polyene macrolide antibiotics. Structures of cholesterol and ergosterol determined by X-ray methods have been used as initial geometries of these molecules for force field calculations. The calculation of steric energy has also been made for conformations which do not appear in the crystal. The latter conformers have different conformations of the side chain as well as different conformations of rings A and D. The rotational barriers around bonds C17–C20 and C20–C22 have also been calculated. The results obtained on differences and similarities in the conformations of cholesterol and ergosterol allow us to postulate a mechanism for differential interaction with the antibiotics. The relatively rigid side chain of ergosterol (stretched molecule) in comparison with the flexible side chain of cholesterol (bent molecule), allows better intermolecular contact of the first sterol molecule with a polyene macrolide and in consequence facilitates complex formation involving Van der Waals forces.

On the possibility of the amphotericin B-sterol complex formation in cholesterol- and ergosterol-containing lipid bilayers: a molecular dynamics study.

Amphotericin B (AmB) is a well-known membrane-active antibiotic that has been used to treat systemic fungal infections for more than 45 years. Therapeutic application of AmB is based on the fact that it is more active against ergosterol-containing membranes of fungal cells than against mammalian membranes with cholesterol. In this paper, we examine the hypothesis according to which the selectivity of the AmBs membrane action originates from its different ability to form the binary complexes with the relevant sterols. To this end, molecular dynamics simulations were performed for systems containing the preformed models of AmB/sterol complexes embedded in lipid bilayers containing either cholesterol or ergosterol. The initial structures of the studied binary associates were selected on the basis of a systematic scan of all possible mutual positions and orientations of the two molecules. The results obtained demonstrate that in general the complexes with ergosterol are more stable on the 100 ns time scale. Furthermore, on the basis of motional correlation analysis, taking into account the effects of lipid environment, we propose that, within the sterol-enriched liquid-ordered membrane phases, AmB molecules exhibit a greater tendency to bind ergosterol than cholesterol. The analysis of the interactions suggests that this affinity difference is of enthalpic origin and may arise from the considerable difference in the energy of the van der Waals interactions between AmB and the two types of sterols. Thus, our current results: (i) support the hypothesis that binary AmB/sterol complexes form within a lipid membrane and (ii) suggest that the higher toxicity may at least partly be attributed to the higher affinity of AmB for ergosterol than for cholesterol within a lipid membrane environment.

Induction of unique structural changes in guanine-rich DNA regions by the triazoloacridone C-1305, a topoisomerase II inhibitor with antitumor activities

We recently reported that the antitumor triazoloacridone, compound C-1305, is a topoisomerase II poison with unusual properties. In this study we characterize the DNA interactions of C-1305 in vitro, in comparison with other topoisomerase II inhibitors. Our results show that C-1305 binds to DNA by intercalation and possesses higher affinity for GC- than AT-DNA as revealed by surface plasmon resonance studies. Chemical probing with DEPC indicated that C-1305 induces structural perturbations in DNA regions with three adjacent guanine residues. Importantly, this effect was highly specific for C-1305 since none of the other 22 DNA interacting drugs tested was able to induce similar structural changes in DNA. Compound C-1305 induced stronger structural changes in guanine triplets at higher pH which suggested that protonation/deprotonation of the drug is important for this drug-specific effect. Molecular modeling analysis predicts that the zwitterionic form of C-1305 intercalates within the guanine triplet, resulting in widening of both DNA grooves and aligning of the triazole ring with the N7 atoms of guanines. Our results show that C-1305 binds to DNA and induces very specific and unusual structural changes in guanine triplets which likely plays an important role in the cytotoxic and antitumor activity of this unique compound.

Conformational properties of amphotericin B amide derivatives – impact on selective toxicity

Even though it is highly toxic, Amphotericin B (AmB), an amphipathic polyene macrolide antibiotic, is used in the treatment of severe systemic fungal infections as a life-saving drug. To examine the influence of conformational factors on selective toxicity of these compounds, we have investigated the conformational properties of five AmB amide derivatives. It was found that the extended conformation with torsional angles (φ,ψ)=(290°,180° ) is a common minimum of the potential energy surfaces (PES) of unsubstituted AmB and its amide derivatives. The extended conformation of the studied compounds allows for the formation of an intermolecular hydrogen bond network between adjacent antibiotic molecules in the open channel configuration. Therefore, the extended conformation is expected to be the dominant conformer in an open AmB (or its amide derivatives) membrane channel. The derivative compounds for calculations were chosen according to their selective toxicity compared to AmB and they had a wide range of selective toxicity. Except for two AmB derivatives, the PES maps of the derivatives reveal that the molecules can coexist in more than one conformer. Taking into account the cumulative conclusions drawn from the earlier MD simulation studies of AmB membrane channel, the results of the potential energy surface maps, and the physical considerations of the molecular structures, we hypothesize a new model of structure-selective toxicity of AmB derivatives. In this proposed model the presence of the extended conformation as the only well defined global conformer for AmB derivatives is taken as the indicator of their higher selective toxicity. This model successfully explains our results. To further test our model, we also investigated an AmB derivative whose selective toxicity has not been experimentally measured before. Our prediction for the selective toxicity of this compound can be tested in experiments to validate or invalidate the proposed model.

The effect of sterols on amphotericin B self-aggregation in a lipid bilayer as revealed by free energy simulations.

Amphotericin B (AmB) is an effective but toxic antifungal drug, known to increase the permeability of the cell membrane, presumably by assembling into transmembrane pores in a sterol-dependent manner. The aggregation of AmB molecules in a phospholipid bilayer is, thus, crucial for the drugs activity. To provide an insight into the molecular nature of this process, here, we report an atomistic molecular dynamics simulation study of AmB head-to-head dimerization in a phospholipid bilayer, a possible early stage of aggregation. To compare the effect of sterols on the thermodynamics of aggregation and the architecture of the resulting AmB-AmB complexes, free energy profiles for the dimerization in ergosterol- or cholesterol-containing and sterol-free membranes are derived from the simulations. These profiles demonstrate that although AmB dimers are formed in all the systems studied, they are significantly less favorable in the bilayer with ergosterol than in the cholesterol-containing or sterol-free ones. We investigate the structural and energetic determinants of this difference and discuss its consequences for the AmB mechanism of action.

MM/PBSA analysis of molecular dynamics simulations of bovine β‐lactoglobulin: Free energy gradients in conformational transitions?

The pH‐driven opening and closure of β‐lactoglobulin EF loop, acting as a lid and closing the internal cavity of the protein, has been studied by molecular dynamics (MD) simulations and free energy calculations based on molecular mechanics/Poisson–Boltzmann (PB) solvent‐accessible surface area (MM/PBSA) methodology. The forms above and below the transition pH differ presumably only in the protonation state of residue Glu89. MM/PBSA calculations are able to reproduce qualitatively the thermodynamics of the transition. The analysis of MD simulations using a combination of MM/PBSA methodology and the colony energy approach is able to highlight the driving forces implied in the transition. The analysis suggests that global rearrangements take place before the equilibrium local conformation is reached. This conclusion may bear general relevance to conformational transitions in all lipocalins and proteins in general. Proteins 2005.

Distribution of electrostatic potential around amphotericin B and its membrane targets

Abstract Amphotericin B is a well-known anti-fungal polyene macrolide antibiotic used widely in the treatment of systemic fungal infections. The molecular mechanism of the anti-fungal action of amphotericin B is still not elucidated and is under study. In the present article molecular electrostatic properties of amphotericin B and its membrane targets (cholesterol and ergosterol) were studied. The electrostatic potential around all three molecules was calculated using Poisson-Boltzmann methods. The physiological environment (water and membrane) was taken into account in these calculations. The distribution of electrostatic potential was presented on molecular surfaces. It was found that both targets have quite different electrostatic patterns that may be important for interaction with the antibiotic. The calculations of electrostatic potential for Amphotericin B also revealed some new features about the hydrophobic/hydrophilic pattern of this drug. The knowledge of how the electrostatic potential distributes around the studied molecules will be useful to build a model of amphotericin B-sterol channel.

Influence of a lipid bilayer on the conformational behavior of amphotericin B derivatives — A molecular dynamics study

Amphotericin B (AmB) is an effective but very toxic antifungal antibiotic. In our laboratory a series of AmB derivatives of improved selectivity of action was synthesized and tested. To understand molecular basis of this improvement, comparative conformational studies of amphotericin B and its two more selective derivatives were carried out in an aqueous solution and in a lipid membrane. These molecular simulation studies revealed that within a membrane environment the conformational behavior of the derivatives differs significantly from the one observed for the parent molecule. Possible reasons for such a difference are analyzed. Furthermore, we hypothesize that the observed conformational transition within the polar head of AmB derivatives may lead to destabilization of antibiotic-induced transmembrane channels. Consequently, the selective toxicity of the derivatives should increase as ergosterol-rich liquid-ordered domains are more rigid and conformationally ordered than their cholesterol-containing counterparts, and as such may better support less stable channel structure.