Essential dynamics sampling study of adenylate kinase: comparison to citrate synthase and implication for the hinge and shear mechanisms of domain motions

Essential dynamics sampling simulations of the domain conformations of unliganded Escherichia coli adenylate kinase have been performed to determine whether the ligand-induced closed-domain conformation is accessible to the open unliganded enzyme. Adenylate kinase is a three- domain protein with a central CORE domain and twoflanking domains, the LID and the NMPbind domains. The sampling simulations were applied to the CORE and NMPbind domain pair and the CORE and LID domain pair separately. One aim is to compare the results to those of a similar study on the enzyme citrate synthase to determine whether a similar domain-locking mechanism operates in adenylate kinase. Although for adenylate kinase the simulations suggest that the closed-domain conformation of the unliganded enzyme is at a slightly higher free energy than the open for both domain pairs, the results are radically different to those found for citrate synthase. In adenylate kinase the targeted domain conformations could always be achieved, whereas this was not the case in citrate synthase due to an apparent free-energy barrier between the open and closed conformations. Adenylate kinase has been classified as a protein that undergoes closure through a hinge mechanism, whereas citrate synthase has been assigned to the shear mechanism. This was quantified here in terms of the change in the number of interdomain contacting atoms upon closure which showed a considerable increase in adenylate kinase. For citrate synthase this number remained largely the same, suggesting that the domain faces slide over each other during closure. This suggests that shear and hinge mechanisms of domain closure may relate to the existence or absence of an appreciable barrier to closure for the unliganded protein, as the latter can hinge comparatively freely, whereas the former must follow a more constrained path. In general though it appears a bias toward keeping the unliganded enzyme in the open-domain conformation may be a common feature of domain enzymes.     

Dynamics of Small Molecules in a Dense Polymer Matrix: Molecular Dynamics Studies

Abstract

The anomalous diffusion regime appearing in the self-diffusion of small molecules in bulk amorphous polymers has been extensively studied by molecular dynamics simulations. A rather long simulation of duration ～ 10 ^?8 s is performed on a polyethylene-like simple polymer model containing either oxygen molecules or helium atoms as a diffusant. Dynamic properties evaluated for these diffusants are the mean-square displacement, the van Hove self correlation function, and the self part of the density autocorrelation. It is first confirmed that the anomalous diffusion regime appears in a few hundred picoseconds for oxygen molecule, while the Einstein relation adopted beyond this regime results in the self-diffusion coefficient of the order of ～ 10^?5 cm²/s. This anomaly is still observed for helium that diffuses much faster than oxygen. In the anomalous diffusion regime, it is found that the correlation functions for the two diffusants show characteristic features and become essentially the same as time is scaled appropriately. These features allow the estimation of the two characteristic spatial scales which are probably dominated by the microstructure of the polymer matrix, namely, the cage size and the distance between adjacent cages. The time dependence of the mean-square displacements of the two diffusants can be well interpreted by these characteristic spatial scales as time is scaled with the self-diffusion coefficients. It is shown that the anomalous diffusion regime arises from the inhomogeneous microstructure of the polymer matrix.     

Protein Surface Dynamics: Interaction with Water and Small Solutes

Previous time resolved measurements had indicated that protons could propagate on the surface of a protein, or a membrane, by a special mechanism that enhances the shuttle of the proton towards a specific site [1]. It was proposed that a proper location of residues on the surface contributes to the proton shuttling function. In the present study, this notion was further investigated using molecular dynamics, with only the mobile charge replaced by Na⁺ and Cl⁻ ions. A molecular dynamics simulation of a small globular protein (the S6 of the bacterial ribosome) was carried out in the presence of explicit water molecules and four pairs of Na⁺ and Cl⁻ ions. A 10 ns simulation indicated that the ions and the protein's surface were in equilibrium, with rapid passage of the ions between the protein's surface and the bulk. Yet it was noted that, close to some domains, the ions extended their duration near the surface, suggesting that the local electrostatic potential prevented them from diffusing to the bulk. During the time frame in which the ions were detained next to the surface, they could rapidly shuttle between various attractor sites located under the electrostatic umbrella. Statistical analysis of molecular dynamics and electrostatic potential/entropy consideration indicated that the detainment state is an energetic compromise between attractive forces and entropy of dilution. The similarity between the motion of free ions next to a protein and the proton transfer on the protein's surface are discussed.     

Stabilization of the integrase‐DNA complex by Mg2+ ions and prediction of key residues for binding HIV‐1 integrase inhibitors

The HIV‐1 integrase is an attractive target for the therapeutics development against AIDS, as no host homologue of this protein has been identified. The integrase strand transfer inhibitors (INSTIs), including raltegravir, specifically target the second catalytic step of the integration process by binding to the DDE motif of the catalytic site and coordinating Mg²⁺ ions. Recent X‐ray crystallographic structures of the integrase/DNA complex from prototype foamy virus allowed to investigate the role of the different partners (integrase, DNA, Mg²⁺ ions, raltegravir) in the complex stability using molecular dynamics (MD) simulations. The presence of Mg²⁺ ions is found to be essential for the stability, whereas the simultaneous presence of raltegravir and Mg²⁺ ions has a destabilizing influence. A homology model of HIV‐1 integrase was built on the basis of the X‐ray crystallographic information, and protein marker residues for the ligand binding were detected by clustering the docking poses of known HIV‐1 integrase inhibitors on the model. Interestingly, we had already identified some of these residues to be involved in HIV‐1 resistance mutations and in the stabilization of the catalytic site during the MD simulations. Classification of protein conformations along MD simulations, as well as of ligand docking poses, was performed by using an original learning method, based on self‐organizing maps. This allows us to perform a more in‐depth investigation of the free‐energy basins populated by the complex in MD simulations on the one hand, and a straightforward classification of ligands according to their binding residues on the other hand. Proteins 2014; 82:466–478. © 2013 Wiley Periodicals, Inc.     

Interaction of RNA with Phospholipid Membranes

RNAs binding with liposomes under near-physiological conditions were obtained by molecular selection. Structural analysis showed that the RNAs could form complexes owing to complementary sequences located in loops. Oligomerization of the RNAs selected was experimentally confirmed. The results and published data testified that formation of high-molecular-weight complexes is a major mechanism increasing the RNA affinity for phospholipid membranes. The role of RNA–membrane interactions in early evolution is discussed in terms of the RNA world hypothesis.     

Interaction of an Amphiphilic Peptide with a Phospholipid Bilayer Surface by Molecular Dynamics Simulation Study

Abstract

Corticotropin-releasing factor (CRF) is the principal neuroregulator of adrenocorticotropic hormone (ACTH) secretion. Previous experiments have demonstrated that CRF binds avidly to the surface of single egg phosphatidylcholine vesicles and its amphiphilic secondary structure might play an important role in the function. In this study, the interaction of the residues 13–41 in human CRF with the surface of a DOPC bilayer was investigated by molecular dynamics (MD) simulation in order to understand the role of the membrane surface in the formation of the amphiphilic α helix as well as to determine the effects of the peptide on the lipid bilayer. The model used included 60 DOPC molecules, 1 helical peptide (CRF_13–41) on the bilayer surface, and explicit waters of solvation in the lipid polar head group regions, together with constant-volume periodic boundary conditions in three dimensions. The MD simulation was carried out for 510 ps. In addition, CRF_13–41, initially in a helical form, was simulated in vacuo as a control. The results indicate that while it was completely unstable in vacuo, the peptide helical form was generally maintained on the bilayer surface, but with distortions near the terminal ends. The peptide was confined to the bilayer headgroup/water region, similar to that reported from neutron diffraction measurement of tripeptides bound to the phosphatidylcholine bilayer surface (Ref 1). The amphiphilicity of the peptide matched that of the bilayer headgroup environment, with the hydrophilic side oriented toward water and the hydrophobic side making contact with the bilayer hydrocarbon core. These results support the hypothesis that the amphiphilic environment of a membrane surface is important in the induction of peptide amphiphilic α-helical secondary structure. Two major effects of the peptide on the lipids were found: the first CH₂ segment in the lipid chains was significantly disordered and the lipid headgroup distribution was broadened towards the water region.     

Study on the interactions between diketo-acid inhibitors and prototype foamy virus integrase-DNA complex via molecular docking and comparative molecular dynamics simulation methods

Human immunodeficiency virus type 1 (HIV-1) integrase (IN) is an important drug target for anti-acquired immune deficiency disease (AIDS) treatment and diketo-acid (DKA) inhibitors are potent and selective inhibitors of HIV-1 IN. Due to lack of three-dimensional structures including detail interactions between HIV-1 IN and its substrate viral DNA, the drug design and screening platform remains incompleteness and deficient. In addition, the action mechanism of DKA inhibitors with HIV-1 IN is not well understood. In view of the high homology between the structure of prototype foamy virus (PFV) IN and that of HIV-1 IN, we used PFV IN as a surrogate model for HIV-1 IN to investigate the inhibitory mechanism of raltegravir (RLV) and the binding modes with a series of DKA inhibitors. Firstly, molecular dynamics simulations of PFV IN, IN-RLV, IN-DNA, and IN-DNA-RLV systems were performed for 10?ns each. The interactions and inhibitory mechanism of RLV to PFV IN were explored through overall dynamics behaviors, catalytic loop conformation distribution, and hydrogen bond network analysis. The results show that the coordinated interactions of RLV with IN and viral DNA slightly reduce the flexibility of catalytic loop region of IN, and remarkably restrict the mobility of the CA end of viral DNA, which may lead to the partial loss of the inhibitory activity of IN. Then, we docked a series of DKA inhibitors into PFV IN-DNA receptor and obtained the IN-DNA-inhibitor complexes. The docking results between PFV IN-DNA and DKA inhibitors agree well with the corresponding complex of HIV-1 IN, which proves the dependability of PFV IN-DNA used for the anti-AIDS drug screening. Our study may help to make clear some theoretical questions and to design anti-AIDS drug based on the structure of IN.     

Investigating the accessibility of the closed domain conformation of citrate synthase using essential dynamics sampling

A molecular dynamics study of pig heart citrate synthase is presented that aims to directly address the question of whether, for this enzyme, the ligand-induced closed domain conformation is accessible to the open unliganded enzyme. The approach utilises the technique of essential dynamics sampling, which is used in two modes. In exploring mode, the enzyme is encouraged to explore domain conformations it might not normally sample in free molecular dynamics simulation. In targeting mode, the enzyme is encouraged to adopt the domain conformation of a target structure. Using both modes extensively, it has been found that when the enzyme is prepared from a crystallographic open-domain structure and is in the unliganded state, it is unable to adopt the crystallographic closed-domain conformation of the liganded enzyme. Likewise, when the enzyme is prepared from the crystallographic closed liganded conformation with the ligands removed, it is unable to adopt the crystallographic open domain conformation. Structural investigations point to a common structural difference that is the source of this energy barrier; namely, the shift of alpha-helix 328-341 along its own axis relative to the large domain. Without this shift, the domains are unable to close or open fully. The charged substrate, oxaloacetate, binds near the base of this helix in the large domain and the interaction of Arg329 at the base of the helix with oxaloacetate is one that is consistent with the shift of this helix in going from the crystallographic open to closed structure. Therefore, the results suggest that without the substrate the enzyme remains in a partially open conformation ready to receive the substrate. In this way, the efficiency of the enzyme should be increased over one that is closed part of the time, with its binding site inaccessible to the substrate.     

A Molecular Dynamics Study of the Stereoselective Interaction between an Enantiomeric Amino-Acid Ester and an l-Histidine-Containing Catalyst in the Bilayer Membrane

Molecular dynamics simulations were performed to investigate the stereoselective interaction between an enantiomeric amino acid substrate (N-acylated-l-phenylalanine-p-nitrophenyl ester) and an l-histidine-containing dipeptide catalyst (N-(N-benzyloxycarbonyl-l-leucyl)-l-histidine) in the bilayer of cationic surfactants (N,N-bisdodecyl-N,N-dimethylammonium chloride). We found that a catalyst–substrate complex, which had an interamide hydrogen bond, was formed spontaneously in vacuum at 500 K. This complex was found to be stable both in vacuum and in the bilayer membrane for 100 ps at 300 K. The distances between the hydrophobic side chains in the complex were consistent with experimental results. The interamide hydrogen bond was retained in the hydrophobic core of the membrane. These results suggest that the catalyst–substrate complex found in this work is relevant to the stereoselective hydrolysis of the l-enantiomer of the substrate.     

HIV-1整合酶核心区野生型和F185K突变型活性和溶解性的比较及分子动力学模拟分析

表达纯化了野生型(WT)及F185K突变型HIV-1整合酶核心区蛋白(IN_C),并对二者的溶解性和活性进行了比较．实验结果表明：F185K 突变后IN_C溶解性显著提高,活性有一定程度降低．对WT和F185K IN_C体系进行了1800 ps的分子动力学模拟．模拟结果表明：F185K IN_C功能loop区柔性和蛋白质整体运动性降低,使蛋白质活性降低,F185K突变后盐桥网络的变化驱动了IN_C局部构象改变,引起IN_C表面的部分疏水残基被包埋,亲水残基暴露,相对亲水溶剂可接近面积增大,同时,突变后IN_C与水之间形成氢键的数量增加,与水之间作用加强,以上变化使IN_C溶解性提高．分子动力学模拟与实验结果相吻合．为理解蛋白质溶解性和对蛋白质进行可溶性改造提供了一定的理论依据．     

Molecular Dynamics Simulations of some Small Organic Molecules

Abstract

Free energy differences between different conformers of D-ribofuranose, L-malic acid and meso-tartaric acid in solution were calculated using Molecular Dynamics simulations. In case of ribose the α → β transition was studied. For the acids attention was focussed on the transitions between the three possible staggered conformers with respect to the central C-C bond. In all cases a thermodynamic integration method was employed to evaluate the free energy difference. The use of an alternative technique, umbrella sampling, for ribose did not give promising results.

It was shown that one needs a fairly accurate picture of the accessible conformational space in case of flexible molecules like the ones considered here before one can determine meaningful free energy differences. Large hysteresis effects between forward and reverse simulated transitions were observed, but contrary to the general belief they are no direct measure of the accuracy of the calculated ΔG values. In all cases the ΔG values resulting from the simulations and from NMR experiments agree within the, considerable, error limits and for the different forms of D-ribose, L-malic acid and L-tartaric acid the relative order of their populations is also correctly reproduced.     

Bilayer Mixing,Fusion, and Lysis Following the Interaction of Populations of Cationic and Anionic Phospholipid Bilayer Vesicles

Cationic, O-alkylphosphatidylcholines, recently developed as DNA transfection agents, form bilayers indistinguishable from those of natural phospholipids and undergo fusion with anionic bilayers. Membrane merging (lipid mixing), contents release, and contents mixing between populations of positive vesicles containing O-ethylphosphatidylcholine (EDOPC) and negative vesicles containing dioleolylphosphatidylglycerol (DOPG) have been determined with standard fluorometric vesicle-population assays. Surface-charge densities were varied from zero to full charge. All interactions depended critically on surface-charge density, as expected from the adhesion-condensation mechanism. Membrane mixing ranged from zero to 100%, with significant mixing (>10 <70%) occurring between cationic vesicles that were fully charged and anionic vesicles that had fractional surface charges as low as 0.1. Such mixing with membranes as weakly charged as cell membranes should be relevant to transfection with cationic lipids. Unexpectedly, lipid mixing was higher at high than at low ionic strength when one lipid dispersion was prepared from EDOPC plus DOPG (in different proportions), especially when the other vesicles were of EDOPC; this may somehow be a consequence of the ability of the former mixture to assume non-lamellar phases. Leakage of aqueous contents was also a strong function of charge, with fully charged vesicles releasing essentially all of their contents less than 1 min after mixing. EDOPC was more active in this regard than was DOPG, which probably reflects stronger intermolecular interactions of DOPG. Fusion, as measured by contents mixing, exhibited maximal values of 10% at intermediate surface charge. Reduced fusion at higher charge is attributed to multiple vesicle interactions leading to rupture. The existence of previously published data on individual interactions of vesicles of the same composition made it possible for the first time to compare pairwise with population interactions, confirming the likelihood of population studies to overestimate rupture and hemifusion and underestimate true vesicle fusion.     

Interaction of Piscidin-1 with zwitterionic versus anionic membranes: a comparative molecular dynamics study

Plasma membrane of each micro-organism has a unique set of lipid composition as a consequence of the environmental adaptation or a response to exposure to antimicrobial peptides (AMPs) as antibiotic agents. Understanding the relationship between lipid composition and action of antimicrobial peptides or considering how different lipid bilayers respond to AMPs may help us design more effective peptide drugs in the future. In this contribution, we intend to elucidate how two currently used membrane models, namely palmitoyl-oleoyl-phosphtidylglycerol (POPG) and 1-palmitoyl-oleoyl-glycero-phosphocholine (POPC), respond to antimicrobial peptide Piscidin-1 (Pis-1).The computed density profile of the peptide as it moves from the bulk solvent toward the membrane core suggests that Pis-1 penetrates into the POPG bilayer less than the POPC membrane. Furthermore, we showed that the two model membranes used in this study have different behavior in the presence of Pis-1. Hence, we suggest that membrane composition could be an important factor in determining lytic ability of peptide drugs to kill a unique bacterial species.

An animated interactive 3D complement (I3DC) is available in Proteopedia at http://proteopedia.org/w/Journal:JBSD:37     

HIV-1整合酶四聚体结构模拟及其活性位点分析

HIV-1复制需要HIV-1整合酶将其环状DNA整合进宿主DNA中,这其中包括2个重要反应,即“3′-加工”和“链转移”,两者均由HIV-1整合酶催化完成.阻断其中的任一反应,都能达到抑制HIV-1复制的目的.因此,了解HIV-1整合酶的完整结构和聚合状态,对深入探讨其作用机理及设计新型抑制剂具有重要的指导作用.然而,迄今为止仅有HIV-1整合酶单独结构域的晶体结构可供参考,而其全酶晶体结构尚未获得解析.本研究利用分子模拟技术,通过蛋白质 蛋白质/DNA分子对接、动力学模拟等方法,构建了全长整合酶四聚体的结构模型、HIV-1 DNA与整合酶复合物的结构模型,进一步从理论上证实HIV-1整合酶是以四聚体形态发挥催化作用,明确“3′-加工”和“链转移”在HIV-1整合酶上的催化位点.同时,通过与作用机理相似的细菌转座子Tn5转座酶等的结构比对,推测HIV-1整合酶的核心结构域中应有第2个Mg^２＋存在,其位置螯合于Asp64与Glu152之间.在HIV-1整合酶结构研究的基础上,有望进一步设计出新的抗艾滋病药物.     

The conformational feasibility for the formation of reaching dimer in ASV and HIV integrase: a molecular dynamics study

Retroviral integrases are reported to form alternate dimer assemblies like the core–core dimer and reaching dimer. The core–core dimer is stabilized predominantly by an extensive interface between two catalytic core domains. The reaching dimer is stabilized by N-terminal domains that reach to form intermolecular interfaces with the other subunit’s core and C-terminal domains (CTD), as well as CTD–CTD interactions. In this study, molecular dynamics (MD), Brownian dynamics (BD) simulations, and free energy analyses, were performed to elucidate determinants for the stability of the reaching dimer forms of full-length Avian Sarcoma Virus (ASV) and Human Immunodeficiency Virus (HIV) IN, and to examine the role of the C-tails (the last ~16–18 residues at the C-termini) in their structural dynamics. The dynamics of an HIV reaching dimer derived from small angle X-ray scattering and protein crosslinking data, was compared with the dynamics of a core–core dimer model derived from combining the crystal structures of two-domain fragments. The results showed that the core domains in the ASV reaching dimer express free dynamics, whereas those in the HIV reaching dimer are highly stable. BD simulations suggest a higher rate of association for the HIV core–core dimer than the reaching dimer. The predicted stability of these dimers was therefore ranked in the following order: ASV reaching dimer < HIV reaching dimer < composite core–core dimer. Analyses of MD trajectories have suggested residues that are critical for intermolecular contacts in each reaching dimer. Tests of these predictions and insights gained from these analyses could reveal a potential pathway for the association and dissociation of full-length IN multimers.     

Structural Basis and Mechanism of Chiral Benzedrine Molecules Interacting With Third Dopamine Receptor

In order to investigate the chiral benzedrine molecules corresponding to their different characteristics in biochemical systems, we studied their interaction with D₃R using the docking method, molecular dynamic simulation, and quantum chemistry. The obtained results indicate that the active residues for R‐benzedrine (RAT) bound with D₃R are Ala132, Asp133, and Tyr55, while Asn57, Asp133, Asp168, Cys172, Gly54, Trp24, and Vall136 act as the active residues for S‐benzedrine (SAT). The different active pockets are observed for ART or SAT because they possess different active residues. The binding energies between RAT and SAT with D₃R were determined to be ?44.0 kJ.mol^?1 and ?71.2 kJ.mol^?1, respectively. These results demonstrate that SAT within the studied pocket of D₃R has a stronger capability of binding with D₃R, while it is more feasible for RAT to leave from the interior positions of D₃R. In addition, the results suggest that the D₃R protein can recognize chiral benzedrine molecules and influence their different addictive and pharmacological effects in biochemical systems. Chirality 28:674–685, 2016. © 2016 Wiley Periodicals, Inc.     

Probing the initial stage of aggregation of the Abeta(10-35)-protein: assessing the propensity for peptide dimerization

Characterization of the early stages of peptide aggregation is of fundamental importance in elucidating the mechanism of the formation of deposits associated with amyloid disease. The initial step in the pathway of aggregation of the Abeta-protein, whose monomeric NMR structure is known, was studied through the simulation of the structure and stability of the peptide dimer in aqueous solution. A protocol based on shape complementarity was used to generate an assortment of possible dimer structures. The structures generated based on shape complementarity were evaluated using rapidly computed estimates of the desolvation and electrostatic interaction energies to identify a putative stable dimer structure. The potential of mean force associated with the dimerization of the peptides in aqueous solution was computed for both the hydrophobic and the electrostatic driven forces using umbrella sampling and classical molecular dynamics simulation at constant temperature and pressure with explicit solvent and periodic boundary conditions. The comparison of the two free energy profiles suggests that the structure of the peptide dimer is determined by the favorable desolvation of the hydrophobic residues at the interface. Molecular dynamics trajectories originating from two putative dimer structures indicate that the peptide dimer is stabilized primarily through hydrophobic interactions, while the conformations of the peptide monomers undergo substantial structural reorganization in the dimerization process. The finding that the phi-dimer may constitute the ensemble of stable Abeta(10-35) dimer has important implications for fibril formation. In particular, the expulsion of water molecules at the interface might be a key event, just as in the oligomerization of Abeta(16-22) fragments. We conjecture that events prior to the nucleation process themselves might involve crossing free energy barriers which depend on the peptide-peptide and peptide-water interactions. Consistent with existing experimental studies, the peptides within the ensemble of aggregated states show no signs of formation of secondary structure.     

Interaction of alginate single-chain polyguluronate segments with mono- and divalent metal cations: a comparative molecular dynamics study

The interaction of a set of monovalent (Na⁺, K⁺) and divalent (Mg²⁺, Ca²⁺) metal cations with single-chain polyguluronate (periodic chain based on a dodecameric repeat unit, 2₁-helical conformation) is investigated using explicit-solvent molecular dynamics simulations (at 300 K and 1 bar). A total of 14 (neutralising) combinations of the different ions are considered (single type of cation or simultaneous presence of two types of cation, either in the presence or absence of chloride anions). The main observations are: (1) the chain conformation and intramolecular hydrogen bonding is insensitive to the counter-ion environment; (2) the binding of the cations is essentially non-specific for all ions considered (counter-ion atmosphere confined within a cylinder of high ionic density, but no well-defined binding sites); (3) the density and tightness of the distributions of the different cations within the counter-ion atmosphere follow the approximate sequence Ca²⁺>Mg²⁺>K⁺～Na⁺; (4) the solvent-separated binding of the cations to the carboxylate groups of the chain is frequent, and its occurrence follows the approximate sequence K⁺>Na⁺>Ca²⁺>Mg²⁺ (contact-binding events as well as the binding of a cation to multiple carboxylate groups are very infrequent); and (5) the counter-ion atmosphere typically leads to a complete screening of the chain charge within 1.0–1.2 nm of the chain axis and, for most systems, to a charge reversal at about 1.5 nm (i.e. the effective chain charge becomes positive at this distance and as high in magnitude as one-quarter of the bare chain charge, before slowly decreasing to zero). These findings agree well (in a qualitative sense) with available experimental data and predictions from simple analytical models, and provide further insight concerning the nature of alginate–cation interactions in aqueous solution.     

A comparison of the dynamic behavior of monomeric and dimeric insulin shows structural rearrangements in the active monomer

Molecular dynamics (MD) simulations (5-10ns in length) and normal mode analyses were performed for the monomer and dimer of native porcine insulin in aqueous solution; both starting structures were obtained from an insulin hexamer. Several simulations were done to confirm that the results obtained are meaningful. The insulin dimer is very stable during the simulation and remains very close to the starting X-ray structure; the RMS fluctuations calculated from the MD simulation agree with the experimental B-factors. Correlated motions were found within each of the two monomers; they can be explained by persistent non-bonded interactions and disulfide bridges. The correlated motions between residues B24 and B26 of the two monomers are due to non-bonded interactions between the side-chains and backbone atoms. For the isolated monomer in solution, the A chain and the helix of the B chain are found to be stable during 5ns and 10ns MD simulations. However, the N-terminal and the C-terminal parts of the B chain are very flexible. The C-terminal part of the B chain moves away from the X-ray conformation after 0.5-2.5ns and exposes the N-terminal residues of the A chain that are thought to be important for the binding of insulin to its receptor. Our results thus support the hypothesis that, when monomeric insulin is released from the hexamer (or the dimer in our study), the C-terminal end of the monomer (residues B25-B30) is rearranged to allow binding to the insulin receptor. The greater flexibility of the C-terminal part of the beta chain in the B24 (Phe-->Gly) mutant is in accord with the NMR results. The details of the backbone and side-chain motions are presented. The transition between the starting conformation and the more dynamic structure of the monomers is characterized by displacements of the backbone of Phe B25 and Tyr B26; of these, Phe B25 has been implicated in insulin activation.     

Stabilization of duplex DNA and RNA by dangling ends studied by free energy simulations

Single unpaired nucleotides at the end of double‐stranded nucleic acids, termed dangling ends, can contribute to duplex stability. Umbrella sampling free energy simulations of dangling cytosine and guanine nucleotides at the end of duplex and single stranded RNA and DNA molecules have been used to investigate the molecular origin of dangling end effects. In unrestraint simulations, the dangling end nucleotides stayed close to placements observed in experimental structures. Calculated free energy contributions associated with the presence of dangling nucleotides were in reasonable agreement with experiment predicting the general trend of a more stabilizing effect of purine vs. pyrimidine dangling ends. In addition, the calculations indicate a more significant stabilizing effect of dangling ends at the 5′‐end vs. 3′‐end in case of DNA and the opposite trend in case of RNA. Both electrostatic and van der Waals interactions contribute to the duplex stabilizing effect of dangling end nucleotides. The free energy simulation scheme could also be used to design dangling end nucleotides that result in enhanced duplex stabilization. © 2013 Wiley Periodicals, Inc. Biopolymers 101: 418–427, 2014.