Permeability of Dermatological Solutes through the Short Periodicity Phase of Human Stratum Corneum Lipid Bilayers.

Determining the permeability of drug-like solutes through the densely packed and heterogeneous stratum corneum lipid layer presents a significant challenge. In this study, we employed umbrella sampling with a periodic weighing function applied to the center of mass of the lipid bilayers. Precise umbrella sampling was conducted with an interframe distance of 0.5 Å, spanning from the bilayer center to the water phase, and each frame was simulated for at least 20 ns. Autocorrelation functions, potential of mean force (PMF), and diffusivity profiles were analyzed for six solutes (testosterone, benzene, caffeine, ethanol, mannitol, and histidine). The results revealed that autocorrelations were dependent solely on the medium, whether water or lipid phase. Diffusivity and PMF profiles along the reaction coordinate were influenced by the hydrophilicity of the solute rather than its size. For hydrophobic solutes, the PMF curves exhibited a minimum at the bilayer center, while for hydrophilic solutes, the PMFs peaked at the bilayer center and lipid tails (where the lipid tails are not interacting with the cholesterol). Diffusivity curves were low at the bilayer center and water phase, with peaks observed at the headgroup or the boundary between fatty acid and cholesterol (1 nm from the bilayer center). The quantitative findings presented in this work hold significance for pharmacists and drug designers.

Probing Disaccharide Binding to Triplatin as Models for Tumor Cell Heparan Sulfate (GAG) Interactions.

In this study, we have used [1H, 15N] NMR spectroscopy to investigate the interactions of the trinuclear platinum anticancer drug triplatin (1) (1,0,1/t,t,t or BBR3464) with site-specific sulfated and carboxylated disaccharides. Specifically, the disaccharides GlcNS(6S)-GlcA (I) and GlcNS(6S)-IdoA(2S) (II) are useful models of longer-chain glycosaminoglycans (GAGs) such as heparan sulfate (HS). For both the reactions of 15N-1 with I and II, equilibrium conditions were achieved more slowly (65 h) compared to the reaction with the monosaccharide GlcNS(6S) (9 h). The data suggest both carboxylate and sulfate binding of disaccharide I to the Pt with the sulfato species accounting for <1% of the total species at equilibrium. The rate constant for sulfate displacement of the aqua ligand (kL2) is 4 times higher than the analogous rate constant for carboxylate displacement (kL1). There are marked differences in the equilibrium concentrations of the chlorido, aqua, and carboxy-bound species for reactions with the two disaccharides, notably a significantly higher concentration of carboxylate-bound species for II, where sulfate-bound species were barely detectable. The trend mirrors that reported for the corresponding dinuclear platinum complex 1,1/t,t, where the rate constant for sulfate displacement of the aqua ligand was 3 times higher than that for acetate. Also similar to what we observed for the reactions of 1,1/t,t with the simple anions, aquation of the sulfato group is rapid, and the rate constant k-L2 is 3 orders of magnitude higher than that for displacement of the carboxylate (k-L1). Molecular dynamics calculations suggest that extra hydrogen-bonding interactions with the more sulfated disaccharide II may prevent or diminish sulfate binding of the triplatin moiety. The overall results suggest that Pt-O donor interactions should be considered in any full description of platinum complex cellular chemistry.

Understanding the Binding of Starch Fragments to Granule-Bound Starch Synthase.

Granule-bound starch synthase (GBSS) plays a major role, that of chain elongation, in the biosynthesis of amylose, a starch component with mostly (1 → 4)-α connected long chains of glucose with a few (1 → 6)-α branch points. Chain-length distributions (CLDs) of amylose affect functional properties, which can be controlled by changing appropriate residues on granule-bound starch synthase (GBSS). Knowing the binding of GBSS and amylose at a molecular level can help better determine the key amino acids on GBSS that affect CLDs of amylose for subsequent use in molecular engineering. Atomistic molecular dynamics simulations with explicit solvent and docking approaches were used in this study to build a model of the binding between rice GBSS and amylose. Amylose fragments containing 3-12 linearly linked glucose units were built to represent the starch fragments. The stability of the complexes, interactions between GBSS and sugars, and difference in structure/conformation of bound and free starch fragments were analyzed. The study found that starch/amylose fragments with 5 or 6 glucose units were suitable for modeling starch binding to GBSS. The removal of an interdomain disulfide on GBSS was found to affect both GBSS and starch stability. Key residues that could affect the binding ability were also indicated. This model can help rationalize the design of mutants and suggest ways to make single-point mutations, which could be used to develop plants producing starches with improved functional properties.

Crystal Structures of Protein-Bound Cyclic Peptides.

Cyclization is an important post-translational modification of peptides and proteins that confers key advantages such as protection from proteolytic degradation, altered solubility, membrane permeability, bioavailability, and especially restricted conformational freedom in water that allows the peptide backbone to adopt the major secondary structure elements found in proteins. Non-ribosomal synthesis in bacteria, fungi, and plants or synthetic chemistry can introduce unnatural amino acids and non-peptidic constraints that modify peptide backbones and side chains to fine-tune cyclic peptide structure. Structures can be potentially altered further upon binding to a protein in biological environments. Here we analyze three-dimensional crystal structures for 211 bioactive cyclic peptides bound to 65 different proteins. The protein-bound cyclic peptides were examined for similarities and differences in bonding modes, for main-chain and side-chain structure, and for the importance of polarity, hydrogen bonds, hydrophobic effects, and water molecules in interactions with proteins. Many protein-bound cyclic peptides show backbone structures like those (strands, sheets, turns, helices, loops, or distorted variations) found at protein-protein binding interfaces. However, the notion of macrocycles simply as privileged scaffolds that primarily project side-chain substituents for complementary interactions with proteins is dispelled here. Unlike small-molecule drugs, the cyclic peptides do not rely mainly upon hydrophobic and van der Waals interactions for protein binding; they also use their main chain and side chains to form polar contacts and hydrogen bonds with proteins. Compared to small-molecule ligands, cyclic peptides can bind across larger, polar, and water-exposed protein surface areas, making many more contacts that can increase affinity, selectivity, biological activity, and ligand-receptor residence time. Cyclic peptides have a greater capacity than small-molecule drugs to modulate protein-protein interfaces that involve large, shallow, dynamic, polar, and water-exposed protein surfaces.

Effect of Binding on Enantioselectivity of Epoxide Hydrolase.

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.

Molecular dynamics simulations reveal structural insights into inhibitor binding modes and functionality in human Group IIA phospholipase A2.

Human Group IIA phospholipase A2 (hGIIA) promotes inflammation in immune-mediated pathologies by regulating the arachidonic acid pathway through both catalysis-dependent and -independent mechanisms. The hGIIA crystal structure, both alone and inhibitor-bound, together with structures of closely related snake-venom-derived secreted phospholipase enzymes has been well described. However, differentiation of biological and nonbiological contacts and the relevance of structures determined from snake venom enzymes to human enzymes are not clear. We employed molecular dynamics (MD) and docking approaches to understand the binding of inhibitors that selectively or nonselectively block the catalysis-independent mechanism of hGIIA. Our results indicate that hGIIA behaves as a monomer in the solution environment rather than a dimer arrangement that is in the asymmetric unit of some crystal structures. The binding mode of a nonselective inhibitor, KH064, was validated by a combination of the experimental electron density and MD simulations. The binding mode of the selective pentapeptide inhibitor FLSYK to hGIIA was stipulated to be different to that of the snake venom phospholipases A2 of Daboia russelli pulchella (svPLA2 ). Our data suggest that the application of MD approaches to crystal structure data is beneficial in evaluating the robustness of conclusions drawn based on crystal structure data alone. Proteins 2017; 85:827-842. © 2016 Wiley Periodicals, Inc.

Learning epistatic interactions from sequence-activity data to predict enantioselectivity.

Enzymes with a high selectivity are desirable for improving economics of chemical synthesis of enantiopure compounds. To improve enzyme selectivity mutations are often introduced near the catalytic active site. In this compact environment epistatic interactions between residues, where contributions to selectivity are non-additive, play a significant role in determining the degree of selectivity. Using support vector machine regression models we map mutations to the experimentally characterised enantioselectivities for a set of 136 variants of the epoxide hydrolase from the fungus Aspergillus niger (AnEH). We investigate whether the influence a mutation has on enzyme selectivity can be accurately predicted through linear models, and whether prediction accuracy can be improved using higher-order counterparts. Comparing linear and polynomial degree = 2 models, mean Pearson coefficients (r) from [Formula: see text]-fold cross-validation increase from 0.84 to 0.91 respectively. Equivalent models tested on interaction-minimised sequences achieve values of [Formula: see text] and [Formula: see text]. As expected, testing on a simulated control data set with no interactions results in no significant improvements from higher-order models. Additional experimentally derived AnEH mutants are tested with linear and polynomial degree = 2 models, with values increasing from [Formula: see text] to [Formula: see text] respectively. The study demonstrates that linear models perform well, however the representation of epistatic interactions in predictive models improves identification of selectivity-enhancing mutations. The improvement is attributed to higher-order kernel functions that represent epistatic interactions between residues.

Elucidating the Spatial Arrangement of Emitter Molecules in Organic Light-Emitting Diode Films.

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.

Binding of Starch Fragments to the Starch Branching Enzyme: Implications for Developing Slower-Digesting Starch.

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.

Testing and validation of the Automated Topology Builder (ATB) version 2.0: prediction of hydration free enthalpies.

To test and validate the Automated force field Topology Builder and Repository (ATB; http://compbio.biosci.uq.edu.au/atb/ ) the hydration free enthalpies for a set of 214 drug-like molecules, including 47 molecules that form part of the SAMPL4 challenge have been estimated using thermodynamic integration and compared to experiment. The calculations were performed using a fully automated protocol that incorporated a dynamic analysis of the convergence and integration error in the selection of intermediate points. The system has been designed and implemented such that hydration free enthalpies can be obtained without manual intervention following the submission of a molecule to the ATB. The overall average unsigned error (AUE) using ATB 2.0 topologies for the complete set of 214 molecules was 6.7 kJ/mol and for molecules within the SAMPL4 7.5 kJ/mol. The root mean square error (RMSE) was 9.5 and 10.0 kJ/mol respectively. However, for molecules containing functional groups that form part of the main GROMOS force field the AUE was 3.4 kJ/mol and the RMSE was 4.0 kJ/mol. This suggests it will be possible to further refine the parameters provided by the ATB based on hydration free enthalpies.

On the Validation of Protein Force Fields Based on Structural Criteria.

Molecular dynamics simulations depend critically on the quality of the force field used to describe the interatomic interactions and the extent to which it has been validated for use in a specific application. Using a curated test set of 52 high-resolution structures, 39 derived from X-ray diffraction and 13 solved using NMR, we consider the extent to which different parameter sets of the GROMOS protein force field can be distinguished based on comparing a range of structural criteria, including the number of backbone hydrogen bonds, the number of native hydrogen bonds, polar and nonpolar solvent-accessible surface area, radius of gyration, the prevalence of secondary structure elements, J-coupling constants, nuclear Overhauser effect (NOE) intensities, positional root-mean-square deviations (RMSD), and the distribution of backbone ϕ and ψ dihedral angles. It is shown that while statistically significant differences between the average values of individual metrics could be detected, these were in general small. Furthermore, improvements in agreement in one metric were often offset by loss of agreement in another. The work establishes a framework and test set against which protein force fields can be validated. It also highlights the danger of inferring the relative quality of a given force field based on a small range of structural properties or small number of proteins.

Development of a model for granule-bound starch synthase I activity using free-energy calculations.

Starch is a branched polymer of glucose with two components, both of which have (1 â†’ 4)-α linear links and (1 â†’ 6)-α branch points: amylopectin, of high molecular weight with many short branches, and amylose, of lower molecular weight and only a few long-chain branches. Granule-bound starch synthase I (GBSSI) is one of the main enzymes controlling amylose synthesis and chain-length distribution. As production of different GBSSI mutants is time-consuming and laborious, molecular dynamics (MD) simulations are used here to predict the binding of different GBSSI mutants to a representative amylose fragment. The simulations were atomistic, with explicit solvent and docking, a method successfully used to understand the binding of wild-type GBSSI to amylose fragments. The binding of GBSSI to G5 (a pentasaccharide amylose fragment) is combined with free-energy calculations employing a thermodynamic integration method to predict the effects of mutations on enzyme activity. Ten GBSSI mutants with different enzyme activities were analyzed to find the structural and energy changes among different single amino-acid mutants and their possible relationship to starch characteristics. Comparing the structural changes and the relative binding free energy of G5 to the wild type GBSSI and GBSSI mutants, it was found that mutants with negative binding energy (lower than -2.0 kcal/mol) are more likely to have higher enzyme activity and amylose content compared to the wild type. This theoretical paper used simulations and robust free energy calculations to interpret in planta data with potential predictions as to what mutants might be generated to give desired properties. This study can be used to help develop grains with improved functional properties.

Molecular dynamics unlocks atomic level self-assembly of the exopolysaccharide matrix of water-treatment granular biofilms.

Biofilm formation, in which bacteria are embedded within an extracellular matrix, is the default form of microbial life in most natural and engineered habitats. In this work, atomistic molecular dynamics simulations were employed to examine the self-assembly of the polysaccharide Granulan to provide insight into the molecular interactions that lead to biofilm formation. Granulan is a major gel forming matrix component of granular microbial biofilms found in used-water treatment systems. Molecular dynamics simulations showed that Granulan forms an antiparallel double helix stabilized by complementary hydrogen bonds between the ß-glucosamine of one strand and the N-acetyl-ß-galactosamine-2-acetoamido-2-deoxy-α-galactopyranuronic pair of the other in both the presence and absence of Ca(2+). It is shown that Ca(2+) binds primarily to the carboxyl group of the terminal hexuronic acid of the sugar branch and that interactions between branches mediated by Ca(2+) suggest a possible mechanism for strengthening gels by facilitating interhelical bridging.

A dynamic pharmacophore drives the interaction between Psalmotoxin-1 and the putative drug target acid-sensing ion channel 1a.

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.

Challenges in the determination of the binding modes of non-standard ligands in X-ray crystal complexes.

Despite its central role in structure based drug design the determination of the binding mode (position, orientation and conformation in addition to protonation and tautomeric states) of small heteromolecular ligands in protein:ligand complexes based on medium resolution X-ray diffraction data is highly challenging. In this perspective we demonstrate how a combination of molecular dynamics simulations and free energy (FE) calculations can be used to correct and identify thermodynamically stable binding modes of ligands in X-ray crystal complexes. The consequences of inappropriate ligand structure, force field and the absence of electrostatics during X-ray refinement are highlighted. The implications of such uncertainties and errors for the validation of virtual screening and fragment-based drug design based on high throughput X-ray crystallography are discussed with possible solutions and guidelines.

Bisubstrate Ether-Linked Uridine-Peptide Conjugates as O-GlcNAc Transferase Inhibitors.

The O-linked ß-N-acetylglucosamine (O-GlcNAc) transferase (OGT) is a master regulator of installing O-GlcNAc onto serine or threonine residues on a multitude of target proteins. Numerous nuclear and cytosolic proteins of varying functional classes, including translational factors, transcription factors, signaling proteins, and kinases are OGT substrates. Aberrant O-GlcNAcylation of proteins is implicated in signaling in metabolic diseases such as diabetes and cancer. Selective and potent OGT inhibitors are valuable tools to study the role of OGT in modulating a wide range of effects on cellular functions. We report linear bisubstrate ether-linked uridine-peptide conjugates as OGT inhibitors with micromolar affinity. In vitro evaluation of the compounds revealed the importance of donor substrate, linker and acceptor substrate in the rational design of bisubstrate analogue inhibitors. Molecular dynamics simulations shed light on the binding of this novel class of inhibitors and rationalized the effect of amino acid truncation of acceptor peptide on OGT inhibition.

Platinum complexes act as shielding agents against virus infection.

We determine that the substitution-inert polynuclear platinum complex (PPC) TriplatinNC is an antiviral agent and protects cells from enterovirus 71 and human metapneumovirus infection. This protection occurs through the formation of adducts with cell-surface glycosaminoglycans. Our detailed mechanistic investigation demonstrates that TriplatinNC blocks viral entry by shielding cells from virus attack, opening new directions for metalloshielding antiviral drug development.

SARM1 is a metabolic sensor activated by an increased NMN/NAD+ ratio to trigger axon degeneration.

Axon degeneration is a central pathological feature of many neurodegenerative diseases. Sterile alpha and Toll/interleukin-1 receptor motif-containing 1 (SARM1) is a nicotinamide adenine dinucleotide (NAD+)-cleaving enzyme whose activation triggers axon destruction. Loss of the biosynthetic enzyme NMNAT2, which converts nicotinamide mononucleotide (NMN) to NAD+, activates SARM1 via an unknown mechanism. Using structural, biochemical, biophysical, and cellular assays, we demonstrate that SARM1 is activated by an increase in the ratio of NMN to NAD+ and show that both metabolites compete for binding to the auto-inhibitory N-terminal armadillo repeat (ARM) domain of SARM1. We report structures of the SARM1 ARM domain bound to NMN and of the homo-octameric SARM1 complex in the absence of ligands. We show that NMN influences the structure of SARM1 and demonstrate via mutagenesis that NMN binding is required for injury-induced SARM1 activation and axon destruction. Hence, SARM1 is a metabolic sensor responding to an increased NMN/NAD+ ratio by cleaving residual NAD+, thereby inducing feedforward metabolic catastrophe and axonal demise.

MyD88 TIR domain higher-order assembly interactions revealed by microcrystal electron diffraction and serial femtosecond crystallography.

MyD88 and MAL are Toll-like receptor (TLR) adaptors that signal to induce pro-inflammatory cytokine production. We previously observed that the TIR domain of MAL (MALTIR) forms filaments in vitro and induces formation of crystalline higher-order assemblies of the MyD88 TIR domain (MyD88TIR). These crystals are too small for conventional X-ray crystallography, but are ideally suited to structure determination by microcrystal electron diffraction (MicroED) and serial femtosecond crystallography (SFX). Here, we present MicroED and SFX structures of the MyD88TIR assembly, which reveal a two-stranded higher-order assembly arrangement of TIR domains analogous to that seen previously for MALTIR. We demonstrate via mutagenesis that the MyD88TIR assembly interfaces are critical for TLR4 signaling in vivo, and we show that MAL promotes unidirectional assembly of MyD88TIR. Collectively, our studies provide structural and mechanistic insight into TLR signal transduction and allow a direct comparison of the MicroED and SFX techniques.

Simplified Novel Muraymycin Analogues; using a Serine Template Strategy for Linking Key Pharmacophores.

The present status of antibiotic research requires the urgent invention of novel agents that act on multidrug-resistant bacteria. The World Health Organization has classified antibiotic-resistant bacteria into critical, high and medium priority according to the urgency of need for new antibiotics. Naturally occurring uridine-derived "nucleoside antibiotics" have shown promising activity against numerous priority resistant organisms by inhibiting the transmembrane protein MraY (translocase I), which is yet to be explored in a clinical context. The catalytic activity of MraY is an essential process for bacterial cell viability and growth including that of priority organisms. Muraymycins are one subclass of naturally occurring MraY inhibitors. Despite having potent antibiotic properties, the structural complexity of muraymycins advocates for simplified analogues as potential lead structures. Herein, we report a systematic structure-activity relationship (SAR) study of serine template-linked, simplified muraymycin-type analogues. This preliminary SAR lead study of serine template analogues successfully revealed that the complex structure of naturally occurring muraymycins could be easily simplified to afford bioactive scaffolds against resistant priority organisms. This study will pave the way for the development of novel antibacterial lead compounds based on a simplified serine template.