Identification of Triazolopyrimidinyl Scaffold SARS-CoV-2 Papain-Like Protease (PL) Inhibitor.

The global impact of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and its companion disease, COVID-19, has reminded us of the importance of basic coronaviral research. In this study, a comprehensive approach using molecular docking, in vitro assays, and molecular dynamics simulations was applied to identify potential inhibitors for SARS-CoV-2 papain-like protease (PL), a key and underexplored viral enzyme target. A focused protease inhibitor library was initially created and molecular docking was performed using CmDock software (v0.

Automatic Assembly and Calibration of Models of Enzymatic Reactions Based on Ordinary Differential Equations.

Enzymatic reactions have been studied for more than a 100 years. Indeed, isolation of enzymes from biological materials is no longer the main source of enzymes today, as they are now largely produced using recombinant technology, or can even be synthesized from scratch. Studies of the three-dimensional structures of enzymes can provide answers to many questions, but the kinetics of enzymatic reactions is the only method that can lead to better understanding of their function.

CHARMM Force-Field Parameters for Morphine, Heroin, and Oliceridine, and Conformational Dynamics of Opioid Drugs.

Opioid drug binding to specialized G protein-coupled receptors (GPCRs) can lead to analgesia upon activation via downstream G protein signaling and to severe side effects via activation of the β-arrestin signaling pathway. Knowledge of how different opioid drugs interact with receptors is essential, as it can inform and guide the design of safer therapeutics. We performed quantum and classical mechanical computations to explore the potential energy landscape of four opioid drugs: morphine and its derivatives heroin and fentanyl and for the unrelated oliceridine.

Molecular Dynamics Simulations Reveal Interactions of an IgG1 Antibody With Selected Fc Receptors.

In a survey of novel interactions between an IgG1 antibody and different Fcγ receptors (FcγR), molecular dynamics simulations were performed of interactions of monoclonal antibody involved complexes with FcγRs. Free energy simulations were also performed of isolated wild-type and substituted Fc regions bound to FcγRs with the aim of assessing their relative binding affinities. Two different free energy calculation methods, Molecular Mechanical/Generalized Born Molecular Volume (MM/GBMV) and Bennett Acceptance Ratio (BAR), were used to evaluate the known effector substitution G236A that is known to selectively increase antibody dependent cellular phagocytosis.

Potential Energy Function for Fentanyl-Based Opioid Pain Killers.

Opioids are molecules whose binding to specialized G-Protein Coupled Receptors (GPCRs) triggers a signaling cascade that leads to the downregulation of pain pathways. Binding of an opioid to the membrane-embedded GPCR occurs when the opioid molecule is protonated, which provides a potential strategy to design nontoxic opioids that are protonated and bind to the GPCR only at the low pH of injured or inflamed tissue. Excellent model systems to study protonation-dependent binding of opioids to GPCRs are fentanyl, which is protonated and binds to the GPCR at both physiological and low pH, and the fluorinated fentanyl derivative NFEPP, which is protonated and binds to the GPCR only at low pH.

Loop Grafting between Similar Local Environments for Fc-Silent Antibodies.

Reduction of the affinity of the fragment crystallizable (Fc) region with immune receptors by substitution of one or a few amino acids, known as Fc-silencing, is an established approach to reduce the immune effector functions of monoclonal antibody therapeutics. This approach to Fc-silencing, however, is problematic as it can lead to instability and immunogenicity of the developed antibodies. We evaluated loop grafting as a novel approach to Fc-silencing in which the Fc loops responsible for immune receptor binding were replaced by loops of up to 20 amino acids from similar local environments in other human and mouse antibodies.

Insight into Inhibitor Binding in the Eukaryotic Proteasome: Computations of the 20S CP.

A combination of molecular dynamics (MD) simulations and computational analyses uncovers structural features that may influence substrate passage and exposure to the active sites within the proteolytic chamber of the 20S proteasome core particle (CP). MD simulations of the CP reveal relaxation dynamics in which the CP slowly contracts over the 54 ns sampling period. MD simulations of the SyringolinA (SylA) inhibitor within the proteolytic B 1 ring chamber of the CP indicate that favorable van der Waals and electrostatic interactions account for the predominant association of the inhibitor with the walls of the proteolytic chamber.

Simulations of NPC1(NTD):NPC2 Protein Complex Reveal Cholesterol Transfer Pathways.

The Niemann Pick type C (NPC) proteins, NPC1 and NPC2, are involved in the lysosomal storage disease, NPC disease. The formation of a NPC1⁻NPC2 protein⁻protein complex is believed to be necessary for the transfer of cholesterol and lipids out of the late endosomal (LE)/lysosomal (Lys) compartments. Mutations in either NPC1 or NPC2 can lead to an accumulation of cholesterol and lipids in the LE/Lys, the primary phenotype of the NPC disease.

Crystal Structures of Diaryliodonium Fluorides and Their Implications for Fluorination Mechanisms.

The radiofluorination of diaryliodonium salts is of value for producing radiotracers for positron emission tomography. We report crystal structures for two diaryliodonium fluorides. Whereas diphenyliodonium fluoride (1 a) exists as a tetramer bridged by four fluoride ions, 2-methylphenyl(phenyl)iodonium fluoride (2 a) forms a fluoride-bridged dimer that is further halogen bonded to two other monomers.

Role of magnesium ions in the reaction mechanism at the interface between Tm1631 protein and its DNA ligand.

A protein, Tm1631 from the hyperthermophilic organism Thermotoga maritima belongs to a domain of unknown function protein family. It was predicted that Tm1631 binds with the DNA and that the Tm1631-DNA complex is an endonuclease repair system with a DNA repair function (Konc et al. PLoS Comput Biol 9(11): e1003341, 2013).

Molecular dynamics to enhance structure-based virtual screening on cathepsin B.

Molecular dynamics (MD) and molecular docking are commonly used to study molecular interactions in drug discovery. Most docking approaches consider proteins as rigid, which can decrease the accuracy of predicted docked poses. Therefore MD simulations can be used prior to docking to add flexibility to proteins.

Rapid identification of atypical tetracyclines using tandem mass spectrometric fragmentation patterns.

Rationale: When applying biosynthetic engineering approaches at the early stages of drug discovery, e.g. aiming to develop novel tetracycline analogues, target compounds are generally produced by engineered microorganisms in low yields.

Structural insight into the unique binding properties of pyridylethanol(phenylethyl)amine inhibitor in human CYP51.

Sterol 14α-demethylase (CYP51) is the main drug target for the treatment of fungal infections. The discovery of new efficient fungal CYP51 inhibitors requires an understanding of the structural requirements for selectivity for the fungal over the human ortholog. In this study, a binding mode of the pyridylethanol(phenylethyl)amine type CYP51 inhibitor to the human ortholog was determined at the atomic level.

Efficiently computing pathway free energies: New approaches based on chain-of-replica and Non-Boltzmann Bennett reweighting schemes.

Background: Accurately modeling condensed phase processes is one of computation's most difficult challenges. Include the possibility that conformational dynamics may be coupled to chemical reactions, where multiscale (i.e.

Common force field thermodynamics of cholesterol.

Four different force fields are examined for dynamic characteristics using cholesterol as a case study. The extent to which various types of internal degrees of freedom become thermodynamically relevant is evaluated by means of principal component analysis. More complex degrees of freedom (angle bending, dihedral rotations) show a trend towards force field independence.

Structure-based function prediction of uncharacterized protein using binding sites comparison.

A challenge in structural genomics is prediction of the function of uncharacterized proteins. When proteins cannot be related to other proteins of known activity, identification of function based on sequence or structural homology is impossible and in such cases it would be useful to assess structurally conserved binding sites in connection with the protein's function. In this paper, we propose the function of a protein of unknown activity, the Tm1631 protein from Thermotoga maritima, by comparing its predicted binding site to a library containing thousands of candidate structures.

Ligand binding site identification by higher dimension molecular dynamics.

We propose a new molecular dynamics (MD) protocol to identify the binding site of a guest within a host. The method utilizes a four spatial (4D) dimension representation of the ligand allowing for rapid and efficient sampling within the receptor. We applied the method to two different model receptors characterized by diverse structural features of the binding site and different ligand binding affinities.

The binding mode of second-generation sulfonamide inhibitors of MurD: clues for rational design of potent MurD inhibitors.

A series of optimized sulfonamide derivatives was recently reported as novel inhibitors of UDP-N-acetylmuramoyl-L-alanine:D-glutamate ligase (MurD). These are based on naphthalene-N-sulfonyl-D-glutamic acid and have the D-glutamic acid replaced with rigidified mimetics. Here we have defined the binding site of these novel ligands to MurD using (1)H/(13)C heteronuclear single quantum correlation.

Correlating protein hot spot surface analysis using ProBiS with simulated free energies of protein-protein interfacial residues.

A protocol was developed for the computational determination of the contribution of interfacial amino acid residues to the free energy of protein-protein binding. Thermodynamic integration, based on molecular dynamics simulation in CHARMM, was used to determine the free energy associated with single point mutations to glycine in a protein-protein interface. The hot spot amino acids found in this way were then correlated to structural similarity scores detected by the ProBiS algorithm for local structural alignment.

Replica exchanging self-guided Langevin dynamics for efficient and accurate conformational sampling.

This work presents a replica exchanging self-guided Langevin dynamics (RXSGLD) simulation method for efficient conformational searching and sampling. Unlike temperature-based replica exchanging simulations, which use high temperatures to accelerate conformational motion, this method uses self-guided Langevin dynamics (SGLD) to enhance conformational searching without the need to elevate temperatures. A RXSGLD simulation includes a series of SGLD simulations, with simulation conditions differing in the guiding effect and/or temperature.

Double locking of an Escherichia coli promoter by two repressors prevents premature colicin expression and cell lysis.

The synthesis of Eschericha coli colicins is lethal to the producing cell and is repressed during normal growth by the LexA transcription factor, which is the master repressor of the SOS system for repair of DNA damage. Following DNA damage, LexA is inactivated and SOS repair genes are induced immediately, but colicin production is delayed and induced only in terminally damaged cells. The cause of this delay is unknown.

Comparison of Three Chain-of-States Methods: Nudged Elastic Band and Replica Path with Restraints or Constraints.

Chain-of-state methods are becoming important tools in studying the chemical reaction mechanisms, especially for biomacromolecules. In this article, three chain-of-state methods, nudged elastic band (NEB) method and the replica path method with restraints or constraints, were tested and compared using three model systems with various sizes and at different levels of theory: alanine dipeptide isomerization, β-alanine intramolecular condensation, and the matrix metalloproteinase 2 inhibition mechanism. The levels of theory used to describe the three model systems include molecular mechanics (MM), quantum mechanics (QM), and combined quantum mechanics and molecular mechanics (QM/MM).

ENZO: a web tool for derivation and evaluation of kinetic models of enzyme catalyzed reactions.

We describe a web tool ENZO (Enzyme Kinetics), a graphical interface for building kinetic models of enzyme catalyzed reactions. ENZO automatically generates the corresponding differential equations from a stipulated enzyme reaction scheme. These differential equations are processed by a numerical solver and a regression algorithm which fits the coefficients of differential equations to experimentally observed time course curves.

MSCALE: A General Utility for Multiscale Modeling.

The combination of theoretical models of macromolecules that exist at different spatial and temporal scales has become increasingly important for addressing complex biochemical problems. This work describes the extension of concurrent multiscale approaches, introduces a general framework for carrying out calculations, and describes its implementation into the CHARMM macromolecular modeling package. This functionality, termed MSCALE, generalizes both the additive and subtractive multiscale scheme (e.

Conformational structure and energetics of 2-methylphenyl(2'-methoxyphenyl)iodonium chloride: evidence for solution clusters.

Diaryliodonium salts allow the efficient incorporation of cyclotron-produced [(18)F]fluoride ions into electron-rich and electron-deficient arenes to provide potential radiotracers for molecular imaging in vivo with positron emission tomography (PET). This process (ArI(+)Ar'+(18)F(-)→Ar(18)F+Ar'I) is still not well understood mechanistically. To better understand this and similar reactions, it would be valuable to understand the structures of diaryliodonium salts in organic media, where the reactions are typically conducted.