Structure of the Human Respiratory Syncytial Virus M2-1 Protein in Complex with a Short Positive-Sense Gene-End RNA

In vitro higher-order oligomeric assembly of the respiratory syncytial virus M2-1 protein with longer RNAs

The respiratory syncytial virus (RSV) M2-1 protein is a transcriptional antitermination factor crucial for efficiently synthesizing multiple full-length viral mRNAs. During RSV infection, M2-1 exists in a complex with mRNA within cytoplasmic compartments called inclusion body-associated granules (IBAGs). Prior studies showed that M2-1 can bind along the entire length of viral mRNAs instead of just gene-end (GE) sequences, suggesting that M2-1 has more sophisticated RNA recognition and binding characteristics. Here, we analyzed the higher oligomeric complexes formed by M2-1 and RNAs in vitro using size exclusion chromatography (SEC), electrophoretic mobility shift assays (EMSA), negative stain electron microscopy (EM), and mutagenesis. We observed that the minimal RNA length for such higher oligomeric assembly is about 14 nucleotides for polyadenine sequences, and longer RNAs exhibit distinct RNA-induced binding modality to M2-1, leading to enhanced particle formation frequency and particle homogeneity as the local RNA concentration increases. We showed that particular cysteine residues of the M2-1 cysteine-cysteine-cystine-histidine (CCCH) zinc-binding motif are essential for higher oligomeric assembly. Furthermore, complexes assembled with long polyadenine sequences remain unaffected when co-incubated with ribonucleases or a zinc chelation agent. Our study provided new insights into the higher oligomeric assembly of M2-1 with longer RNA.IMPORTANCERespiratory syncytial virus (RSV) causes significant respiratory infections in infants, the elderly, and immunocompromised individuals. The virus forms specialized compartments to produce genetic material, with the M2-1 protein playing a pivotal role. M2-1 acts as an anti-terminator in viral transcription, ensuring the creation of complete viral mRNA and associating with both viral and cellular mRNA. Our research focuses on understanding M2-1's function in viral mRNA synthesis by modeling interactions in a controlled environment. This approach is crucial due to the challenges of studying these compartments in vivo. Reconstructing the system in vitro uncovers structural and biochemical aspects and reveals the potential functions of M2-1 and its homologs in related viruses. Our work may contribute to identifying targets for antiviral inhibitors and advancing RSV infection treatment.

Flexible Gaussian Accelerated Molecular Dynamics to Enhance Biological Sampling

Molecular dynamics simulations often struggle to obtain sufficient sampling to study complex molecular events due to high energy barriers separating the minima of interest. Multiple enhanced sampling techniques have been developed and improved over the years to tackle this issue. Gaussian accelerated molecular dynamics (GaMD) is a recently developed enhanced sampling technique that works by adding a biasing potential, lifting the energy landscape up, and decreasing the height of its barriers. GaMD allows one to increase the sampling of events of interest without the need of a priori knowledge of the system or the relevant coordinates. All required acceleration parameters can be obtained from a previous search run. Upon its development, several improvements for the methodology have been proposed, among them selective GaMD in which the boosting potential is selectively applied to the region of interest. There are currently four selective GaMD methods that have shown promising results. However, all of these methods are constrained on the number, location, and scenarios in which this selective boosting potential can be applied to ligands, peptides, or protein-protein interactions. In this work, we showcase a GROMOS implementation of the GaMD methodology with a fully flexible selective GaMD approach that allows the user to define, in a straightforward way, multiple boosting potentials for as many regions as desired. We show and analyze the advantages of this flexible selective approach on two previously used test systems, the alanine dipeptide and the chignolin peptide, and extend these examples to study its applicability and potential to study conformational changes of glycans and glycosylated proteins.

Hardening of Respiratory Syncytial Virus Inclusion Bodies by Cyclopamine Proceeds through Perturbation of the Interactions of the M2-1 Protein with RNA and the P Protein

Respiratory syncytial virus (RSV) RNA synthesis takes place in cytoplasmic viral factories also called inclusion bodies (IBs), which are membrane-less organelles concentrating the viral RNA polymerase complex. The assembly of IBs is driven by liquid-liquid phase separation promoted by interactions between the viral nucleoprotein N and the phosphoprotein P. We recently demonstrated that cyclopamine (CPM) inhibits RSV multiplication by disorganizing and hardening IBs. Although a single mutation in the viral transcription factor M2-1 induced resistance to CPM, the mechanism of action of CPM still remains to be characterized. Here, using FRAP experiments on reconstituted pseudo-IBs both in cellula and in vitro, we first demonstrated that CPM activity depends on the presence of M2-1 together with N and P. We showed that CPM impairs the competition between P and RNA binding to M2-1. As mutations on both P and M2-1 induced resistance against CPM activity, we suggest that CPM may affect the dynamics of the M2-1-P interaction, thereby affecting the relative mobility of the proteins contained in RSV IBs. Overall, our results reveal that stabilizing viral protein-protein interactions is an attractive new antiviral approach. They pave the way for the rational chemical optimization of new specific anti-RSV molecules.

Screening antivirals with a mCherry-expressing recombinant bovine respiratory syncytial virus: a proof of concept using cyclopamine

Bovine respiratory syncytial virus (BRSV) is a pathogenic pneumovirus and a major cause of acute respiratory infections in calves. Although different vaccines are available against BRSV, their efficiency remains limited, and no efficient and large-scale treatment exists. Here, we developed a new reverse genetics system for BRSV expressing the red fluorescent protein mCherry, based on a field strain isolated from a sick calf in Sweden. Although this recombinant fluorescent virus replicated slightly less efficiently compared to the wild type virus, both viruses were shown to be sensitive to the natural steroidal alkaloid cyclopamine, which was previously shown to inhibit human RSV replication. Our data thus point to the potential of this recombinant fluorescent BRSV as a powerful tool in preclinical drug discovery to enable high throughput compound screening.

Sigmoid Accelerated Molecular Dynamics: An Efficient Enhanced Sampling Method for Biosystems

Gaussian accelerated molecular dynamics (GaMD) is recognized as a popular enhanced sampling method for tackling long-standing challenges in biomolecular simulations. Inspired by GaMD, Sigmoid accelerated molecular dynamics (SaMD) is proposed in this work by adding a Sigmoid boost potential to improve the balance between the highest acceleration and accurate reweighting. Compared with GaMD, SaMD extends the accessible time scale and improves the computational efficiency as tested in three tasks. In the alanine dipeptide task, SaMD can produce the free energy landscape with better accuracy and efficiency. In the chignolin folding task, the estimated Gibbs free energy difference can converge to the experimental value ∼30% faster. In the protein-ligand binding task, the bound conformations are closer to the crystal structure with a minimal ligand root-mean-square deviation of 1.7 Å. The binding of the ligand XK263 to the HIV protease is reproduced by SaMD in ∼60% less simulation time.

The NEK family of serine/threonine kinases as a biomarker for cancer

Cancer is defined by unrestrained cell proliferation due to impaired protein activity. Cell cycle-related proteins are likely to play a role in human cancers, including proliferation, invasion, and therapeutic resistance. The serine/threonine NEK kinases are the part of Never In Mitosis A Kinases (NIMA) family, which are less explored kinase family involved in the cell cycle, checkpoint regulation, and cilia biology. They comprise of eleven members, namely NEK1, NEK2, NEK3, NEK4, NEK5, NEK6, NEK7, NEK8, NEK9, NEK10, and NEK11, located in different cellular regions. Recent research has shown the role of NEK family in various cancers by perversely expressing. Therefore, this review aimed to provide a systematic account of our understanding of NEK kinases; structural details; and its role in the cell cycle regulation. Furthermore, we have comprehensively reviewed the NEK kinases in terms of their expression and regulation in different cancers. Lastly, we have emphasized on some of the potential NEK inhibitors reported so far.

Microscale thermophoresis as a powerful growing analytical technique for the investigation of biomolecular interaction and the determination of binding parameters

The in vitro panel of technologies to address biomolecular interactions are in play, however microscale thermophoresis is continuously increasing in use to represent a key player in this arena. This review highlights the usefulness of microscale thermophoresis in the determination of molecular and biomolecular affinity interactions. This work reviews the literature from January 2016 to January 2022 about microscale thermophoresis. It gives a summarized overview about both the state-of the art and the development in the field of microscale thermophoresis. The principle of microscale thermophoresis is also described supported with self-created illustrations. Moreover, some recent advances are mentioned that showing application of the technique in investigating biomolecular interactions in different fields. Finally, advantages as well as drawbacks of the technique in comparison with other competing techniques are summarized.

Binding Analysis Using Accelerated Molecular Dynamics Simulations and Future Perspectives

Biomolecular recognition such as binding of small molecules, nucleic acids, peptides and proteins to their target receptors plays key roles in cellular function and has been targeted for therapeutic drug design. Molecular dynamics (MD) is a computational approach to analyze these binding processes at an atomistic level, which provides valuable understandings of the mechanisms of biomolecular recognition. However, the rather slow biomolecular binding events often present challenges for conventional MD (cMD), due to limited simulation timescales (typically over hundreds of nanoseconds to tens of microseconds). In this regard, enhanced sampling methods, particularly accelerated MD (aMD), have proven useful to bridge the gap and enable all-atom simulations of biomolecular binding events. Here, we will review the recent method developments of Gaussian aMD (GaMD), ligand GaMD (LiGaMD) and peptide GaMD (Pep-GaMD), which have greatly expanded our capabilities to simulate biomolecular binding processes. Spontaneous binding of various biomolecules to their receptors has been successfully simulated by GaMD. Microsecond LiGaMD and Pep-GaMD simulations have captured repetitive binding and dissociation of small-molecule ligands and highly flexible peptides, and thus enabled ligand/peptide binding thermodynamics and kinetics calculations. We will also present relevant application studies in simulations of important drug targets and future perspectives for rational computer-aided drug design.

Mechanism of RNA recognition by a Musashi RNA-binding protein

The Musashi RNA-binding proteins (RBPs) regulate translation of target mRNAs and maintenance of cell stemness and tumorigenesis. Musashi-1 (MSI1), long considered as an intestinal and neural stem cell marker, has been more recently found to be over expressed in many cancers. It has served as an important drug target for treating acute myeloid leukemia and solid tumors such as ovarian, colorectal and bladder cancer. One of the reported binding targets of MSI1 is Numb, a negative regulator of the Notch signaling. However, the dynamic mechanism of Numb RNA binding to MSI1 remains unknown, largely hindering effective drug design targeting this critical interaction. Here, we have performed extensive all-atom microsecond-timescale simulations using a robust Gaussian accelerated molecular dynamics (GaMD) method, which successfully captured multiple times of spontaneous and highly accurate binding of the Numb RNA from bulk solvent to the MSI1 protein target site. GaMD simulations revealed that Numb RNA binding to MSI1 involved largely induced fit in both the RNA and protein. The simulations also identified important low-energy intermediate conformational states during RNA binding, in which Numb interacted mainly with the β2-β3 loop and C terminus of MSI1. The mechanistic understanding of RNA binding obtained from our GaMD simulations is expected to facilitate rational structure-based drug design targeting MSI1 and other RBPs.

Gaussian accelerated molecular dynamics (GaMD): principles and applications

Gaussian accelerated molecular dynamics (GaMD) is a robust computational method for simultaneous unconstrained enhanced sampling and free energy calculations of biomolecules. It works by adding a harmonic boost potential to smooth biomolecular potential energy surface and reduce energy barriers. GaMD greatly accelerates biomolecular simulations by orders of magnitude. Without the need to set predefined reaction coordinates or collective variables, GaMD provides unconstrained enhanced sampling and is advantageous for simulating complex biological processes. The GaMD boost potential exhibits a Gaussian distribution, thereby allowing for energetic reweighting via cumulant expansion to the second order (i.e., "Gaussian approximation"). This leads to accurate reconstruction of free energy landscapes of biomolecules. Hybrid schemes with other enhanced sampling methods, such as the replica exchange GaMD (rex-GaMD) and replica exchange umbrella sampling GaMD (GaREUS), have also been introduced, further improving sampling and free energy calculations. Recently, new "selective GaMD" algorithms including the ligand GaMD (LiGaMD) and peptide GaMD (Pep-GaMD) enabled microsecond simulations to capture repetitive dissociation and binding of small-molecule ligands and highly flexible peptides. The simulations then allowed highly efficient quantitative characterization of the ligand/peptide binding thermodynamics and kinetics. Taken together, GaMD and its innovative variants are applicable to simulate a wide variety of biomolecular dynamics, including protein folding, conformational changes and allostery, ligand binding, peptide binding, protein-protein/nucleic acid/carbohydrate interactions, and carbohydrate/nucleic acid interactions. In this review, we present principles of the GaMD algorithms and recent applications in biomolecular simulations and drug design.

Structural Insights into the Respiratory Syncytial Virus RNA Synthesis Complexes

RNA synthesis in respiratory syncytial virus (RSV), a negative-sense (-) nonsegmented RNA virus, consists of viral gene transcription and genome replication. Gene transcription includes the positive-sense (+) viral mRNA synthesis, 5'-RNA capping and methylation, and 3' end polyadenylation. Genome replication includes (+) RNA antigenome and (-) RNA genome synthesis. RSV executes the viral RNA synthesis using an RNA synthesis ribonucleoprotein (RNP) complex, comprising four proteins, the nucleoprotein (N), the large protein (L), the phosphoprotein (P), and the M2-1 protein. We provide an overview of the RSV RNA synthesis and the structural insights into the RSV gene transcription and genome replication process. We propose a model of how the essential four proteins coordinate their activities in different RNA synthesis processes.

Respiratory syncytial virus M2-1 protein associates non-specifically with viral messenger RNA and with specific cellular messenger RNA transcripts

Respiratory syncytial virus (RSV) is a major cause of respiratory disease in infants and the elderly. RSV is a non-segmented negative strand RNA virus. The viral M2-1 protein plays a key role in viral transcription, serving as an elongation factor to enable synthesis of full-length mRNAs. M2-1 contains an unusual CCCH zinc-finger motif that is conserved in the related human metapneumovirus M2-1 protein and filovirus VP30 proteins. Previous biochemical studies have suggested that RSV M2-1 might bind to specific virus RNA sequences, such as the transcription gene end signals or poly A tails, but there was no clear consensus on what RSV sequences it binds. To determine if M2-1 binds to specific RSV RNA sequences during infection, we mapped points of M2-1:RNA interactions in RSV-infected cells at 8 and 18 hours post infection using crosslinking immunoprecipitation with RNA sequencing (CLIP-Seq). This analysis revealed that M2-1 interacts specifically with positive sense RSV RNA, but not negative sense genome RNA. It also showed that M2-1 makes contacts along the length of each viral mRNA, indicating that M2-1 functions as a component of the transcriptase complex, transiently associating with nascent mRNA being extruded from the polymerase. In addition, we found that M2-1 binds specific cellular mRNAs. In contrast to the situation with RSV mRNA, M2-1 binds discrete sites within cellular mRNAs, with a preference for A/U rich sequences. These results suggest that in addition to its previously described role in transcription elongation, M2-1 might have an additional role involving cellular RNA interactions.

Efficient purification and assembly of ribonucleoprotein complex for interaction analysis by MST assay coupled with GaMD simulations

Here, we describe a generic protocol for monitoring protein-RNA interaction using a cleavable GFP fusion of a recombinant RNA-binding protein. We detail each expression and purification step, including high salt and heparin column for contaminant RNA removal. After the assembly of RNA into the ribonucleoprotein complex, the MicroScale Thermophoresis assay enables the binding affinity to be obtained quickly with a small amount of sample. Further Gaussian accelerated molecular dynamics simulations allow us to analyze protein:RNA interactions in detail.

For complete details on the use and execution of this protocol, please refer to Gao et al. (2020).

•

Using cleavable GFP fusion to monitor RNA-binding protein expression and purification

•

High salt and heparin column to remove contamination RNA from RNA-binding proteins

•

Use the MicroScale Thermophoresis (MST) assay to obtain the binding affinity (K_d)

•

Simulate protein:RNA interactions with Gaussian accelerated molecular dynamics (GaMD)

Biophysical and Dynamic Characterization of Fine-Tuned Binding of the Human Respiratory Syncytial Virus M2-1 Core Domain to Long RNAs

The human respiratory syncytial virus (hRSV) M2-1 protein functions as a processivity and antitermination factor of the viral polymerase complex. Here, the first evidence that the hRSV M2-1 core domain (cdM2-1) alone has an unfolding activity for long RNAs is presented and the biophysical and dynamic characterization of the cdM2-1/RNA complex is provided. The main contact region of cdM2-1 with RNA was the α1-α2-α5-α6 helix bundle, which suffered local conformational changes and promoted the RNA unfolding activity. This activity may be triggered by base-pairing recognition. RNA molecules wrap around the whole cdM2-1, protruding their termini over the domain. The α2-α3 and α3-α4 loops of cdM2-1 were marked by an increase in picosecond internal motions upon RNA binding, even though they are not directly involved in the interaction. The results revealed that the cdM2-1/RNA complex originates from a fine-tuned binding, contributing to the unraveling interaction aspects necessary for M2-1 activity.IMPORTANCE The main outcome is the molecular description of the fine-tuned binding of the cdM2-1/RNA complex and the provision of evidence that the domain alone has unfolding activity for long RNAs. This binding mode is essential in the understanding of the function in the full-length protein. Human respiratory syncytial virus (hRSV), an orthopneumovirus, stands out for the unique role of its M2-1 protein as a transcriptional antitermination factor able to increase RNA polymerase processivity.

RSV M2-1 Protein in Complex with RNA: Old Questions Are Answered and a New One Emerges

The respiratory syncytial virus (RSV) M2-1 protein is essential for virus transcription. In this issue of Structure, Gao et al. (2020) describe the crystal structure of M2-1 in complex with RNA. This structure provides new insight into the RSV transcription mechanism but also raises an intriguing question regarding M2-1's function.