Inhibitor binding influences the protonation states of histidines in SARS-CoV-2 main protease

The main protease (M ^pro ) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is an attractive target for antiviral therapeutics. Recently, many high-resolution apo and inhibitor-bound structures of M ^pro , a cysteine protease, have been determined, facilitating structure-based drug design. M ^pro plays a central role in the viral life cycle by catalyzing the cleavage of SARS-CoV-2 polyproteins. In addition to the catalytic dyad His41-Cys145, M ^pro contains multiple histidines including His163, His164, and His172. The protonation states of these histidines and the catalytic nu-cleophile Cys145 have been debated in previous studies of SARS-CoV M ^pro , but have yet to be investigated for SARS-CoV-2. In this work we have used molecular dynamics simulations to determine the structural stability of SARS-CoV-2 M ^pro as a function of the protonation assignments for these residues. We simulated both the apo and inhibitor-bound enzyme and found that the conformational stability of the binding site, bound inhibitors, and the hydrogen bond networks of M ^pro are highly sensitive to these assignments. Additionally, the two inhibitors studied, the peptidomimetic N3 and an α -ketoamide, display distinct His41/His164 protonation-state-dependent stabilities. While the apo and the N3-bound systems favored N _δ (HD) and N _ϵ (HE) protonation of His41 and His164, respectively, the α -ketoamide was not stably bound in this state. Our results illustrate the importance of using appropriate histidine protonation states to accurately model the structure and dynamics of SARS-CoV-2 M ^pro in both the apo and inhibitor-bound states, a necessary prerequisite for drug-design efforts.

Supercomputer-Based Ensemble Docking Drug Discovery Pipeline with Application to Covid-19

We present a supercomputer-driven pipeline for in-silico drug discovery using enhanced sampling molecular dynamics (MD) and ensemble docking. We also describe preliminary results obtained for 23 systems involving eight protein targets of the proteome of SARS CoV-2. THe MD performed is temperature replica-exchange enhanced sampling, making use of the massively parallel supercomputing on the SUMMIT supercomputer at Oak Ridge National Laboratory, with which more than 1ms of enhanced sampling MD can be generated per day. We have ensemble docked repurposing databases to ten configurations of each of the 23 SARS CoV-2 systems using AutoDock Vina. We also demonstrate that using Autodock-GPU on SUMMIT, it is possible to perform exhaustive docking of one billion compounds in under 24 hours. Finally, we discuss preliminary results and planned improvements to the pipeline, including the use of quantum mechanical (QM), machine learning, and AI methods to cluster MD trajectories and rescore docking poses.

GPU-Accelerated Drug Discovery with Docking on the Summit Supercomputer: Porting, Optimization, and Application to COVID-19 Research

Protein-ligand docking is an in silico tool used to screen potential drug compounds for their ability to bind to a given protein receptor within a drug-discovery campaign. Experimental drug screening is expensive and time consuming, and it is desirable to carry out large scale docking calculations in a high-throughput manner to narrow the experimental search space. Few of the existing computational docking tools were designed with high performance computing in mind. Therefore, optimizations to maximize use of high-performance computational resources available at leadership-class computing facilities enables these facilities to be leveraged for drug discovery. Here we present the porting, optimization, and validation of the AutoDock-GPU program for the Summit supercomputer, and its application to initial compound screening efforts to target proteins of the SARS-CoV-2 virus responsible for the current COVID-19 pandemic.

Carotenoids promote lateral packing and condensation of lipid membranes

Carotenoids are pigment molecules that protect biomembranes against degradation and may be involved in the formation of functional bacterial membrane microdomains. Little is known on whether different types of carotenoids have different effects on the membrane or if there is any concentration dependence of these effects. In this work, we present results from molecular dynamics simulations of phospholipid bilayers containing different amounts of either β-carotene or zeaxanthin. Both β-carotene and zeaxanthin show the ability to laterally condense the membrane lipids and reduce their inter-leaflet interactions. With increasing concentrations, both carotenoids increase the bilayer thickness and rigidity. The results reveal that carotenoids have similar effects to cholesterol on regulating the behavior of fluid-phase membranes, suggesting that they could function as sterol substitutes and confirming their potential role in the formation of functional membrane domains.

Capturing Deuteration Effects in a Molecular Mechanics Force Field: Deuterated THF and the THF-Water Miscibility Gap

Deuteration is a common chemical modification used in conjunction with experiments such as neutron scattering, NMR, and Fourier-transform infrared for the study of molecular systems. Under the Born-Oppenheimer (BO) approximation, while the underlying potential energy surface remains unchanged by isotopic substitutions, isotopic substitution still alters intramolecular vibrations, which in turn may alter intermolecular interactions. Molecular mechanics (MM) force fields used in classical molecular dynamics simulations are assumed to represent local approximations of the BO potential energy surfaces, and hence, MD simulations using simple isotopic mass substitutions should capture BO-compatible isotope effects. However, standard MM force-field parameterizations do not directly fit to the local harmonic quantum mechanical (QM) Hessian that describes the BO surface, but rather to QM normal-modes and/or mass-dependent internal-coordinate derived distortion energies. Here, using tetrahydrofuran (THF)-water mixtures as our model system, we show that not only does a simple mass-substitution approach fail to capture an experimentally characterized deuteration effect (the loss of the closed-loop miscibility gap associated with the complete deuteration of THF) but also it is necessary to generate new MM force-field parameters that correctly describe isotopic dependent vibrations to capture the experimental deuteration effect. We show that the origin of this failure is a result of using mass-dependent features to fit the THF MM force field, which unintentionally biases the bonded terms of the force field to represent only the isotopologue used during the original force-field parameterization. In addition, we make use of our isotopologue-corrected force field for D₈THF to examine the molecular origins of the isotope-dependent loss of the THF-water miscibility gap.

A Minimal Membrane Metal Transport System: Dynamics and Energetics of mer Proteins

The mer operon in bacteria encodes a set of proteins and enzymes that impart resistance to environmental mercury toxicity by importing Hg2+ and reducing it to volatile Hg(0). Because the reduction occurs in the cytoplasm, mercuric ions must first be transported across the cytoplasmic membrane by one of a few known transporters. MerF is the smallest of these, containing only two transmembrane helices and two pairs of vicinal cysteines that coordinate mercuric ions. In this work, we use molecular dynamics simulations to characterize the dynamics of MerF in its apo and Hg2+ -bound states. We find that the apo state positions one of the cysteine pairs closer to the periplasmic side of the membrane, while in the bound state the same pair approaches the cytoplasmic side. This finding is consistent with the functional requirement of accepting Hg2+ from the periplasmic space, sequestering it on acceptance, and transferring it to the cytoplasm. Conformational changes in the TM helices facilitate the functional interaction of the two cysteine pairs. Free-energy calculations provide a barrier of 16 kcal/mol for the association of the periplasmic Hg2+ -bound protein MerP with MerF and 7 kcal/mol for the subsequent association of MerF's two cysteine pairs. Despite the significant conformational changes required to move the binding site across the membrane, coarse-grained simulations of multiple copies of MerF support the expectation that it functions as a monomer. Our results demonstrate how conformational changes and binding thermodynamics could lead to such a small membrane protein acting as an ion transporter. Published 2019. This article is a U.S. Government work and is in the public domain in the USA.

Structure-based group A streptococcal vaccine design: Helical wheel homology predicts antibody cross-reactivity among streptococcal M protein-derived peptides

Group A streptococcus (Strep A) surface M protein, an α-helical coiled-coil dimer, is a vaccine target and a major determinant of streptococcal virulence. The sequence-variable N-terminal region of the M protein defines the M type and also contains epitopes that promote opsonophagocytic killing of streptococci. Recent reports have reported considerable cross-reactivity among different M types, suggesting the prospect of identifying cross-protective epitopes that would constitute a broadly protective multivalent vaccine against Strep A isolates. Here, we have used a combination of immunological assays, structural biology, and cheminformatics to construct a recombinant M protein-based vaccine that included six Strep A M peptides that were predicted to elicit antisera that would cross-react with an additional 15 nonvaccine M types of Strep A. Rabbit antisera against this recombinant vaccine cross-reacted with 10 of the 15 nonvaccine M peptides. Two of the five nonvaccine M peptides that did not cross-react shared high sequence identity (≥50%) with the vaccine peptides, implying that high sequence identity alone was insufficient for cross-reactivity among the M peptides. Additional structural analyses revealed that the sequence identity at corresponding polar helical-wheel heptad sites between vaccine and nonvaccine peptides accurately distinguishes cross-reactive from non-cross-reactive peptides. On the basis of these observations, we developed a scoring algorithm based on the sequence identity at polar heptad sites. When applied to all epidemiologically important M types, this algorithm should enable the selection of a minimal number of M peptide-based vaccine candidates that elicit broadly protective immunity against Strep A.

Discovery of multidrug efflux pump inhibitors with a novel chemical scaffold

Multidrug efflux is a major contributor to antibiotic resistance in Gram-negative bacterial pathogens. Inhibition of multidrug efflux pumps is a promising approach for reviving the efficacy of existing antibiotics. Previously, inhibitors targeting both the efflux transporter AcrB and the membrane fusion protein AcrA in the Escherichia coli AcrAB-TolC efflux pump were identified. Here we use existing physicochemical property guidelines to generate a filtered library of compounds for computational docking. We then experimentally test the top candidate coumpounds using in vitro binding assays and in vivo potentiation assays in bacterial strains with controllable permeability barriers. We thus identify a new class of inhibitors of E. coli AcrAB-TolC. Six molecules with a shared scaffold were found to potentiate the antimicrobial activity of erythromycin and novobiocin in hyperporinated E. coli cells. Importantly, these six molecules were also active in wild-type strains of both Acinetobacter baumannii and Klebsiella pneumoniae, potentiating the activity of erythromycin and novobiocin up to 8-fold.

Conformational Dynamics of AcrA Govern Multidrug Efflux Pump Assembly

Multidrug efflux pumps of pathogenic, Gram-negative bacteria comprise an innate resistance mechanism and are key contributors to the emerging global pandemic of antibiotic resistance. Several increasingly detailed cryo-electron microscopy maps have been resolved of an entire efflux pump complex, AcrAB-TolC, resulting in atomistic structural models. Using a recent model, we have carried out nearly 40 μs of molecular dynamics simulations to study one of the key components of the protein complex AcrA, the membrane fusion protein that connects the inner-membrane-bound AcrB to the outer-membrane-bound TolC. We determined a three-dimensional potential of mean force (PMF) for AcrA, which displays two main conformational basins representing assembly competent and incompetent states. Corresponding experiments show that stabilizing mutations at an interdomain interface shift the dynamic equilibrium between these states to the incompetent one, disrupting pump assembly and function and resensitizing bacteria to existing antibiotics. The modulation of AcrA dynamics through pharmacological intervention therefore presents a promising route for the development of new antibiotics.

The AQUA-MER databases and aqueous speciation server: A web resource for multiscale modeling of mercury speciation

To assess the chemical reactivity, toxicity, and mobility of pollutants in the environment, knowledge of their species distributions is critical. Because their direct measurement is often infeasible, speciation modeling is widely adopted. Mercury (Hg) is a representative pollutant for which study of its speciation benefits from modeling. However, Hg speciation modeling is often hindered by a lack of reliable thermodynamic constants. Although computational chemistry (e.g., density functional theory [DFT]) can generate these constants, methods for directly coupling DFT and speciation modeling are not available. Here, we combine computational chemistry and continuum-scale modeling with curated online databases to ameliorate the problem of unreliable inputs to Hg speciation modeling. Our AQUA-MER databases and web server (https://aquamer.ornl.gov) provides direct speciation results by combining web-based interfaces to a speciation calculator, databases of thermodynamic constants, and a computational chemistry toolkit to estimate missing constants. Although Hg is presented as a concrete use case, AQUA-MER can also be readily applied to other elements. © 2019 Wiley Periodicals, Inc.

Horizontal transfer of a pathway for coumarate catabolism unexpectedly inhibits purine nucleotide biosynthesis

A microbe's ecological niche and biotechnological utility are determined by its specific set of co-evolved metabolic pathways. The acquisition of new pathways, through horizontal gene transfer or genetic engineering, can have unpredictable consequences. Here we show that two different pathways for coumarate catabolism failed to function when initially transferred into Escherichia coli. Using laboratory evolution, we elucidated the factors limiting activity of the newly acquired pathways and the modifications required to overcome these limitations. Both pathways required host mutations to enable effective growth with coumarate, but the necessary mutations differed. In one case, a pathway intermediate inhibited purine nucleotide biosynthesis, and this inhibition was relieved by single amino acid replacements in IMP dehydrogenase. A strain that natively contains this coumarate catabolism pathway, Acinetobacter baumannii, is resistant to inhibition by the relevant intermediate, suggesting that natural pathway transfers have faced and overcome similar challenges. Molecular dynamics simulation of the wild type and a representative single-residue mutant provide insight into the structural and dynamic changes that relieve inhibition. These results demonstrate how deleterious interactions can limit pathway transfer, that these interactions can be traced to specific molecular interactions between host and pathway, and how evolution or engineering can alleviate these limitations.

Ligand-Dependent Sodium Ion Dynamics within the A2A Adenosine Receptor: A Molecular Dynamics Study

Sodium ions have long been known to reduce the binding of agonists in many class-A GPCRs while having little effect on antagonist binding. Here, using long-time scale classical all-atom molecular dynamics simulations, we explore, in atomic detail, the motion of sodium ions within the ligand-binding pocket of the A_2A adenosine receptor (A2A-AR) both in the presence and absence of ligands and in the active and inactive state. We identify novel secondary ion binding sites within the pocket and find that the types of ion motions within the pocket are highly dependent on the presence and type of ligand within the pocket. Our results provide a first step toward developing a molecular understanding of the impact of sodium ions on class-A GPCRs.

Using Small-Angle Scattering Data and Parametric Machine Learning to Optimize Force Field Parameters for Intrinsically Disordered Proteins

Intrinsically disordered proteins (IDPs) and proteins with intrinsically disordered regions (IDRs) play important roles in many aspects of normal cell physiology, such as signal transduction and transcription, as well as pathological states, including Alzheimer's, Parkinson's, and Huntington's disease. Unlike their globular counterparts that are defined by a few structures and free energy minima, IDP/IDR comprise a large ensemble of rapidly interconverting structures and a corresponding free energy landscape characterized by multiple minima. This aspect has precluded the use of structural biological techniques, such as X-ray crystallography and nuclear magnetic resonance (NMR) for resolving their structures. Instead, low-resolution techniques, such as small-angle X-ray or neutron scattering (SAXS/SANS), have become a mainstay in characterizing coarse features of the ensemble of structures. These are typically complemented with NMR data if possible or computational techniques, such as atomistic molecular dynamics, to further resolve the underlying ensemble of structures. However, over the past 10–15 years, it has become evident that the classical, pairwise-additive force fields that have enjoyed a high degree of success for globular proteins have been somewhat limited in modeling IDP/IDR structures that agree with experiment. There has thus been a significant effort to rehabilitate these models to obtain better agreement with experiment, typically done by optimizing parameters in a piecewise fashion. In this work, we take a different approach by optimizing a set of force field parameters simultaneously, using machine learning to adapt force field parameters to experimental SAXS scattering profiles. We demonstrate our approach in modeling three biologically IDP ensembles based on experimental SAXS profiles and show that our optimization approach significantly improve force field parameters that generate ensembles in better agreement with experiment.

A Multifunctional Cosolvent Pair Reveals Molecular Principles of Biomass Deconstruction

The complex structure of plant cell walls resists chemical or biological degradation, challenging the breakdown of lignocellulosic biomass into renewable chemical precursors that could form the basis of future production of green chemicals and transportation fuels. Here, experimental and computational results reveal that the effect of the tetrahydrofuran (THF)-water cosolvents on the structure of lignin and on its interactions with cellulose in the cell wall drives multiple synergistic mechanisms leading to the efficient breakdown and fractionation of biomass into valuable chemical precursors. Molecular simulations show that THF-water is an excellent "theta" solvent, such that lignin dissociates from itself and from cellulose and expands to form a random coil. The expansion of the lignin molecules exposes interunit linkages, rendering them more susceptible to depolymerization by acid-catalyzed cleavage of aryl-ether bonds. Nanoscale infrared sensors confirm cosolvent-mediated molecular rearrangement of lignin in the cell wall of micrometer-thick hardwood slices and track the disappearance of lignin. At bulk scale, adding dilute acid to the cosolvent mixture liberates the majority of the hemicellulose and lignin from biomass, allowing unfettered access of cellulolytic enzymes to the remaining cellulose-rich material, allowing them to sustain high rates of hydrolysis to glucose without enzyme deactivation. Through this multiscale analysis, synergistic mechanisms for biomass deconstruction are identified, portending a paradigm shift toward first-principles design and evaluation of other cosolvent methods to realize low cost fuels and bioproducts.

Abstract 4584: Malignant pancreatic cancer cells respond to immune selection pressure to foster immunosuppression

Background:

Host immunity is relevant to pancreatic cancer (PC). Pancreatic tumors are infiltrated with a high number of immune cells, and they play a potential important role in disease evolution. However, there is no FDA approved immunotherapy for PC yet, suggesting that a deeper understanding of immune-cancer interaction is needed. The aim of this study was to investigate the effects of immune selection pressure on malignant pancreatic epithelial cells and on subsequent cancer immune evasion mechanisms.

Materials and Methods:

1 X105 mT3-2D murine pancreatic cancer cells (Kras+/G12D, p53+/-R172H (Boj et al. 2015)) were injected subcutaneously into syngeneic C57BL/6J (WT), B6.CB17-Prkdcscid/SzJ (SCID) and T cell-depleted WT mice (10 mice/group). Tumor growth was monitored over time. 1 cm3 tumors then were harvested, processed and studied. To consider the complex biology of pancreatic cancer immunology, we analyzed tumor specimens in their entirety by histopathological evaluation, multiple color flow cytometry, NanoString nCounter, and global proteomics profiling. Furthermore, to evaluate the influence of immune selection pressure on neoplastic cells, we selectively analyzed malignant epithelial cell gene expression in the cell line, and in cells isolated from tumors grown in WT and SCID mice, respectively, by RNAseq on mT3-2D-GFP+ fluorescence-activated cell sorter (FACS)-sorted cells. Also, whole exome sequencing was performed on WT tumors.

Results

Here we demonstrate that mT3-2D tumors engrafted in WT mice grow more slowly as compared with syngeneic SCID and T cell-depleted tumors, demonstrating that the reduced mT3-2D tumor growth rate in WT mice is T cell dependent. Putative cancer neoantigens were identified from whole exome sequences with pVACtools. Supporting evidence of tumor immunogenicity came from the observation that WT tumors are moderately infiltrated with T cells, by IHC and flow cytometry. Proteomics and NanoString nCounter analysis demonstrated the selective presence of myeloid suppressive cells in WT tumors, supporting their role in immune evasion. Notably, RNAseq showed selective upregulation of malignant epithelial cell-derived genes known to stimulate this immunosuppressive phenotype, including myeloid- and complement-related genes.

Conclusion

Despite the immunogenicity of the WT tumors and the infiltration of T cells in this PC model, T cell immunity incompletely controls tumor growth. Malignant epithelial cells mediate this cancer immunosuppressive phenotype in response to immune selection pressure, suggesting novel translational opportunities.

References:

Boj SF, Hwang C-I, Baker LA, Chio IIC, Engle DD, Corbo V, Jager M, Ponz-Sarvise M, Tiriac H, Spector MS, et al. 2015. Organoid models of human and mouse ductal pancreatic cancer. Cell 160:324–338.

Citation Format: Reham Ajina, Annie Zuo, Maha Moussa, Sarah J. Cooper, Yue Shen, Quentin R. Johnson, Jerry M. Parks, Jeremy C. Smith, Elana Fertig, Marta Catalfamo, Sandra A. Jablonski, Louis M. Weiner. Malignant pancreatic cancer cells respond to immune selection pressure to foster immunosuppression [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2019; 2019 Mar 29-Apr 3; Atlanta, GA. Philadelphia (PA): AACR; Cancer Res 2019;79(13 Suppl):Abstract nr 4584.

Environmental Mercury Chemistry - In Silico

Mercury (Hg) is a global environmental contaminant. Major anthropogenic sources of Hg emission include gold mining and the burning of fossil fuels. Once deposited in aquatic environments, Hg can undergo redox reactions, form complexes with ligands, and adsorb onto particles. It can also be methylated by microorganisms. Mercury, especially its methylated form methylmercury, can be taken up by organisms, where it bioaccumulates and biomagnifies in the food chain, leading to detrimental effects on ecosystem and human health. In support of the recently enforced Minamata Convention on Mercury, a legally binding international convention aimed at reducing the anthropogenic emission of-and human exposure to-Hg, its global biogeochemical cycle must be understood. Thus, a detailed understanding of the molecular-level interactions of Hg is crucial. The ongoing rapid development of hardware and methods has brought computational chemistry to a point that it can usefully inform environmental science. This is particularly true for Hg, which is difficult to handle experimentally due to its ultratrace concentrations in the environment and its toxicity. The current account provides a synopsis of the application of computational chemistry to filling several major knowledge gaps in environmental Hg chemistry that have not been adequately addressed experimentally. Environmental Hg chemistry requires defining the factors that determine the relative affinities of different ligands for Hg species, as they are critical for understanding its speciation, transformation and bioaccumulation in the environment. Formation constants and the nature of bonding have been determined computationally for environmentally relevant Hg(II) complexes such as chlorides, hydroxides, sulfides and selenides, in various physical phases. Quantum chemistry has been used to determine the driving forces behind the speciation of Hg with hydrochalcogenide and halide ligands. Of particular importance is the detailed characterization of solvation effects. Indeed, the aqueous phase reverses trends in affinities found computationally in the gas phase. Computation has also been used to investigate complexes of methylmercury with (seleno)amino acids, providing a molecular-level understanding of the toxicological antagonism between Hg and selenium (Se). Furthermore, evidence is emerging that ice surfaces play an important role in Hg transport and transformation in polar and alpine regions. Therefore, the diffusion of Hg and its ions through an idealized ice surface has been characterized. Microorganisms are major players in environmental mercury cycling. Some methylate inorganic Hg species, whereas others demethylate methylmercury. Quantum chemistry has been used to investigate catalytic mechanisms of enzymatic Hg methylation and demethylation. The complex interplay between the myriad chemical reactions and transport properties both in and outside microbial cells determines net biogeochemical cycling. Prospects for scaling up molecular work to obtain a mechanistic understanding of Hg cycling with comprehensive multiscale biogeochemical modeling are also discussed.

Identification of Binding Sites for Efflux Pump Inhibitors of the AcrAB-TolC Component AcrA

The overexpression of multidrug efflux pumps is an important mechanism of clinical resistance in Gram-negative bacteria. Recently, four small molecules were discovered that inhibit efflux in Escherichia coli and interact with the AcrAB-TolC efflux pump component AcrA. However, the binding site(s) for these molecules was not determined. Here, we combine ensemble docking and molecular dynamics simulations with tryptophan fluorescence spectroscopy, site-directed mutagenesis, and antibiotic susceptibility assays to probe binding sites and effects of binding of these molecules. We conclude that clorobiocin and SLU-258 likely bind at a site located between the lipoyl and β-barrel domains of AcrA.

Hydration-mediated stiffening of collective membrane dynamics by cholesterol

Hydration water governs the cholesterol-induced changes in collective headgroup dynamics in lipid bilayers.

Neutron scattering in the biological sciences: progress and prospects

The scattering of neutrons can be used to provide information on the structure and dynamics of biological systems on multiple length and time scales. Pursuant to a National Science Foundation-funded workshop in February 2018, recent developments in this field are reviewed here, as well as future prospects that can be expected given recent advances in sources, instrumentation and computational power and methods. Crystallography, solution scattering, dynamics, membranes, labeling and imaging are examined. For the extraction of maximum information, the incorporation of judicious specific deuterium labeling, the integration of several types of experiment, and interpretation using high-performance computer simulation models are often found to be particularly powerful.

Quantum Mechanical/Molecular Mechanical Analysis of the Catalytic Mechanism of Phosphoserine Phosphatase

Phosphoserine phosphatase (PSP), a member of the haloacid dehalogenase (HAD) superfamily that comprises the vast majority of phosphotransferases, is likely a steady-state regulator of the level of d-serine in the brain. The proposed catalytic cycle of PSP consists of a two-step mechanism: formation of a phospho-enzyme intermediate by phosphate transfer to Asp11 and its subsequent hydrolysis. Our combined quantum mechanical/molecular mechanical (QM/MM) calculations of the reaction pathways favour a dissociative mechanism of nucleophilic substitution via a trigonal-planar metaphosphate-like configuration for both steps, associated with proton transfer to the leaving group or from the nucleophile. This proton transfer is facilitated by active site residue Asp13 that acts as both a general base and a general acid. Free energy calculation on the reaction pathways further support the structural role of the enzymatic environment and the active site architecture. The choice of a proper reaction coordinate along which to bias the free energy calculations can be guided by a projection of the canonical reaction coordinate obtained from a chain-of-state optimisation onto important internal coordinates.

The importance of the membrane interface as the reference state for membrane protein stability

The insertion of nascent polypeptide chains into lipid bilayer membranes and the stability of membrane proteins crucially depend on the equilibrium partitioning of polypeptides. For this, the transfer of full sequences of amino-acid residues into the bilayer, rather than individual amino acids, must be understood. Earlier studies have revealed that the most likely reference state for partitioning very hydrophobic sequences is the membrane interface. We have used μs-scale simulations to calculate the interface-to-transmembrane partitioning free energies ΔG_S→TM for two hydrophobic carrier sequences in order to estimate the insertion free energy for all 20 amino acid residues when bonded to the center of a partitioning hydrophobic peptide. Our results show that prior single-residue scales likely overestimate the partitioning free energies of polypeptides. The correlation of ΔG_S→TM with experimental full-peptide translocon insertion data is high, suggesting an important role for the membrane interface in translocon-based insertion. The choice of carrier sequence greatly modulates the contribution of each single-residue mutation to the overall partitioning free energy. Our results demonstrate the importance of quantifying the observed full-peptide partitioning equilibrium, which is between membrane interface and transmembrane inserted, rather than combining individual water-to-membrane amino acid transfer free energies.

Jörg Langowski: his scientific legacy and the future it promises

Background

With the passing of Jörg Langowski 6 May 2017 in a sailplane accident, the scientific community was deprived of a strident and effective voice for DNA and chromatin molecular and computational biophysics, for open access publishing and for the creation of effective scientific research networks.

Methods

Here, after reviewing some of Jörg's key research contributions and ideas, we offer through the personal remembrance of his closest collaborators, a deep analysis of the major results of his research and the future directions they have engendered.

Conclusions

The legacy of Jörg Langowski has been to propel a way of viewing biological function that considers living systems as dynamic and in three dimensions. This physical view of biology that he pioneered is now, finally, becoming established also because of his great effort.

GPCR6A Is a Molecular Target for the Natural Products Gallate and EGCG in Green Tea

SCOPE

The molecular mechanisms whereby gallates in green tea exert metabolic effects are poorly understood.

METHODS AND RESULTS

We found that GPRC6A, a multi-ligand-sensing G-protein-coupled receptor that regulates energy metabolism, sex hormone production, and prostate cancer progression, is a target for gallates. Sodium gallate (SG), gallic acid (GA) > ethyl gallate (EG) > octyl gallate (OG) dose dependently activated ERK in HEK-293 cells transfected with GPRC6A but not in non-transfected controls. SG also stimulated insulin secretion in β-cells isolated from wild-type mice similar to the endogenous GPRC6A ligands, osteocalcin (Ocn) and testosterone (T). Side-chain additions to create OG resulted in loss of GPRC6A agonist activity. Another component of green tea, epigallocatechin 3-gallate (EGCG), dose-dependently inhibited Ocn activation of GPRC6A in HEK-293 cells transfected with GPRC6A and blocked the effect of Ocn in stimulating glucose production in CH10T1/2 cells. Using structural models of the venus fly trap (VFT) and 7-transmembrane (7-TM) domains of GPRC6A, calculations suggest that l-amino acids and GA bind to the VFT, whereas EGCG is calculated to bind to sites in both the VFT and 7-TM.

CONCLUSION

GA and EGCG have offsetting agonist and antagonist effects on GPRC6A that may account for the variable metabolic effect of green tea consumption.

Ensemble Docking in Drug Discovery

Ensemble docking corresponds to the generation of an "ensemble" of drug target conformations in computational structure-based drug discovery, often obtained by using molecular dynamics simulation, that is used in docking candidate ligands. This approach is now well established in the field of early-stage drug discovery. This review gives a historical account of the development of ensemble docking and discusses some pertinent methodological advances in conformational sampling.

Determination of Dynamical Heterogeneity from Dynamic Neutron Scattering of Proteins

Motional displacements of hydrogen (H) in proteins can be measured using incoherent neutron-scattering methods. These displacements can also be calculated numerically using data from molecular dynamics simulations. An enormous amount of data on the average mean-square motional displacement (MSD) of H as a function of protein temperature, hydration, and other conditions has been collected. H resides in a wide spectrum of sites in a protein. Some H are tightly bound to molecular chains, and the H motion is dictated by that of the chain. Other H are quite independent. As a result, there is a distribution of motions and MSDs of H within a protein that is denoted dynamical heterogeneity. The goal of this paper is to incorporate a distribution of MSDs into models of the H incoherent intermediate scattering function, I(Q,t), that is calculated and observed. The aim is to contribute information on the distribution as well as on the average MSD from comparison of the models with simulations and experiment. For example, we find that simulations of I(Q,t) in lysozyme are well reproduced if the distribution of MSDs is bimodal with two broad peaks rather than a single broad peak.

Dynamic Neutron Scattering by Biological Systems

Dynamic neutron scattering directly probes motions in biological systems on femtosecond to microsecond timescales. When combined with molecular dynamics simulation and normal mode analysis, detailed descriptions of the forms and frequencies of motions can be derived. We examine vibrations in proteins, the temperature dependence of protein motions, and concepts describing the rich variety of motions detectable using neutrons in biological systems at physiological temperatures. New techniques for deriving information on collective motions using coherent scattering are also reviewed.

Catalysis of Ground State cis[Formula: see text] trans Isomerization of Bacteriorhodopsin's Retinal Chromophore by a Hydrogen-Bond Network

For the photocycle of the membrane protein bacteriorhodopsin to proceed efficiently, the thermal 13-cis to all-trans back-isomerization of the retinal chromophore must return the protein to its resting state on a time-scale of milliseconds. Here, we report on quantum mechanical/molecular mechanical energy calculations examining the structural and energetic determinants of the retinal cis-trans isomerization in the protein environment. The results suggest that a hydrogen-bonded network consisting of the retinal Schiff base, active site amino acid residues, and water molecules can stabilize the twisted retinal, thus reducing the intrinsic energy cost of the cis-trans thermal isomerization barrier.

Double J Stents in the Treatment of Gynaecological Injury to the Ureter

Ureteric injury is a recognized complication of hysterectomy and may present with obstruction or fistula.

Between 1987 and 1989 in Oxford nine patients with 10 injured ureters underwent attempted retrograde placement of double J stents. Three patients had successful outcomes and one patient with bilateral ureteric obstruction required reimplantation of the right ureter after successful stenting of the left ureter. One patient required removal of a stent due to irritation but her fistula eventually closed. In three patients placement was unsuccessful and in one patient injury to the bladder base prevented the ureteric orifices from being seen and hence stenting was not possible. Thus five of these 10 injured ureters were managed successfully with double J stents.

We advocate the initial use of double J stents in gynaecological ureteric injury. This approach is simple and may cure the fistula. If it is unsuccessful, subsequent reimplantation is not hindered.

Effects of carotenoids on lipid bilayers

Carotenoids have been found to be important in improving the integrity of biomembranes in eukaryotes. However, the molecular details of how carotenoids modulate the physical properties of biomembranes are unknown. To this end, we have conducted a series of molecular dynamics simulations of different biologically-relevant membranes in the presence of carotenoids. The carotenoid effect on the membrane was found to be specific to the identity of the carotenoid and the composition of the membrane itself. Therefore, different classes of carotenoids produce a different effect on the membrane, and different membrane phases are affected differently by carotenoids. It is apparent from our data that carotenoids do trigger the bilayer to become thinner. The mechanism by which this occurs depends on two competing factors, the ability of the lipid tails of opposing monolayers to either (1) compress or (2) interdigitate as the bilayer condenses. Indeed, carotenoids directly influence the physical properties via these two mechanisms, thus compacting the bilayer. However, the degree to which these competing mechanisms are utilized depends on the bilayer phase and the carotenoid identity.

"To Be or Not to Be" Protonated: Atomic Details of Human Carbonic Anhydrase-Clinical Drug Complexes by Neutron Crystallography and Simulation

Human carbonic anhydrases (hCAs) play various roles in cells, and have been drug targets for decades. Sequence similarities of hCA isoforms necessitate designing specific inhibitors, which requires detailed structural information for hCA-inhibitor complexes. We present room temperature neutron structures of hCA II in complex with three clinical drugs that provide in-depth analysis of drug binding, including protonation states of the inhibitors, hydration water structure, and direct visualization of hydrogen-bonding networks in the enzyme's active site. All sulfonamide inhibitors studied bind to the Zn metal center in the deprotonated, anionic, form. Other chemical groups of the drugs can remain neutral or be protonated when bound to hCA II. MD simulations have shown that flexible functional groups of the inhibitors may alter their conformations at room temperature and occupy different sub-sites. This study offers insights into the design of specific drugs to target cancer-related hCA isoform IX.

Ensemble docking to difficult targets in early-stage drug discovery: Methodology and application to fibroblast growth factor 23

Ensemble docking is now commonly used in early-stage in silico drug discovery and can be used to attack difficult problems such as finding lead compounds which can disrupt protein-protein interactions. We give an example of this methodology here, as applied to fibroblast growth factor 23 (FGF23), a protein hormone that is responsible for regulating phosphate homeostasis. The first small-molecule antagonists of FGF23 were recently discovered by combining ensemble docking with extensive experimental target validation data (Science Signaling, 9, 2016, ra113). Here, we provide a detailed account of how ensemble-based high-throughput virtual screening was used to identify the antagonist compounds discovered in reference (Science Signaling, 9, 2016, ra113). Moreover, we perform further calculations, redocking those antagonist compounds identified in reference (Science Signaling, 9, 2016, ra113) that performed well on drug-likeness filters, to predict possible binding regions. These predicted binding modes are rescored with the molecular mechanics Poisson-Boltzmann surface area (MM/PBSA) approach to calculate the most likely binding site. Our findings suggest that the antagonist compounds antagonize FGF23 through the disruption of protein-protein interactions between FGF23 and fibroblast growth factor receptor (FGFR).

Dynamics of the lignin glass transition

Despite lignin being a heterogenous polyphenolic, its glass transition obeys well-established polymer theory concepts.

Relationship between lignocellulosic biomass dissolution and physicochemical properties of ionic liquids composed of 3-methylimidazolium cations and carboxylate anions

Experimental and simulation studies identify 1-allyl-3-methylimidazolium formate as an efficient biomass solvent, mainly due to strong interactions with hemicellulose.

Polycystin-1 interacts with TAZ to stimulate osteoblastogenesis and inhibit adipogenesis

The molecular mechanisms that transduce the osteoblast response to physical forces in the bone microenvironment are poorly understood. Here, we used genetic and pharmacological experiments to determine whether the polycystins PC1 and PC2 (encoded by Pkd1 and Pkd2) and the transcriptional coactivator TAZ form a mechanosensing complex in osteoblasts. Compound-heterozygous mice lacking 1 copy of Pkd1 and Taz exhibited additive decrements in bone mass, impaired osteoblast-mediated bone formation, and enhanced bone marrow fat accumulation. Bone marrow stromal cells and osteoblasts derived from these mice showed impaired osteoblastogenesis and enhanced adipogenesis. Increased extracellular matrix stiffness and application of mechanical stretch to multipotent mesenchymal cells stimulated the nuclear translocation of the PC1 C-terminal tail/TAZ (PC1-CTT/TAZ) complex, leading to increased runt-related transcription factor 2-mediated (Runx2-mediated) osteogenic and decreased PPARγ-dependent adipogenic gene expression. Using structure-based virtual screening, we identified a compound predicted to bind to PC2 in the PC1:PC2 C-terminal tail region with helix:helix interaction. This molecule stimulated polycystin- and TAZ-dependent osteoblastogenesis and inhibited adipogenesis. Thus, we show that polycystins and TAZ integrate at the molecular level to reciprocally regulate osteoblast and adipocyte differentiation, indicating that the polycystins/TAZ complex may be a potential therapeutic target to increase bone mass.

Organosolv-Water Cosolvent Phase Separation on Cellulose and its Influence on the Physical Deconstruction of Cellulose: A Molecular Dynamics Analysis

Deconstruction of cellulose is crucial for the chemical conversion of lignocellulose into fuel/bioproduct precursors. Recently, a water-organosolv cosolvent system (THF-water) has been shown to both phase-separate on cellulose surfaces and partially deconstruct Avicel  (cellulose) in the absence of acid. Here we employ molecular dynamics simulations to determine whether other common water-organosolv cosolvent systems (acetone, ethanol, and γ-valerolactone) exhibit phase separation at cellulose surface and whether this alters a purely physical cellulose dissociation pathway. Despite finding varied degrees of phase-separation of organosolv on cellulose surfaces, physical dissociation is not enhanced. Interestingly, however, the total amount the median water-cellulose contact lifetimes increases for the cosolvent systems in the order of THF > acetone > ethanol > γ-valerolactone. Together our results indicate two points: a purely physical process for deconstruction of cellulose is unlikely for these cosolvents, and in THF-water, unlike γ-valerolactone- (and some concentrations of acetone and ethanol) water cosolvents, a significant fraction of surface water is slowed. This slowing may be of importance in enhancing chemical deconstruction of cellulose, as it permits an increase in potential THF-water-cellulose reactions, even while the amount of water near cellulose is decreased.

Dynamics at a Peptide-TiO2 Anatase (101) Interface

The interface between biological matter and inorganic materials is a widely investigated research topic due to possible applications in biomedicine and nanotechnology. In this context, the molecular level adsorption mechanism that drives specific recognition between small peptide sequences and inorganic surfaces represents an important topic likely to provide much information useful for designing bioderived materials. Here, we investigate the dynamics at the interface between a Ti-binding peptide sequence (AMRKLPDAPGMHC) and a TiO₂ anatase surface by using molecular dynamics (MD) simulations. In the simulations the adsorption mechanism is characterized by diffusion of the peptide from the bulk water phase toward the TiO₂ surface, followed by the anchoring of the peptide to the surface. The anchoring is mediated by the interfacial water layers by means of the charged groups of the side chains of the peptide. The peptide samples anchored and dissociated states from the surface and its conformation is not affected by the surface when anchored.

Modeling of the Passive Permeation of Mercury and Methylmercury Complexes Through a Bacterial Cytoplasmic Membrane

Cellular uptake and export are important steps in the biotransformation of mercury (Hg) by microorganisms. However, the mechanisms of transport across biological membranes remain unclear. Membrane-bound transporters are known to be relevant, but passive permeation may also be involved. Inorganic Hg^II and methylmercury ([CH₃Hg^II]⁺) are commonly complexed with thiolate ligands. Here, we have performed extensive molecular dynamics simulations of the passive permeation of Hg^II and [CH₃Hg^II]⁺ complexes with thiolate ligands through a model bacterial cytoplasmic membrane. We find that the differences in free energy between the individual complexes in bulk water and at their most favorable position within the membrane are ∼2 kcal mol^-1. We provide a detailed description of the molecular interactions that drive the membrane crossing process. Favorable interactions with carbonyl and tail groups of phospholipids stabilize Hg-containing solutes in the tail-head interface region of the membrane. The calculated permeability coefficients for the neutral compounds CH₃S-Hg^II-SCH₃ and CH₃Hg^II-SCH₃ are on the order of 10^-5 cm s^-1. We conclude that small, nonionized Hg-containing species can permeate readily through cytoplasmic membranes.

Dynamical Transition of Collective Motions in Dry Proteins

Water is widely assumed to be essential for protein dynamics and function. In particular, the well-documented "dynamical" transition at ∼200  K, at which the protein changes from a rigid, nonfunctional form to a flexible, functional state, as detected in hydrogenated protein by incoherent neutron scattering, requires hydration. Here, we report on coherent neutron scattering experiments on perdeuterated proteins and reveal that a transition occurs in dry proteins at the same temperature resulting primarily from the collective heavy-atom motions. The dynamical transition discovered is intrinsic to the energy landscape of dry proteins.

Community detection in sequence similarity networks based on attribute clustering

Networks are powerful tools for the presentation and analysis of interactions in multi-component systems. A commonly studied mesoscopic feature of networks is their community structure, which arises from grouping together similar nodes into one community and dissimilar nodes into separate communities. Here, the community structure of protein sequence similarity networks is determined with a new method: Attribute Clustering Dependent Communities (ACDC). Sequence similarity has hitherto typically been quantified by the alignment score or its expectation value. However, pair alignments with the same score or expectation value cannot thus be differentiated. To overcome this deficiency, the method constructs, for pair alignments, an extended alignment metric, the link attribute vector, which includes the score and other alignment characteristics. Rescaling components of the attribute vectors qualitatively identifies a systematic variation of sequence similarity within protein superfamilies. The problem of community detection is then mapped to clustering the link attribute vectors, selection of an optimal subset of links and community structure refinement based on the partition density of the network. ACDC-predicted communities are found to be in good agreement with gold standard sequence databases for which the "ground truth" community structures (or families) are known. ACDC is therefore a community detection method for sequence similarity networks based entirely on pair similarity information. A serial implementation of ACDC is available from https://cmb.ornl.gov/resources/developments.