AlphaFold accelerated discovery of psychotropic agonists targeting the trace amine-associated receptor 1.

Artificial intelligence is revolutionizing protein structure prediction, providing unprecedented opportunities for drug design. To assess the potential impact on ligand discovery, we compared virtual screens using protein structures generated by the AlphaFold machine learning method and traditional homology modeling. More than 16 million compounds were docked to models of the trace amine-associated receptor 1 (TAAR1), a G protein-coupled receptor of unknown structure and target for treating neuropsychiatric disorders.

Identification of 14 Known Drugs as Inhibitors of the Main Protease of SARS-CoV-2.

A consensus virtual screening protocol has been applied to ca. 2000 approved drugs to seek inhibitors of the main protease (M) of SARS-CoV-2, the virus responsible for COVID-19. 42 drugs emerged as top candidates, and after visual analyses of the predicted structures of their complexes with M, 17 were chosen for evaluation in a kinetic assay for M inhibition.

Identification of 14 Known Drugs as Inhibitors of the Main Protease of SARS-CoV-2.

A consensus virtual screening protocol has been applied to ca. 2000 approved drugs to seek inhibitors of the main protease (M ) of SARS-CoV-2, the virus responsible for COVID-19. 42 drugs emerged as top candidates, and after visual analyses of the predicted structures of their complexes with M , 17 were chosen for evaluation in a kinetic assay for M inhibition.

Absolute Free Energy of Binding Calculations for Macrophage Migration Inhibitory Factor in Complex with a Druglike Inhibitor.

Calculation of the absolute free energy of binding (Δ) for a complex in solution is challenging owing to the need for adequate configurational sampling and an accurate energetic description, typically with a force field (FF). In this study, Monte Carlo (MC) simulations with improved side-chain and backbone sampling are used to assess Δ for the complex of a druglike inhibitor (MIF180) with the protein macrophage migration inhibitory factor (MIF) using free energy perturbation (FEP) calculations. For comparison, molecular dynamics (MD) simulations were employed as an alternative sampling method for the same system.

Computation of protein-ligand binding free energies using quantum mechanical bespoke force fields.

A quantum mechanical bespoke molecular mechanics force field is derived for the L99A mutant of T4 lysozyme and used to compute absolute binding free energies of six benzene analogs to the protein. Promising agreement between theory and experiment highlights the potential for future use of system-specific force fields in computer-aided drug design.

Robust Free Energy Perturbation Protocols for Creating Molecules in Solution.

Accurate methods to estimate free energies play an important role for studying diverse condensed-phase problems in chemistry and biochemistry. The most common methods used in conjunction with molecular dynamics (MD) and Monte Carlo statistical mechanics (MC) simulations are free energy perturbation (FEP) and thermodynamic integration (TI). For common applications featuring the conversion of one molecule to another, simulations are run in stages or multiple "λ-windows" to promote convergence of the results.

Enhanced Monte Carlo Methods for Modeling Proteins Including Computation of Absolute Free Energies of Binding.

The generation of a complete ensemble of geometrical configurations is required to obtain reliable estimations of absolute binding free energies by alchemical free energy methods. Molecular dynamics (MD) is the most popular sampling method, but the representation of large biomolecular systems may be incomplete owing to energetic barriers that impede efficient sampling of the configurational space. Monte Carlo (MC) methods can possibly overcome this issue by adapting the attempted movement sizes to facilitate transitions between alternative local-energy minima.

LigParGen web server: an automatic OPLS-AA parameter generator for organic ligands.

The accurate calculation of protein/nucleic acid-ligand interactions or condensed phase properties by force field-based methods require a precise description of the energetics of intermolecular interactions. Despite the progress made in force fields, small molecule parameterization remains an open problem due to the magnitude of the chemical space; the most critical issue is the estimation of a balanced set of atomic charges with the ability to reproduce experimental properties. The LigParGen web server provides an intuitive interface for generating OPLS-AA/1.

Ecoupling server: A tool to compute and analyze electronic couplings.

Electron transfer processes are often studied through the evaluation and analysis of the electronic coupling (EC). Since most standard QM codes do not provide readily such a measure, additional, and user-friendly tools to compute and analyze electronic coupling from external wave functions will be of high value. The first server to provide a friendly interface for evaluation and analysis of electronic couplings under two different approximations (FDC and GMH) is presented in this communication.

New Monte Carlo Based Technique To Study DNA-Ligand Interactions.

We present a new all-atom Monte Carlo technique capable of performing quick and accurate DNA-ligand conformational sampling. In particular, and using the PELE software as a frame, we have introduced an additional force field, an implicit solvent, and an anisotropic network model to effectively map the DNA energy landscape. With these additions, we successfully generated DNA conformations for a test set composed of six DNA fragments of A-DNA and B-DNA.

Porphyrin Binding to Gun4 Protein, Facilitated by a Flexible Loop, Controls Metabolite Flow through the Chlorophyll Biosynthetic Pathway.

In oxygenic phototrophs, chlorophylls, hemes, and bilins are synthesized by a common branched pathway. Given the phototoxic nature of tetrapyrroles, this pathway must be tightly regulated, and an important regulatory role is attributed to magnesium chelatase enzyme at the branching between the heme and chlorophyll pathway. Gun4 is a porphyrin-binding protein known to stimulate in vitro the magnesium chelatase activity, but how the Gun4-porphyrin complex acts in the cell was unknown.

A theoretical multiscale treatment of protein-protein electron transfer: The ferredoxin/ferredoxin-NADP(+) reductase and flavodoxin/ferredoxin-NADP(+) reductase systems.

In the photosynthetic electron transfer (ET) chain, two electrons transfer from photosystem I to the flavin-dependent ferredoxin-NADP(+) reductase (FNR) via two sequential independent ferredoxin (Fd) electron carriers. In some algae and cyanobacteria (as Anabaena), under low iron conditions, flavodoxin (Fld) replaces Fd as single electron carrier. Extensive mutational studies have characterized the protein-protein interaction in FNR/Fd and FNR/Fld complexes.

Direct Measurement of the Nanomechanical Stability of a Redox Protein Active Site and Its Dependence upon Metal Binding.

The structural basis of the low reorganization energy of cupredoxins has long been debated. These proteins reconcile a conformationally heterogeneous and exposed metal-chelating site with the highly rigid copper center required for efficient electron transfer. Here we combine single-molecule mechanical unfolding experiments with statistical analysis and computer simulations to show that the metal-binding region of apo-azurin is mechanically flexible and that high mechanical stability is imparted by copper binding.

Conformational response to ligand binding in phosphomannomutase2: insights into inborn glycosylation disorder.

The most common glycosylation disorder is caused by mutations in the gene encoding phosphomannomutase2, producing a disease still without a cure. Phosphomannomutase2, a homodimer in which each chain is composed of two domains, requires a bisphosphate sugar (either mannose or glucose) as activator, opening a possible drug design path for therapeutic purposes. The crystal structure of human phosphomannomutase2, however, lacks bound substrate and a key active site loop.

Atomic level rendering of DNA-drug encounter.

Computer simulations have been demonstrated to be important for unraveling atomic mechanisms in biological systems. In this study, we show how combining unbiased molecular dynamic simulations with appropriate analysis tools can successfully describe metal-based drug interactions with DNA. To elucidate the noncovalent affinity of cisplatin's family to DNA, we performed extensive all-atom molecular dynamics simulations (3.

Tryptogalinin is a tick Kunitz serine protease inhibitor with a unique intrinsic disorder.

Background: A salivary proteome-transcriptome project on the hard tick Ixodes scapularis revealed that Kunitz peptides are the most abundant salivary proteins. Ticks use Kunitz peptides (among other salivary proteins) to combat host defense mechanisms and to obtain a blood meal. Most of these Kunitz peptides, however, remain functionally uncharacterized, thus limiting our knowledge about their biochemical interactions.

H-bond network optimization in protein-protein complexes: are all-atom force field scores enough?

Structural prediction of protein-protein complexes given the structures of the two interacting compounds in their unbound state is a key problem in biophysics. In addition to the problem of sampling of near-native orientations, one of the modeling main difficulties is to discriminate true from false positives. Here, we present a hierarchical protocol for docking refinement able to discriminate near native poses from a group of docking candidates.