Search

Background and Purpose The spider‐venom peptide PcTx1 is the most potent and selective inhibitor of acid‐sensing ion channel (ASIC) 1a. It has centrally acting analgesic activity and is neuroprotective in rodent models of ischaemic stroke. Understanding the molecular details of the PcTx1 : ASIC1a interaction should facilitate development of therapeutically useful ASIC1a modulators. Previously, we showed that several key pharmacophore residues of PcTx1 reside in a dynamic β‐hairpin loop; conclusions confirmed by recent crystal structures of the complex formed between PcTx1 and chicken ASIC1 (cASIC1). Numerous peptide : channel contacts were observed in these crystal structures, but it remains unclear which of these are functionally important.

The multidrug transporter P-glycoprotein (P-gp) is central to the development of multidrug resistance in cancer. While residues essential for transport and binding have been identified, the location, composition, and specificity of potential drug binding sites are uncertain. Here molecular dynamics simulations are used to calculate the free energy profile for the binding of morphine and nicardipine to P-gp. We show that morphine and nicardipine primarily interact with key residues implicated in binding and transport from mutational studies, binding at different but overlapping sites within the transmembrane pore. Their permeation pathways were distinct but involved overlapping sets of residues. The results indicate that the binding location and permeation pathways of morphine and nicardipine are not well separated and cannot be considered as unique. This has important implications for our understanding of substrate uptake and transport by P-gp. Our results are independent of the choice of starting structure and consistent with a range of experimental studies.

As of late 2014, the Protein Data Bank (PDB) contained over 105,000 structures, of which over 26,000 X-ray crystal structures contained at least one of >18,000 organic ligand molecules. The number ...

Cyclopropane fatty acids are widespread in bacteria. As their concentration increases on exposure to hostile environments, they have been proposed to protect membranes. Here, the effect of cyclopropane and unsaturated fatty acids, both in cis and trans configurations, on the packing, order, and fluidity of lipid bilayers is explored using molecular dynamics simulations. It is shown that cyclopropane fatty acids disrupt lipid packing, favor the occurrence of gauche defects in the chains, and increase the lipid lateral diffusion, suggesting that they enhance fluidity. At the same time, they generally induce a greater degree of order than unsaturated fatty acids of the same configuration and limit the rotation about the bonds surrounding the cyclopropane ring. This indicates that cyclopropane fatty acids may fulfill a dual function: stabilizing membranes against adverse conditions while simultaneously promoting their fluidity. Marked differences in the effect of cis- and trans-monocyclopropanated fatty acids were also observed, suggesting that they may play alternative roles in membranes.

A general method for parametrizing atomic interaction functions is presented. The method is based on an analysis of surfaces corresponding to the difference between calculated and target data as a function of alternative combinations of parameters (parameter space mapping). The consideration of surfaces in parameter space as opposed to local values or gradients leads to a better understanding of the relationships between the parameters being optimized and a given set of target data. This in turn enables for a range of target data from multiple molecules to be combined in a robust manner and for the optimal region of parameter space to be trivially identified. The effectiveness of the approach is illustrated by using the method to refine the chlorine 6-12 Lennard-Jones parameters against experimental solvation free enthalpies in water and hexane as well as the density and heat of vaporization of the liquid at atmospheric pressure for a set of 10 aromatic-chloro compounds simultaneously. Single-step perturbation is used to efficiently calculate solvation free enthalpies for a wide range of parameter combinations. The capacity of this approach to parametrize accurate and transferrable force fields is discussed.

Finding and enumerating common molecular substructures is an important task in cheminformatics, where small molecules are often modeled as molecular graphs. We introduce the problem of enumerating all maximal k-common molecular fragments of a pair of molecular graphs. A k-common fragment is a common connected induced subgraph that consists of a common core and a common k-neighborhood. It is thus a generalization of the NP-hard task to enumerate all maximal common connected induced subgraphs (MCCIS) of two graphs, which corresponds to the k = 0 case. We extend the MCCIS enumeration algorithm by Ina Koch and apply algorithm engineering techniques to solve practical instances fast for the general k > 0 case, which is relevant, for example, for automatically generating force field topologies for molecular dynamics (MD) simulations. We find that our methods achieve good performance on a real-world benchmark of all-against-all comparisons of 255 molecules. Our software is available under the LGPL open source license at https://github.com/enitram/mogli .

Atomistic nonequilibrium molecular dynamics simulations have been used to model the induction of molecular orientation anisotropy within the emission layer of an organic light-emitting diode (OLED) formed by vapor deposition. Two emitter species were compared: racemic fac-tris(2-phenylpyridine)iridium(III) (Ir(ppy)3) and trans-bis(2-phenylpyridine)(acetylacetonate)iridium(III) (Ir(ppy)2(acac)). The simulations show that the molecular symmetry axes of both emitters preferentially align perpendicular to the surface during deposition. The molecular arrangement formed on deposition combined with consideration of the transition dipole moments provides insight into experimental reports that Ir(ppy)3 shows isotropic emission, while Ir(ppy)2(acac) displays improved efficiency due to an apparent preferential alignment of the transition dipole vectors parallel to the substrate. The simulations indicate that this difference is not due to differences in the extent of emitter alignment, but rather differences in the d...

Computer simulation of molecular systems enables structure-energy-function relationships of molecular processes to be described at the sub-atomic, atomic, supra-atomic, or supra-molecular level. To interpret results of such simulations appropriately, the quality of the calculated properties must be evaluated. This depends on the way the simulations are performed and on the way they are validated by comparison to values Qexp of experimentally observable quantities Q. One must consider 1) the accuracy of Qexp , 2) the accuracy of the function Q(rN ) used to calculate a Q-value based on a molecular configuration rN of N particles, 3) the sensitivity of the function Q(rN ) to the configuration rN , 4) the relative time scales of the simulation and experiment, 5) the degree to which the calculated and experimental properties are equivalent, and 6) the degree to which the system simulated matches the experimental conditions. Experimental data is limited in scope and generally corresponds to averages over both time and space. A critical analysis of the various factors influencing the apparent degree of (dis)agreement between simulations and experiment is presented and illustrated using examples from the literature. What can be done to enhance the validation of molecular simulation is also discussed.

Sterols and hopanoids have been suggested to reinforce membranes and protect against unfavorable environmental conditions. In particular, hopanoids are found in high concentrations in membranes of thermotolerant and thermophilic bacteria. However, the mechanism whereby sterols and hopanoids stabilize membranes at elevated temperatures is poorly understood. Here, the effect of temperature on the ordering of lipids in bilayers containing cholesterol or the hopanoids bacteriohopanetetrol and diplopterol was explored using molecular dynamics simulations. It is shown that cholesterol induces a high level of ordering over a wide range of temperatures. Bacteriohopanetetrol promotes order within the lipid tails but enhances fluid-like properties of the head groups at high temperatures. In contrast, diplopterol partitions in the midplane of the bilayer. This suggests that individual hopanoids fulfill distinct functions in membranes, with the ordering properties of bacteriohopanetetrol being particularly well suited to maintain the integrity of membranes at temperatures preferred by thermotolerant and thermophilic bacteria.

The peptides AM1 and Lac21E self-organize into switchable films at an air-water interface. In an earlier study, it was proposed that both AM1 and Lac21E formed monolayers of alpha-helical peptides based on consistency with neutron reflectivity data. In this article, molecular dynamics simulations of assemblies of helical and nonhelical AM1 and Lac21E at an air-water interface suggest some tendency for the peptides to spontaneously adopt an alpha-helical conformation. However, irrespective of the structure of the peptides, the simulations reproduced not only the structural properties of the films (thickness and distribution of the hydrophobic and hydrophilic amino acids) but also the experimental neutron reflectivity measurements at different contrast variations. This suggests that neutron reflectometry alone cannot be used to determine the structure of the peptides in this case. However, together with molecular dynamics simulations, it is possible to obtain a detailed understanding of peptide films at an atomic level.

Phox-homology (PX) domains target proteins to the organelles of the secretary and endocytic systems by binding to phosphatidylinositol phospholipids (PIPs). Among all the structures of PX domains that have been solved, only three have been solved in a complex with the main physiological ligand: PtdIns3P. In this work, molecular dynamic simulations have been used to explore the structure and dynamics of the p40(phox)-PX domain and the SNX17-PX domain and their interaction with membrane-bound PtdIns3P. In the simulations, both PX domains associated spontaneously with the membrane-bound PtdIns3P and formed stable complexes. The interaction between the p40(phox)-PX domain and PtdIns3P in the membrane was found to be similar to the crystal structure of the p40(phox)-PX-PtdIns3P complex that is available. The interaction between the SNX17-PX domain and PtdIns3P was similar to that observed in the p40(phox)-PX-PtdIns3P complex; however, some residues adopted different orientations. The simulations also showed that nonspecific interactions between the beta 1-beta 2 loop and the membrane play an important role in the interaction of membrane bound PtdIns3P and different PX domains. The behaviour of unbound PtdIns3P within a 2-oleoyl-1-palmitoyl-sn-glycero-3-phosphocholine (POPC) membrane environment was also examined and compared to the available experimental data and simulation studies of related molecules. (C) 2014 Wiley Periodicals, Inc.

Phosphorescent emissive materials in organic light-emitting diodes (OLEDs) manufactured using evaporation are usually blended with host materials at a concentration of 3–15 wt% to avoid concentration quenching of the luminescence. Here, experimental measurements of hole mobility and photoluminescence are related to the atomic level morphology of films created using atomistic nonequilibrium molecular dynamics simulations mimicking the evaporation process with similar guest concentrations as those used in operational test devices. For blends of fac-tris[2-phenylpyridinato-C2,N]iridium(III) [Ir(ppy)] in tris(4-carbazoyl-9-ylphenyl)amine (TCTA), it is found that clustering of the Ir(ppy) (surface of the molecules within ≈0.4 nm) in the simulated films is directly relatable to the experimentally-measured hole mobility. Films containing 1–10 wt% of Ir(ppy) in TCTA have a mobility of up to two orders of magnitude lower (≈10 cm V s) than the neat TCTA film, which is consistent with the Ir(ppy) molecules acting as hole traps due to their smaller ionization potential. Comparison of the simulated film morphologies with the measured photoluminescence properties shows that for luminescence quenching to occur, the Ir(ppy) molecules have to have their ligands partially overlapping. Thus, the results show that the effect of guest interactions on charge transport and luminescence are markedly different for OLED light-emitting layers.

Molecular simulation techniques are increasingly being used to study biomolecular systems at an atomic level. Such simulations rely on empirical force fields to represent the intermolecular interactions. There are many different force fields available—each based on a different set of assumptions and thus requiring different parametrization procedures. Recently, efforts have been made to fully automate the assignment of force-field parameters, including atomic partial charges, for novel molecules. In this work, we focus on a problem arising in the automated parametrization of molecules for use in combination with the gromos family of force fields: namely, the assignment of atoms to charge groups such that for every charge group the sum of the partial charges is ideally equal to its formal charge. In addition, charge groups are required to have size at most k. We show \documentclass{aastex}\usepackage{amsbsy}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{bm}\usepackage{mathrsfs}\usepackage{pifont}\usepackage{stmaryrd}\usepackage{textcomp}\usepackage{portland, xspace}\usepackage{amsmath, amsxtra}\pagestyle{empty}\DeclareMathSizes{10}{9}{7}{6}\begin{document}$${ \cal N P}$$\end{document}-hardness and give an exact algorithm that solves practical problem instances to provable optimality in a fraction of a second.

The antibiotic vancomycin targets lipid II, blocking cell wall synthesis in Gram-positive bacteria. Despite extensive study, questions remain regarding how it recognizes its primary ligand and what is the most biologically relevant form of vancomycin. In this study, molecular dynamics simulation techniques have been used to examine the process of ligand binding and dimerization of vancomycin. Starting from one or more vancomycin monomers in solution, together with different peptide ligands derived from lipid II, the simulations predict the structures of the ligated monomeric and dimeric complexes to within 0.1 nm rmsd of the structures determined experimentally. The simulations reproduce the conformation transitions observed by NMR and suggest that proposed differences between the crystal structure and the solution structure are an artifact of the way the NMR data has been interpreted in terms of a structural model. The spontaneous formation of both back-to-back and face-to-face dimers was observed in the simulations. This has allowed a detailed analysis of the origin of the cooperatively between ligand binding and dimerization and suggests that the formation of face-to-face dimers could be functionally significant. The work also highlights the possible role of structural water in stabilizing the vancomycin ligand complex and its role in the manifestation of vancomycin resistance.

During the past half century, the number and accuracy of experimental techniques that can deliver values of observables for biomolecular systems have been steadily increasing. The conversion of a measured value Qexp of an observable quantity Q into structural information is, however, a task beset with theoretical and practical problems: (i) insufficient or inaccurate values of Qexp, (ii) inaccuracies in the function Q(~r) used to relate the quantity Q to structure ~r, (iii) how to account for the averaging inherent in the measurement of Qexp, (iv) how to handle the possible multiple-valuedness of the inverse ~r(Q) of the function Q(~r), to mention a few. These apply to a variety of observable quantities Q and measurement techniques such as X-ray and neutron diffraction, small-angle and wide-angle X-ray scattering, free-electron laser imaging, cryo-electron microscopy, nuclear magnetic resonance, electron paramagnetic resonance, infrared and Raman spectroscopy, circular dichroism, Förster resonance energy transfer, atomic force microscopy and ion-mobility mass spectrometry. The process of deriving structural information from measured data is reviewed with an eye to non-experts and newcomers in the field using examples from the literature of the effect of the various choices and approximations involved in the process. A list of choices to be avoided is provided.

The effect of varying the emitter concentration on the structural properties of an archetypal phosphorescent blend consisting of 4,4′-bis(N-carbazolyl)biphenyl and tris(2-phenylpyridyl)iridium(III) has been investigated using non-equilibrium molecular dynamics (MD) simulations that mimic the process of vacuum deposition. By comparison with reflectometry measurements, we show that the simulations provide an accurate model of the average density of such films. The emitter molecules were found not to be evenly distributed throughout film, but rather they can form networks that provide charge and/or energy migration pathways, even at emitter concentrations as low as ≈5 weight percent. At slightly higher concentrations, percolated networks form that span the entire system. While such networks would give improved charge transport, they could also lead to more non-radiative pathways for the emissive state and a resultant loss of efficiency.

Plants use intracellular immunity receptors, known as nucleotide-binding oligomerization domain-like receptors (NLRs), to recognize specific pathogen effector proteins and induce immune responses. These proteins provide resistance to many of the world’s most destructive plant pathogens, yet we have a limited understanding of the molecular mechanisms that lead to defense signaling. We examined the wheat NLR protein, Sr33, which is responsible for strain-specific resistance to the wheat stem rust pathogen, Puccinia graminis f. sp. tritici. We present the solution structure of a coiled-coil (CC) fragment from Sr33, which adopts a four-helix bundle conformation. Unexpectedly, this structure differs from the published dimeric crystal structure of the equivalent region from the orthologous barley powdery mildew resistance protein, MLA10, but is similar to the structure of the distantly related potato NLR protein, Rx. We demonstrate that these regions are, in fact, largely monomeric and adopt similar folds in solution in all three proteins, suggesting that the CC domains from plant NLRs adopt a conserved fold. However, larger C-terminal fragments of Sr33 and MLA10 can self-associate both in vitro and in planta, and this self-association correlates with their cell death signaling activity. The minimal region of the CC domain required for both cell death signaling and self-association extends to amino acid 142, thus including 22 residues absent from previous biochemical and structural protein studies. These data suggest that self-association of the minimal CC domain is necessary for signaling but is likely to involve a different structural basis than previously suggested by the MLA10 crystallographic dimer.

Crystallographic studies of ligands bound to biological macromolecules (proteins and nucleic acids) represent\ud an important source of information concerning drug-target interactions, providing atomic level insights\ud into the physical chemistry of complex formation between macromolecules and ligands. Of the\ud more than 115,000 entries extant in the Protein Data Bank (PDB) archive, ~75% include at least one non-polymeric\ud ligand. Ligand geometrical and stereochemical quality, the suitability of ligand models for in silico drug\ud discovery and design, and the goodness-of-fit of ligand models to electron-density maps vary widely across\ud the archive. We describe the proceedings and conclusions from the first Worldwide PDB/Cambridge Crystallographic\ud Data Center/Drug Design Data Resource (wwPDB/CCDC/D3R) Ligand Validation Workshop\ud held at the Research Collaboratory for Structural Bioinformatics at Rutgers University on July 30–31, 2015.\ud Experts in protein crystallography from academe and industry came together with non-profit and for-profit\ud software providers for crystallography and with experts in computational chemistry and data archiving to\ud discuss and make recommendations on best practices, as framed by a series of questions central to structural\ud studies of macromolecule-ligand complexes. What data concerning bound ligands should be archived\ud in the PDB? How should the ligands be best represented? How should structural models of macromoleculeligand\ud complexes be validated? What supplementary information should accompany publications of structural\ud studies of biological macromolecules? Consensus recommendations on best practices developed in\ud response to each of these questions are provided, together with some details regarding implementation.\ud Important issues addressed but not resolved at the workshop are also enumerated.

The Type-I Cytokine Receptors regulate diverse cellular functions that encompass immune responses, growth, metabolism, lactation and red blood cell production. A sub-group of receptors that are models for receptor activation includes the growth hormone receptor (GHR), prolactin receptor (PRLR) and erythropoietin receptor (EPOR). While the cellular effects of receptor activation have been well characterised, the atomic details of how these receptors mechanically couple cytokine binding to the extracellular domains through the plasma membrane has remained elusive. Using fully atomistic molecular dynamics simulations, we have shown the extracellular domains of the GHR and PRLR homodimers undergo a relative rigid body rotation upon activation. Recently, we have tested if the rotation-based mechanism could be extended to erythropoietin (EPO) induced activation of the EPOR homodimer along with agonistic and antagonistic EPO mimetic peptides. In the simulations, the EPOR protomers rotate as rigid bodies similarly to the GHR and PRLR upon cytokine activation. However, analysis of the mimetic-bound EPOR could not distinguish between agonised or antagonised conformations of the dimer. Simulations of the GHR, PRLR and EPOR embedded in membranes exhibited different behavior in cholesterol rich membranes, suggesting a key role of cholesterol in receptor activation. The simulations demonstrate a shared cytokine induced rotational-based mechanism of activation of Type-I Cytokine Receptors. The atomic resolution of the mechanical coupling of the receptors through the plasma membrane provides clearer understanding of how the signal of cytokine binding is transmitted from the extracellular to the cytosolic domains.

The growth hormone receptor (GHR) has been the archetype for the class I cytokine receptor family. The original paradigm for how growth hormone (GH) binding to the GHR led to activation of the associated JAK2 was based on biophysical and structural studies over 20 years ago which showed that growth hormone binding to the extracellular domain of the receptor occurred in a sequential fashion, first binding to a high affinity site on a single receptor (site 1) and then binding to a lower affinity site (site 2) on a second receptor. We propose a new model for JAK2 activation by the growth hormone receptor. We utilised multiple approaches to define this model including FRET to monitor positioning of the JAK2 binding motif, leucine-zipper receptor constructs to control receptor transmembrane helix position, atomistic modelling of TM helix interactions and docking of JAK2 kinase and inhibitory pseudokinase crystal structures with an opposing pair in trans. Surprisingly, we identified that receptor activation induces a separation of its JAK2 binding motifs which leads to removal of the pseudokinase domain from the kinase domain of the partner JAK2 and pairing of the two kinase domains, facilitating trans-activation. In addition, we have recently shown that the receptor intracellular domain is intrinsically disordered and has specific interactions with the lipids of the inner membrane. This model of receptor activation may represent a common mechanism to other class I cytokine receptors.

Free to read at publisher's site. ProTx-II is a disulfide-rich peptide toxin from tarantula venom able to inhibit the human voltage-gated sodium channel 1.7 (hNaV1.7), a channel reported to be involved in nociception, and thus it might have potential as a pain therapeutic. ProTx-II acts by binding to the membrane-embedded voltage sensor domain of hNaV1.7, but the precise peptide channel-binding site and the importance of membrane binding on the inhibitory activity of ProTx-II remain unknown. In this study, we examined the structure and membrane-binding properties of ProTx-II and several analogues using NMR spectroscopy, surface plasmon resonance, fluorescence spectroscopy, and molecular dynamics simulations. Our results show a direct correlation between ProTx-II membrane binding affinity and its potency as an hNaV1.7 channel inhibitor. The data support a model whereby a hydrophobic patch on the ProTx-II surface anchors the molecule at the cell surface in a position that optimizes interaction of the peptide with the binding site on the voltage sensor domain. This is the first study to demonstrate that binding of ProTx-II to the lipid membrane is directly linked to its potency as an hNaV1.7 channel inhibitor.

The sensitivity of the structure and dynamics of a fully hydrated pure bilayer of 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) in molecular dynamics simulations to changes in force-field and simulation parameters has been assessed. Three related force fields (the Gromos 54A7 force field, a Gromos 53A6-derived parameter set and a variant of the Berger parameters) in combination with either particle-mesh Ewald (PME) or a reaction field (RF) were compared. Structural properties such as the area per lipid, carbon-deuterium order parameters, electron density profile and bilayer thicknesses, are reproduced by all the parameter sets within the uncertainty of the available experimental data. However, there are clear differences in the ordering of the glycerol backbone and choline headgroup, and the orientation of the headgroup dipole. In some cases, the degree of ordering was reminiscent of a liquid-ordered phase. It is also shown that, although the lateral diffusion of the lipids in the plane of the bilayer is often used to validate lipid force fields, because of the uncertainty in the experimental measurements and the fact that the lateral diffusion is dependent on the choice of the simulation conditions, it should not be employed as a measure of quality. Finally, the simulations show that the effect of small changes in force-field parameters on the structure and dynamics of a bilayer is more significant than the treatment of the long-range electrostatic interactions using RF or PME. Overall, the Gromos 54A7 best reproduced the range of experimental data examined.

The efficiency of using a variant of Hamiltonian replica-exchange molecular dynamics (Chaperone H-replica-exchange molecular dynamics [CH-REMD]) for the refinement of protein structural models generated de novo is investigated. In CH-REMD, the interaction between the protein and its environment, specifically, the electrostatic interaction between the protein and the solvating water, is varied leading to cycles of partial unfolding and refolding mimicking some aspects of folding chaperones. In 10 of the 15 cases examined, the CH-REMD approach sampled structures in which the root-mean-square deviation (RMSD) of secondary structure elements (SSE-RMSD) with respect to the experimental structure was more than 1.0 A lower than the initial de novo model. In 14 of the 15 cases, the improvement was more than 0.5 A. The ability of three different statistical potentials to identify near-native conformations was also examined. Little correlation between the SSE-RMSD of the sampled structures with respect to the experimental structure and any of the scoring functions tested was found. The most effective scoring function tested was the DFIRE potential. Using the DFIRE potential, the SSE-RMSD of the best scoring structures was on average 0.3 A lower than the initial model. Overall the work demonstrates that targeted enhanced-sampling techniques such as CH-REMD can lead to the systematic refinement of protein structural models generated de novo but that improved potentials for the identification of near-native structures are still needed.

The aim of fragment-based drug design (FBDD) is to identify molecular fragments that bind to alternate subsites within a given binding pocket leading to cooperative binding when linked. In this study, the binding of fragments to human phenylethanolamine N-methyltransferase is used to illustrate how (a) current protocols may fail to detect fragments that bind cooperatively, (b) theoretical approaches can be used to validate potential hits, and (c) apparent false positives obtained when screening against cocktails of fragments may in fact indicate promising leads.

Acid-sensing ion channel 1a (ASIC1a) is a primary acid sensor in the peripheral and central nervous system. It has been implicated as a novel therapeutic target for a broad range of pathophysiological conditions including pain, ischemic stroke, depression, and autoimmune diseases such as multiple sclerosis. The only known selective blocker of ASIC1a is π-TRTX-Pc1a (PcTx1), a disulfide-rich 40-residue peptide isolated from spider venom. π-TRTX-Pc1a is an effective analgesic in rodent models of acute pain and it provides neuroprotection in a mouse model of ischemic stroke. Thus, understanding the molecular basis of the π-TRTX-Pc1a-ASIC1a interaction should facilitate development of therapeutically useful ASIC1a blockers. We therefore developed an efficient bacterial expression system to produce a panel of π-TRTX-Pc1a mutants for probing structure-activity relationships as well as isotopically labeled toxin for determination of its solution structure and dynamics. We demonstrate that the toxin pharmacophore resides in a β-hairpin loop that was revealed to be mobile over a wide range of time scales using molecular dynamics simulations in combination with NMR spin relaxation and relaxation dispersion measurements. The toxin-receptor interaction was modeled by in silico docking of the toxin structure onto a homology model of rat ASIC1a in a restraints-driven approach that was designed to take account of the dynamics of the toxin pharmacophore and the consequent remodeling of side-chain conformations upon receptor binding. The resulting model reveals new insights into the mechanism of action of π-TRTX-Pc1a and provides an experimentally validated template for the rational design of therapeutically useful π-TRTX-Pc1a mimetics.

A fundamental challenge in computational drug design is the availability of reliable and validated experimental binding and structural data against which theoretical calculations can be compared. In this work a combination of molecular dynamics (MD) simulations and free energy calculations has been used to analyze the structural and thermodynamic basis of ligand recognition by phenylethanolamine N-methyltransferase (PNMT) in an attempt to resolve uncertainties in the available binding and structural data. PNMT catalyzes the conversion of norepinephrine into epinephrine (adrenaline), and inhibitors of PNMT are of potential therapeutic importance in Alzheimer's and Parkinson's disease. Excellent agreement between the calculated and recently revised relative binding free energies to human PNMT was obtained with the average deviation between the calculated and the experimentally determined values being only 0.8 kJ/mol. In this case, the variation in the experimental data over time is much greater than the uncertainties in the theoretical estimates. The calculations have also enabled the refinement of structure-activity relationships in this system, to understand the basis of enantiomeric selectivity of substitution at position three of tetrahydroisoquinoline and to identify the role of specific structural waters. Finally, the calculations suggest that the preferred binding mode of trans-(1S,2S)-2-amino-1-tetralol is similar to that of its epimer cis-(1R,2S)-2-amino-1-tetralol and that the ligand does not adopt the novel binding mode proposed in the pdb entry 2AN5 . The work demonstrates how MD simulations and free energy calculations can be used to resolve uncertainties in experimental binding affinities, binding modes, and other aspects related to X-ray refinement and computational drug design.

Short cationic antimicrobial peptides (AMPs) are believed to act either by inducing transmembrane pores or disrupting membranes in a detergent-like manner. For example, the antimicrobial peptides aurein 1.2, citropin 1.1, maculatin 1.1 and caerin 1.1, despite being closely related, appear to act by fundamentally different mechanisms depending on their length. Using molecular dynamics simulations, the structural properties of these four peptides have been examined in solution as well as in a variety of membrane environments. It is shown that each of the peptides has a strong preference for binding to regions of high membrane curvature and that the structure of the peptides is dependent on the degree of local curvature. This suggests that the shorter peptides aurein 1.2 and citropin 1.1 act via a detergent-like mechanism because they can induce high local, but not long-range curvature, whereas the longer peptides maculatin 1.1 and caerin 1.1 require longer range curvature to fold and thus bind to and stabilize transmembrane pores.

The article: “Consciousness in the universe: A review of the ‘OrchOR’ theory” by Hameroff and Penrose reviews work concerning a proposal for the origin of consciousness by Penrose. This model postulates that microtubules can sustain long-lived quantum states and support quantum information processing associated with the effects of quantum gravity (“OrchOR”)...

Enhanced hydrostatic pressure can induce phase transitions in hydrated lipid bilayers especially those composed of saturated phospholipids. In this work, the phase behavior of fully hydrated DPPC (1,2-dipalmitoyl-sn-glycero-3-phosphocholine) and POPC (2-oleoyl-1-palmitoyl-sn-glycero-3-phosphocholine) bilayers as a function of pressure up to 3000 atm has been examined in atomic detail on time scales of up to 1.0 μs, using the molecular dynamics simulation technique. DPPC bilayers formed a rippled gel-like phase comprising a minor disordered fluid-like region and a major ordered gel-like region at 1000 atm, a partially interdigitated gel-like phase at 2000 atm, and a gel-like phase with most of the lipid acyl chains tilted with respect to the plane of the bilayer at 3000 atm. POPC bilayers formed a rippled gel-like phase at 1800, 2400, and 3000 atm. The phase behavior observed for both DPPC and POPC bilayers is in agreement with experiment. The simulations provide insight into the structural changes of DPPC and POPC bilayers as a function of pressure and demonstrate the ability to model biologically relevant lipid systems under high hydrostatic pressure.

The possibility of estimating equilibrium free-energy profiles from multiple non-equilibrium simulations using the fluctuation-dissipation theory or the relation proposed by Jarzynski has attracted much attention. Although the Jarzynski estimator has poor convergence properties for simulations far from equilibrium, corrections have been derived for cases in which the work is Gaussian distributed. Here, we examine the utility of corrections proposed by Gore and collaborators using a simple dissipative system as a test case. The system consists of a single methane-like particle in explicit water. The Jarzynski equality is used to estimate the change in free energy associated with pulling the methane particle a distance of 3.9 nm at rates ranging from similar to 0.1 to 100 ms(-1). It is shown that although the corrections proposed by Gore and collaborators have excellent numerical performance, the profiles still converge slowly. Even when the corrections are applied in an ideal case where the work distribution is necessarily Gaussian, performing simulations under quasi-equilibrium conditions is still most efficient. Furthermore, it is shown that even for a single methane molecule in water, pulling rates as low as 1 ms(-1) can be problematic. The implications of this finding for studies in which small molecules or even large biomolecules are pulled through inhomogeneous environments at similar pulling rates are discussed.

Atomistic molecular dynamics simulations have been used to investigate the conformational changes associated with the binding of human growth hormone (hGH) to the extracellular domains (ECD) of the human growth hormone receptor (hGHR), thereby shedding light on the mechanism of activation. It is shown that the removal of hGH from the hormone-bound receptor complex results in a counter-clockwise rotation of the twosubunits relative to each other by 30°–64° (average 45° ± 14°), in close agreement in terms of both the magnitude and direction of the rotation with that proposed based on mutagenesis experiments. In addition to providing evidence to support a rotational activation mechanism, the simulations have enabled the nature of the interaction interfaces in both the cytokine-bound and unliganded hGHR states to be analyzed in detail. Proteins 2010. © 2009 Wiley-Liss, Inc.

A new force field for the simulation of dipalmitoylphosphatidylcholine (DPPC) in the liquid-crystalline, fluid phase at zero surface tension is presented. The structure of the bilayer with the area per lipid (0.629 nm(2); experiment 0.629-0.64 nm(2)), the volume per lipid (1.226 nm(3); experiment 1.229-1.232 nm(3)), and the ordering of the palmitoyl chains (order parameters) are all in very good agreement with experiment. Experimental electron density profiles are well reproduced in particular with regard to the penetration of water into the bilayer. The force field was further validated by simulating the spontaneous assembly of DPPC into a bilayer in water. Notably, the timescale on which membrane sealing was observed using this model appears closer to the timescales for membrane resealing suggested by electroporation experiments than previous simulations using existing models.

The free-energy landscape of a small protein, the FBP 28 WW domain, has been explored using molecular dynamics (MD) simulations with alternative descriptions of the molecule. The molecular models used range from coarse-grained to all-atom with either an implicit or explicit treatment of the solvent. Sampling of conformation space was performed using both conventional and temperature-replica exchange MD simulations. Experimental chemical shifts and NOEs were used to validate the simulations, and experimental phi values both for validation and as restraints. This combination of different approaches has provided insight into the free energy landscape and barriers encountered by the protein during folding and enabled the characterization of native, denatured and transition states which are compatible with the available experimental data. All the molecular models used stabilize well defined native and denatured basins; however, the degree of agreement with the available experimental data varies. While the most detailed, explicit solvent model predicts the data reasonably accurately, it does not fold despite a simulation time 10 times that of the experimental folding time. The less detailed models performed poorly relative to the explicit solvent model: an implicit solvent model stabilizes a ground state which differs from the experimental native state, and a structure-based model underestimates the size of the barrier between the two states. The use of experimental phi values both as restraints, and to extract structures from unfolding simulations, result in conformations which, although not necessarily true transition states, appear to share the geometrical characteristics of transition state structures. In addition to characterizing the native, transition and denatured states of this particular system in this work, the advantages and limitations of using varying levels of representation are discussed. (C) 2008 Wiley Periodicals, Inc. J Comput Chem 30: 1059-1068, 2009

When estimating binding affinities of a ligand, which can exists in multiple forms, for a target molecule, one must consider all possible competing equilibria. Here, a method is presented that estimates the contribution of the protonation equilibria of a ligand in solution to the measured or calculated binding affinity. The method yields a correction to binding constants that are based on the total concentration of inhibitor (the sum of all ionized forms of the inhibitor in solution) to account for the complexed form of the inhibitor only. The method is applied to the calculation of the difference in binding affinity of two inhibitors, 2-phosphoglycolate (PGA) and its phoshonate analog 3-phosphonopropionate (3PP), for the glycolytic enzyme triosephosphate isomerase. Both inhibitors have three titrating sites and exist in solution as a mixture of different forms. In this case the form that actually binds to the enzyme is present at relative low concentrations. The contributions of the alternative forms to the difference in binding energies is estimated by means of molecular dynamics simulations and corrections. The inhibitors undergo a pKa shift upon binding that is estimated by ab initio calculations. An interesting finding is that the affinity difference of the two inhibitors is not due to different interactions in the active site of the enzyme, but rather due to the difference in the solvation properties of the inhibitors. Protein 2009. © 2008 Wiley-Liss, Inc.

A large variety of antimicrobial peptides have been shown to act, at least in vitro, by potation of the lipid membrane. The nanometre size of these pores, however, complicates their structural characterization by experimental techniques. Here we use molecular dynamics simulations, to study the interaction of a specific class of antimicrobial peptides, melittin, with a dipalmitoylphosphatidylcholine bilayer in atomic detail. We show that transmembrane pores spontaneously form above a critical peptide to lipid ratio. The lipid molecules bend inwards to form a toroidally shaped pore but with only one or two peptides lining the pore. This is in strong contrast to the traditional models of toroidal pores in which the peptides are assumed to adopt a transmembrane orientation. We find that peptide aggregation, either prior or after binding to the membrane surface. is a prerequisite to pore formation. The presence of a stable helical secondary structure of the peptide, however is not. Furthermore, results obtained with modified peptides point to the importance of electrostatic interactions in the potation process. Removing the charges of the basic amino-acid residues of melittin prevents pore formation. It was also found that in the absence of counter ions pores not only form more rapidly but lead to membrane rupture. The rupture process occurs via a novel recursive potation pathway, which we coin the Droste mechanism. (c) 2008 Elsevier B.V. All rights reserved.

A method is presented to enhance the efficiency of simulations of lipid vesicles. The method increases computational speed by eliminating water molecules that either surround the vesicle or reside in the interior of the vesicle, without altering the properties of the water at the membrane interface. Specifically, mean field force approximation (MFFA) boundary potentials are used to replace both the internal and external excess bulk solvent. In addition to reducing the cost of simulating preformed vesicles, the molding effect of the boundary potentials also enhances the formation and equilibration of vesicles from random solutions of lipid in water. Vesicles with diameters in the range from 20 to 60 nm were obtained on a nanosecond time scale, without any noticeable effect of the boundary potentials on their structure.

Electrophoresis is widely used to determine the electrostatic potential of colloidal particles. Oil droplets in pure water show negative or positive electrophoretic mobilities depending on the pH. This is commonly attributed to the adsorption of hydroxyl or hydronium ions, resulting in a negative or positive surface charge, respectively. This explanation, however, is not in agreement with the difference in isoelectric point and point of zero charge observed in experiment. Here we present molecular dynamics simulations of oil droplets in water in the presence of an external electric field but in the absence of any ions. The simulations reproduce the negative sign and the order of magnitude of the oil droplet mobilities at the point of zero charge in experiment. The electrostatic potential in the oil with respect to the water phase, induced by anisotropic dipole orientation in the interface, is positive. Our results suggest that electrophoretic mobility does not always reflect the net charge or electrostatic potential of a suspended liquid droplet and, thus, the interpretation of electrophoresis in terms of purely continuum effects may need to be reevaluated. (C) 2007 Elsevier Inc. All rights reserved.

Class I hydrophobins are fungal proteins that self-assemble into robust amphipathic rodlet monolayers on the surface of aerial structures such as spores and fruiting bodies. These layers share many structural characteristics with amyloid fibrils and belong to the growing family of functional amyloid-like materials produced by microorganisms. Although the three-dimensional structure of the soluble monomeric form of a class I hydrophobin has been determined, little is known about the molecular structure of the rodlets or their assembly mechanism. Several models have been proposed, some of which suggest that the Cys3–Cys4 loop has a critical role in the initiation of assembly or in the polymeric structure. In order to provide insight into the relationship between hydrophobin sequence and rodlet assembly, we investigated the role of the Cys3–Cys4 loop in EAS, a class I hydrophobin from Neurospora crassa. Remarkably, deletion of up to 15 residues from this 25-residue loop does not impair rodlet formation or reduce the surface activity of the protein, and the physicochemical properties of rodlets formed by this mutant are indistinguishable from those of its full-length counterpart. In addition, the core structure of the truncation mutant is essentially unchanged. Molecular dynamics simulations carried out on the full-length protein and this truncation mutant binding to an air–water interface show that, although it is hydrophobic, the loop does not play a role in positioning the protein at the surface. These results demonstrate that the Cys3–Cys4 loop does not have an integral role in the formation or structure of the rodlets and that the major determinant of the unique properties of these proteins is the amphipathic core structure, which is likely to be preserved in all hydrophobins despite the high degree of sequence variation across the family.

A series of molecular dynamics simulations in explicit solvent were started from nine structural models of the transition state of the SH3 domain of alpha-spectrin, which were generated by Lindorff Larsen et al. (Nat Struct Mol Biol 2004;11:443-449) using molecular dynamics simulations in which experimental (Phi-values were incorporated as restraints. Two of the nine models were simulated 10 times for 200 ns and the remaining models simulated two times for 200 us. Complete folding was observed in one case, while in the other simulations partial folding or unfolding even ts were observed, which were characterized by a regularization of elements of secondary structure. These results are consistent with recent experimental evidence that the folding of SH3 domains involves low populated intermediate states. (C) 2007 Wiley-Liss, Inc.

The Death Receptor 5 (DR5), a member of tumor necrosis factor receptor (TNFR) superfamily of receptors, triggers apoptosis (programmed cell death) when stimulated by its tridentate ligand TRAIL. Until recently it was generally assumed that the activation of DR5 resulted from the recruitment of three independent receptor units, leading to the trimerization of intracellular domains. However, there is mounting evidence to suggest that, in the absence of ligand, such cytokine receptors primarily reside as preformed complexes. In this work, molecular dynamics simulations of the TRALL-DR5 complex, the unbound receptor trimer and individual receptor monomers are compared to gain insight in the mechanism of activation. The results suggest that, in the absence of TRAIL, DR5 has a strong propensity to self-associate and that this is primarily mediated through interactions of the membrane proximal domains. The association of the free receptors leads to a loss of the threefold symmetry found within the receptor-ligand complex. The simulations suggest that the primary role of TRAIL is to induce threefold-symmetry within the DR5 complex and to constrain the receptor to a specific conformation. The implications of this in terms of the mechanism by which the receptor switches from an inactive to an active state are discussed.

The question of why plants release isoprene when heat stressed has been continuously debated for more than half a century. In this work we use molecular dynamics simulation techniques to directly investigate the interaction between isoprene and a model phospholipid membrane in atomic detail. It is found that isoprene partitions preferentially in the center of the membrane and in a dose dependent manner enhances the order within the membrane without significantly changing the dynamical properties of the system. At a concentration of 20 mol% isoprene (16 isoprene molecules per 64 lipid molecules) the effect of the addition of isoprene on the membrane order is equivalent to a reduction in temperature of 10 K, rising to a reduction of 30 K at 43 mol% isoprene. The significance of the work is that it provides for the first time direct evidence that isoprene stabilizes lipid membranes and reduces the likelihood of a phospholipid membrane undergoing a heat induced phase transition. Furthermore it provides a clear mechanistic picture as to why plants specifically utilize isoprene for this purpose.

Molecular dynamics simulations of fully hydrated pure bilayers of four widely studied phospholipids, 1,2-dilauroyl-sn-glycero-3-phosphocholine (DLPC), 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC), 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC), and 2-oleoyl-1-palmitoyl-sn-glycero-3-phosphocholine (POPC) using a recent revision of the GROMOS96 force field are reported. It is shown that the force field reproduces the structure and the hydration of bilayers formed by each of the four lipids with high accuracy. Specifically, the solvation and the orientation of the dipole of the phosphocholine headgroup and of the ester carbonyls show that the structure of the primary hydration shell in the simulations closely matches experimental findings. This work highlights the need to reproduce a broad range of properties beyond the area per lipid, which is poorly defined experimentally, and to consider the effect of system size and sampling times well beyond those commonly used.

The Automated force field Topology Builder (ATB, http://compbio.biosci.uq.edu.au/atb ) is a Web-accessible server that can provide topologies and parameters for a wide range of molecules appropriate for use in molecular simulations, computational drug design, and X-ray refinement. The ATB has three primary functions: (1) to act as a repository for molecules that have been parametrized as part of the GROMOS family of force fields, (2) to act as a repository for pre-equilibrated systems for use as starting configurations in molecular dynamics simulations (solvent mixtures, lipid systems pre-equilibrated to adopt a specific phase, etc.), and (3) to generate force field descriptions of novel molecules compatible with the GROMOS family of force fields in a variety of formats (GROMOS, GROMACS, and CNS). Force field descriptions of novel molecules are derived using a multistep process in which results from quantum mechanical (QM) calculations are combined with a knowledge-based approach to ensure compatibility (as far as possible) with a specific parameter set of the GROMOS force field. The ATB has several unique features: (1) It requires that the user stipulate the protonation and tautomeric states of the molecule. (2) The symmetry of the molecule is analyzed to ensure that equivalent atoms are assigned identical parameters. (3) Charge groups are assigned automatically. (4) Where the assignment of a given parameter is ambiguous, a range of possible alternatives is provided. The ATB also provides several validation tools to assist the user to assess the degree to which the topology generated may be appropriate for a given task. In addition to detailing the steps involved in generating a force field topology compatible with a specific GROMOS parameter set (GROMOS 53A6), the challenges involved in the automatic generation of force field parameters for atomic simulations in general are discussed.

The stability of the crystal structure of the multidrug transporter P-glycoprotein proposed by Aller et al. (PDBid 3G5U ) has been examined under different environmental conditions using molecular dynamics. We show that in the presence of the detergent cholate, the structure of P-glycoprotein solved at pH 7.5 is stable. However, when incorporated into a cholesterol-enriched POPC membrane in the presence of 150 mM NaCl, the structure rapidly deforms. Only when the simulation conditions closely matched the experimental conditions under which P-glycoprotein is transport active was a stable conformation obtained. Specifically, the presence of Mg(2+), which bound to distinct sites in the nucleotide binding domains (NBDs), and the double protonation of the catalytic histidines (His583 and His1228) and His149 were required. While the structure obtained in a membrane environment under these conditions is very similar to the crystal structure of Aller et al., there are several key differences. The NBDs are in direct contact, reminiscent of the open state of MalK. The angle between the transmembrane domains is also increased, resulting in an outward motion of the intracellular loops. Notably, the structures obtained from the simulations provide a better match to a range of experimental cross-linking data than does the original 3G5U-a crystal structure. This work highlights the effect small changes in environmental conditions can have of the conformation of a membrane protein and the importance of representing the experimental conditions appropriately in modeling studies.

Fatty acids containing a terminal cyclic group such as cyclohexyl and cycloheptyl are commonly found in prokaryotic membranes, especially in those of thermo-acidophilic bacteria. These so-called ω-alicyclic fatty acids have been proposed to stabilize the membranes of bacteria by reducing the fluidity in membranes and increasing lipid packing and lipid chain order. In this article, molecular dynamics simulations are used to examine the effect of 3- to 7-membered cycloalkyl saturated and unsaturated (cyclopent-2-enyl and phenyl) rings in ω-alicyclic fatty acyl chains on the structure (lipid packing, lipid chain order, and fraction of gauche defects in the chains) and dynamics (lateral lipid diffusion) of a model lipid bilayer. It was found that ω-alicyclic chains in which the ring was saturated reduced lipid condensation and lowered chain order which would be associated with enhanced fluidity. However, this effect was limited. The lateral diffusion of the lipids diminished as the ring size increased. In particular, ω-cyclohexyl and ω-cycloheptyl acyl tails led to a decrease in lipid diffusion. In contrast, ω-alicyclic acyl chains that contain an unsaturated ring promoted membrane fluidity both in terms of changes in membrane structure and lipid diffusion. This may indicate that saturated and unsaturated terminal rings in ω-alicyclic fatty acids fulfill alternative functions within membranes. Overall, the simulations suggest that ω-alicyclic fatty acids in which the terminal ring is saturated might protect the membrane of thermo-acidophilic bacteria from high-temperature and low-pH conditions through a "dynamical barrier" that would limit lipid diffusion and transmembrane diffusion of undesired ions and molecules.

Molecular weight distributions of starch branches affect functional properties, which can be controlled by engineering starch branching enzymes (SBEs). Molecular dynamics and docking approaches are used to examine interactions between SBE and starch fragments. In the native protein, three residues formed stable interactions with starch fragments in the central binding region; these residues may play key roles in substrate recognition. Fragments containing 5-12 glucose units interacted more tightly with SBE than smaller fragments, suggesting a minimal functional fragment size of 5, agreeing with experiment. Effects of three different point mutations on interactions with maltopentaose in the central binding region correlated well with experiment. Simulations indicate that SBE may template formation of the crystalline lamellae characteristic of native starch, consistent with the observation that crystalline lamellae formed by starch in a plant, are not necessarily the state of lowest free energy. The methodology will help develop starches with optimized functional properties.

Hydrophobins are small ( approximately 100 aa) proteins that have an important role in the growth and development of mycelial fungi. They are surface active and, after secretion by the fungi, self-assemble into amphipathic membranes at hydrophobic/hydrophilic interfaces, reversing the hydrophobicity of the surface. In this study, molecular dynamics simulation techniques have been used to model the process by which a specific class I hydrophobin, SC3, binds to a range of hydrophobic/hydrophilic interfaces. The structure of SC3 used in this investigation was modeled based on the crystal structure of the class II hydrophobin HFBII using the assumption that the disulfide pairings of the eight conserved cysteine residues are maintained. The proposed model for SC3 in aqueous solution is compact and globular containing primarily beta-strand and coil structures. The behavior of this model of SC3 was investigated at an air/water, an oil/water, and a hydrophobic solid/water interface. It was found that SC3 preferentially binds to the interfaces via the loop region between the third and fourth cysteine residues and that binding is associated with an increase in alpha-helix formation in qualitative agreement with experiment. Based on a combination of the available experiment data and the current simulation studies, we propose a possible model for SC3 self-assembly on a hydrophobic solid/water interface.