There is an urgent need to repurpose drugs against severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Recent computational-experimental screenings have identified several existing drugs that could serve as effective inhibitors of the virus’ main protease, Mpro, which is involved in gene expression and replication. Among these, ebselen (2-phenyl-1,2-benzoselenazol-3-one) appears to be particularly promising. Here, we examine, at a molecular level, the potential of ebselen to decrease Mpro activity. We find that it exhibits a distinct affinity for the catalytic region. Our results reveal a higher-affinity, previously unknown binding site localized between the II and III domains of the protein. A detailed strain analysis indicates that, on such a site, ebselen exerts a pronounced allosteric effect that regulates catalytic site access through surface-loop interactions, thereby inducing a reconfiguration of water hotspots. Together, these findings highlight the promise of ebselen as a repurposed drug against SARS-CoV-2.

Recent efforts to repurpose drugs to combat COVID-19 have identified Remdesivir as a candidate. It acts on the RNA-dependent, RNA polymerase (RdRp) of the SARS-CoV-2 virus, a protein complex responsible for mediating replication of the virus’s genome. However, its exact action mechanism, and that of other nucleotide analogue inhibitors, is not known. In this study, we examine at the molecular level the interaction of this drug and that of similar nucleotide analogue inhibitors, ribavirin and favilavir, by relying on atomistic molecular simulations and advanced sampling. By analyzing the binding free energies of these different drugs, it is found that all of them bind strongly at the active site. Surprisingly, however, ribavirin and favilavir do not bind the nucleotide on the complementary strand as effectively and seem to act by a different mechanism than remdesivir. Remdesivir exhibits similar binding interactions to the natural base adenine. Moreover, by analyzing remdesivir at downstream positions of the RNA, we also find that, consistent with a “delayed” termination mechanism, additional nucleotides can be incorporated after remdesivir is added, and its highly polar 1′-cyano group induces a set of conformational changes that can affect the normal RdRp complex function. By analyzing the fluctuations of residues that are altered by remdesivir binding, and comparing them to those induced by lethal point mutations, we find a possible secondary mechanism in which remdesivir destabilizes the protein complex and its interactions with the RNA strands.

Molecular dynamics simulations are used to unravel the SARS-CoV-2 NSP13 behavior in interaction with natural ligands and inhibitors.

The multiscale organizational structure of chromatin in eukaryotic cells is instrumental to DNA transcription, replication, and repair. At mesoscopic length scales, nucleosomes pack in a manner that serves to regulate gene expression through condensation and expansion of the genome. The particular structures that arise and their respective thermodynamic stabilities, however, have yet to be fully resolved. In this study, we combine molecular modeling using the 1CPN mesoscale model of chromatin with nonlinear manifold learning to identify and characterize the structure and free energy of metastable states of short chromatin segments comprising between 4- and 16-nucleosomes. Our results reveal the formation of two previously characterized tetranucleosomal conformations, the “α-tetrahedron” and the “β-rhombus”, which have been suggested to play an important role in the accessibility of DNA and, respectively, induce local chromatin compaction or elongation. The spontaneous formation of these motifs is potentially responsible for the slow nucleosome dynamics observed in experimental studies. Increases of the nucleosome repeat length are accompanied by more pronounced structural irregularity and flexibility and, ultimately, a dynamic liquid-like behavior that allows for frequent structural reorganization. Our findings indicate that tetranucleosome motifs are intrinsically stable structural states, driven by local internucleosomal interactions, and support a mechanistic picture of chromatin packing, dynamics, and accessibility that is strongly influenced by emergent local mesoscale structure.

All-atom molecular dynamics simulations of butyrylcholinesterase (BChE) sans inhibitor and in complex with each of fifteen dialkyl phenyl phosphate derivatives were conducted to characterize inhibitor binding modes and strengths. Each system was sampled on the 250 ns timescale in explicit ionic solvent, for a total of over 4 μs of simulation time. A K-means algorithm was used to cluster the resulting structures into distinct binding modes, which were further characterized based on atomic-level contacts between inhibitor chemical groups and active site residues. Comparison of experimentally observed inhibition constants (K I) with the resulting contact tables provides structural explanations for relative binding coefficients and highlights several notable interaction motifs. These include ubiquitous contact between glycines in the oxyanion hole and the inhibitor phosphate group; a sterically-driven binding preference for positional isomers that extend aromaticity; a stereochemical binding preference for choline-containing inhibitors, which mimic natural BChE substrates; and the mechanically-induced opening of the omega loop region to fully expose the active site gorge in the presence of choline-containing inhibitors. Taken together, these observations can greatly inform future design of BChE inhibitors, and the approach reported herein is generalizable to other enzyme-inhibitor systems and similar complexes that depend on non-covalent molecular recognition.