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

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.

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.

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

In this study, we have used [H, N] NMR spectroscopy to investigate the interactions of the trinuclear platinum anticancer drug triplatin () (1,0,1 or BBR3464) with site-specific sulfated and carboxylated disaccharides. Specifically, the disaccharides GlcNS(6)-GlcA () and GlcNS(6)-IdoA(2S) () are useful models of longer-chain glycosaminoglycans (GAGs) such as heparan sulfate (HS). For both the reactions of N with and , equilibrium conditions were achieved more slowly (65 h) compared to the reaction with the monosaccharide GlcNS(6S) (9 h).

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.

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.

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 (MAL) forms filaments in vitro and induces formation of crystalline higher-order assemblies of the MyD88 TIR domain (MyD88). 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).

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.

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.

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.

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.

Automated Topology Builder Version 3.0: Prediction of Solvation Free Enthalpies in Water and Hexane.

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 (Δ G) and solvation free enthalpies in the apolar solvent hexane (Δ G) is presented.

Predicting the Prevalence of Alternative Warfarin Tautomers in Solution.

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.

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.

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.

Optimization of Empirical Force Fields by Parameter Space Mapping: A Single-Step Perturbation Approach.

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.

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.

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

Human Group IIA phospholipase A (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.

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.

The characterization of modified starch branching enzymes: toward the control of starch chain-length distributions.

Starch is a complex branched glucose polymer whose branch molecular weight distribution (the chain-length distribution, CLD) influences nutritionally important properties such as digestion rate. Chain-stopping in starch biosynthesis is by starch branching enzyme (SBE). Site-directed mutagenesis was used to modify SBEIIa from Zea mays (mSBEIIa) to produce mutants, each differing in a single conserved amino-acid residue.

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.

Comparative Occupancy Analysis (CoOAn) - A Straightforward and Directly Applicable 3D-QSAR Formalism to Extract Molecular Features Obligatory for Designing Potent Leads.

A simple and directly applicable 3D-QSAR method, termed Comparative Occupancy Analysis (CoOAn), has been developed. The method is based on the comparison of local occupancies of fragments of an aligned set of molecules in a 3D-grid space. The formalism commendably extracts the crucial position-specific molecular features and correlates them quantitatively to their biological endpoints.

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.

Missing fragments: detecting cooperative binding in fragment-based drug design.

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.

An Automated Force Field Topology Builder (ATB) and Repository: Version 1.0.

The Automated force field Topology Builder (ATB, http://compbio.biosci.uq.