Performance of Relative Binding Free Energy Calculations on an Automatically Generated Dataset of Halogen-Deshalogen Matched Molecular Pairs.

In this study, we generated a matched molecular pair dataset of halogen/deshalogen compounds with reliable binding affinity data and structural binding mode information from public databases. The workflow includes automated system preparation and setup of free energy perturbation relative binding free energy calculations. We demonstrate the suitability of these datasets to investigate the performance of molecular mechanics force fields and molecular simulation algorithms for the purpose of in silico affinity predictions in lead optimization.

Large-Scale Assessment of Binding Free Energy Calculations in Active Drug Discovery Projects.

Accurate ranking of compounds with regards to their binding affinity to a protein using computational methods is of great interest to pharmaceutical research. Physics-based free energy calculations are regarded as the most rigorous way to estimate binding affinity. In recent years, many retrospective studies carried out both in academia and industry have demonstrated its potential.

Large-scale in vitro functional testing and novel variant scoring via protein modeling provide insights into alkaline phosphatase activity in hypophosphatasia.

Hypophosphatasia (HPP) is a rare metabolic disorder characterized by low tissue-nonspecific alkaline phosphatase (TNSALP) typically caused by ALPL gene mutations. HPP is heterogeneous, with clinical presentation correlating with residual TNSALP activity and/or dominant-negative effects (DNE). We measured residual activity and DNE for 155 ALPL variants by transient transfection and TNSALP enzymatic activity measurement.

Thermodynamic integration network approach to ion transport through protein channels: Perspectives and limits.

We present a molecular dynamics simulation study of alkali metal cation transport through the double-helical and the head-to-head conformers of the gramicidin ion channel. Our approach is based on a thermodynamic integration network, which consists of a sequence of transport reactions, absolute free energies of solvation and cycles of alchemical transmutations of the ions. In this manner, we can reliably estimate free energies and their statistical errors via a least-squares method without imposing external forces on the system.

Best of Two Worlds? How MD Simulations of Amphiphilic Helical Peptides in Membranes Can Complement Data from Oriented Solid-State NMR.

The membrane alignment of helical amphiphilic peptides in oriented phospholipid bilayers can be obtained as ensemble and time averages from solid state H NMR by fitting the quadrupolar splittings to ideal α-helices. At the same time, molecular dynamics (MD) simulations can provide atomistic insight into peptide-membrane systems. Here, we evaluate the potential of MD simulations to complement the experimental NMR data that is available on three exemplary systems: the natural antimicrobial peptide PGLa and the two designer-made peptides MSI-103 and KIA14, whose sequences were derived from PGLa.

Predicting resistance of clinical Abl mutations to targeted kinase inhibitors using alchemical free-energy calculations.

The therapeutic effect of targeted kinase inhibitors can be significantly reduced by intrinsic or acquired resistance mutations that modulate the affinity of the drug for the kinase. In cancer, the majority of missense mutations are rare, making it difficult to predict their impact on inhibitor affinity. This complicates the practice of precision medicine, pairing of patients with clinical trials, and development of next-generation inhibitors.

Prediction of Accurate Binding Modes Using Combination of Classical and Accelerated Molecular Dynamics and Free-Energy Perturbation Calculations: An Application to Toxicity Studies.

Estimating the correct binding modes of ligands in protein-ligand complexes is crucial not only in the drug discovery process but also for elucidating potential toxicity mechanisms. In the current paper, we propose a computational modeling workflow using the combination of docking, classical molecular dynamics (cMD), accelerated molecular dynamics (aMD) and free-energy perturbation (FEP+ protocol) for identification of possible ligand binding modes. It was applied for investigation of selected perfluorocarboxyl acids (PFCAs) in the PPARγ nuclear receptor.

Validation of tautomeric and protomeric binding modes by free energy calculations. A case study for the structure based optimization of D-amino acid oxidase inhibitors.

Optimization of fragment size D-amino acid oxidase (DAAO) inhibitors was investigated using a combination of computational and experimental methods. Retrospective free energy perturbation (FEP) calculations were performed for benzo[d]isoxazole derivatives, a series of known inhibitors with two potential binding modes derived from X-ray structures of other DAAO inhibitors. The good agreement between experimental and computed binding free energies in only one of the hypothesized binding modes strongly support this bioactive conformation.

A Water-Bridged Cysteine-Cysteine Redox Regulation Mechanism in Bacterial Protein Tyrosine Phosphatases.

The emergence of multidrug-resistant () strains highlights the need to develop more efficacious and potent drugs. However, this goal is dependent on a comprehensive understanding of virulence protein effectors at the molecular level. Here, we used a post-expression cysteine (Cys)-to-dehydrolanine (Dha) chemical editing strategy to identify a water-mediated motif that modulates accessibility of the protein tyrosine phosphatase A (PtpA) catalytic pocket.

Automated Transition State Search and Its Application to Diverse Types of Organic Reactions.

Transition state search is at the center of multiple types of computational chemical predictions related to mechanistic investigations, reactivity and regioselectivity predictions, and catalyst design. The process of finding transition states in practice is, however, a laborious multistep operation that requires significant user involvement. Here, we report a highly automated workflow designed to locate transition states for a given elementary reaction with minimal setup overhead.

Free Energy Perturbation Calculations of the Thermodynamics of Protein Side-Chain Mutations.

Protein side-chain mutation is fundamental both to natural evolutionary processes and to the engineering of protein therapeutics, which constitute an increasing fraction of important medications. Molecular simulation enables the prediction of the effects of mutation on properties such as binding affinity, secondary and tertiary structure, conformational dynamics, and thermal stability. A number of widely differing approaches have been applied to these predictions, including sequence-based algorithms, knowledge-based potential functions, and all-atom molecular mechanics calculations.

Acylguanidine Beta Secretase 1 Inhibitors: A Combined Experimental and Free Energy Perturbation Study.

A series of acylguanidine beta secretase 1 (BACE1) inhibitors with modified scaffold and P3 pocket substituent was synthesized and studied with free energy perturbation (FEP) calculations. The resulting molecules showed potencies in enzymatic BACE1 inhibition assays up to 1 nM. The correlation between the predicted activity from the FEP calculations and the experimental activity was good for the P3 pocket substituents.

Predicting the Effect of Amino Acid Single-Point Mutations on Protein Stability-Large-Scale Validation of MD-Based Relative Free Energy Calculations.

The stability of folded proteins is critical to their biological function and for the efficacy of protein therapeutics. Predicting the energetic effects of protein mutations can improve our fundamental understanding of structural biology, the molecular basis of diseases, and possible routes to addressing those diseases with biological drugs. Identifying the effect of single amino acid point mutations on the thermodynamic equilibrium between the folded and unfolded states of a protein can pinpoint residues of critical importance that should be avoided in the process of improving other properties (affinity, solubility, viscosity, etc.

A Comprehensive Computational Study of the Interaction between Human Serum Albumin and Fullerenes.

Human serum albumin (HSA) is the most abundant blood plasma protein, which transports fatty acids, hormones, and drugs. We consider nanoparticle-HSA interactions by investigating the binding of HSA with three fullerene analogs. Long MD simulations, quantum mechanical (fragment molecular orbital, energy decomposition analysis, atoms-in-molecules), and free energy methods elucidated the binding mechanism in these complexes.

Accurate Binding Free Energy Predictions in Fragment Optimization.

Predicting protein-ligand binding free energies is a central aim of computational structure-based drug design (SBDD)--improved accuracy in binding free energy predictions could significantly reduce costs and accelerate project timelines in lead discovery and optimization. The recent development and validation of advanced free energy calculation methods represents a major step toward this goal. Accurately predicting the relative binding free energy changes of modifications to ligands is especially valuable in the field of fragment-based drug design, since fragment screens tend to deliver initial hits of low binding affinity that require multiple rounds of synthesis to gain the requisite potency for a project.

Stereoselective synthesis and molecular modeling of chiral cyclopentanes.

The reaction of 3-methyseleno-2-methylselenomethyl-propene with benzyl 2,3-anhydro-4-O-triflyl-β-L-ribopyranoside provides a major convenient enantiomeric product of 1-methylene-(benzyl3,4-dideoxy-α-D-arabinopyranoso)-[3,4-c]-cyclopentane, with benzyl-2,3-anhydro-4-deoxy-4-C-(2-methyl- propen-3-yl)-α-D-lyxopyranoside as a minor product. While the reaction of 3-methyseleno-2-[methylselenomethyl]-propene with benzyl 2,3-anhydro-4-O-triflyl-α-D-ribopyranoside produces a good yield of benzyl-2,3-anhydro-4-deoxy-4-C-(2-methylpropen-3-yl)-α-D-lyxo-pyranoside. Molecular modeling and molecular dynamics simulations indicate that the intermediate in the reaction of the β-L sugar frequently occupies an optimal conformation that leads to the formation of cyclopentane, while the intermediate in the reaction of the α-D sugar has a very small probability.

Pyrido[4,3-e][1,2,4]triazolo[4,3-a]pyrazines as Selective, Brain Penetrant Phosphodiesterase 2 (PDE2) Inhibitors.

A novel series of pyrido[4,3-e][1,2,4]triazolo[4,3-a]pyrazines is reported as potent PDE2/PDE10 inhibitors with drug-like properties. Selectivity for PDE2 was obtained by introducing a linear, lipophilic moiety on the meta-position of the phenyl ring pending from the triazole. The SAR and protein flexibility were explored with free energy perturbation calculations.

Accurate and reliable prediction of relative ligand binding potency in prospective drug discovery by way of a modern free-energy calculation protocol and force field.

Designing tight-binding ligands is a primary objective of small-molecule drug discovery. Over the past few decades, free-energy calculations have benefited from improved force fields and sampling algorithms, as well as the advent of low-cost parallel computing. However, it has proven to be challenging to reliably achieve the level of accuracy that would be needed to guide lead optimization (∼5× in binding affinity) for a wide range of ligands and protein targets.

3D hydrophobic moment vectors as a tool to characterize the surface polarity of amphiphilic peptides.

The interaction of membranes with peptides and proteins is largely determined by their amphiphilic character. Hydrophobic moments of helical segments are commonly derived from their two-dimensional helical wheel projections, and the same is true for β-sheets. However, to the best of our knowledge, there exists no method to describe structures in three dimensions or molecules with irregular shape.

Ligand photo-isomerization triggers conformational changes in iGluR2 ligand binding domain.

Neurological glutamate receptors bind a variety of artificial ligands, both agonistic and antagonistic, in addition to glutamate. Studying their small molecule binding properties increases our understanding of the central nervous system and a variety of associated pathologies. The large, oligomeric multidomain membrane protein contains a large and flexible ligand binding domains which undergoes large conformational changes upon binding different ligands.

Fragment orbital based description of charge transfer in peptides including backbone orbitals.

Charge transfer in peptides and proteins can occur on different pathways, depending on the energetic landscape as well as the coupling between the involved orbitals. Since details of the mechanism and pathways are difficult to access experimentally, different modeling strategies have been successfully applied to study these processes in the past. These can be based on a simple empirical pathway model, efficient tight binding type atomic orbital Hamiltonians or ab initio and density functional calculations.

Controlling biological activity with light: diarylethene-containing cyclic peptidomimetics.

Photobiological processes in nature are usually triggered by nonpeptidic chromophores or by modified side chains. A system is presented in which the polypeptide backbone itself can be conformationally switched by light. An amino acid analogue was designed and synthesized based on a reversibly photoisomerizable diarylethene scaffold.

Ground state spectroscopy of hydroxyquinolines: evidence for the formation of protonated species in water-rich dioxane-water mixtures.

We have recently used 6-, 7-, and 8-hydroxyquinolines (HQs) as fluorescent probes to study the binding mechanism in one of the drug binding sites of human serum albumin. In the present work we study the absorption spectra of the HQ molecules in neat and binary mixtures of dioxane and water in order to identify the different tautomeric species in the ground state. This study should help in identifying the environment in nanocavities of macromolecules when HQs are used as local reporters.

Charge transfer in E. coli DNA photolyase: understanding polarization and stabilization effects via QM/MM simulations.

We study fast hole transfer events in E. coli DNA photolyase, a key step in the photoactivation process, using a multiscale computational method that combines nonadiabatic propagation schemes and linear-scaling quantum chemical methods with molecular mechanics force fields. This scheme allows us to follow the time-dependent evolution of the electron hole in an unbiased fashion; that is, no assumptions about hole wave function localization, time scale separation, or adiabaticity of the process have to be made beforehand.