A computationally designed binding mode flip leads to a novel class of potent tri-vector cyclophilin inhibitors

Sire: An interoperability engine for prototyping algorithms and exchanging information between molecular simulation programs

Sire is a Python/C++ library that is used both to prototype new algorithms and as an interoperability engine for exchanging information between molecular simulation programs. It provides a collection of file parsers and information converters that together make it easier to combine and leverage the functionality of many other programs and libraries. This empowers researchers to use sire to write a single script that can, for example, load a molecule from a PDBx/mmCIF file via Gemmi, perform SMARTS searches via RDKit, parameterize molecules using BioSimSpace, run GPU-accelerated molecular dynamics via OpenMM, and then display the resulting dynamics trajectory in a NGLView Jupyter notebook 3D molecular viewer. This functionality is built on by BioSimSpace, which uses sire's molecular information engine to interconvert with programs such as GROMACS, NAMD, Amber, and AmberTools for automated molecular parameterization and the running of molecular dynamics, metadynamics, and alchemical free energy workflows. Sire comes complete with a powerful molecular information search engine, plus trajectory loading and editing, analysis, and energy evaluation engines. This, when combined with an in-built computer algebra system, gives substantial flexibility to researchers to load, search for, edit, and combine molecular information from multiple sources and use that to drive novel algorithms by combining functionality from other programs. Sire is open source (GPL3) and is available via conda and at a free Jupyter notebook server at https://try.openbiosim.org. Sire is supported by the not-for-profit OpenBioSim community interest company.

Alchemical Free Energy Workflows for the Computation of Protein-Ligand Binding Affinities

Alchemical free energy methods can be used for the efficient computation of relative binding free energies during preclinical drug discovery stages. In recent years, this has been facilitated further by the implementation of workflows that enable non-experts to quickly and consistently set up the required simulations. Given the correct input structures, workflows handle the difficult aspects of setting up perturbations, including consistently defining the perturbable molecule, its atom mapping and topology generation, perturbation network generation, running of the simulations via different sampling methods, and analysis of the results. Different academic and commercial workflows are discussed, including FEW, FESetup, FEPrepare, CHARMM-GUI, Transformato, PMX, QLigFEP, TIES, ProFESSA, PyAutoFEP, BioSimSpace, FEP+, Flare, and Orion. These workflows differ in various aspects, such as mapping algorithms or enhanced sampling methods. Some workflows can accommodate more than one molecular dynamics (MD) engine and use external libraries for tasks. Differences between workflows can present advantages for different use cases, however a lack of interoperability of the workflows' components hinders systematic comparisons.

Cyclophilin D in Mitochondrial Dysfunction: A Key Player in Neurodegeneration?

Mitochondrial dysfunction plays a pivotal role in numerous complex diseases. Understanding the molecular mechanisms by which the "powerhouse of the cell" turns into the "factory of death" is an exciting yet challenging task that can unveil new therapeutic targets. The mitochondrial matrix protein CyPD is a peptidylprolyl cis-trans isomerase involved in the regulation of the permeability transition pore (mPTP). The mPTP is a multi-conductance channel in the inner mitochondrial membrane whose dysregulated opening can ultimately lead to cell death and whose involvement in pathology has been extensively documented over the past few decades. Moreover, several mPTP-independent CyPD interactions have been identified, indicating that CyPD could be involved in the fine regulation of several biochemical pathways. To further enrich the picture, CyPD undergoes several post-translational modifications that regulate both its activity and interaction with its clients. Here, we will dissect what is currently known about CyPD and critically review the most recent literature about its involvement in neurodegenerative disorders, focusing on Alzheimer's Disease and Parkinson's Disease, supporting the notion that CyPD could serve as a promising therapeutic target for the treatment of such conditions. Notably, significant efforts have been made to develop CyPD-specific inhibitors, which hold promise for the treatment of such complex disorders.

Comparison of Receptor-Ligand Restraint Schemes for Alchemical Absolute Binding Free Energy Calculations

Alchemical absolute binding free energy calculations are of increasing interest in drug discovery. These calculations require restraints between the receptor and ligand to restrict their relative positions and, optionally, orientations. Boresch restraints are commonly used, but they must be carefully selected in order to sufficiently restrain the ligand and to avoid inherent instabilities. Applying multiple distance restraints between anchor points in the receptor and ligand provides an alternative framework without inherent instabilities which may provide convergence benefits by more strongly restricting the relative movements of the receptor and ligand. However, there is no simple method to calculate the free energy of releasing these restraints due to the coupling of the internal and external degrees of freedom of the receptor and ligand. Here, a method to rigorously calculate free energies of binding with multiple distance restraints by imposing intramolecular restraints on the anchor points is proposed. Absolute binding free energies for the human macrophage migration inhibitory factor/MIF180, system obtained using a variety of Boresch restraints and rigorous and nonrigorous implementations of multiple distance restraints are compared. It is shown that several multiple distance restraint schemes produce estimates in good agreement with Boresch restraints. In contrast, calculations without orientational restraints produce erroneously favorable free energies of binding by up to approximately 4 kcal mol-1. These approaches offer new options for the deployment of alchemical absolute binding free energy calculations.

Discovery and molecular basis of subtype-selective cyclophilin inhibitors

Although cyclophilins are attractive targets for probing biology and therapeutic intervention, no subtype-selective cyclophilin inhibitors have been described. We discovered novel cyclophilin inhibitors from the in vitro selection of a DNA-templated library of 256,000 drug-like macrocycles for cyclophilin D (CypD) affinity. Iterated macrocycle engineering guided by ten X-ray co-crystal structures yielded potent and selective inhibitors (half maximal inhibitory concentration (IC50) = 10 nM) that bind the active site of CypD and also make novel interactions with non-conserved residues in the S2 pocket, an adjacent exo-site. The resulting macrocycles inhibit CypD activity with 21- to >10,000-fold selectivity over other cyclophilins and inhibit mitochondrial permeability transition pore opening in isolated mitochondria. We further exploited S2 pocket interactions to develop the first cyclophilin E (CypE)-selective inhibitor, which forms a reversible covalent bond with a CypE S2 pocket lysine, and exhibits 30- to >4,000-fold selectivity over other cyclophilins. These findings reveal a strategy to generate isoform-selective small-molecule cyclophilin modulators, advancing their suitability as targets for biological investigation and therapeutic development.

Non-Immunosuppressive Cyclophilin Inhibitors

Cyclophilins, enzymes with peptidyl-prolyl cis/trans isomerase activity, are relevant to a large variety of biological processes. The most abundant member of this enzyme family, cyclophilin A, is the cellular receptor of the immunosuppressive drug cyclosporine A (CsA). As a consequence of the pathophysiological role of cyclophilins, particularly in viral infections, there is a broad interest in cyclophilin inhibition devoid of immunosuppressive activity. This Review first gives an introduction into the physiological and pathophysiological roles of cyclophilins. The presentation of non-immunosuppressive cyclophilin inhibitors will commence with drugs based on chemical modifications of CsA. The naturally occurring macrocyclic sanglifehrins have become other lead structures for cyclophilin-inhibiting drugs. Finally, de novo designed compounds, whose structures are not derived from or inspired by natural products, will be presented. Relevant synthetic concepts will be discussed, but the focus will also be on biochemical studies, structure-activity relationships, and clinical studies.

Small-molecule inhibitors of cyclophilin D as potential therapeutics in mitochondria-related diseases

Cyclophilin D (CypD) is a key regulator of mitochondrial permeability transition pore (mPTP) opening. This pathophysiological phenomenon is associated with the development of several human diseases, including ischemia-reperfusion injury and neurodegeneration. Blocking mPTP opening through CypD inhibition could be a novel and promising therapeutic approach for these conditions. While numerous CypD inhibitors have been discovered to date, none have been introduced into clinical practice, mostly owing to their high toxicity, unfavorable pharmacokinetics, and low selectivity for CypD over other cyclophilins. This review summarizes current knowledge of CypD inhibitors, with a particular focus on small-molecule compounds with regard to their in vitro activity, their selectivity for CypD, and their binding mode within the enzyme's active site. Finally, approaches for improving the molecular design of CypD inhibitors are discussed.

Repurposing of cyclophilin A inhibitors as broad-spectrum antiviral agents

Cyclophilin A (CypA) is linked to diverse human diseases including viral infections. With the worldwide emergence of severe acute respiratory coronavirus 2 (SARS-CoV-2), drug repurposing has been highlighted as a strategy with the potential to speed up antiviral development. Because CypA acts as a proviral component in hepatitis C virus, coronavirus and HIV, its inhibitors have been suggested as potential treatments for these infections. Here, we review the structure of cyclosporin A and sanglifehrin A analogs as well as synthetic micromolecules inhibiting CypA; and we discuss their broad-spectrum antiviral efficacy in the context of the virus lifecycle.

Implementation of the QUBE Force Field in SOMD for High-Throughput Alchemical Free-Energy Calculations

The quantum mechanical bespoke (QUBE) force-field approach has been developed to facilitate the automated derivation of potential energy function parameters for modeling protein-ligand binding. To date, the approach has been validated in the context of Monte Carlo simulations of protein-ligand complexes. We describe here the implementation of the QUBE force field in the alchemical free-energy calculation molecular dynamics simulation package SOMD. The implementation is validated by demonstrating the reproducibility of absolute hydration free energies computed with the QUBE force field across the SOMD and GROMACS software packages. We further demonstrate, by way of a case study involving two series of non-nucleoside inhibitors of HIV-1 reverse transcriptase, that the availability of QUBE in a modern simulation package that makes efficient use of graphics processing unit acceleration will facilitate high-throughput alchemical free-energy calculations.

Hybrid Alchemical Free Energy/Machine-Learning Methodology for the Computation of Hydration Free Energies

A methodology that combines alchemical free energy calculations (FEP) with machine learning (ML) has been developed to compute accurate absolute hydration free energies. The hybrid FEP/ML methodology was trained on a subset of the FreeSolv database and retrospectively shown to outperform most submissions from the SAMPL4 competition. Compared to pure machine-learning approaches, FEP/ML yields more precise estimates of free energies of hydration and requires a fraction of the training set size to outperform standalone FEP calculations. The ML-derived correction terms are further shown to be transferable to a range of related FEP simulation protocols. The approach may be used to inexpensively improve the accuracy of FEP calculations and to flag molecules which will benefit the most from bespoke force field parametrization efforts.

Sampling Conformational Changes of Bound Ligands Using Nonequilibrium Candidate Monte Carlo and Molecular Dynamics

Flexible ligands often have multiple binding modes or bound conformations that differ by rotation of a portion of the molecule around internal rotatable bonds. Knowledge of these binding modes is important for understanding the interactions stabilizing the ligand in the binding pocket, and other studies indicate it is important for calculating accurate binding affinities. In this work, we use a hybrid molecular dynamics (MD)/nonequilibrium candidate Monte Carlo (NCMC) method to sample the different binding modes of several flexible ligands and also to estimate the population distribution of the modes. The NCMC move proposal is divided into three parts. The flexible part of the ligand is alchemically turned off by decreasing the electrostatics and steric interactions gradually, followed by rotating the rotatable bond by a random angle and then slowly turning the ligand back on to its fully interacting state. The alchemical steps prior to and after the move proposal help the surrounding protein and water atoms in the binding pocket relax around the proposed ligand conformation and increase move acceptance rates. The protein-ligand system is propagated using classical MD in between the NCMC proposals. Using this MD/NCMC method, we were able to correctly reproduce the different binding modes of inhibitors binding to two kinase targets-c-Jun N-terminal kinase-1 and cyclin-dependent kinase 2-at a much lower computational cost compared to conventional MD and umbrella sampling. This method is available as a part of the BLUES software package.

Dynamic design: manipulation of millisecond timescale motions on the energy landscape of cyclophilin A

Proteins need to interconvert between many conformations in order to function, many of which are formed transiently, and sparsely populated. Particularly when the lifetimes of these states approach the millisecond timescale, identifying the relevant structures and the mechanism by which they interconvert remains a tremendous challenge. Here we introduce a novel combination of accelerated MD (aMD) simulations and Markov state modelling (MSM) to explore these 'excited' conformational states. Applying this to the highly dynamic protein CypA, a protein involved in immune response and associated with HIV infection, we identify five principally populated conformational states and the atomistic mechanism by which they interconvert. A rational design strategy predicted that the mutant D66A should stabilise the minor conformations and substantially alter the dynamics, whereas the similar mutant H70A should leave the landscape broadly unchanged. These predictions are confirmed using CPMG and R1ρ solution state NMR measurements. By efficiently exploring functionally relevant, but sparsely populated conformations with millisecond lifetimes in silico, our aMD/MSM method has tremendous promise for the design of dynamic protein free energy landscapes for both protein engineering and drug discovery.

Molecular design, synthesis, and biological evaluation of bisamide derivatives as cyclophilin A inhibitors for HCV treatment

Hepatitis C virus (HCV) is a major cause of end-stage liver diseases. Direct-acting antivirals (DAAs), including inhibitors of nonstructural proteins (NS3/4A protease, NS5A, and NS5B polymerase), represent key components of anti-HCV treatment. However, some DAAs are associated with increased drug resistance and undesired side effects. Previous reports have shown that bisamides could be a novel class of cyclophilin A (CypA) inhibitors for treating HCV as a member of combinational therapies. To fully elucidate structure-activity relationships of bisamide derivatives and find a better hit compound with diverse binding modes, 16 biamides were designed with the help of docking program. They were then synthesized using one-pot four-component Ugi reaction. 7e with selectivity index of more than 18.9 (50% effective concentration of 5.3 μM, but no cytotoxicity at 100 μM) and unique binding mode that could be dived into gatekeeper pocket was selected as a new hit compound. Surface plasmon resonance experiments revealed that 7e is able to bind to CypA with a K_D of 3.66 μM. Taken together, these results suggest that 7e as a CypA inhibitor could be used as an alternative anti-HCV agent in combinational therapy in the future.