Polymerization of a vitelline envelope (VE)-like structure produced by using fractionated VE extracts from carp eggs

Fractionation of vitelline envelope (VE) extracts from carp eggs made possible the efficient polymerization of a VE-like structure. The structure corresponded to the fourth layer of the VE or fertilization envelope (FE), and its organization was achieved by reassembly in vitro after solubilization of the sheets composed of filamentous substances or network-like aggregates which were induced by a cortical alveolus sialoglycoprotein or thrombin. The sialoglycoprotein was a serine proteinase and immunolocalized only in the structure at the periphery of cortical alveoli, not in the VE and yolk granules. Ultrastructural features of the VE-like structure suggested that reassembly in vitro occurred via several intermediates in the process of polymerization. A polyclonal antibody produced against one of the assembled VE components, a 64 kDa protein, more intensely immunostained the outer periphery of the VEs than other areas, and immunoelectron microscopy showed that immunogold particles specifically labeled reassembled VE-like structures and major skeletons of the networks or network-like sheets. The protein with a molecular weight of 64 kDa was found to be a DNase. Thus, these results suggest a new approach to investigating not only the FE assembly process in vitro but also the organizing relationship between the major skeleton of the VE or FE and other additional constituents.     

Computational modeling of human coreceptor CCR5 antagonist as a HIV-1 entry inhibitor: using an integrated homology modeling,docking, and membrane molecular dynamics simulation analysis approach

Chemokine receptor 5 (CCR5) is an integral membrane protein that is utilized during human immunodeficiency virus type-1 entry into host cells. CCR5 is a G-protein coupled receptor that contains seven transmembrane (TM) helices. However, the crystal structure of CCR5 has not been reported. A homology model of CCR5 was developed based on the recently reported CXCR4 structure as template. Automated docking of the most potent (14), medium potent (37), and least potent (25) CCR5 antagonists was performed using the CCR5 model. To characterize the mechanism responsible for the interactions between ligands (14, 25, and 37) and CCR5, membrane molecular dynamic (MD) simulations were performed. The position and orientation of ligands (14, 25, and 37) were found to be changed after MD simulations, which demonstrated the ability of this technique to identify binding modes. Furthermore, at the end of simulation, it was found that residues identified by docking were changed and some new residues were introduced in the proximity of ligands. Our results are in line with the majority of previous mutational reports. These results show that hydrophobicity is the determining factor of CCR5 antagonism. In addition, salt bridging and hydrogen bond contacts between ligands (14, 25, and 37) and CCR5 are also crucial for inhibitory activity. The residues newly identified by MD simulation are Ser160, Phe166, Ser180, His181, and Trp190, and so far no site-directed mutagenesis studies have been reported. To determine the contributions made by these residues, additional mutational studies are suggested. We propose a general binding mode for these derivatives based on the MD simulation results of higher (14), medium (37), and lower (25) potent inhibitors. Interestingly, we found some trend for these inhibitors such as, salt bridge interaction between basic nitrogen of ligand and acidic Glu283 seemed necessary for inhibitory activity. Also, two aromatic pockets (pocket I – TM1-3 and pocket II – TM3-6) were linked by the central polar region in TM7, and the simulated inhibitors show important interactions with the Trp86, Tyr89, Tyr108, Phe112, Ile198, Tyr251, Leu255, and Gln280 and Glu283 residues. These results shed light on the usage of MD simulation to identify more stable, optimal binding modes of the inhibitors.     

Functional and direct interaction between the RNA binding protein HuD and active Akt1

The RNA binding protein HuD plays essential roles in neuronal development and plasticity. We have previously shown that HuD stimulates translation. Key for this enhancer function is the linker region and the poly(A) binding domain of HuD that are also critical for its function in neurite outgrowth. Here, we further explored the underlying molecular interactions and found that HuD but not the ubiquitously expressed HuR interacts directly with active Akt1. We identify that the linker region of HuD is required for this interaction. We also show by using chimeric mutants of HuD and HuR, which contain the reciprocal linker between RNA-binding domain 2 (RBD2) and RBD3, respectively, and by overexpressing a dominant negative mutant of Akt1 that the HuD-Akt1 interaction is functionally important, as it is required for the induction of neurite outgrowth in PC12 cells. These results suggest the model whereby RNA-bound HuD functions as an adapter to recruit Akt1 to trigger neurite outgrowth. These data might also help to explain how HuD enhances translation of mRNAs that encode proteins involved in neuronal development.     

Study on the binding mode of a pyrrolotriazin derivative with JAK2 by docking and MD simulation

The pyrrolotriazin derivative 2-(4-(4-((7-(3-(N-methylmethylsulfonamido)phenyl)pyrrolo [2,1-f][1,2,4]triazin-2-yl)amino)phenyl)piperidin-1-yl)acetamide (PPA) is a potential Janus kinase 2 (JAK2) inhibitor. The binding mode between PPA and JAK2 was investigated by using a combined method of docking, molecular dynamics (MD) simulation and binding free-energy calculation. The docking calculations preliminarily indicated that there were two possible binding modes 1 and 2; MD simulations and binding free-energy calculations identified that binding mode 1 was more stable and favourable, with the lower MM-PBSA binding free energy of ?34.00?±?0.17?kcal/mol. Moreover, some valuable binding information is revealed as follows: the inhibitor PPA is suitably located at the ATP-binding site of JAK2 and the hydrophobic interaction plays an essential role. PPA not only interacts with residues Leu855, Val863, Ala880, Tyr931, Leu932 and Leu983 via hydrophobic interaction but also interacts with Ser936 and Asp994 by hydrogen bonds. These two factors are advantageous for PPA to strongly bind to JAK2. These results help to understand the action mechanisms and designing new compounds with a higher affinity to JAK2.     

Molecular dynamics simulations reveal how the reticulon-homology domain of the autophagy receptor RETREG1/FAM134B remodels membranes for efficient selective reticulophagy

ABSTRACT

The autophagy receptor for selective reticulophagy, RETREG1/FAM134B is essential for ER maintenance, and its dysfunction is associated with neuronal disorders, vascular dementia, or viral infections. The protein consists of the reticulon-homology domain (RHD) that is flanked at the N- and C-termini by an intrinsically disordered protein region (IDPR), where the C terminal IDPR carries the indispensable LC3-interacting region (LIR) motif for the interaction with LC3. The RHD of RETREG1 is presumed to play a role in membrane remodeling, but the absence of a known 3D structure of this domain so far prevented researchers from gaining mechanistic insights into how the RETREG1 RHD curves membranes, and thereby facilities reticulophagy. The recent study by Bhaskara et al., which is described in this editor’s corner article, used molecular dynamics (MD) simulations to create a structural model of the RETREG1 RHD. MD simulations along with in vitro liposome remodeling experiments reveal how the RHD domain acts on the ER membrane and, in concert with the C terminal IDPR, executes the function of RETREG1 in selective reticulophagy.

Abbreviations: ER, endoplasmic reticulum; IDPR, intrinsically disordered protein region; LIR, LC3-interacting region; MD, molecular dynamics; RHD, reticulon-homology domain; TM, transmembrane     

Crystal contact-free conformation of an intrinsically flexible loop in protein crystal: Tim21 as the case study

BackgroundIn protein crystals, flexible loops are frequently deformed by crystal contacts, whereas in solution, the large motions result in the poor convergence of such flexible loops in NMR structure determinations. We need an experimental technique to characterize the structural and dynamic properties of intrinsically flexible loops of protein molecules.MethodsWe designed an intended crystal contact-free space (CCFS) in protein crystals, and arranged the flexible loop of interest in the CCFS. The yeast Tim 21 protein was chosen as the model protein, because one of the loops (loop 2) is distorted by crystal contacts in the conventional crystal.ResultsYeast Tim21 was fused to the MBP protein by a rigid α-helical linker. The space created between the two proteins was used as the CCFS. The linker length provides adjustable freedom to arrange loop 2 in the CCFS. We re-determined the NMR structure of yeast Tim21, and conducted MD simulations for comparison. Multidimensional scaling was used to visualize the conformational similarity of loop 2. We found that the crystal contact-free conformation of loop 2 is located close to the center of the ensembles of the loop 2 conformations in the NMR and MD structures.ConclusionsLoop 2 of yeast Tim21 in the CCFS adopts a representative, dominant conformation in solution.General significanceNo single powerful technique is available for the characterization of flexible structures in protein molecules. NMR analyses and MD simulations provide useful, but incomplete information. CCFS crystallography offers a third route to this goal.     

Synthesis of New Pyrrolo[1,2-a]quinoxaline Derivatives as Potential Inhibitors of Akt Kinase

Akt kinases are attractive targets for small molecule drug discovery because of their key role in tumor cell survival/proliferation and their overexpression/activation in many human cancers. Recent efforts in the development and biological evaluation of small molecule inhibitors of Akt have led to the identification of novel Akt kinase inhibitors, based on a quinoxaline or pyrazinone scaffold. A series of new substituted pyrrolo[1,2-a]quinoxaline derivatives, structural analogues of these active quinoxaline or pyrazinone pharmacophores, was synthesized from various substituted 2-nitroanilines or 1,2-phenylenediamine via multistep heterocyclization process. These new compounds were tested for their in vitro ability to inhibit the proliferation of the human leukemic cell lines K562, U937 and HL60, and the breast cancer cell line MCF7. Three of these human cell lines (K562, U937 and MCF7) exhibited an active phosphorylated Akt form. The most promising active pyrroloquinoxalines were found to be 1a that inhibited K562 cell line proliferation with an IC₅₀ of 4.5 μM, and 1h that inhibited U937 and MCF7 cell lines with IC₅₀ of 5 and 8 μM, respectively. These two candidates exhibited more potent activities than the reference inhibitor A6730.     

壳聚糖与细菌细胞膜相互作用的分子动力学研究

目的 壳聚糖（chitosan,CS）是一种天然的广谱抗菌活性物质。现有研究表明，壳聚糖与细菌细胞膜的相互作用是其发挥抗菌功能的关键。受限于传统实验技术的表征能力，壳聚糖与细菌细胞膜相互作用的具体机制仍有待研究。本文旨在研究壳聚糖与细菌细胞膜相互作用的分子机制。方法 本研究利用全原子分子动力学模拟技术主要探究了完全脱乙酰化的不同聚合度壳聚糖（八聚糖、十二聚糖和十六聚糖）与革兰氏阴性菌外膜（outer membrane,OM）和革兰氏阳性菌质膜（cytoplasmic membrane,CM）相互作用的动态过程。结果 壳聚糖主要依靠其氨基、碳6位羟基和碳3位羟基与OM和CM的头部极性区发生快速结合。随后壳聚糖末端糖基单元倾向于插入OM内部，深度约1 nm，并与脂质分子脂肪酸链上的羰基形成稳定的氢键相互作用。与之相比，壳聚糖分子难以稳定地插入CM内部。壳聚糖结合对膜结构性质产生影响，主要表现在降低OM和CM的单分子脂质面积，显著减少OM和CM极性区的Ca2+和Na+数目，破坏阳离子介导的脂质间相互作用。结论 本研究证明，壳聚糖带正电的氨基基团是介导其与膜相互作用的关键，并破环脂质间的相互作...     

A CHARMm Based Force Field for Carbohydrates Using the CHEAT Approach: Carbohydrate Hydroxyl Groups Represented by Extended Atoms

Abstract

Computational studies of carbohydrates that do not consider explicit solvent molecules suffer from the strong tendency of the carbohydrate pendant hydroxyl groups to form intramolecular hydrogen bonds that are unlikely to be present in protic media. In this paper a novel approach towards molecular modelling of carbohydrates is described. The average effect of intra- and intermolecular hydrogen bonding is introduced into the potential energy function by adding a new (extended) atom type representing a carbohydrate hydroxyl group to the CHARMm force field; we coin the name CHEAT (Carbohydrate Hydroxyls represented by Extended AToms) for the resulting force field. As a training set for the parametrisation of CHEAT we used ethylene glycol, 10 cyclohexanols, 5 inositols, and 12 glycopyranoses for which in total 64 conformational energy differences were estimated using a set of steric interaction energies between hydroxyl and/or methyl groups on six-membered ring compounds as derived by Angyal (Angew. Chem., 8, 172-182, (1969)). The root-mean-square deviation between the estimated energy differences and the corresponding values obtained by our CHEAT approach amounts to 0.37 kcal/mol (n = 64). The CHEAT approach, which is claimed to calculate aqueous state compatible energetical and conformational properties of carbohydrates, is computationally very efficient and facilitates the calculation of nanosecond range MD trajectories as well as systematic conformational searches of oligosaccharides.     

On the Occurrence of Three-Center Hydrogen Bonds in Cyclodextrins in Crystalline Form and in Aqueous Solution: Comparison of Neutron Diffraction and Molecular Dynamics Results

Abstract

Three-center (bifurcated) hydrogen bonds may play a role by serving as an intermediate state between different dynamically changing hydrogen bonding patterns. Hydrogen bonding configurations can be studied experimentally by neutron diffraction and theoretically by computer simulation techniques. Here, both methods are used to analyse the occurrence of three-center hydrogen bonds in crystals of cyclodextrins.

Almost all experimentally observed three-center hydrogen bonds in the crystal are reproduced in the molecular dynamics (MD) simulations, even as far as the detailed asymmetric geometry is concerned. On the basis of this result a MD simulation of cyclodextrin in aqueous solution is searched for the occurrence of three-center hydrogen bonds. Significant differences are found. In solution more different three-center hydrogen bonds per α-cyclodextrin molecule are observed than in the crystal but the population (existence as percent of the simulation period) of each three-center hydrogen bond is lower in solution than in crystal. These may indeed serve as intermediate states in the process of changing one hydrogen bonding pattern into another.     

Molecular dynamics simulation of proteins under high pressure: Structure,function and thermodynamics

BackgroundMolecular dynamics (MD) simulation is well-recognized as a powerful tool to investigate protein structure, function, and thermodynamics. MD simulation is also used to investigate high pressure effects on proteins. For conducting better MD simulation under high pressure, the main issues to be addressed are: (i) protein force fields and water models were originally developed to reproduce experimental properties obtained at ambient pressure; and (ii) the timescale to observe the pressure effect is often much longer than that of conventional MD simulations.Scope of reviewFirst, we describe recent developments in MD simulation methodologies for studying the high-pressure structure and dynamics of protein molecules. These developments include force fields for proteins and water molecules, and enhanced simulation techniques. Then, we summarize recent studies of MD simulations of proteins in water under high pressure.Major conclusionsRecent MD simulations of proteins in solution under pressure have reproduced various phenomena identified by experiments using high pressure, such as hydration, water penetration, conformational change, helix stabilization, and molecular stiffening.General significanceMD simulations demonstrate differences in the properties of proteins and water molecules between ambient and high-pressure conditions. Comparing the results obtained by MD calculations with those obtained experimentally could reveal the mechanism by which biological molecular machines work well in collaboration with water molecules.     

Effect of Calcium ion on synaptotagmin-like protein during pre-fusion of vesicle for exocytosis in blood-brain barrier

BackgroundCalcium signaling and membrane fusion play key roles in exocytosis of drug-containing vesicles through the blood-brain barrier (BBB). Identifying the role of synaptotagmin-like protein4-a (Slp4-a) in the presence of Ca²⁺ ions, at the pre-fusion stage of a vesicle with the basolateral membrane of endothelial cell, can reveal brain drug transportation across BBB.MethodsWe utilized molecular dynamics (MD) simulations with a coarse-grained PACE force field to investigate the behaviors of Slp4-a with vesicular and endothelial membranes at the pre-fusion stage of exocytosis since all-atom MD simulation or experiments are more time-consuming and expensive to capture these behaviors.ResultsThe Slp4-a pulls lipid membranes (vesicular and endothelial) into close proximity and disorganizes lipid arrangement at contact points, which are predictors for initiation of fusion. Our MD results also indicate that Slp4-a needs Ca²⁺ to bind with weakly-charged POPE lipids (phosphatidylethanolamine).ConclusionsSlp4-a is an important trigger for membrane fusion in BBB exocytosis. It binds to lipid membranes at multiple binding sites and triggers membrane disruption for fusion in calcium-dependent case.General significanceUnderstanding the prefusion process of the vesicle will help to design better drug delivery mechanisms to the brain through formidable BBB.     

Identifying potential GluN2B subunit containing N-Methyl-D-aspartate receptor inhibitors: an integrative in silico and molecular modeling approach

Abstract

N-methyl-D-aspartate receptors (NMDARs), a class of ligand-gated ion channels, are involved in non-selective cation transport across the membrane. These are contained in glutamatergic synapse and produce excitatory effects leading to synaptic plasticity and memory function. GluN1-GluN2B, a subtype of NMDAR(s), has significant role in neurodegeneration, amyloid β (Aβ) induced synaptic dysfunction and loss. Thus, targeting and inhibiting GluN1-GluN2B may be effective in the management of neurodegenerative diseases including Alzheimer’s disease. In the present study, ligand and structure-based approaches were tried to identify the inhibitors. The pharmacophore, developed from co-crystallised ifenprodil, afforded virtual hits, which were further subjected through drug likeliness and PAINS filters to remove interfering compounds. Further comprehensive docking studies, free energy calculations and ADMET studies resulted in two virtual leads. The leads, ZINC257261614 and ZINC95977857 displayed good docking scores of ?12.90 and ?12.20?Kcal/mol and free binding energies of ?60.83 and ?61.83?Kcal/mol, respectively. The compounds were having acceptable predicted ADMET profiles and were subjected to molecular dynamic (MD) studies. The MD simulation produced stable complexes of these ligands with GluN1-GluN2B subunit having protein and ligand RMSD in acceptable limit. Abbreviations AD Alzheimer's disease

ADME Absorption distribution metabolism and excretion

ATD Amino terminal domain

BBB Blood-brain barrier

CNS Central nervous system

CREB cAMP response element binding protein

CTD Carboxy-terminal domain

Glu Glutamate

GMQE Global model quality estimation

HTVS High throughput virtual screening

HIA Human intestinal absorption

LGA Lamarckian genetic algorithm

MD Molecular dynamics

MM-GBSA Molecular mechanics, the Generalised Born model for Solvent Accessibility

NMDAR N-methyl-D-aspartate receptors

PAINS Pan assay interference compounds

RMSD Root-mean square deviation

RMSF Root-mean-square fluctuation

SMARTS SMILES arbitrary target specification

SP standard precision

XP extra precision

Communicated by Ramaswamy H. Sarma     

Pharmacoinformatics analysis to identify inhibitors of Mtb-ASADH

Aspartate-semialdehyde dehydrogenase (ASADH; EC 1.2.1.11) is a key enzyme in the biosynthesis of essential amino acids in prokaryotes and fungi, inhibition of ASADH leads to the development of novel antitubercular agents. In the present work, a combined structure and ligand-based pharmacophore modeling, molecular docking, and molecular dynamics (MD) approaches were employed to identify potent inhibitors of mycobacterium tuberculosis (Mtb)-ASADH. The structure-based pharmacophore hypothesis consists of three hydrogen bond acceptor (HBA), two negatively ionizable, and one positively ionizable center, while ligand-based pharmacophore consists of additional one HBA and one hydrogen bond donor features. The validated pharmacophore models were used to screen the chemical databases (ZINC and NCI). The screened hits were subjected to ADME and toxicity filters, and subsequently to the molecular docking analysis. Best-docked 25 compounds carry the characteristics of highly electronegative functional groups (–COOH and –NO₂) on both sides and exhibited the H-bonding interactions with highly conserved residues Arg99, Arg249, and His256. For further validation of docking results, MD simulation studies were carried out on two representative compounds NSC51108 and ZINC04203124. Both the compounds remain bound to the key active residues of Mtb-ASADH during the MD simulations. These identified hits can be further used for lead optimization and in the design more potent inhibitors against Mtb-ASADH.     

A computational study of binding between 3-(4-fluorophenyl)-N-((4-fluorophenyl)sulphonyl)acrylamide and tubulin

3-(4-Fluorophenyl)-N-((4-fluorophenyl)sulphonyl)acrylamide (FFSA) is a potential tubulin polymerisation inhibitor. In this article, a theoretical study of the binding between FFSA and tubulin in colchicine site was carried out by molecular docking, molecular dynamics (MD) simulation and binding free energy calculations. The docking calculations preliminarily indicate that there are three possible binding modes 1, 2 and 3; MD simulations and binding free energy calculations identify that binding mode 2 is the most favourable, with the lowest binding free energy of ? 29.54 kcal/mol. Moreover, our valuable results for the binding are as follows: the inhibitor FFSA is suitably located at the colchicine site of tubulin, where it not only interacts with residues Leu248β, Lys254β, Leu255β, Lys352β, Met259β and Val181a by hydrophilic interaction, but also interacts with Val181α and Thr179α by hydrogen bond interaction. These two factors are both essential for FFSA strongly binding to tubulin. These theoretical results help understanding the action mechanism and designing new compounds with higher affinity to tubulin.     

Membrane raft redox signalosomes in endothelial cells

Abstract

Membrane rafts (MRs) are specialized microdomains in the cell membrane with an altered lipid composition. Upon various stimulations, MRs can be clustered to aggregate or recruit NADPH oxidase sub-units and related proteins to form MR redox signalosomes in the membrane of cells like vascular endothelial cells (ECs). Multiple protein complexes, like MR redox signalosomes, are now considered to play a crucial role in the regulation of cell function and in the development of different cell dysfunctions. To form such redox signalosomes, ceramide will be generated from the hydrolysis of sphingomyelin by lysosomal acid sphingomyelinase that has been translocated via lysosome fusion to the MR area. In this brief review, current information is provided to help understand the occurrence and function of MR redox signalosomes. This may increase enthusiasm of the scientific community for further studies on the molecular mechanisms and the functional significance of forming such MR redox signalosomes.     

Molecular Dynamics Study on Protein and it's Water Structure at High Pressure

Abstract

We have performed NPT molecular dynamics simulations (Langevin Piston Method) on two types of solvated proteins-‘denaturation-unfavorable’ protein (insulin) and ‘denaturation-favorable protein’ (ribonuclease A) at high pressure (from 1 bar up to 20 kbar). The method is based on the extended system formalism introduced by Andersen, where the deterministic equations of motion for the piston degree of freedom are replaced by Langevin equation. We report the structural changes of proteins (ribonuclease A and insulin) and water molecules through radius of gyration, solvent accessible surface area, hydrogen bond pattern, and the topology of water clusters connected by the hydrogen bonded circular network. The solvent accessibility of ribonuclease A is mainly decreased by hydrophilic residues rather than hydrophobic residues under high pressure. From the results of hydrogen bond analysis, we have found that α-helix is more stable than β-sheet under high pressure. In addition, from the analysis of the water cluster, we have observed that for ribonuclease A, 5-membered ring structure is more favorable than 6-membered ring at higher pressure. However, for insulin, the ratio of 5 to 6-ring is constant over the pressure ranges for which we have performed MD simulation. This indicates that the water structure around insulin does not change under high pressure.     

Network pharmacology and molecular docking analyses on Lianhua Qingwen capsule indicate Akt1 is a potential target to treat and prevent COVID‐19

ObjectivesCoronavirus disease 2019 (COVID‐19) is rapidly spreading worldwide. Lianhua Qingwen capsule (LQC) has shown therapeutic effects in patients with COVID‐19. This study is aimed to discover its molecular mechanism and provide potential drug targets.Materials and MethodsAn LQC target and COVID‐19–related gene set was established using the Traditional Chinese Medicine Systems Pharmacology database and seven disease‐gene databases. Gene ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis and protein‐protein interaction (PPI) network were performed to discover the potential mechanism. Molecular docking was performed to visualize the patterns of interactions between the effective molecule and targeted protein.ResultsA gene set of 65 genes was generated. We then constructed a compound‐target network that contained 234 nodes of active compounds and 916 edges of compound‐target pairs. The GO and KEGG indicated that LQC can act by regulating immune response, apoptosis and virus infection. PPI network and subnetworks identified nine hub genes. The molecular docking was conducted on the most significant gene Akt1, which is involved in lung injury, lung fibrogenesis and virus infection. Six active compounds of LQC can enter the active pocket of Akt1, namely beta‐carotene, kaempferol, luteolin, naringenin, quercetin and wogonin, thereby exerting potential therapeutic effects in COVID‐19.ConclusionsThe network pharmacological strategy integrates molecular docking to unravel the molecular mechanism of LQC. Akt1 is a promising drug target to reduce tissue damage and help eliminate virus infection.     

Fluorescent α-Anomeric 1,N(6)Etheno- Deoxyadenosine in DNA Duplexes. the α-εdA / dG Pair

Abstract

Structural properties of the fluorescent α-anomeric 1,N(6)ethenodeoxyadenosine residue placed in opposition to all four canonical deoxynucleotide units within 11-mer DNA duplexes have been studied. The duplex with α-εedA / dG pairing is most thermodynamically stable while the α-edA / dC one is the least stable. Fluorescence measurements confirm the thermodynamic data and indicate base-pair dependent stacking properties of α-edA within duplex structures. Results of molecular dynamics (MD) simulations in aqueous solution for the most stable duplex point to the presence of different conformational states of the α-1,N(6)etheno-deoxyadenosine residue, including formation of a hydrogen bonded pair with the dG and possible occurrence of severe kinking in the duplex.