ATOMDANCE: Kernel-based denoising and choreographic analysis for protein dynamic comparison.

Comparative methods in molecular evolution and structural biology rely heavily upon the site-wise analysis of DNA sequence and protein structure, both static forms of information. However, it is widely accepted that protein function results from nanoscale nonrandom machine-like motions induced by evolutionarily conserved molecular interactions. Comparisons of molecular dynamics (MD) simulations conducted between homologous sites representative of different functional or mutational states can potentially identify local effects on binding interaction and protein evolution. In addition, comparisons of different (i.e., nonhomologous) sites within MD simulations could be employed to identify functional shifts in local time-coordinated dynamics indicative of logic gating within proteins. However, comparative MD analysis is challenged by the large fraction of protein motion caused by random thermal noise in the surrounding solvent. Therefore, properly denoised MD comparisons could reveal functional sites involving these machine-like dynamics with good accuracy. Here, we introduce ATOMDANCE, a user-interfaced suite of comparative machine learning-based denoising tools designed for identifying functional sites and the patterns of coordinated motion they can create within MD simulations. ATOMDANCE-maxDemon4.0 employs Gaussian kernel functions to compute site-wise maximum mean discrepancy between learned features of motion, thereby assessing denoised differences in the nonrandom motions between functional or evolutionary states (e.g., ligand bound versus unbound, wild-type versus mutant). ATOMDANCE-maxDemon4.0 also employs maximum mean discrepancy to analyze potential random amino acid replacements allowing for a site-wise test of neutral versus nonneutral evolution on the divergence of dynamic function in protein homologs. Finally, ATOMDANCE-Choreograph2.0 employs mixed-model analysis of variance and graph network to detect regions where time-synchronized shifts in dynamics occur. Here, we demonstrate ATOMDANCE's utility for identifying key sites involved in dynamic responses during functional binding interactions involving DNA, small-molecule drugs, and virus-host recognition, as well as understanding shifts in global and local site coordination occurring during allosteric activation of a pathogenic protease.

Curriculum Guidelines for Graduate and Undergraduate Virology Courses.

Curriculum guidelines for virology are needed to best guide student learning due to the continuous and ever-increasing volume of virology information, the need to ensure that undergraduate and graduate students have a foundational understanding of key virology concepts, and the importance in being able to communicate that understanding to both other virologists and nonvirologists. Such guidelines, developed by virology educators and the American Society for Virology Education and Career Development Committee, are described herein.

Evolution of drug resistance drives destabilization of flap region dynamics in HIV-1 protease.

The HIV-1 protease is one of several common key targets of combination drug therapies for human immunodeficiency virus infection and acquired immunodeficiency syndrome. During the progression of the disease, some individual patients acquire drug resistance due to mutational hotspots on the viral proteins targeted by combination drug therapies. It has recently been discovered that drug-resistant mutations accumulate on the "flap region" of the HIV-1 protease, which is a critical dynamic region involved in nonspecific polypeptide binding during invasion and infection of the host cell. In this study, we utilize machine learning-assisted comparative molecular dynamics, conducted at single amino acid site resolution, to investigate the dynamic changes that occur during functional dimerization and drug binding of wild-type and common drug-resistant versions of the main protease. We also use a multiagent machine learning model to identify conserved dynamics of the HIV-1 main protease that are preserved across simian and feline protease orthologs. We find that a key conserved functional site in the flap region, a solvent-exposed isoleucine (Ile50) that controls flap dynamics is functionally targeted by drug resistance mutations, leading to amplified molecular dynamics affecting the functional ability of the flap region to hold the drugs. We conclude that better long-term patient outcomes may be achieved by designing drugs that target protease regions that are less dependent upon single sites with large functional binding effects.

Identifying vaccine escape sites via statistical comparisons of short-term molecular dynamics.

The identification of viral mutations that confer escape from antibodies is crucial for understanding the interplay between immunity and viral evolution. We describe a molecular dynamics (MD)-based approach that goes beyond contact mapping, scales well to a desktop computer with a modern graphics processor, and enables the user to identify functional protein sites that are prone to vaccine escape in a viral antigen. We first implement our MD pipeline to employ site-wise calculation of Kullback-Leibler divergence in atom fluctuation over replicate sets of short-term MD production runs thus enabling a statistical comparison of the rapid motion of influenza hemagglutinin (HA) in both the presence and absence of three well-known neutralizing antibodies. Using this simple comparative method applied to motions of viral proteins, we successfully identified in silico all previously empirically confirmed sites of escape in influenza HA, predetermined via selection experiments and neutralization assays. Upon the validation of our computational approach, we then surveyed potential hotspot residues in the receptor binding domain of the SARS-CoV-2 virus in the presence of COVOX-222 and S2H97 antibodies. We identified many single sites in the antigen-antibody interface that are similarly prone to potential antibody escape and that match many of the known sites of mutations arising in the SARS-CoV-2 variants of concern. In the Omicron variant, we find only minimal adaptive evolutionary shifts in the functional binding profiles of both antibodies. In summary, we provide an inexpensive and accurate computational method to monitor hotspots of functional evolution in antibody binding footprints.

Computational prediction of intracellular targets of wild-type or mutant vesicular stomatitis matrix protein.

The matrix (M) protein of vesicular stomatitis virus (VSV) has a complex role in infection and immune evasion, particularly with respect to suppression of Type I interferon (IFN). Viral strains bearing the wild-type (wt) M protein are able to suppress Type I IFN responses. We recently reported that the 22-25 strain of VSV encodes a wt M protein, however its sister plaque isolate, strain 22-20, carries a M[MD52G] mutation that perturbs the ability of the M protein to block NFκB, but not M-mediated inhibition of host transcription. Therefore, although NFκB is activated in 22-20 infected murine L929 cells infected, no IFN mRNA or protein is produced. To investigate the impact of the M[D52G] mutation on immune evasion by VSV, we used transcriptomic data from L929 cells infected with wt, 22-25, or 22-20 to define parameters in a family of executable logical models with the aim of discovering direct targets of viruses encoding a wt or mutant M protein. After several generations of pruning or fixing hypothetical regulatory interactions, we identified specific predicted targets of each strain. We predict that wt and 22-25 VSV both have direct inhibitory actions on key elements of the NFκB signaling pathway, while 22-20 fails to inhibit this pathway.

Modular Synthesis of Peptide-Based Single- and Multimodal Targeted Molecular Imaging Agents.

A practical, modular synthesis of targeted molecular imaging agents (TMIAs) containing near-infrared dyes for optical molecular imaging (OMI) or chelated metals for magnetic resonance imaging (MRI) and single-photon emission correlation tomography (SPECT) or positron emission tomography (PET) has been developed. In the method, imaging modules are formed early in the synthesis by attaching imaging agents to the side chain of protected lysines. These modules may be assembled to provide a given set of single- or dual-modal imaging agents, which may be conjugated in the last steps of the synthesis under mild conditions to linkers and targeting groups. A key discovery was the ability of a metal such as gadolinium, useful in MRI, to serve as a protecting group for the chelator, 1,4,7,10-tetraazacyclododecane-1,4,7,10-tetraacetic acid (DOTA). It was further discovered that two lanthanide metals, La and Ce, can double as protecting groups and placeholder metals, which may be transmetalated under mild conditions by metals used for PET in the final step. The modular method enabled the synthesis of discrete targeted probes with two of the same or different dyes, two same or different metals, or mixtures of dyes and metals. The approach was exemplified by the synthesis of single- or dual-modal imaging modules for MRI-OMI, PET-OMI, and PET-MRI, followed by conjugation to the integrin-seeking peptide, c(RGDyK). For Gd modules, their efficacy for MRI was verified by measuring the NMR spin-lattice relaxivity. To validate functional imaging of TMIAs, dual-modal agents containing Cy5.5 were shown to target A549 cancer cells by confocal fluorescence microscopy.

The Propagation, Quantification, and Storage of Vesicular Stomatitis Virus.

Vesicular stomatitis virus (VSV) is the prototypical member of the Rhabdoviridae family of negative-sense single-stranded RNA viruses. This virus has been used as a powerful model system for decades and is currently being used as a vaccine platform and an oncolytic agent. Here, we present methods to propagate, quantitate, and store VSV. We also review the proper safety protocol for the handling of VSV, which is classified as a Biosafety Level 2 pathogen by the United States Centers for Disease Control and Prevention. © 2020 Wiley Periodicals LLC. Basic Protocol 1: Generation, purification, and storage of vesicular stomatitis virus stocks Basic Protocol 2: Quantification of vesicular stomatitis virus by plaque assay Support Protocol: Propagation of Vero cells.

The VSV matrix protein inhibits NF-κB and the interferon response independently in mouse L929 cells.

The matrix (M) protein of vesicular stomatitis virus (VSV) plays a key role in immune evasion. While VSV has been thought to suppress the interferon (IFN) response primarily by inhibiting host cell transcription and translation, our recent findings indicate that the M protein also targets NF-κB activation. Therefore, the M protein may utilize two distinct mechanisms to limit expression of antiviral genes, inhibiting both host gene expression and NF-κB activation. Here we characterize a recently reported mutation in the M protein [M(D52G)] of VSV isolate 22-20, which suppressed IFN mRNA and protein production despite activating NF-κB. 22-20 inhibited reporter gene expression from multiple promoters, suggesting that 22-20 suppressed the IFN response via M-mediated inhibition of host cell transcription. We propose that suppression of the IFN response and regulation of NF-κB are independent, genetically separable functions of the VSV M protein.

Near-Complete Genome Sequences of Vesicular Stomatitis Virus Indiana Laboratory Strains HR and T1026R1 and Plaque Isolates 22-20 and 22-25.

We report four near-complete genome sequences of vesicular stomatitis virus (VSV) Indiana obtained with Sanger and Illumina next-generation sequencing, namely, laboratory strains HR (heat resistant) and T1026R1 and isolates 22-20 and 22-25. Previously, only the M gene of these viruses had been sequenced, and these sequences were not deposited in GenBank.

The vesicular stomatitis virus matrix protein inhibits NF-κB activation in mouse L929 cells.

A previous study found that NF-κB activation is delayed in L929 cells infected with wild-type (wt) strains of VSV, while activation occurred earlier in cells infected with mutant strain T1026R1 (R1) that encodes a mutation in the cytotoxic matrix (M) protein. The integrity of the other R1 proteins is unknown; therefore our goal was to identify the viral component responsible for preventing NF-κB activation in L929 cells. We found that the M protein inhibits viral-mediated activation of NF-κB in the context of viral infection and when expressed alone via transfection, and that the M51R mutation in M abrogates this function. Addition of an IκB kinase (IKK) inhibitor blocked NF-κB activation and interferon-ß mRNA expression in cells infected with viruses encoding the M51R mutation in M. These results indicate that the VSV M protein inhibits activation of NF-κB by targeting an event upstream of IKK in the canonical pathway.

Repression of HPV16 early region transcription by the E2 protein.

HPV16 DNA is often integrated in cancers, disrupting the E1 or E2 genes. E2 can repress the E6/E7 promoter, but other models have been proposed to explain why integration promotes malignant progression. E1 and E2 are required for viral replication, and so genetic analysis of their role in transcriptional regulation is complex. Therefore, we developed an extrachromosomal vector containing HPV16 to undertake a genetic analysis of the E1 and E2 genes. We demonstrate that the E2 protein is primarily a transcriptional repressor when expressed from the virus. Furthermore, repression requires both the transactivation function of E2 and specific binding of E2 to the LCR. We find no evidence that the E1 protein directly modulates HPV16 gene expression. However, certain E1 mutations modulated transcription indirectly by altering splicing of E2 mRNA species. These data provide important insight into which E1 and E2 functions are optimal targets for anti-viral therapies.