Molecular aspects of the interaction between amphotericin B and a phospholipid bilayer: molecular dynamics studies

Amphotericin B (AmB) is a polyene macrolide antibiotic used to treat systemic fungal infections. The molecular mechanism of AmB action is still only partly characterized. AmB interacts with cell-membrane components and forms membrane channels that eventually lead to cell death. The interaction between AmB and the membrane surface can be regarded as the first (presumably crucial) step on the way to channel formation. In this study molecular dynamics simulations were performed for an AmB–lipid bilayer model in order to characterize the molecular aspects of AmB–membrane interactions. The system studied contained a box of 200 dimyristoylphosphatidylcholine (DMPC) molecules, a single AmB molecule placed on the surface of the lipid bilayer and 8,065 water molecules. Two molecular dynamics simulations (NVT ensemble), each lasting 1 ns, were performed for the model studied. Two different programs, CHARMM and NAMD2, were used in order to test simulation conditions. The analysis of MD trajectories brought interesting information concerning interactions between polar groups of AmB and both DMPC and water molecules. Our studies show that AmB preferentially took a vertical position, perpendicular to the membrane surface, with no propensity to enter the membrane. Our finding may suggest that a single AmB molecule entering the membrane is very unlikely.Figure The figure presents the whole structure of the system simulated—starting point. AmB is presented as a space-filling model, DMPC molecules—green sticks, water molecules—red sticks     

Toward the understanding of the structure and dynamics of protein-carbohydrate interactions: molecular dynamics studies of the complexes between hevein and oligosaccharidic ligands

Herein we study, through all atom molecular dynamics simulations, the complex between hevein and two N-acetylated chitin oligomers, namely N,N(')-diacetylchitobiose and N,N('),N(")-triacetylchitotriose. The results of the simulations for two disaccharide complexes and one trisaccharide complex show that a carbohydrate oligomer is able to move on the surface of the relatively flat binding pocket of hevein, therefore occupying different binding subpockets. Statistical analysis methods were also applied in order to define the principal overall motions in the complexes, showing how the different ligands in the simulations modulate the protein motions. The oligosaccharide binding can be considered as defined by a subtle balance between enthalpic (formation of intermolecular interactions between the ligand and the receptor) and entropic (due mainly to the possibility for the sugar to move on the surface of the protein domain) effects, determining multiple binding conformations. This structural and dynamical view could parallel the results obtained by regularly used restrained MD simulations based on NOE NMR data that provide a well defined structure for both the disaccharide and trisaccharide complexes, and agrees with the observations for longer oligosaccharide chains.     

Cellular uptake of transportan 10 and its analogs in live cells: Selectivity and structure-activity relationship studies

Transportan 10 (TP10) is an amphipathic cell-penetrating peptide with high translocation ability. In order to obtain more details of structure-activity relationship of TP10, we evaluated the effects of structure and charge on its translocation ability. Our results demonstrated that disrupting the helical structure or Arg substitution could remarkably decrease the cellular uptake of TP10. However, increasing the number of positive charge was an effective strategy to enhance translocation ability of TP10. Furthermore, the molecular dynamics simulation supported the results derived from experiments, suggesting that higher membrane disturbance leads to higher cellular uptake of peptides. In addition, our study also demonstrated TP10 and its analogs preferentially entered cancer cells rather than normal cells. The uptake selectivity toward cancer cells makes TP10 and its analogs as potent CPPs for drug delivery.     

The challenges of understanding glycolipid functions: An open outlook based on molecular simulations

Glycolipids are the most complex lipid type in cell membranes, characterized by a great diversity of different structures and functions. The underlying atomistic/molecular interactions and mechanisms associated with these functions are not well understood. Here we discuss how atomistic and molecular simulations can be used to shed light on the role of glycolipids in membrane structure and dynamics, receptor function, and other phenomena related to emergence of diseases such as Parkinson's. The cases we discuss highlight the challenge to understand how glycolipids function in cell membranes, and the significant added value that one would gain by bridging molecular simulations with experiments. This article is part of a Special Issue entitled Tools to study lipid functions.     

On the importance of van der Waals interaction in the groove binding of DNA with ligands: restrained molecular dynamics study

Comparison of interaction energy between an oligonucleotide and a DNA-binding ligand in the minor and major groove modes was made by use of restrained molecular dynamics. Distortion in DNA was found for the major groove mode whereas less significant changes for both ligand and DNA were detected for the minor groove binding after molecular dynamics simulation. The conformation of the ligand obtained from the major groove mode resembles that computed with the ligand soaked in water. The van der Waals contact energy was found to be as significant as electrostatic energy and more important for difference in binding energy between these two binding modes. The importance of van der Waals force in groove binding was supported by computations on the complex formed by the repressor peptide fragment from the bacteriophage 434 and its operator oligonucleotide.     

Iron-iron interaction through an ethanedithiolate ligand: a magnetic and theoretical study

The compound ([CpFe(dppe)]₂[μ-SCH₂CH₂SS,S′])(PF₆)₂ ([1][PF₆]₂) has been synthesized and its magnetic properties have been investigated by susceptometer quantum interface device (SQUID) measurements in the temperature range 5-300 K. The d⁵-d⁵1²⁺ complex exhibits intramolecular antiferromagnetic behavior, with a magnetic coupling constant of −6.4 cm⁻¹. Density functional theory (DFT) calculations on a model of 1²⁺ (as well as on models of 1⁺ and 1) allow to determine its molecular structure and analyse its bonding and magnetic properties. The computed spin density exhibits significant localization on both the Fe and S centers. Replacing the heteroatoms of the bridging ligand by CH₂ groups leads to a relocalization of the spin density on the metal atoms and favors ferromagnetic coupling.     

KcsA closed and open: modelling and simulation studies

Bacterial homologues of mammalian potassium channels provide structures of two states of a gated K channel. Thus, the crystal structure of KcsA represents a closed state whilst that of MthK represents an open state. Using homology modelling and molecular dynamics simulations we have built a model of the transmembrane domain of KcsA in an open state and have compared its conformational stability with that of the same domain of KcsA in a closed state. Approximate Born energy calculations of monovalent cations within the two KcsA channel states suggest that the intracellular hydrophobic gate in the closed state provides a barrier of height ~5 kT to ion permeation, whilst in the open state the barrier is absent. Simulations (10 ns duration) in an octane slab (a simple membrane mimetic) suggest that closed- and open-state models are of comparable conformational stability, both exhibiting conformational drifts of ~3.3 Å C RMSD relative to the respective starting models. Substantial conformational fluctuations are observed in the intracellular gate region during both simulations (closed state and open state). In the simulation of open-state KcsA, rapid (<5 ns) exit of all three K⁺ ions occurs through the intracellular mouth of the channel. Helix kink and swivel motion is observed at the molecular hinge formed by residue G99 of the M2 helix. This motion is more substantial for the open- than for the closed-state model of the channel.     

Long-range Electrostatic Interactions in Molecular Dynamics: An Endothelin-1 Case Study

Abstract

An extensive conformational search in explicit solvent was performed in order to compare the influence of different long-range electrostatic interaction treatments in molecular dynamics. The short peptide endothelin-1 was selected as the subject of molecular dynamics studies that started from both X-ray and NMR obtained structures. Electrostatic interactions were treated using two of the most common methods—residue-based cutoff and particle mesh Ewald (PME). Analyses of free energy calculations (MM-PBSA method used), secondary structure elements and hydrogen bonds were performed, and there suggested that there is no unambiguous conclusion about which of the two methods of long-range electrostatics treatment should be used in MD simulations in this case. The most reliable data was provided by a trajectory that started with the NMR structure and used the cutoff method to treat electrostatic interactions. This leads to a recommendation that the choice of electrostatics treatment should be made carefully and not automatically by choosing the PME method simply because it is the most widely used.     

Direct evidence of vinculin tail-lipid membrane interaction in beta-sheet conformation

The focal adhesion protein vinculin (1066 residues) plays an important role in cell adhesion and migration. The interaction between vinculin and lipid membranes is necessary to ensure these processes. There are three putative lipid-membrane interaction sites located at the vinculin tail domain two that form amphipathic alpha-helices (residues 935-978 and 1020-1040) and one that remains unstructured (residues 1052-1066) during crystallization. In this work, the structural and biochemical properties of the last 21 residues of the vinculin tail domain were investigated. Differential scanning calorimetry was performed in the presence of lipid vesicles consisting of dimyristoyl-l-α-phosphatidylcholine and dimyristoyl-l-α-phosphatidylglycerol at various molar ratios. The results demonstrate that this peptide inserts into lipid vesicle membranes. Examining the secondary structure of this peptide by molecular dynamics simulations and circular dichroism spectroscopy, we show that it adopts an antiparallel beta sheet backbone geometry that could ensure the association with lipid vesicles.     

DFT: a dynamic study of the interaction of ethanol and methanol with platinum

The interaction between alcohol molecules and platinum (Pt) was studied using molecular dynamics (MD; Born-Oppenheimer method). Alcohol molecules like ethanol and methanol present a similar molecular structure, with a methyl group (CH₃) at one end and a fragment of hydroxyl (OH) at the other. This fact generates two orientations that are considered in the interaction with Pt. The MD calculation results for these two orientations indicate a preferential orientation due to energy interactions. A plausible reaction mechanism that takes into account the interaction between Pt and alcohol is presented. The charge transference obtained from the Pt–alcohol interaction was also analyzed. The energy for the two orientations was calculated by indicating the preferential orientation. The methyl and hydroxyl groups are involved in heterolytic breakage of hydrogen bonds, joined to a carbon atom in the former and to an oxygen atom in the latter; however, the methyl group reaction seems to be the most important.     

Dissecting substrate specificity of two rice BADH isoforms: Enzyme kinetics, docking and molecular dynamics simulation studies

Fragrance rice (Oryza sativa) contains two isoforms of BADH, named OsBADH1 and OsBADH2. OsBADH1 is implicated in acetaldehyde oxidation in rice plant peroxisomes, while the non-functional OsBADH2 is believed to be involved in the accumulation of 2-acetyl-1-pyrroline, the major compound of aroma in fragrance rice. In the present study, site-directed mutagenesis, molecular docking and molecular dynamics simulation studies were used to investigate the substrate specificity towards Bet-ald and GAB-ald. Consistent with our previous study, kinetics data indicated that the enzymes catalyze the oxidation of GAB-ald more efficiently than Bet-ald and the OsBADH1 W172F and OsBADH2 W170F mutants displayed a higher catalytic efficiency towards GAB-ald. Molecular docking analysis and molecular dynamics simulations for the first time provided models for aldehyde substrate-bound complexes of OsBADHs. The amino acid residues, E262, L263, C296 and W461 of OsBADH1 and E260, L261, C294 and W459 of OsBADH2 located within 5 Å of the OsBADH active site mainly interacted with GAB-ald forming strong hydrogen bonds in both OsBADH isoforms. Residues W163, N164, Q294, C296 and F397 of OsBADH1–Bet-ald and Y163, M167, W170, E260, S295 and C453 of OsBADH2–Bet-ald formed the main interaction sites while E260 showed an interaction energy of −14.21 kcal/mol. Unconserved A290 in OsBADH1 and W288 in OsBADH2 appeared to be important for substrate recognition similar to that observed in PsAMADHs. Overall, the results here help to explain how two homologous rice BADHs recognize the aldehyde substrate differently, a key property to their biological role.     

Versatile strategy for the synthesis of biotin-labelled glycans,their immobilization to establish a bioactive surface and interaction studies with a lectin on a biochip

The emerging role of glycans as versatile biochemical signals in diverse aspects of cellular sociology calls for establishment of sensitive methods to monitor carbohydrate recognition by receptors such as lectins. Most of these techniques involve the immobilization of one of the binding partners on a surface, e.g. atomic force microscopy, glycan array and Surface Plasmon Resonance (SPR), hereby simulating cell surface presentation. Here, we report the synthesis of fluorescent glycoconjugates, with a functionalization strategy which avoids the frequently occurring ring opening at the reducing end for further immobilization on a surface or derivatization with biotin. In order to improve the versatility of these derivatized glycans for biological studies, a new approach for the synthesis of biotinylated and fluorescent glycans has also been realized. Finally, to illustrate their usefulness the neoglycoconjugates were immobilized on different surfaces, and the interaction analysis with a model lectin, the toxin from mistletoe, proved them to act as potent ligands, underscoring the merit of the presented synthetic approach. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Helix-capping interaction in λ cro protein: A free energy simulation analysis

The stability mutant Tyr-26 → Asp was studied in the Cro protein from bacteriophage λ using free energy molecular dynamics simulations. The mutant was calculated to be more stable than the wild type by 3.0 ± 1.7 kcal/mol/monomer, in reasonable agreement with experiment (1.4 kcal/mol/monomer). Moreover, the aspartic acid in the mutant was found to form a capping interation with the amino terminus of the third α-helix of Cro. The simulations were analyzed to understand better the source of the stability of this helix-capping interaction and to examine the results in light of previous explanations of stabilizing helix caps-namely, a model of local unsatisfied hydrogen bonds at the helix termini and the helix macro dipole model. Analysis of the simulations shows that the stabilizing effect of this charged helical cap is due both to favorable hydrogen bonds with backbone NH groups at the helix terminus and to favorable electrostatic interactions (but not hydrogen bonds) with their carbonyls (effectively the next row of local dipoles in the helix). However, electrostatic interactions are weak or negligible with backbone dipolar groups in the helix further away from the terminus. Moreover, the importance of other local electrostatic interactions with polar side chains near the helix terminus, which are neglected in most treatments of this effect, are shown to be important. Thus, the results support a model that is intermediate between the two previous explanations: both unsatisfied hydrogen bonds at the helix terminus and other, local preoriented dipolar groups stabilize the helix cap. These findings suggest that similar interactions with preoriented dipolar groups may be important for cooperativity in other charge–dipole interactions and may be employed to advantage for molecular design. © 1994 Wiley-Liss, Inc.     

Rabbit indolethylamine N-methyltransferase three-dimensional structure prediction: a model approach to bridge sequence to function in pharmacogenomic studies

Pharmacogenomics is the study of the genetic basis for individual variation in response to drugs and other xenobiotics. Successful prediction of effects of genetic variations that change encoded amino acid sequences on protein function and their consequent biomedical implications depends on three-dimensional (3D) structures of the encoded amino acid sequences. To bridge sequence to function, thus facilitating an in-depth pharmacogenomic study, we tested the feasibility of the use of a semi-computational approach to predict 3D structures of rabbit and human indolethylamine N-methyltransferases (INMTs) from their amino acid sequences, which share less than 26% sequence identity with known protein 3D structures. Herein, we report 3D models of INMTs predicted by using the crystal structure of rat catechol O-methyltransferase as a template, testing of the models both computationally and experimentally, and successful use of the models in retrospective prediction of the effects of genetic polymorphisms and in identification of residues that contribute to observed species-specific differences in substrate affinity. The results encourage the use of the semi-computational approach to predict 3D protein structures for use in pharmacogenomic studies when de novo prediction of protein 3D structures from their amino acid sequences is still not feasible and X-ray crystallography and/or solution nuclear magnetic resonance spectroscopy can only determine 3D structures for a small number of known amino acid sequences.Electronic Supplementary Material available.     

Effects of surfactant, salt and solvent on the structure and activity of adenosine deaminase: molecular dynamic and spectrophotometric studies

Effects of sodium dodecyl sulfate, dodecyltrimethylammonium bromide, sodium chloride, sodium sulfate, methanol and ethanol, on the structure and activity of adenosine deaminase (ADA) were investigated by UV-Vis, circular dichroism spectrophotometry and molecular dynamics (MDs) studies. Relative activity, experimental and computational helix content, total accessible surface area (ASA) and exposed charged surface area (ECSA) were obtained. The relative activity of ADA in the absence and the presence of denaturants were compared with structural results. It was shown that an increase in the surface area and a decrease in the amount of helicity are associated with a decrease in the activity of ADA.     

Low-mass molecular dynamics simulation for configurational sampling enhancement: More evidence and theoretical explanation

It has been reported recently that classical, isothermal–isobaric molecular dynamics (NTP MD) simulations at a time step of 1.00 fs of the standard-mass time (Δt=1.00 fs^smt) and a temperature of ≤340 K using uniformly reduced atomic masses by tenfold offers better configurational sampling than standard-mass NTP MD simulations at the same time step. However, it has long been reported that atomic masses can also be increased to improve configurational sampling because higher atomic masses permit the use of a longer time step. It is worth investigating whether standard-mass NTP MD simulations at Δt=2.00 or 3.16 fs^smt can offer better or comparable configurational sampling than low-mass NTP MD simulations at Δt=1.00 fs^smt. This article reports folding simulations of two β-hairpins showing that the configurational sampling efficiency of NTP MD simulations using atomic masses uniformly reduced by tenfold at Δt=1.00 fs^smt is statistically equivalent to and better than those using standard masses at Δt=3.16 and 2.00 fs^smt, respectively. The results confirm that, relative to those using standard masses at routine Δt=2.00 fs^smt, the low-mass NTP MD simulations at Δt=1.00 fs^smt are a simple and generic technique to enhance configurational sampling at temperatures of ≤340 K.     

New developments in first-principles excited-state dynamics simulations: unveiling the solvent specificity of excited anionic cluster relaxation and electron solvation

Charge-transfer-to-solvent excited iodide–polar solvent molecule clusters, [I^? (Solv)_n]^*, have attracted substantial interest over the past 20 years as they can undergo intriguing relaxation processes leading ultimately to the formation of gas-phase molecular analogues of the solvated electron. In this review article, we present a comprehensive overview of the development and application of state-of-the-art first-principles molecular dynamics simulation approaches to understand and interpret the results of femtosecond photoelectron spectroscopy experiments on [I^? (Solv)_n]^* relaxation, which point to a high degree of solvent specificity in the electron solvation dynamics. The intricate molecular details of the [I^? (Solv)_n]^* relaxation process are presented, and by contrasting the relaxation mechanisms of clusters with several different solvents (water, methanol and acetonitrile), the molecular basis of the solvent specificity of electron solvation in [I^? (Solv)_n]^* is uncovered, leading to a more refined view of the manifestation of electron solvation in small gas-phase clusters.     

Effects of interaction with sulphur compounds and free volume in imidazolium-based ionic liquid on desulphurisation: a molecular dynamics study

Interaction energy with sulphur compounds and free volume in imidazolium-based ionic liquid were calculated by molecular dynamics (MD) simulations to examine their effects on desulphurisation. From microstructure analysis and energy contribution calculation, it is found that an increasing fractional free volume in ionic liquid and an enhancement of interaction with solute by tuning the structure of ionic liquid or oxidising sulphur compounds are favourable for desulphurisation, which allows more efficient packing of sulphur compounds in ionic liquids and more easily extraction of sulphur compounds from fuel. The MD results are in good agreement with experimental desulphurisation performance.     

Hemoglobins of Kiefferulus,sister genus of Chironomus (Diptera: Insecta): Evolution of the Hb VIIB cluster

A genomic clone containing hemoglobin genes was isolated from a species of the chironomid genus Kiefferulus. Eight genes, including an apparent pseudogene, were sequenced and the amino acid sequences of the putative proteins were determined. By comparison to the previously described hemoglobins in the sister-genus Chironomus, they were identified as members of the dimeric Hb VIIB group. The results indicate that the existence of clusters of hemoglobin genes may be a common feature in chironomids and not just confined to Chironomus. The Kiefferulus genes show greatest similarity of amino acid sequence to Hb VIIB-7 from the Chironomus cluster. The results suggest that the ancestral cluster contained at least two gene types, one of which gave rise to VIIB-7 and the Kiefferulus genes while the other gave rise to the other Chironomus VIIB genes. Both clusters appear to have increased in size by duplication or unequal crossing over since the separation of the genera. It also appears that an unrelated gene present in the Chironomus cluster, Hb-Y, arose from a completely independent origin with no apparent equivalent gene anywhere in the genome of Kiefferulus or some other Chironomus species. Correspondence to: J. Martin