Identification of 1H-indene-(1,3,5,6)-tetrol derivatives as potent pancreatic lipase inhibitors using molecular docking and molecular dynamics approach.

Pancreatic lipase is a potential therapeutic target to treat diet-induced obesity in humans, as obesity-related diseases continue to be a global problem. Despite intensive research on finding potential inhibitors, very few compounds have been introduced to clinical studies. In this work, new chemical scaffold 1H-indene-(1,3,5,6)-tetrol was proposed using knowledge-based approach, and 36 inhibitors were derived by modifying its functional groups at different positions in scaffold. To explore binding affinity and interactions of ligands with protein, CDOCKER and AutoDock programs were used for molecular docking studies. Analyzing results of rigid and flexible docking algorithms, inhibitors C_12, C_24, and C_36 were selected based on different properties and high predicted binding affinities for further analysis. These three inhibitors have different moieties placed at different functional groups in scaffold, and to characterize structural rationales for inhibitory activities of compounds, molecular dynamics simulations were performed (500 nSec). It has been shown through simulations that two structural fragments (indene and indole) in inhibitor can be treated as isosteric structures and their position at binding cleft can be replaced by each other. Taking into account these information, two lines of inhibitors can further be developed, each line based on a different core scaffold, that is, indene/indole.

Interaction of cisplatin and two potential antitumoral platinum(II) complexes with a model lipid membrane: a combined NMR and MD study.

In this study, the interaction of cisplatin (1) and two potential antitumoral Pt(II) complexes (2 and 3) with a model DMPC bilayer was investigated by multinuclear NMR spectroscopy and MD simulations in order to understand its implication for the different antitumoral properties shown by the three complexes. In particular, (31)P, (13)C and (2)H solid state NMR experiments were performed to obtain information on the phase structure, phase transitions and structural and dynamic changes in the phospholipid bilayer upon interaction with the platinum complexes. On the other hand, MD calculations yielded free energy profiles for the different complexes across the bilayer; the results were analysed to obtain MD predictions on complex distribution with respect to the bilayer, as well as to establish their effects on the conformational equilibrium of the DMPC acyl chains. The combination of NMR and MD approaches highlighted that, whereas the more hydrophilic cisplatin tends to remain in the polar head group region causing a decrease in flexibility of the bilayer, the two new complexes enter into the bilayer. In particular, complex 2 is preferentially located relatively close to the surface, only slightly affecting the bilayer structure and mobility, while complex 3 penetrates more deeply, strongly perturbing the bilayer and giving rise to lateral phase separation.

Structural, functional, and stability change predictions in human telomerase upon specific point mutations.

Overexpression of telomerase is one of the hallmarks of human cancer. Telomerase is important for maintaining the integrity of the ends of chromosomes, which are called telomeres. A growing number of human disease syndromes are associated with organ failure caused by mutations in telomerase (hTERT or hTR). Mutations in telomerase lead to telomere shortening by decreasing the stability of the telomerase complex, reducing its accumulation, or directly affecting its enzymatic activity. In this work, potential human telomerase mutations were identified by a systematic computational approach. Moreover, molecular docking methods were used to predict the effects of these mutations on the affinity of certain ligands (C_9i, C_9k, 16A, and NSC749234). The C_9k inhibitor had the best binding affinity for wild-type (WT) telomerase. Moreover, C_9i and C_9k had improved interactions with human telomerase in most of the mutant models. The R631 and Y717 residues of WT telomerase formed interactions with all studied ligands and these interactions were also commonly found in most of the mutant models. Residues forming stable interactions with ligands in molecular dynamics (MD) were traced, and the MD simulations showed that the C_9k ligand formed different conformations with WT telomerase than the C_9i ligand.

Structural and dynamic insights on the EmrE protein with TPP+ and related substrates through molecular dynamics simulations.

EmrE is a bacterial transporter protein that forms an anti-parallel homodimer with four transmembrane helices in each monomer. EmrE transports positively charged aromatic compounds, such as TPP+ and its derivatives. We performed molecular dynamics (MD) simulations of EmrE in complex with TPP+, MeTPP+, and MBTPP+ embedded in a membrane. The detailed molecular properties and interactions were analysed for all EmrE-ligand complexes. Our MD results identified that Lys22, Tyr40, Phe44, Trp45, and Trp63 formed potential π interactions with all three ligands and further confirmed the essential role of Glu14. Moreover, distance analysis and structural changes in the EmrE translocation pathway suggest that ligand recognition and protein conformational changes depend on the structural properties of the substrate. Analysis of the movement of the ligand in the protein binding site and rotation of the ligand's aromatic rings confirm that substrates with aromatic moieties, such as MBTPP+, exhibit relatively stable binding to EmrE. Interestingly, the aromatic rings of Tyr40, Phe44, Trp45, and Trp63 underwent parallel movements with the aromatic rings of TPP+. Based on the MD results, we propose that π interactions, as well as the mutual rotation of the aromatic rings in the protein and ligand, can be regarded as sources of ligand movement, and thus, the whole complex may work as a "molecular propeller".

Electrostatic and non-electrostatic contributions to the binding free energies of anthracycline antibiotics to DNA.

The knowledge about molecular factors driving simple ligand-DNA interactions is still limited. The aim of the present study was to investigate the electrostatic and non-electrostatic contributions to the binding free energies of anthracycline compounds with DNA. Theoretical calculations based on continuum methods (Poisson-Boltzmann and solvent accessible surface area) were performed to estimate the binding free energies of five selected anthracycline ligands (daunomycin, adriamycin, 9-deoxyadriamycin, hydroxyrubicin, and adriamycinone) to DNA. The free energy calculations also took into account the conformational change that DNA undergoes upon ligand binding. This conformational change appeared to be very important for estimating absolute free energies of binding. Our studies revealed that the absolute values of all computed contributions to the binding free energy were quite large compared to the total free energy of binding. However, the sum of these large positive and negative values produced a small negative value of the free energy around -10 kcal/mol. This value is in good agreement with experimental data. Experimental values for relative binding free energies were also reproduced for charged ligands by our calculations. Together, it was found that the driving force for ligand-DNA complex formation is the non-polar interaction between the ligand and DNA even if the ligand is positively charged.

The role of amphotericin B amino group basicity in its antifungal action. A theoretical approach.

The role of basicity of the amine group of amphotericin B in the molecular mechanism of antifungal activity of this antibiotic has been investigated by AM1 and MNDO quantum chemistry methods. Calculations of proton affinity of the amine group, as a measure of its basicity, for appropriate models of free amphotericin B and its N-alkyl derivatives were carried out. These studies were preceded by a critical examination of the usefulness and reliability of both methods to predict the proton affinities of several aliphatic amines. It has been concluded that the diminution of protonability of the substituted amine group of amphotericin B correlates with the decrease of antifungal activity of the appropriate derivatives of antibiotic. It was experimentally demonstrated (A. Czerwinski et al., J. Antibiot. 44 (1991) 979) that the introduction of additional amine groups in such a derivative restores antifungal activity of the compound. In our studies it was evidenced, using theoretical methods, that the proton affinity of this additional amine group is similar to that in free amphotericin B.

Conformational analysis of Amphotericin B.

Within a theoretical approach to the problem of antifungal action of Amphotericin B (AmB), a conformational analysis of the neutral and zwitterionic form of this antibiotic in vacuo was performed by the MM2P and AM1 methods. The analysis was carried out with regard to the mutual orientation of the macrolidic and glycosidic fragments of the molecule, which is defined by the phi and psi steric angles. This orientation defines the overall shape of the molecule and is postulated to be important for the antifungal action of the drug. As a result of the MM2P calculations, phi, psi steric energy and population maps were prepared. Several conformers were found on these maps but only two of them (one each for the zwitterionic and the neutral forms of the antibiotic) were previously observed experimentally for isolated molecules. Our other calculated conformers were not observed experimentally but we propose that they may also appear in the AmB channel structure. The results of our conformational analysis were compared with experimental NMR data (nuclear Overhauser effects between selected hydrogen atoms) obtained previously. New structural information obtained for AmB in the present work will be useful for building a molecular model of AmB-target interactions as well as for designing new derivatives of AmB.

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.

Thermodynamic and electrostatic properties of ternary Oxytricha nova TEBP-DNA complex.

Telomeres constitute the nucleoprotein ends of eukaryotic chromosomes which are essential for their proper function. Telomere end binding protein (TEBP) from Oxytricha nova was among the first telomeric proteins, which were well characterized biologically. TEBP consists of two protein subunits (alpha, beta) and forms a ternary complex with single stranded telomeric DNA containing tandem repeats TTTTGGGG. This work presents the characterization of the thermodynamic and electrostatic properties of this complex by computational chemistry methods (continuum Poisson-Boltzmann and solvent accessible surface calculations). Our calculations give a new insight into molecular properties of studied system. Based on the thermodynamic analysis we provide a rationale for the experimental observation that alpha and ssDNA forms a binary complex and the beta subunit joins alpha:ssDNA complex only after the latter is formed. Calculations of distribution of the molecular electrostatic potential for protein subunits alone and for all possible binary complexes revealed the important role of the "guiding funnel" potential generated by alpha:ssDNA complex. This potential may help the beta subunit to dock to the already formed alpha:DNA intermediate in highly steric and electrostatic favorable manner. Our pK(a) calculations of TEBP are able to explain the experimental mobility shifts of the complex in electrophoretic non-denaturating gels.

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 Waal's forces.

Molecular properties of amphotericin B membrane channel: a molecular dynamics simulation.

Amphotericin B is a powerful but toxic antifungal antibiotic that is used to treat systemic infections. It forms ionic membrane channels in fungal cells. These antibiotic/sterol channels are responsible for the leakage of ions, which causes cell destruction. The detailed molecular properties and structure of amphotericin B channels are still unknown. In the current study, two molecular dynamic simulations were performed of a particular model of amphotericin B/cholesterol channel. The water and phospholipid environment were included in our simulations, and the results obtained were compared with available experimental data. It was found that it is mainly the hydrogen bonding interactions that keep the channel stable in its open form. Our study also revealed the important role of the intermolecular interactions among the hydroxyl, amino, and carboxyl groups of the channel-forming molecules; in particular, some hydroxyl groups stand out as new "hot spots" that are potentially useful for chemotherapeutic investigations. Our results also help to clarify why certain antibiotic derivatives, with a blocked amino group, are less active. We present a hypothesis for the role of membrane lipids and cholesterol in the channel.

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 (phi,psi)=(290 degrees,180 degrees) 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.