Search

The fusion of the viral and target cell membranes is a key step in the life cycle of all enveloped viruses. Here, a range of structural data is used to generate an evidence-based model of the active conformation of an archetypical type-I fusion protein, the Ebola glycoprotein 2 (GP2). The stability of the trimeric complex is demonstrated using molecular dynamics and validated by simulating the interaction of the complex with a lipid bilayer. In this model, the fusion peptides project away from the central helix bundle parallel to the target membrane. This maximizes contact with the host membrane, enhances lateral stability, and would explain why, when activated, viral fusion proteins are trimeric.

α-Helical membrane-active antimicrobial peptides (AMPs) are known to act via a range of mechanisms, including the formation of barrel-stave and toroidal pores and the micellisation of the membrane (carpet mechanism). Different mechanisms imply that the peptides adopt different 3D structures when bound at the water-membrane interface, a highly amphipathic environment. Here, an evolutionary algorithm is used to predict the 3D structure of a range of α-helical membrane-active AMPs at the water-membrane interface by optimising amphipathicity. This amphipathic structure prediction (ASP) is capable of distinguishing between curved and linear peptides solved experimentally, potentially allowing the activity and mechanism of action of different membrane-active AMPs to be predicted. The ASP algorithm is accessible via a web interface at http://atb.uq.edu.au/asp/.

Atomistic nonequilibrium molecular dynamics simulations have been used to model the orientation of a side-chain substituted dicyanovinyl oligothiophene, DCV4T-Et2, in thin films formed by vacuum deposition as used in organic photovoltaics. The orientation of the DCV4T-Et2 molecules was analyzed in neat layers deposited onto smooth or rough C60 substrates. The average angle between the long axis of DCV4T-Et2 and the horizontal was 21 +/- 1 degrees in the layer deposited on smooth C-60, in agreement with experimental measurements of layers deposited on silicon. In the layer deposited on rough C-60, the average angle was 25 +/- 2 degrees, which could lead to a decrease of 6% in the efficiency of light absorption compared to the layer on the smooth substrate. Thus, the simulations suggest that surface roughness must be considered in the design of model systems in order to maintain relevance to devices. Importantly, the simulations provide the ability to analyze the orientation of each individual molecule as opposed to experiments that provide an average across a film. The orientation of the thiophene rings of DCV4T-Et2 and the variation in orientation with respect to distance from the substrate were also analyzed and found to be similar for films on smooth and rough substrates. In addition, we show that X-ray diffraction patterns of the layers can be constructed from the simulations.

Warfarin, a widely used oral anticoagulant, is prescribed as a racemic mixture. Each enantiomer of neutral Warfarin can exist in 20 possible tautomeric states leading to complex pharmacokinetics and uncertainty as to the relevant species under different conditions. Here, the ability of alternative computational approaches to predict the preferred tautomeric form(s) of neutral Warfarin in different solvents is examined. It is shown that varying the method used to estimate the heat of formation in vacuum (direct or via homodesmic reactions), whether entropic corrections were included, and the method used to estimate the free enthalpy of solvation (i.e., PCM, COSMO, or SMD implicit models or explicit solvent) lead to large differences in the predicted rank and relative populations of the tautomers. In this case, only a combination of the enthalpy of formation using homodesmic reactions and explicit solvent to estimate the free enthalpy of solvation yielded results compatible with the available experimental data. The work also suggests that a small but significant subset of the possible Warfarin tautomers are likely to be physiologically relevant.

E-cadherin is a transmembrane glycoprotein that facilitates inter-cellular adhesion in the epithelium. The ectodomain of the native structure is comprised of five repeated immunoglobulin-like domains. All E-cadherin crystal structures show the protein in one of three alternative conformations: a monomer, a strand-swapped trans homodimer and the so-called X-dimer, which is proposed to be a kinetic intermediate to forming the strand-swapped trans homodimer. However, previous studies have indicated that even once the trans strand-swapped dimer is formed, the complex is highly dynamic and the E-cadherin monomers may reorient relative to each other. Here, molecular dynamics simulations have been used to investigate the stability and conformational flexibility of the human E-cadherin trans strand-swapped dimer. In four independent, 100 ns simulations, the dimer moved away from the starting structure and converged to a previously unreported structure, which we call the Y-dimer. The Y-dimer was present for over 90% of the combined simulation time, suggesting that it represents a stable conformation of the E-cadherin dimer in solution. The Y-dimer conformation is stabilised by interactions present in both the trans strand-swapped dimer and X-dimer crystal structures, as well as additional interactions not found in any E-cadherin dimer crystal structures. The Y-dimer represents a previously unreported, stable conformation of the human E-cadherin trans strand-swapped dimer and suggests that the available crystal structures do not fully capture the conformations that the human E-cadherin trans homodimer adopts in solution.

Based on differences between the x-ray crystal structures of ligand-bound and unbound forms, the activation of the erythropoietin receptor (EPOR) was initially proposed to involve a cross-action scissorlike motion. However, the validity of the motions involved in the scissorlike model has been recently challenged. Here, atomistic molecular dynamics simulations are used to examine the structure of the extracellular domain of the EPOR dimer in the presence and absence of erythropoietin and a series of agonistic or antagonistic mimetic peptides free in solution. The simulations suggest that in the absence of crystal packing effects, the EPOR chains in the different dimers adopt very similar conformations with no clear distinction between the agonist and antagonist-bound complexes. This questions whether the available x-ray crystal structures of EPOR truly represent active or inactive conformations. The study demonstrates the difficulty in using such structures to infer a mechanism of action, especially in the case of membrane receptors where just part of the structure has been considered in addition to potential confounding effects that arise from the comparison of structures in a crystal as opposed to a membrane environment. The work highlights the danger of assigning functional significance to small differences between structures of proteins bound to different ligands in a crystal environment without consideration of the effects of the crystal lattice and thermal motion.

Biobutanol production by fermentation is potentially a sustainable alternative to butanol production from fossil fuels. However, the toxicity of butanol to fermentative bacteria, resulting largely from cell membrane fluidization, limits production titers and is a major factor limiting the uptake of the technology. Here, studies were undertaken, in vitro and in silico, on the butanol effects on a representative bacterial (i.e. Escherichia coli) inner cell membrane. A critical butanol : lipid ratio for stability of 2 : 1 was observed, computationally, consistent with complete interdigitation. However, at this ratio the bilayer was ∼20% thicker than for full interdigitation. Furthermore, butanol intercalation induced acyl chain bending and increased disorder, measured as a 27% lateral diffusivity increase experimentally in a supported lipid bilayer. There was also a monophasic Tm reduction in butanol-treated large unilamellar vesicles. Both behaviours are inconsistent with an interdigitated gel. Butanol thus causes only partial interdigitation at physiological temperatures, due to butanol accumulating at the phospholipid headgroups. Acyl tail disordering (i.e. splaying and bending) fills the subsequent voids. Finally, butanol short-circuits the bilayer and creates a coupled system where interdigitated and splayed phospholipids coexist. These findings will inform the design of strategies targeting bilayer stability for increasing biobutanol production titers. Ministry of Education (MOE) Nanyang Technological University National Research Foundation (NRF) Accepted version The computational work for this article was performed on resources of the National Supercomputing Centre, Singapore (https://www.nscc.sg). This work was supported by the Ministry of Education, Singapore (MOE) Grants M4360005 and Tier 1 RG146/17. SCELSE is funded by Singapore’s Ministry of Education, National Research Federation, Nanyang Technological University (NTU), and National University of Singapore (NUS) and hosted by NTU in partnership with NUS.

Biological membranes display a great diversity in lipid composition and lateral structure that is crucial in a variety of cellular functions. Simulations of membranes have contributed significantly to the understanding of the properties, functions and behaviour of membranes and membrane-protein assemblies. This success relies on the ability of the force field used to describe lipid-lipid and lipid-environment interactions accurately, reproducibly and realistically. In this review, we present some recent progress in lipid force-field development and validation strategies. In particular, we highlight how a range of properties obtained from various experimental techniques on lipid bilayers and membranes, can be used to assess the quality of a force field. We discuss the limitations and assumptions that are inherent to both computational and experimental approaches and how these can influence the comparison between simulations and experimental data. This article is part of a Special Issue entitled: Membrane Proteins edited by J.C. Gumbart and Sergei Noskov.

The apparent activity of the multidrug transporter P-glycoprotein (P-gp) is enhanced by the presence of cholesterol. Whether this is due to the direct effect of cholesterol on the activity of P-gp, its effect on the local concentration of substrate in the membrane, or its effect on the rate of entry of the drug into the cell, is unknown. In this study, molecular dynamics simulation techniques coupled with potential of mean force calculations have been used to investigate the role of cholesterol in the movement of four P-gp substrates across a POPC bilayer in the presence or absence of 10% cholesterol. The simulations suggest that the presence of cholesterol lowers the free energy associated with entering the middle of the bilayer in a substrate-specific manner. These findings suggest that P-gp substrates may preferentially accumulate in cholesterol-rich regions of the membrane, which may explain its enhanced transport activity.

Free to read at publisher's site. Many venom peptides are potent and selective inhibitors of voltage-gated ion channels, including channels that are validated therapeutic targets for treatment of a wide range of human diseases. However, the development of novel venom-peptide-based therapeutics requires an understanding of their mechanism of action. In the case of voltage-gated ion channels, venom peptides act either as pore blockers that bind to the extracellular side of the channel pore or gating modifiers that bind to one or more of the membrane-embedded voltage sensor domains. In the case of gating modifiers, it has been debated whether the peptide must partition into the membrane to reach its binding site. In this study, we used surface plasmon resonance, fluorescence spectroscopy and molecular dynamics to directly compare the lipid-binding properties of two gating modifiers (μ-TRTX-Hd1a and ProTx-I) and two pore blockers (ShK and KIIIA). Only ProTx-I was found to bind to model membranes. Our results provide further evidence that the ability to insert into the lipid bilayer is not a requirement to be a gating modifier. In addition, we characterised the surface of ProTx-I that mediates its interaction with neutral and anionic phospholipid membranes and show that it preferentially interacts with anionic lipids.

The human multidrug transporter P-glycoprotein (P-gp) transports over 200 chemically diverse substrates, influencing their bioavailability and tissue distribution. Pharmacological studies have identified both competitive and noncompetitive P-gp substrates, but neither the precise location of the substrate binding sites, nor the basis of competitive and noncompetitive interactions has been fully characterized. Here, potential of mean force (PMF) calculations are used to identify the transport-competent minimum free energy binding locations of five compounds, Hoechst 33342, Rhodamine 123, paclitaxel, tariquidar, and verapamil to P-gp. Unrestrained molecular dynamics simulations were also performed to confirm the substrates were stable in the energy wells determined using the PMF calculations. All compounds had energy minima within the P-gp transmembrane (TM) pore. For Hoechst 33342 and Rhodamine 123, a second minimum outside the TM pore was also identified. Based on this and previous studies of nicardipine and morphine [ Subramanian et al. J. Chem. Inf. Model. 2015 , 55 , 1202 ], a general scheme that accounts for the observed noncompetitive and competitive substrate interactions with P-gp is proposed.

The ability of atomic interaction parameters generated using the Automated Topology Builder and Repository version 3.0 (ATB3.0) to predict experimental hydration free enthalpies (ΔGwater) and solvation free enthalpies in the apolar solvent hexane (ΔGhexane) is presented. For a validation set of 685 molecules the average unsigned error (AUE) between ΔGwater values calculated using the ATB3.0 and experiment is 3.8 kJ·mol-1. The slope of the line of best fit is 1.00, the intercept -1.0 kJ·mol-1, and the R2 0.90. For the more restricted set of 239 molecules used to validate OPLS3 (J. Chem. Theory Comput. 2016, 12, 281-296, DOI: 10.1021/acs.jctc.5b00864) the AUE using the ATB3.0 is just 2.7 kJ·mol-1 and the R2 0.93. A roadmap for further improvement of the ATB parameters is presented together with a discussion of the challenges of validating force fields against the available experimental data.

Molecular dynamics simulations and free energy calculations have been used to investigate the effect of ligand binding on the enantioselectivity of an epoxide hydrolase (EH) from Aspergillus niger. Despite sharing a common mechanism, a wide range of alternative mechanisms have been proposed to explain the origin of enantiomeric selectivity in EHs. By comparing the interactions of ( R)- and ( S)-glycidyl phenyl ether (GPE) with both the wild type (WT, E = 3) and a mutant showing enhanced enantioselectivity to GPE (LW202, E = 193), we have examined whether enantioselectivity is due to differences in the binding pose, the affinity for the ( R)- or ( S)- enantiomers, or a kinetic effect. The two enantiomers were easily accommodated within the binding pockets of the WT enzyme and LW202. Free energy calculations suggested that neither enzyme had a preference for a given enantiomer. The two substrates sampled a wide variety of conformations in the simulations with the sterically hindered and unhindered carbon atoms of the GPE epoxide ring both coming in close proximity to the nucleophilic aspartic acid residue. This suggests that alternative pathways could lead to the formation of a ( S)- and ( R)-diol product. Together, the calculations suggest that the enantioselectivity is due to kinetic rather than thermodynamic effects and that the assumption that one substrate results in one product when interpreting the available experimental data and deriving E-values may be inappropriate in the case of EHs.

Despite decades of research, the mechanism of action of the ABC multidrug transporter P-glycoprotein (P-gp) remains elusive. Due to experimental limitations, many researchers have turned to molecular dynamics simulation studies in order to investigate different aspects of P-gp function. However, such studies are challenging and caution is required when interpreting the results. P-gp is highly flexible and the time scale on which it can be simulated is limited. There is also uncertainty regarding the accuracy of the various crystal structures available, let alone the structure of the protein in a physiologically relevant environment. In this study, three alternative structural models of mouse P-gp (3G5U, 4KSB, 4M1M), all resolved to 3.8 Å, were used to initiate sets of simulations of P-gp in a membrane environment in order to determine: a) the sensitivity of the results to differences in the starting configuration; and b) the extent to which converged results could be expected on the times scales commonly simulated for this system. The simulations suggest that the arrangement of the nucleotide binding domains (NBDs) observed in the crystal structures is not stable in a membrane environment. In all simulations, the NBDs rapidly associated (within 10 ns) and changes within the transmembrane helices were observed. The secondary structure within the transmembrane domain was best preserved in the 4M1M model under the simulation conditions used. However, the extent to which replicate simulations diverged on a 100 to 200 ns timescale meant that it was not possible to draw definitive conclusions as to which structure overall was most stable, or to obtain converged and reliable results for any of the properties examined. The work brings into question the reliability of conclusions made in regard to the nature of specific interactions inferred from previous simulation studies on this system involving similar sampling times. It also highlights the need to demonstrate the statistical significance of any results obtained in simulations of large flexible proteins, especially where the initial structure is uncertain.

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.

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...

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.

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.

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.

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.

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 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.