Ebola virus VP35 interacts non-covalently with ubiquitin chains to promote viral replication.

Ebolavirus (EBOV) belongs to a family of highly pathogenic viruses that cause severe hemorrhagic fever in humans. EBOV replication requires the activity of the viral polymerase complex, which includes the cofactor and Interferon antagonist VP35. We previously showed that the covalent ubiquitination of VP35 promotes virus replication by regulating interactions with the polymerase complex. In addition, VP35 can also interact non-covalently with ubiquitin (Ub); however, the function of this interaction is unknown. Here, we report that VP35 interacts with free (unanchored) K63-linked polyUb chains. Ectopic expression of Isopeptidase T (USP5), which is known to degrade unanchored polyUb chains, reduced VP35 association with Ub and correlated with diminished polymerase activity in a minigenome assay. Using computational methods, we modeled the VP35-Ub non-covalent interacting complex, identified the VP35-Ub interacting surface, and tested mutations to validate the interface. Docking simulations identified chemical compounds that can block VP35-Ub interactions leading to reduced viral polymerase activity. Treatment with the compounds reduced replication of infectious EBOV in cells and in vivo in a mouse model. In conclusion, we identified a novel role of unanchored polyUb in regulating Ebola virus polymerase function and discovered compounds that have promising anti-Ebola virus activity.

DockOpt: A Tool for Automatic Optimization of Docking Models.

Molecular docking is a widely used technique for leveraging protein structure for ligand discovery, but it remains difficult to utilize due to limitations that have not been adequately addressed. Despite some progress toward automation, docking still requires expert guidance, hindering its adoption by a broader range of investigators. To make docking more accessible, we developed a new utility called DockOpt, which automates the creation, evaluation, and optimization of docking models prior to their deployment in large-scale prospective screens. DockOpt outperforms our previous automated pipeline across all 43 targets in the DUDE-Z benchmark data set, and the generated models for 84% of targets demonstrate sufficient enrichment to warrant their use in prospective screens, with normalized LogAUC values of at least 15%. DockOpt is available as part of the Python package Pydock3 included in the UCSF DOCK 3.8 distribution, which is available for free to academic researchers at https://dock.compbio.ucsf.edu and free for everyone upon registration at https://tldr.docking.org.

Ebola Virus VP35 Interacts Non-Covalently with Ubiquitin Chains to Promote Viral Replication Creating New Therapeutic Opportunities.

Ebolavirus (EBOV) belongs to a family of highly pathogenic viruses that cause severe hemorrhagic fever in humans. EBOV replication requires the activity of the viral polymerase complex, which includes the co-factor and Interferon antagonist VP35. We previously showed that the covalent ubiquitination of VP35 promotes virus replication by regulating interactions with the polymerase complex. In addition, VP35 can also interact non-covalently with ubiquitin (Ub); however, the function of this interaction is unknown. Here, we report that VP35 interacts with free (unanchored) K63-linked polyUb chains. Ectopic expression of Isopeptidase T (USP5), which is known to degrade unanchored polyUb chains, reduced VP35 association with Ub and correlated with diminished polymerase activity in a minigenome assay. Using computational methods, we modeled the VP35-Ub non-covalent interacting complex, identified the VP35-Ub interacting surface and tested mutations to validate the interface. Docking simulations identified chemical compounds that can block VP35-Ub interactions leading to reduced viral polymerase activity that correlated with reduced replication of infectious EBOV. In conclusion, we identified a novel role of unanchored polyUb in regulating Ebola virus polymerase function and discovered compounds that have promising anti-Ebola virus activity.

The DynaSig-ML Python package: automated learning of biomolecular dynamics-function relationships.

The DynaSig-ML ('Dynamical Signatures-Machine Learning') Python package allows the efficient, user-friendly exploration of 3D dynamics-function relationships in biomolecules, using datasets of experimental measures from large numbers of sequence variants. It does so by predicting 3D structural dynamics for every variant using the Elastic Network Contact Model (ENCoM), a sequence-sensitive coarse-grained normal mode analysis model. Dynamical Signatures represent the fluctuation at every position in the biomolecule and are used as features fed into machine learning models of the user's choice. Once trained, these models can be used to predict experimental outcomes for theoretical variants. The whole pipeline can be run with just a few lines of Python and modest computational resources. The compute-intensive steps are easily parallelized in the case of either large biomolecules or vast amounts of sequence variants. As an example application, we use the DynaSig-ML package to predict the maturation efficiency of human microRNA miR-125a variants from high-throughput enzymatic assays. AVAILABILITY AND IMPLEMENTATION: DynaSig-ML is open-source software available at https://github.com/gregorpatof/dynasigml_package.

Sequence-sensitive elastic network captures dynamical features necessary for miR-125a maturation.

The Elastic Network Contact Model (ENCoM) is a coarse-grained normal mode analysis (NMA) model unique in its all-atom sensitivity to the sequence of the studied macromolecule and thus to the effect of mutations. We adapted ENCoM to simulate the dynamics of ribonucleic acid (RNA) molecules, benchmarked its performance against other popular NMA models and used it to study the 3D structural dynamics of human microRNA miR-125a, leveraging high-throughput experimental maturation efficiency data of over 26 000 sequence variants. We also introduce a novel way of using dynamical information from NMA to train multivariate linear regression models, with the purpose of highlighting the most salient contributions of dynamics to function. ENCoM has a similar performance profile on RNA than on proteins when compared to the Anisotropic Network Model (ANM), the most widely used coarse-grained NMA model; it has the advantage on predicting large-scale motions while ANM performs better on B-factors prediction. A stringent benchmark from the miR-125a maturation dataset, in which the training set contains no sequence information in common with the testing set, reveals that ENCoM is the only tested model able to capture signal beyond the sequence. This ability translates to better predictive power on a second benchmark in which sequence features are shared between the train and test sets. When training the linear regression model using all available data, the dynamical features identified as necessary for miR-125a maturation point to known patterns but also offer new insights into the biogenesis of microRNAs. Our novel approach combining NMA with multivariate linear regression is generalizable to any macromolecule for which relatively high-throughput mutational data is available.

Modelling conformational state dynamics and its role on infection for SARS-CoV-2 Spike protein variants.

The SARS-CoV-2 Spike protein needs to be in an open-state conformation to interact with ACE2 to initiate viral entry. We utilise coarse-grained normal mode analysis to model the dynamics of Spike and calculate transition probabilities between states for 17081 variants including experimentally observed variants. Our results correctly model an increase in open-state occupancy for the more infectious D614G via an increase in flexibility of the closed-state and decrease of flexibility of the open-state. We predict the same effect for several mutations on glycine residues (404, 416, 504, 252) as well as residues K417, D467 and N501, including the N501Y mutation recently observed within the B.1.1.7, 501.V2 and P1 strains. This is, to our knowledge, the first use of normal mode analysis to model conformational state transitions and the effect of mutations on such transitions. The specific mutations of Spike identified here may guide future studies to increase our understanding of SARS-CoV-2 infection mechanisms and guide public health in their surveillance efforts.

The NRGTEN Python package: an extensible toolkit for coarse-grained normal mode analysis of proteins, nucleic acids, small molecules and their complexes.

SUMMARY: The Najmanovich Research Group Toolkit for Elastic Networks (NRGTEN) is a Python toolkit that implements four different NMA models in addition to popular and novel metrics to benchmark and measure properties from these models. Furthermore, the toolkit is available as a public Python package and is easily extensible for the development or implementation of additional normal mode analysis models. The inclusion of the Elastic Network Contact Model developed in our group within NRGTEN is noteworthy, owing to its account for the specific chemical nature of atomic interactions. AVAILABILITY AND IMPLEMENTATION: https://github.com/gregorpatof/nrgten_package/.

RNA-Puzzles Round IV: 3D structure predictions of four ribozymes and two aptamers.

RNA-Puzzles is a collective endeavor dedicated to the advancement and improvement of RNA 3D structure prediction. With agreement from crystallographers, the RNA structures are predicted by various groups before the publication of the crystal structures. We now report the prediction of 3D structures for six RNA sequences: four nucleolytic ribozymes and two riboswitches. Systematic protocols for comparing models and crystal structures are described and analyzed. In these six puzzles, we discuss (i) the comparison between the automated web servers and human experts; (ii) the prediction of coaxial stacking; (iii) the prediction of structural details and ligand binding; (iv) the development of novel prediction methods; and (v) the potential improvements to be made. We show that correct prediction of coaxial stacking and tertiary contacts is essential for the prediction of RNA architecture, while ligand binding modes can only be predicted with low resolution and simultaneous prediction of RNA structure with accurate ligand binding still remains out of reach. All the predicted models are available for the future development of force field parameters and the improvement of comparison and assessment tools.

A three-way inter-molecular network accounts for the CaVα2δ1-induced functional modulation of the pore-forming CaV1.2 subunit.

L-type CaV1.2 channels are essential for the excitation-contraction coupling in cardiomyocytes and are hetero-oligomers of a pore-forming CaVα1C assembled with CaVß and CaVα2δ1 subunits. A direct interaction between CaVα2δ1 and Asp-181 in the first extracellular loop of CaVα1 reproduces the native properties of the channel. A 3D model of the von Willebrand factor type A (VWA) domain of CaVα2δ1 complexed with the voltage sensor domain of CaVα1C suggests that Ser-261 and Ser-263 residues in the metal ion-dependent adhesion site (MIDAS) motif are determinant in this interaction, but this hypothesis is untested. Here, coimmunoprecipitation assays and patch-clamp experiments of single-substitution variants revealed that CaVα2δ1 Asp-259 and Ser-261 are the two most important residues in regard to protein interactions and modulation of CaV1.2 currents. In contrast, mutating the side chains of CaVα2δ1 Ser-263, Thr-331, and Asp-363 with alanine did not completely prevent channel function. Molecular dynamics simulations indicated that the carboxylate side chain of CaVα2δ1 Asp-259 coordinates the divalent cation that is further stabilized by the oxygen atoms from the hydroxyl side chain of CaVα2δ1 Ser-261 and the carboxylate group of CaVα1C Asp-181. In return, the hydrogen atoms contributed by the side chain of Ser-261 and the main chain of Ser-263 bonded the oxygen atoms of CaV1.2 Asp-181. We propose that CaVα2δ1 Asp-259 promotes Ca2+ binding necessary to produce the conformation of the VWA domain that locks CaVα2δ1 Ser-261 and Ser-263 within atomic distance of CaVα1C Asp-181. This three-way network appears to account for the CaVα2δ1-induced modulation of CaV1.2 currents.