A proof-of-concept study of the secondary structure of influenza A, B M2 and MERS- and SARS-CoV E transmembrane peptides using folding molecular dynamics simulations in a membrane mimetic solvent.

Here, we have carried out a proof-of-concept molecular dynamics (MD) simulation with adaptive tempering in a membrane mimetic environment to study the folding of single-pass membrane peptides. We tested the influenza A M2 viroporin, influenza B M2 viroporin, and protein E from coronaviruses MERS-Cov-2 and SARS-CoV-2 peptides with known experimental secondary structures in membrane bilayers. The two influenza-derived peptides are significantly different in the peptide sequence and secondary structure and more polar than the two coronavirus-derived peptides. Through a total of more than 50 µs of simulation time that could be accomplished in trifluoroethanol (TFE), as a membrane model, we characterized comparatively the folding behavior, helical stability, and helical propensity of these transmembrane peptides that match perfectly their experimental secondary structures, and we identified common motifs that reflect their quaternary organization and known (or not) biochemical function. We showed that BM2 is organized into two structurally distinct parts: a significantly more stable N-terminal half, and a fast-converting C-terminal half that continuously folds and unfolds between α-helical structures and non-canonical structures, which are mostly turns. In AM2, both the N-terminal half and C-terminal half are very flexible. In contrast, the two coronavirus-derived transmembrane peptides are much more stable and fast helix-formers when compared with the influenza ones. In particular, the SARS-derived peptide E appears to be the fastest and most stable helix-former of all the four viral peptides studied, with a helical structure that persists almost without disruption for the whole of its 10 µs simulation. By comparing the results with experimental observations, we benchmarked TFE in studying the conformation of membrane and hydrophobic peptides. This work provided accurate results suggesting a methodology to run long MD simulations and predict structural properties of biologically important membrane peptides.

Distinct conformational states of SARS-CoV-2 spike protein.

Intervention strategies are urgently needed to control the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic. The trimeric viral spike (S) protein catalyzes fusion between viral and target cell membranes to initiate infection. Here, we report two cryo-electron microscopy structures derived from a preparation of the full-length S protein, representing its prefusion (2.9-angstrom resolution) and postfusion (3.0-angstrom resolution) conformations, respectively. The spontaneous transition to the postfusion state is independent of target cells. The prefusion trimer has three receptor-binding domains clamped down by a segment adjacent to the fusion peptide. The postfusion structure is strategically decorated by N-linked glycans, suggesting possible protective roles against host immune responses and harsh external conditions. These findings advance our understanding of SARS-CoV-2 entry and may guide the development of vaccines and therapeutics.

Plasticity of the 340-Loop in Influenza Neuraminidase Offers New Insight for Antiviral Drug Development.

The recently discovered 340-cavity in influenza neuraminidase (NA) N6 and N7 subtypes has introduced new possibilities for rational structure-based drug design. However, the plasticity of the 340-loop (residues 342-347) and the role of the 340-loop in NA activity and substrate binding have not been deeply exploited. Here, we investigate the mechanism of 340-cavity formation and demonstrate for the first time that seven of nine NA subtypes are able to adopt an open 340-cavity over 1.8 µs total molecular dynamics simulation time. The finding that the 340-loop plays a role in the sialic acid binding pathway suggests that the 340-cavity can function as a druggable pocket. Comparing the open and closed conformations of the 340-loop, the side chain orientation of residue 344 was found to govern the formation of the 340-cavity. Additionally, the conserved calcium ion was found to substantially influence the stability of the 340-loop. Our study provides dynamical evidence supporting the 340-cavity as a druggable hotspot at the atomic level and offers new structural insight in designing antiviral drugs.

A So-Far Overlooked Secondary Conformation State in the Binding Mode of SARS-CoV-2 Spike Protein to Human ACE2 and Its Conversion Rate Are Crucial for Estimating Infectivity Efficacy of the Underlying Virus Variant.

Since its outbreak in 2019, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has spread with high transmission efficiency across the world, putting health care as well as economic systems under pressure. During the course of the pandemic, the originally identified SARS-CoV-2 variant has been multiple times replaced by various mutant versions, which showed enhanced fitness due to increased infection and transmission rates. In order to find an explanation for why SARS-CoV-2 and its emerging mutated versions showed enhanced transmission efficiency compared with SARS-CoV (2002), an enhanced binding affinity of the spike protein to human angiotensin converting enzyme 2 (hACE2) has been proposed by crystal structure analysis and was identified in cell culture models. Kinetic analysis of the interaction of various spike protein constructs with hACE2 was considered to be best described by a Langmuir-based 1:1 stoichiometric interaction. However, we demonstrate in this report that the SARS-CoV-2 spike protein interaction with hACE2 is best described by a two-step interaction, which is defined by an initial binding event followed by a slower secondary rate transition that enhances the stability of the complex by a factor of ~190 (primary versus secondary state) with an overall equilibrium dissociation constant (KD) of 0.20 nM. In addition, we show that the secondary rate transition is not only present in SARS-CoV-2 wild type ("wt"; Wuhan strain) but also found in the B.1.1.7 variant, where its transition rate is 5-fold increased. IMPORTANCE The current SARS-CoV-2 pandemic is characterized by the high infectivity of SARS-CoV-2 and its derived variants of concern (VOCs). It has been widely assumed that the reason for its increased cell entry compared with SARS-CoV (2002) is due to alterations in the viral spike protein, where single amino acid residue substitutions can increase affinity for hACE2. So far, the interaction of a single unit of the CoV-2 spike protein has been described using the 1:1 Langmuir interaction kinetic. However, we demonstrate here that there is a secondary state binding step that may be essential for novel VOCs in order to further increase their infectivity. These findings are important for quantitatively understanding the infection process of SARS-CoV-2 and characterization of emerging SARS-CoV-2 variants of spike proteins. Thus, they provide a tool for predicting the potential infectivity of the respective viral variants based on secondary rate transition and secondary complex stability.

Recognition of Protein Network for Bioinformatics Knowledge Analysis Using Support Vector Machine.

Protein is the material foundation of living things, and it directly takes part in and runs the process of living things itself. Predicting protein complexes helps us understand the structure and function of complexes, and it is an important foundation for studying how cells work. Genome-wide protein interaction (PPI) data is growing as high-throughput experiments become more common. The aim of this research is that it provides a dual-tree complex wavelet transform which is used to find out about the structure of proteins. It also identifies the secondary structure of protein network. Many computer-based methods for predicting protein complexes have also been developed in the field. Identifying the secondary structure of a protein is very important when you are studying protein characteristics and properties. This is how the protein sequence is added to the distance matrix. The scope of this research is that it can confidently predict certain protein complexes rapidly, which compensates for shortcomings in biological research. The three-dimensional coordinates of C atom are used to do this. According to the texture information in the distance matrix, the matrix is broken down into four levels by the double-tree complex wavelet transform because it has four levels. The subband energy and standard deviation in different directions are taken, and then, the two-dimensional feature vector is used to show the secondary structure features of the protein in a way that is easy to understand. Then, the KNN and SVM classifiers are used to classify the features that were found. Experiments show that a new feature called a dual-tree complex wavelet can improve the texture granularity and directionality of the traditional feature extraction method, which is called secondary structure.

Peptide-Based Dual HIV and Coronavirus Entry Inhibitors.

The continued HIV/AIDS epidemic worldwide and the battle against emerging infectious diseases caused by coronaviruses underscore the need for the development of an ever-expanding repertoire of antiviral drugs. Entry inhibitors are of particular interest because of their potential to be used as therapeutic or prophylactic treatments for blocking viral invasion. HIV and coronaviruses utilize class I fusion proteins to facilitate their entry and membrane fusion. Discovery of a common hexameric coiled-coil fusion complex resulting from the packing of three C-terminal heptad repeat region from the fusion-mediating subunit of viral fusion proteins against trimeric coiled-coil made up by their N-terminal heptad repeat prompted the search for peptides mimicking the heptad repeat regions that could potentially inhibit viral entry. This has led to the development of effective peptides that are specific to the virus that is developed for. In this review, we focus on peptide-based entry dual inhibitors that block fusion process not only of HIV but also coronaviruses through interrupting their fusogenic six-helical bundle core and which hopefully will help to gain insight into the α-helical secondary structure- and coiled-coil superstructure-based strategies to design entry inhibitors with broad-spectrum antiviral activity against enveloped viruses with class I fusion proteins.

Omicron and Delta variant of SARS-CoV-2: A comparative computational study of spike protein.

Emerging severe acute respiratory syndrome coronavirus type 2 (SARS-CoV-2) variants, especially those of concern, may have an impact on the virus's transmissibility and pathogenicity, as well as diagnostic equipment performance and vaccine effectiveness. Even though the SARS-CoV-2 Delta variant (B.1.617.2) emerged during India's second wave of infections, Delta variants have grown dominant internationally and are still evolving. On November 26, 2021, World Health Organization identified the variant B.1.1.529 as a variant of concern, naming it Omicron, based on evidence that Omicron contains numerous mutations that may influence its behavior. However, the mode of transmission and severity of the Omicron variant remains unknown. We used computational studies to examine the Delta and Omicron variants in this study and found that the Omicron variant had a higher affinity for human angiotensin-converting enzyme 2 (ACE2) than the Delta variant due to a significant number of mutations in the SARS-CoV-2 receptor-binding domain (RBD), indicating a higher potential for transmission. Based on docking studies, the Q493R, N501Y, S371L, S373P, S375F, Q498R, and T478K mutations contribute significantly to high binding affinity with human ACE2. In comparison to the Delta variant, both the entire spike protein and the RBD in Omicron include a high proportion of hydrophobic amino acids such as leucine and phenylalanine. These amino acids are located within the protein's core and are required for structural stability. We observed a disorder-order transition in the Omicron variant between spike protein RBD regions 468-473, and it may be significant in the influence of disordered residues/regions on spike protein stability and binding to ACE2. A future study might investigate the epidemiological and biological consequences of the Omicron variant.

Structural polymorphism of coiled-coils from the stalk domain of SARS-CoV-2 spike protein.

Spike trimer plays a key role in SARS-CoV-2 infection and vaccine development. It consists of a globular head and a flexible stalk domain that anchors the protein into the viral membrane. While the head domain has been extensively studied, the properties of the adjoining stalk are poorly understood. Here, we characterize the coiled-coil formation and thermodynamic stability of the stalk domain and its segments. We find that the N-terminal segment of the stalk does not form coiled-coils and remains disordered in solution. The C-terminal stalk segment forms a trimeric coiled-coil in solution, which becomes significantly stabilized in the context of the full-length stalk. Its crystal structure reveals a novel antiparallel tetramer coiled-coil with an unusual combination of a-d and e-a-d hydrophobic core packing. Structural analysis shows that a subset of hydrophobic residues stabilizes different coiled-coil structures: trimer, tetramer, and heterohexamer, underscoring a highly polymorphic nature of the SARS-CoV-2 stalk sequence.

Secondary Structures of the Transmembrane Domain of SARS-CoV-2 Spike Protein in Detergent Micelles.

Spike protein of SARS-CoV-2 contains a single-span transmembrane (TM) domain and plays roles in receptor binding, viral attachment and viral entry to the host cells. The TM domain of spike protein is critical for viral infectivity. Herein, the TM domain of spike protein of SARS-CoV-2 was reconstituted in detergent micelles and subjected to structural analysis using solution NMR spectroscopy. The results demonstrate that the TM domain of the protein forms a helical structure in detergent micelles. An unstructured linker is identified between the TM helix and heptapeptide repeat 2 region. The linker is due to the proline residue at position 1213. Side chains of the three tryptophan residues preceding to and within the TM helix important for the function of S-protein might adopt multiple conformations which may be critical for their function. The side chain of W1212 was shown to be exposed to solvent and the side chains of residues W1214 and W1217 are buried in micelles. Relaxation study shows that the TM helix is rigid in solution while several residues have exchanges. The secondary structure and dynamics of the TM domain in this study provide insights into the function of the TM domain of spike protein.

Computational Design of Miniproteins as SARS-CoV-2 Therapeutic Inhibitors.

A rational therapeutic strategy is urgently needed for combating SARS-CoV-2 infection. Viral infection initiates when the SARS-CoV-2 receptor-binding domain (RBD) binds to the ACE2 receptor, and thus, inhibiting RBD is a promising therapeutic for blocking viral entry. In this study, the structure of lead antiviral candidate binder (LCB1), which has three alpha-helices (H1, H2, and H3), is used as a template to design and simulate several miniprotein RBD inhibitors. LCB1 undergoes two modifications: structural modification by truncation of the H3 to reduce its size, followed by single and double amino acid substitutions to enhance its binding with RBD. We use molecular dynamics (MD) simulations supported by ab initio density functional theory (DFT) calculations. Complete binding profiles of all miniproteins with RBD have been determined. The MD investigations reveal that the H3 truncation results in a small inhibitor with a -1.5 kcal/mol tighter binding to RBD than original LCB1, while the best miniprotein with higher binding affinity involves D17R or E11V + D17R mutation. DFT calculations provide atomic-scale details on the role of hydrogen bonding and partial charge distribution in stabilizing the minibinder:RBD complex. This study provides insights into general principles for designing potential therapeutics for SARS-CoV-2.

The atomic portrait of SARS-CoV-2 as captured by cryo-electron microscopy.

Transmission electron microscopy has historically been indispensable for virology research, as it offers unique insight into virus function. In the past decade, as cryo-electron microscopy (cryo-EM) has matured and become more accessible, we have been able to peer into the structure of viruses at the atomic level and understand how they interact with the host cell, with drugs or with antibodies. Perhaps, there was no time in recent history where cryo-EM was more needed, as SARS-CoV-2 has spread around the globe, causing millions of deaths and almost unquantifiable economic devastation. In this concise review, we aim to mark the most important contributions of cryo-EM to understanding the structure and function of SARS-CoV-2 proteins, from surface spikes to the virus core and from virus-receptor interactions to antibody binding.

1H, 13C and 15N resonance assignment of the SARS-CoV-2 full-length nsp1 protein and its mutants reveals its unique secondary structure features in solution.

Structural characterization of the SARS-CoV-2 full length nsp1 protein will be an essential tool for developing new target-directed antiviral drugs against SARS-CoV-2 and for further understanding of intra- and intermolecular interactions of this protein. As a first step in the NMR studies of the protein, we report the 1H, 13C and 15N resonance backbone assignment as well as the Cß of the apo form of the full-lengthSARS-CoV-2 nsp1 including the folded domain together with the flaking N- and C- terminal intrinsically disordered fragments. The 19.8 kD protein was characterized by high-resolution NMR. Validation of assignment have been done by using two different mutants, H81P and K129E/D48E as well as by amino acid specific experiments. According to the obtained assignment, the secondary structure of the folded domain in solution was almost identical to its previously published X-ray structure as well as another published secondary structure obtained by NMR, but some discrepancies have been detected. In the solution SARS-CoV-2 nsp1 exhibited disordered, flexible N- and C-termini with different dynamic characteristics. The short peptide in the beginning of the disordered C-terminal domain adopted two different conformations distinguishable on the NMR time scale. We propose that the disordered and folded nsp1 domains are not fully independent units but are rather involved in intramolecular interactions. Studies of the structure and dynamics of the SARS-CoV-2 mutant in solution are on-going and will provide important insights into the molecular mechanisms underlying these interactions.

Design of functionalised circular tandem repeat proteins with longer repeat topologies and enhanced subunit contact surfaces.

Circular tandem repeat proteins ('cTRPs') are de novo designed protein scaffolds (in this and prior studies, based on antiparallel two-helix bundles) that contain repeated protein sequences and structural motifs and form closed circular structures. They can display significant stability and solubility, a wide range of sizes, and are useful as protein display particles for biotechnology applications. However, cTRPs also demonstrate inefficient self-assembly from smaller subunits. In this study, we describe a new generation of cTRPs, with longer repeats and increased interaction surfaces, which enhanced the self-assembly of two significantly different sizes of homotrimeric constructs. Finally, we demonstrated functionalization of these constructs with (1) a hexameric array of peptide-binding SH2 domains, and (2) a trimeric array of anti-SARS CoV-2 VHH domains. The latter proved capable of sub-nanomolar binding affinities towards the viral receptor binding domain and potent viral neutralization function.

A generic force field for simulating native protein structures using dissipative particle dynamics.

A coarse-grained force field for molecular dynamics simulations of the native structures of proteins in a dissipative particle dynamics (DPD) framework is developed. The parameters for bonded interactions are derived by mapping the bonds and angles for 20 amino acids onto target distributions obtained from fully atomistic simulations in explicit solvent. A dual-basin potential is introduced for stabilizing backbone angles, to cover a wide spectrum of protein secondary structures. The backbone dihedral potential enables folding of the protein from an unfolded initial state to the folded native structure. The proposed force field is validated by evaluating the structural properties of several model peptides and proteins including the SARS-CoV-2 fusion peptide, consisting of α-helices, ß-sheets, loops and turns. Detailed comparisons with fully atomistic simulations are carried out to assess the ability of the proposed force field to stabilize the different secondary structures present in proteins. The compact conformations of the native states were evident from the radius of gyration and the high intensity peaks of the root mean square deviation histograms, which were found to be within 0.4 nm. The Ramachandran-like energy landscape on the phase space of backbone angles (θ) and dihedrals (ϕ) effectively captured the conformational phase space of α-helices at ∼(ϕ = 50°,θ = 90°) and ß-strands at ∼(ϕ = ±180°,θ = 90-120°). Furthermore, the residue-residue native contacts were also well reproduced by the proposed DPD model. The applicability of the model to multidomain complexes was assessed using lysozyme and a large α-helical bacterial pore-forming toxin, cytolysin A. Our study illustrates that the proposed force field is generic, and can potentially be extended for efficient in silico investigations of membrane bound polypeptides and proteins using DPD simulations.

1H, 13C and 15N Backbone chemical shift assignments of the n-terminal and central intrinsically disordered domains of SARS-CoV-2 nucleoprotein.

The nucleoprotein (N) from SARS-CoV-2 is an essential cofactor of the viral replication transcription complex and as such represents an important target for viral inhibition. It has also been shown to colocalize to the transcriptase-replicase complex, where many copies of N decorate the viral genome, thereby protecting it from the host immune system. N has also been shown to phase separate upon interaction with viral RNA. N is a 419 amino acid multidomain protein, comprising two folded, RNA-binding and dimerization domains spanning residues 45-175 and 264-365 respectively. The remaining 164 amino acids are predicted to be intrinsically disordered, but there is currently no atomic resolution information describing their behaviour. Here we assign the backbone resonances of the first two intrinsically disordered domains (N1, spanning residues 1-44 and N3, spanning residues 176-263). Our assignment provides the basis for the identification of inhibitors and functional and interaction studies of this essential protein.

Functional and structural features of proteins associated with alternative splicing.

The alternative splicing is a mechanism increasing the number of expressed proteins and a variety of these functions. We uncovered the protein domains most frequently lacked or occurred in the splice variants. Proteins presented by several isoforms participate in such processes as transcription regulation, immune response, etc. Our results displayed the association of alternative splicing with branched regulatory pathways. By considering the published data on the protein proteins encoded by the 18th human chromosome, we noted that alternative products display the differences in several functional features, such as phosphorylation, subcellular location, ligand specificity, protein-protein interactions, etc. The investigation of alternative variants referred to the protein kinase domain was performed by comparing the alternative sequences with 3D structures. It was shown that large enough insertions/deletions could be compatible with the kinase fold if they match between the conserved secondary structures. Using the 3D data on human proteins, we showed that conformational flexibility could accommodate fold alterations in splice variants. The investigations of structural and functional differences in splice isoforms are required to understand how to distinguish the isoforms expressed as functioning proteins from the non-realized transcripts. These studies allow filling the gap between genomic and proteomic data.

Intrinsic bias at non-canonical, ß-arrestin-coupled seven transmembrane receptors.

G-protein-coupled receptors (GPCRs), also known as seven transmembrane receptors (7TMRs), typically interact with two distinct signal-transducers, i.e., G proteins and ß-arrestins (ßarrs). Interestingly, there are some non-canonical 7TMRs that lack G protein coupling but interact with ßarrs, although an understanding of their transducer coupling preference, downstream signaling, and structural mechanism remains elusive. Here, we characterize two such non-canonical 7TMRs, namely, the decoy D6 receptor (D6R) and the complement C5a receptor subtype 2 (C5aR2), in parallel with their canonical GPCR counterparts. We discover that D6R and C5aR2 efficiently couple to ßarrs, exhibit distinct engagement of GPCR kinases (GRKs), and activate non-canonical downstream signaling pathways. We also observe that ßarrs adopt distinct conformations for D6R and C5aR2, compared to their canonical GPCR counterparts, in response to common natural agonists. Our study establishes D6R and C5aR2 as ßarr-coupled 7TMRs and provides key insights into their regulation and signaling with direct implication for biased agonism.

Destruction mechanisms of ozone over SARS-CoV-2.

In this pandemic SARS-CoV-2 crisis, any attempt to contain and eliminate the virus will also stop its spread and consequently decrease the risk of severe illness and death. While ozone treatment has been suggested as an effective disinfection process, no precise mechanism of action has been previously reported. This study aimed to further investigate the effect of ozone treatment on SARS-CoV-2. Therefore, virus collected from nasopharyngeal and oropharyngeal swab and sputum samples from symptomatic patients was exposed to ozone for different exposure times. The virus morphology and structure were monitored and analyzed through Atomic Force Microscopy (AFM), Transmission Electron Microscopy (TEM), Atomic Absorption Spectroscopy (AAS), and ATR-FTIR. The obtained results showed that ozone treatment not only unsettles the virus morphology but also alters the virus proteins' structure and conformation through amino acid disturbance and Zn ion release from the virus non-structural proteins. These results could provide a clearer pathway for virus elimination and therapeutics preparation.

RPocket: an intuitive database of RNA pocket topology information with RNA-ligand data resources.

BACKGROUND: RNA regulates a variety of biological functions by interacting with other molecules. The ligand often binds in the RNA pocket to trigger structural changes or functions. Thus, it is essential to explore and visualize the RNA pocket to elucidate the structural and recognition mechanism for the RNA-ligand complex formation. RESULTS: In this work, we developed one user-friendly bioinformatics tool, RPocket. This database provides geometrical size, centroid, shape, secondary structure element for RNA pocket, RNA-ligand interaction information, and functional sites. We extracted 240 RNA pockets from 94 non-redundant RNA-ligand complex structures. We developed RPDescriptor to calculate the pocket geometrical property quantitatively. The geometrical information was then subjected to RNA-ligand binding analysis by incorporating the sequence, secondary structure, and geometrical combinations. This new approach takes advantage of both the atom-level precision of the structure and the nucleotide-level tertiary interactions. The results show that the higher-level topological pattern indeed improves the tertiary structure prediction. We also proposed a potential mechanism for RNA-ligand complex formation. The electrostatic interactions are responsible for long-range recognition, while the Van der Waals and hydrophobic contacts for short-range binding and optimization. These interaction pairs can be considered as distance constraints to guide complex structural modeling and drug design. CONCLUSION: RPocket database would facilitate RNA-ligand engineering to regulate the complex formation for biological or medical applications. RPocket is available at http://zhaoserver.com.cn/RPocket/RPocket.html .

The highly flexible disordered regions of the SARS-CoV-2 nucleocapsid N protein within the 1-248 residue construct: sequence-specific resonance assignments through NMR.

The nucleocapsid protein N from SARS-CoV-2 is one of the most highly expressed proteins by the virus and plays a number of important roles in the transcription and assembly of the virion within the infected host cell. It is expected to be characterized by a highly dynamic and heterogeneous structure as can be inferred by bioinformatics analyses as well as from the data available for the homologous protein from SARS-CoV. The two globular domains of the protein (NTD and CTD) have been investigated while no high-resolution information is available yet for the flexible regions of the protein. We focus here on the 1-248 construct which comprises two disordered fragments (IDR1 and IDR2) in addition to the N-terminal globular domain (NTD) and report the sequence-specific assignment of the two disordered regions, a step forward towards the complete characterization of the whole protein.