Shedding Light on Dark Chemical Matter: The Discovery of a SARS-CoV-2 Mpro Main Protease Inhibitor through Intensive Virtual Screening and In Vitro Evaluation.

The development of specific antiviral therapies targeting SARS-CoV-2 remains fundamental because of the continued high incidence of COVID-19 and limited accessibility to antivirals in some countries. In this context, dark chemical matter (DCM), a set of drug-like compounds with outstanding selectivity profiles that have never shown bioactivity despite being extensively assayed, appears to be an excellent starting point for drug development. Accordingly, in this study, we performed a high-throughput screening to identify inhibitors of the SARS-CoV-2 main protease (Mpro) using DCM compounds as ligands. Multiple receptors and two different docking scoring functions were employed to identify the best molecular docking poses. The selected structures were subjected to extensive conventional and Gaussian accelerated molecular dynamics. From the results, four compounds with the best molecular behavior and binding energy were selected for experimental testing, one of which presented inhibitory activity with a Ki value of 48 ± 5 µM. Through virtual screening, we identified a significant starting point for drug development, shedding new light on DCM compounds.

Fragment dissolved molecular dynamics: a systematic and efficient method to locate binding sites.

Diverse computational methods to support fragment-based drug discovery (FBDD) are available in the literature. Despite their demonstrated efficacy in supporting FBDD campaigns, they exhibit some drawbacks such as protein denaturation or ligand aggregation that have not yet been clearly overcome in the framework of biomolecular simulations. In the present work, we discuss a systematic semi-automatic novel computational procedure, designed to surpass these difficulties. The method, named fragment dissolved Molecular Dynamics (fdMD), utilizes simulation boxes of solvated small fragments, adding a repulsive Lennard-Jones potential term to avoid aggregation, which can be easily used to solvate the targets of interest. This method has the advantage of solvating the target with a low number of ligands, thus preventing the denaturation of the target, while simultaneously generating a database of ligand-solvated boxes that can be used in further studies. A number of scripts are made available to analyze the results and obtain the descriptors proposed as a means to trustfully discard spurious binding sites. To test our method, four test cases of different complexity have been solvated with ligand boxes and four molecular dynamics runs of 200 ns length have been run for each system, which have been extended up to 1 µs when needed. The reported results point out that the selected number of replicas are enough to identify the correct binding sites irrespective of the initial structure, even in the case of proteins having several close binding sites for the same ligand. We also propose a set of descriptors to analyze the results, among which the average MMGBSA and the average KDEEP energies have emerged as the most robust ones.

Molecular Determinants for the Activation/Inhibition of Bak Protein by BH3 Peptides.

Apoptosis is a key cell death pathway in mammalian cells. Understanding this process and its regulation has been a subject of study in the last three decades. Members of the Bcl-2 family of proteins are involved in the regulation of apoptosis through mitochondrial poration with the subsequent initiation of apoptosis. Deregulation of proapoptotic proteins contributes to the progression of many tumor processes. Understanding how these pore-forming Bcl-2 proteins Bak and Bax are activated is key to find new anticancer treatments. As no drug capable of activating Bak has been disclosed yet, the study of the structural features of BH3 peptides-known as Bak activators-relevant for binding along with its binding energy decomposition analysis, becomes essential for designing novel small-molecule mimics of BH3. Interestingly, a BH3 Bim analogue-inactivating Bak has recently been discovered, opening a question on the molecular features that determine the functions of BH3 peptides. Therefore, the present work is aimed at understanding the way BH3 peptides activate or inactivate Bak in order to identify differential structural features that can be used in drug design. For this purpose, complexes of Bak with an activator and an inhibitor have been subjected to a molecular dynamics study. Structural differences were assessed by means of the fluctuations of the corresponding principal components. Moreover, the MMPB/GBSA approach was used to compute the binding free energy of the diverse complexes to identify those residues of the BH3 peptide that exhibit the larger contributions to complex formation. The results obtained in this work show differences between activators and inhibitors, both in structural and energetic terms, which can be used in the design of new molecules that can activate or inactivate proapoptotic Bak.

Inspecting the Electronic Architecture and Semiconducting Properties of a Rosette-Like Supramolecular Columnar Liquid Crystal.

Computational and experimental studies unravel the structural and electronic properties of a novel supramolecular liquid crystal built through a hierarchical assembly process resulting in an H-bonded melamine rosette decorated with peripheral triphenylenes. The six-fold symmetry of the mesogen facilitates the formation of a highly organized hexagonal columnar mesophase stable at room temperature. X-ray diffraction and electron density maps confirm additional intra- and intercolumn segregation of functional subunits, and this paves the way for 1D charge transport. Indeed, hole mobility has been measured and found to be higher than for related mesogens. DFT calculations of HOMO and LUMO levels and parameters such as reorganization energy and transfer integral of the rosette structure have been achieved, and not only validate the columnar organization but also establish the way it translates into a favorable electronic architecture and molecular orbital interactions to promote charge carrier mobility.

DFT study of the effect of fluorine atoms on the crystal structure and semiconducting properties of poly(arylene-ethynylene) derivatives.

The effect of fluorine substitution on the molecular structure, crystal packing, and n-type semiconducting properties of a set of poly(arylene-ethynylene) polymers based on alternating thiadiazole and phenyl units linked through ethynylene groups has been studied by means of Density Functional Theory. As a result, an enlargement in the interplanar distance between cofacial polymer chains, as well as a decrease of the electronic coupling and electron mobility is predicted. On the other hand, fluorination could facilitate electron injection into the material. A polymer containing both alkoxy pendant chains and fluorine atoms is proposed as a compromise solution between efficiency of electron injection and charge transport within the material.

Effect of five-membered ring and heteroatom substitution on charge transport properties of perylene discotic derivatives: A theoretical approach.

Density functional theory calculations were carried out to investigate the evolvement of charge transport properties of a set of new discotic systems as a function of ring and heteroatom (B, Si, S, and Se) substitution on the basic structure of perylene. The replacement of six-membered rings by five-membered rings in the reference compound has shown a prominent effect on the electron reorganization energy that decreases ∼0.2 eV from perylene to the new carbon five-membered ring derivative. Heteroatom substitution with boron also revealed to lower the LUMO energy level and increase the electron affinity, therefore lowering the electron injection barrier compared to perylene. Since the rate of the charge transfer between two molecules in columnar discotic systems is strongly dependent on the orientation of the stacked cores, the total energy and transfer integral of a dimer as a disc is rotated with respect to the other along the stacking axis have been predicted. Aimed at obtaining a more realistic approach to the bulk structure, the molecular geometry of clusters made up of five discs was fully optimized, and charge transfer rate and mobilities were estimated for charge transport along a one dimensional pathway. Heteroatom substitution with selenium yields electron transfer integral values ∼0.3 eV with a relative disc orientation of 25°, which is the preferred angle according to the dimer energy profile. All the results indicate that the tetraselenium-substituted derivative, not synthetized so far, could be a promising candidate among those studied in this work for the fabrication of n-type semiconductors based on columnar discotic liquid crystals materials.

A DFT approach to the charge transport related properties in columnar stacked π-conjugated N-heterocycle cores including electron donor and acceptor units.

We present a density functional theory (DFT) study on charge-transport related properties in a series of discotic systems based on 1,3,5-triazine and tris[1,2,4]triazolo[1,3,5]triazine central cores as electron acceptor units, and phenyl-thiophene and N-carbazolyl-thiophene segments as electron donor units. The presence of both electron donor and acceptor moieties in the π-conjugated core could lead to new discotic liquid crystal (DLC) materials which are predicted to display ambipolar charge transport behavior in such a way that electrons could move through the central part of the next cores while holes mainly do through the peripheral groups. A significant increase in hole mobility when N-carbazolyl is present as an electron donor unit in the peripheral region is predicted. In addition, a detailed topological analysis of the electron charge density within the framework provided by Quantum Theory of Atoms in Molecules (QTAIM) has been performed in order to characterize intra- and intermolecular interactions in terms of hydrogen bonds and/or π···π stacking which contribute to the stabilization of the columnar stack and the helical self-assembly at the molecular scale.

Theoretical estimation of the optical bandgap in a series of poly(aryl-ethynylene)s: a DFT study.

Aimed to optimize the ratio accuracy/computational cost, in this work we study the performance of three different theoretical methodologies in the calculation of the optical bandgap for a test set made of a number of poly(aryl-ethynylene)s related polymers. Infinite, ideal polymer chains were first optimized by means of periodic calculations. Different length oligomers were afterward generated by direct replication of the corresponding periodic structure and their optical bandgaps were calculated by means of different time dependent-density functional theory (TD-DFT) methodologies. These results were fitted to an exponential function for each oligomer family in order to get a theoretical estimation of the optical bandgap for each polymer to be compared to the experimental reported values. The best result was obtained for TD-M06-2X yielding an average deviation of 3.4% with respect to the experimental values.

The Discovery of inhibitors of the SARS-CoV-2 S protein through computational drug repurposing.

SARS-CoV-2 must bind its principal receptor, ACE2, on the target cell to initiate infection. This interaction is largely driven by the receptor binding domain (RBD) of the viral Spike (S) protein. Accordingly, antiviral compounds that can block RBD/ACE2 interactions can constitute promising antiviral agents. To identify such molecules, we performed a virtual screening of the Selleck FDA approved drugs and the Selleck database of Natural Products using a multistep computational procedure. An initial set of candidates was identified from an ensemble docking process using representative structures determined from the analysis of four 3 µ s molecular dynamics trajectories of the RBD/ACE2 complex. Two procedures were used to construct an initial set of candidates including a standard and a pharmacophore guided docking procedure. The initial set was subsequently subjected to a multistep sieving process to reduce the number of candidates to be tested experimentally, using increasingly demanding computational procedures, including the calculation of the binding free energy computed using the MMPBSA and MMGBSA methods. After the sieving process, a final list of 10 candidates was proposed, compounds which were subsequently purchased and tested ex-vivo. The results identified estradiol cypionate and telmisartan as inhibitors of SARS-CoV-2 entry into cells. Our findings demonstrate that the methodology presented here enables the discovery of inhibitors targeting viruses for which high-resolution structures are available.

Crystal structure and charge transport properties of poly(arylene-ethynylene) derivatives: a DFT approach.

In the present study, a series of crystalline poly(arylene-ethynylene) copolymers containing phenylethynylene and 2,5-dialkoxy-phenylethynylene units together with 1,3,4-thiadiazole rings has been modeled by means of periodic calculations. Optimized three-dimensional polymeric structures show interchain distances that are consistent with the experimental values reported for a related polymer. It has also been observed that the presence of pendant alkoxy chains brings on both a further flattening and a separation of the coplanar chains. This fact is linked to a decrease of the interchain cofacial distance. The electron transport character of the polymer crystal structures was assessed through Marcus theory. Electronic coupling between neighboring polymer chains is most influenced by the presence of alkoxy chains giving rise to an expectable enhancement of the electron hopping mobility.