In silico discovery and biological validation of ligands of FAD synthase, a promising new antimicrobial target.

New treatments for diseases caused by antimicrobial-resistant microorganisms can be developed by identifying unexplored therapeutic targets and by designing efficient drug screening protocols. In this study, we have screened a library of compounds to find ligands for the flavin-adenine dinucleotide synthase (FADS) -a potential target for drug design against tuberculosis and pneumonia- by implementing a new and efficient virtual screening protocol. The protocol has been developed for the in silico search of ligands of unexplored therapeutic targets, for which limited information about ligands or ligand-receptor structures is available. It implements an integrative funnel-like strategy with filtering layers that increase in computational accuracy. The protocol starts with a pharmacophore-based virtual screening strategy that uses ligand-free receptor conformations from molecular dynamics (MD) simulations. Then, it performs a molecular docking stage using several docking programs and an exponential consensus ranking strategy. The last filter, samples the conformations of compounds bound to the target using MD simulations. The MD conformations are scored using several traditional scoring functions in combination with a newly-proposed score that takes into account the fluctuations of the molecule with a Morse-based potential. The protocol was optimized and validated using a compound library with known ligands of the Corynebacterium ammoniagenes FADS. Then, it was used to find new FADS ligands from a compound library of 14,000 molecules. A small set of 17 in silico filtered molecules were tested experimentally. We identified five inhibitors of the activity of the flavin adenylyl transferase module of the FADS, and some of them were able to inhibit growth of three bacterial species: C. ammoniagenes, Mycobacterium tuberculosis, and Streptococcus pneumoniae, where the last two are human pathogens. Overall, the results show that the integrative VS protocol is a cost-effective solution for the discovery of ligands of unexplored therapeutic targets.

Flexi-pharma: a molecule-ranking strategy for virtual screening using pharmacophores from ligand-free conformational ensembles.

Computer-aided strategies are useful for reducing the costs and increasing the success-rate in drug discovery. Among these strategies, methods based on pharmacophores (an ensemble of electronic and steric features representing the target active site) are efficient to implement over large compound libraries. However, traditional pharmacophore-based methods require knowledge of active compounds or ligand-receptor structures, and only few ones account for target flexibility. Here, we developed a pharmacophore-based virtual screening protocol, Flexi-pharma, that overcomes these limitations. The protocol uses molecular dynamics (MD) simulations to explore receptor flexibility, and performs a pharmacophore-based virtual screening over a set of MD conformations without requiring prior knowledge about known ligands or ligand-receptor structures for building the pharmacophores. The results from the different receptor conformations are combined using a "voting" approach, where a vote is given to each molecule that matches at least one pharmacophore from each MD conformation. Contrarily to other approaches that reduce the pharmacophore ensemble to some representative models and score according to the matching models or molecule conformers, the Flexi-pharma approach takes directly into account the receptor flexibility by scoring in regards to the receptor conformations. We tested the method over twenty systems, finding an enrichment of the dataset for 19 of them. Flexi-pharma is computationally efficient allowing for the screening of thousands of compounds in minutes on a single CPU core. Moreover, the ranking of molecules by vote is a general strategy that can be applied with any pharmacophore-filtering program.

Optimal Reaction Coordinates and Kinetic Rates from the Projected Dynamics of Transition Paths.

Finding optimal reaction coordinates and predicting accurate kinetic rates for activated processes are two of the foremost challenges of molecular simulations. We introduce an algorithm that tackles the two problems at once: starting from a limited number of reactive molecular dynamics trajectories (transition paths), we automatically generate with a Monte Carlo approach a sequence of different reaction coordinates that progressively reduce the kinetic rate of their projected effective dynamics. Based on a variational principle, the minimal rate accurately approximates the exact one, and it corresponds to the optimal reaction coordinate. After benchmarking the method on an analytic double-well system, we apply it to complex atomistic systems: the interaction of carbon nanoparticles of different sizes in water.

Free Energy Landscapes, Diffusion Coefficients, and Kinetic Rates from Transition Paths.

We address the problem of constructing accurate mathematical models of the dynamics of complex systems projected on a collective variable. To this aim we introduce a conceptually simple yet effective algorithm for estimating the parameters of Langevin and Fokker-Planck equations from a set of short, possibly out-of-equilibrium molecular dynamics trajectories, obtained for instance from transition path sampling or as relaxation from high free-energy configurations. The approach maximizes the model likelihood based on any explicit expression of the short-time propagator, hence it can be applied to different evolution equations. We demonstrate the numerical efficiency and robustness of the algorithm on model systems, and we apply it to reconstruct the projected dynamics of pairs of C60 and C240 fullerene molecules in explicit water. Our methodology allows reconstructing the accurate thermodynamics and kinetics of activated processes, namely free energy landscapes, diffusion coefficients, and kinetic rates. Compared to existing enhanced sampling methods, we directly exploit short unbiased trajectories at a competitive computational cost.

Transition Rates and Efficiency of Collective Variables from Time-Dependent Biased Simulations.

Simulations with adaptive time-dependent bias enable an efficient exploration of the conformational space of a system. However, the dynamic information is altered by the bias. Infrequent metadynamics recovers the transition rate of crossing a barrier, if the collective variables are ideal and there is no bias deposition near the transition state. Unfortunately, these conditions are not always fulfilled. To overcome these limitations, and inspired by single-molecule force spectroscopy, we use Kramers' theory for calculating the barrier-crossing rate when a time-dependent bias is added to the system. We assess the efficiency of collective variables parameter by measuring how efficiently the bias accelerates the transitions. We present approximate analytical expressions of the survival probability, reproducing the barrier-crossing time statistics and enabling the extraction of the unbiased transition rate even for challenging cases. We explore the limits of our method and provide convergence criteria to assess its validity.

dockECR: Open consensus docking and ranking protocol for virtual screening of small molecules.

The development of open computational pipelines to accelerate the discovery of treatments for emerging diseases allows finding novel solutions in shorter periods of time. Consensus molecular docking is one of these approaches, and its main purpose is to increase the detection of real actives within virtual screening campaigns. Here we present dockECR, an open consensus docking and ranking protocol that implements the exponential consensus ranking method to prioritize molecular candidates. The protocol uses four open source molecular docking programs: AutoDock Vina, Smina, LeDock and rDock, to rank the molecules. In addition, we introduce a scoring strategy based on the average RMSD obtained from comparing the best poses from each single program to complement the consensus ranking with information about the predicted poses. The protocol was benchmarked using 15 relevant protein targets with known actives and decoys, and applied using the main protease of the SARS-CoV-2 virus. For the application, different crystal structures of the protease, and frames obtained from molecular dynamics simulations were used to dock a library of 79 molecules derived from previously co-crystallized fragments. The ranking obtained with dockECR was used to prioritize eight candidates, which were evaluated in terms of the interactions generated with key residues from the protease. The protocol can be implemented in any virtual screening campaign involving proteins as molecular targets. The dockECR code is publicly available at: https://github.com/rochoa85/dockECR.

A Bayesian approach to extracting free-energy profiles from cryo-electron microscopy experiments.

Cryo-electron microscopy (cryo-EM) extracts single-particle density projections of individual biomolecules. Although cryo-EM is widely used for 3D reconstruction, due to its single-particle nature it has the potential to provide information about a biomolecule's conformational variability and underlying free-energy landscape. However, treating cryo-EM as a single-molecule technique is challenging because of the low signal-to-noise ratio (SNR) in individual particles. In this work, we propose the cryo-BIFE method (cryo-EM Bayesian Inference of Free-Energy profiles), which uses a path collective variable to extract free-energy profiles and their uncertainties from cryo-EM images. We test the framework on several synthetic systems where the imaging parameters and conditions were controlled. We found that for realistic cryo-EM environments and relevant biomolecular systems, it is possible to recover the underlying free energy, with the pose accuracy and SNR as crucial determinants. We then use the method to study the conformational transitions of a calcium-activated channel with real cryo-EM particles. Interestingly, we recover not only the most probable conformation (used to generate a high-resolution reconstruction of the calcium-bound state) but also a metastable state that corresponds to the calcium-unbound conformation. As expected for turnover transitions within the same sample, the activation barriers are on the order of [Formula: see text]. We expect our tool for extracting free-energy profiles from cryo-EM images to enable more complete characterization of the thermodynamic ensemble of biomolecules.

Exponential consensus ranking improves the outcome in docking and receptor ensemble docking.

Consensus-scoring methods are commonly used with molecular docking in virtual screening campaigns to filter potential ligands for a protein target. Traditional consensus methods combine results from different docking programs by averaging the score or rank of each molecule obtained from individual programs. Unfortunately, these methods fail if one of the docking programs has poor performance, which is likely to occur due to training-set dependencies and scoring-function parameterization. In this work, we introduce a novel consensus method that overcomes these limitations. We combine the results from individual docking programs using a sum of exponential distributions as a function of the molecule rank for each program. We test the method over several benchmark systems using individual and ensembles of target structures from diverse protein families with challenging decoy/ligand datasets. The results demonstrate that the novel method outperforms the best traditional consensus strategies over a wide range of systems. Moreover, because the novel method is based on the rank rather than the score, it is independent of the score units, scales and offsets, which can hinder the combination of results from different structures or programs. Our method is simple and robust, providing a theoretical basis not only for molecular docking but also for any consensus strategy in general.

E5 M7+ (E=C-Pb, M=Li-Cs): A Source of Viable Star-Shaped Clusters.

Herein we report the systematic exploration of the potential energy surfaces of a series of clusters with formula E5 M7+ (E=C-Pb and M=Li-Cs). Fifteen of these combinations adopt a D5h three-dimensional seven-pointed star-like structure in a singlet state, where M atoms interact electrostatically with the E5 ring. The determining factors in the relative preference of having the D5h structure over the most competitive isomer or vice-versa are analyzed. These star-shaped systems satisfy the 4n+2 Hückel's rule and exhibit a strong diatropic (σ and π) response to an external magnetic field.