Coevolutionary forces shaping the fitness of SARS-CoV-2 spike glycoprotein against human receptor ACE2

The coevolutionary landscape of drug resistance in epidermal growth factor receptor: A cancer perspective

Epidermal growth factor receptor (EGFR), the first receptor tyrosine kinase, plays a critical role in neoplastic metastasis, angiogenesis, tumor invasion, and apoptosis, making it a prime target for treating non-small cell lung cancer (NSCLC). Although tyrosine kinase inhibitors (TKIs) have shown high efficacy and promise for cancer patients, resistance to these drugs often develops within a year due to alterations. The present study investigates the compensatory alterations in EGFR to understand the evolutionary process behind drug resistance. Our findings reveal that coevolutionary alterations expand the drug-binding pocket; leading to reduced drug efficacy and suggested that such changes significantly influence the structural adaptation of the EGFR against these drugs. Analysis such as root mean square deviation (RMSD), root mean square fluctuation (RMSF), solvent accessible surface area (SASA), principal component analysis (PCA), and free energy landscape (FEL) demonstrated that structures of wild EGFR docked with gefitinib are more stable which suggests its susceptibility towards drug than coevolution dependent double mutant. The findings were supported by MM-GBSA binding affinity analysis. The insights from this study highlighted the evolution-induced structural changes which contributes to drug resistance in EGFR and may certainly aid in designing more effective drugs.

Selection among site-dependent structurally constrained substitution models of protein evolution by approximate Bayesian computation

MOTIVATION

The selection among substitution models of molecular evolution is fundamental for obtaining accurate phylogenetic inferences. At the protein level, evolutionary analyses are traditionally based on empirical substitution models but these models make unrealistic assumptions and are being surpassed by structurally constrained substitution (SCS) models. The SCS models often consider site-dependent evolution, a process that provides realism but complicates their implementation into likelihood functions that are commonly used for substitution model selection.

RESULTS

We present a method to perform selection among site-dependent SCS models, also among empirical and site-dependent SCS models, based on the approximate Bayesian computation (ABC) approach and its implementation into the computational framework ProteinModelerABC. The framework implements ABC with and without regression adjustments and includes diverse empirical and site-dependent SCS models of protein evolution. Using extensive simulated data, we found that it provides selection among SCS and empirical models with acceptable accuracy. As illustrative examples, we applied the framework to analyze a variety of protein families observing that SCS models fit them better than the corresponding best-fitting empirical substitution models.

AVAILABILITY AND IMPLEMENTATION

ProteinModelerABC is freely available from https://github.com/DavidFerreiro/ProteinModelerABC, can run in parallel and includes a graphical user interface. The framework is distributed with detailed documentation and ready-to-use examples.

Antiviral Drug Target Identification and Ligand Discovery

This chapter intends to provide a general overview of web-based resources available for antiviral drug discovery studies. First, we explain how the structure for a potential viral protein target can be obtained and then highlight some of the main considerations in preparing for the application of receptor-based molecular docking techniques. Thereafter, we discuss the resources to search for potential drug candidates (ligands) against this target protein receptor, how to screen them, and preparing their analogue library. We make specific reference to free, online, open-source tools and resources which can be applied for antiviral drug discovery studies.

Coevolution-based computational approach to detect resistance mechanism of epidermal growth factor receptor

Tyrosine kinase epidermal growth factor receptor (EGFR) correlates the neoplastic cell metastasis, angiogenesis, neoplastic incursion, and apoptosis. Due to the involvement of EGFR in these biological processes, it becomes a most potent target for treating non-small cell lung cancer (NSCLC). The tyrosine kinase inhibitors (TKI) have endorsed high efficacy and anticipation to patients but unfortunately, within a year of treatment, drug targets develop resistance due to mutations. The present study detected the compensatory mutations in EGFR to know the evolutionary mechanism of drug resistance. The results of this study demonstrate that compensatory mutations enlarge the drug-binding pocket which may lead to the altered orientation of the ligand (gefitinib and erlotinib) causing drug resistance. This indicates that coevolutionary forces play a significant role in fine-tuning the structure of EGFR protein against the drugs. The analysis provides insight into the evolution-induced structural aspects of drug resistance changes in EGFR which in turn be useful in designing drugs with better efficacy.

Motifs in SARS-CoV-2 evolution

We present a novel framework enhancing the prediction of whether novel lineage poses the threat of eventually dominating the viral population. The framework is based purely on genomic sequence data, without requiring prior established biological analysis. Its building blocks are sets of coevolving sites in the alignment (motifs), identified via coevolutionary signals. The collection of such motifs forms a relational structure over the polymorphic sites. Motifs are constructed using distances quantifying the coevolutionary coupling of pairs and manifest as coevolving clusters of sites. We present an approach to genomic surveillance based on this notion of relational structure. Our system will issue an alert regarding a lineage, based on its contribution to drastic changes in the relational structure. We then conduct a comprehensive retrospective analysis of the COVID-19 pandemic based on SARS-CoV-2 genomic sequence data in GISAID from October 2020 to September 2022, across 21 lineages and 27 countries with weekly resolution. We investigate the performance of this surveillance system in terms of its accuracy, timeliness, and robustness. Lastly, we study how well each lineage is classified by such a system.

Coevolution based immunoinformatics approach considering variability of epitopes to combat different strains: A case study using spike protein of SARS-CoV-2

In the recent past several vaccines were developed to combat the COVID-19 disease. Unfortunately, the protective efficacy of the current vaccines has been reduced due to the high mutation rate in SARS-CoV-2. Here, we successfully implemented a coevolution based immunoinformatics approach to design an epitope-based peptide vaccine considering variability in spike protein of SARS-CoV-2. The spike glycoprotein was investigated for B- and T-cell epitope prediction. Identified T-cell epitopes were mapped on previously reported coevolving amino acids in the spike protein to introduce mutation. The non-mutated and mutated vaccine components were constructed by selecting epitopes showing overlapping with the predicted B-cell epitopes and highest antigenicity. Selected epitopes were linked with the help of a linker to construct a single vaccine component. Non-mutated and mutated vaccine component sequences were modelled and validated. The in-silico expression level of the vaccine constructs (non-mutated and mutated) in E. coli K12 shows promising results. The molecular docking analysis of vaccine components with toll-like receptor 5 (TLR5) demonstrated strong binding affinity. The time series calculations including root mean square deviation (RMSD), radius of gyration (RGYR), and energy of the system over 100 ns trajectory obtained from all atom molecular dynamics simulation showed stability of the system. The combined coevolutionary and immunoinformatics approach used in this study will certainly help to design an effective peptide vaccine that may work against different strains of SARS-CoV-2. Moreover, the strategy used in this study can be implemented on other pathogens.

Molecular Evolution of SARS-CoV-2 during the COVID-19 Pandemic

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) produced diverse molecular variants during its recent expansion in humans that caused different transmissibility and severity of the associated disease as well as resistance to monoclonal antibodies and polyclonal sera, among other treatments. In order to understand the causes and consequences of the observed SARS-CoV-2 molecular diversity, a variety of recent studies investigated the molecular evolution of this virus during its expansion in humans. In general, this virus evolves with a moderate rate of evolution, in the order of 10-3-10-4 substitutions per site and per year, which presents continuous fluctuations over time. Despite its origin being frequently associated with recombination events between related coronaviruses, little evidence of recombination was detected, and it was mostly located in the spike coding region. Molecular adaptation is heterogeneous among SARS-CoV-2 genes. Although most of the genes evolved under purifying selection, several genes showed genetic signatures of diversifying selection, including a number of positively selected sites that affect proteins relevant for the virus replication. Here, we review current knowledge about the molecular evolution of SARS-CoV-2 in humans, including the emergence and establishment of variants of concern. We also clarify relationships between the nomenclatures of SARS-CoV-2 lineages. We conclude that the molecular evolution of this virus should be monitored over time for predicting relevant phenotypic consequences and designing future efficient treatments.

Severe Acute Respiratory Syndrome Type 2-Causing Coronavirus: Variants and Preventive Strategies

COVID-19 vaccines have constituted a substantial scientific leap in countering severe acute respiratory syndrome type 2-causing coronavirus (SARS-CoV-2), and worldwide implementation of vaccination programs has significantly contributed to the global pandemic effort by saving many lives. However, the continuous evolution of the SARS-CoV-2 viral genome has resulted in different variants with a diverse range of mutations, some with enhanced virulence compared with previous lineages. Such variants are still a great concern as they have the potential to reduce vaccine efficacy and increase the viral transmission rate. This review summarizes the significant variants of SARS-CoV-2 encountered to date (December 2021) and discusses a spectrum of possible preventive strategies, with an emphasis on physical and materials science.

Underlying selection for the diversity of spike protein sequences of SARS-CoV-2

The global spread of SARS-CoV-2 is fast moving and has caused a worldwide public health crisis. In the present article, we analyzed spike protein sequences of SARS-CoV-2 genomes to assess the impact of mutational diversity. We observed from amino acid usage patterns that spike proteins are associated with a diversity of mutational changes and most important underlying cause of variation of amino acid usage is the changes in hydrophobicity of spike proteins. The changing patterns of hydrophobicity of spike proteins over time and its influence on the receptor binding affinity provides crucial information on the SARS-CoV-2 interaction with human receptor. Our results also show that spike proteins have evolved to prefer more hydrophobic residues over time. The present study provides a comprehensive analysis of molecular sequence data to consider that mutational variants might play a crucial role in modulating the virulence and spread of the virus and has immediate implications for therapeutic strategies.

Genomic Variation and Diversification in Begomovirus Genome in Implication to Host and Vector Adaptation

Begomoviruses (family Geminiviridae, genus Begomovirus) are DNA viruses transmitted in a circulative, persistent manner by the whitefly Bemisia tabaci (Gennadius). As revealed by their wide host range (more than 420 plant species), worldwide distribution, and effective vector transmission, begomoviruses are highly adaptive. Still, the genetic factors that facilitate their adaptation to a diverse array of hosts and vectors remain poorly understood. Mutations in the virus genome may confer a selective advantage for essential functions, such as transmission, replication, evading host responses, and movement within the host. Therefore, genetic variation is vital to virus evolution and, in response to selection pressure, is demonstrated as the emergence of new strains and species adapted to diverse hosts or with unique pathogenicity. The combination of variation and selection forms a genetic imprint on the genome. This review focuses on factors that contribute to the evolution of Begomovirus and their global spread, for which an unforeseen diversity and dispersal has been recognized and continues to expand.

COVID-19: breaking down a global health crisis

Coronavirus disease 2019 (COVID-19) is the second pandemic of the twenty-first century, with over one-hundred million infections and over two million deaths to date. It is a novel strain from the Coronaviridae family, named Severe Acute Respiratory Distress Syndrome Coronavirus-2 (SARS-CoV-2); the 7th known member of the coronavirus family to cause disease in humans, notably following the Middle East Respiratory syndrome (MERS), and Severe Acute Respiratory Distress Syndrome (SARS). The most characteristic feature of this single-stranded RNA molecule includes the spike glycoprotein on its surface. Most patients with COVID-19, of which the elderly and immunocompromised are most at risk, complain of flu-like symptoms, including dry cough and headache. The most common complications include pneumonia, acute respiratory distress syndrome, septic shock, and cardiovascular manifestations. Transmission of SARS-CoV-2 is mainly via respiratory droplets, either directly from the air when an infected patient coughs or sneezes, or in the form of fomites on surfaces. Maintaining hand-hygiene, social distancing, and personal protective equipment (i.e., masks) remain the most effective precautions. Patient management includes supportive care and anticoagulative measures, with a focus on maintaining respiratory function. Therapy with dexamethasone, remdesivir, and tocilizumab appear to be most promising to date, with hydroxychloroquine, lopinavir, ritonavir, and interferons falling out of favour. Additionally, accelerated vaccination efforts have taken place internationally, with several promising vaccinations being mass deployed. In response to the COVID-19 pandemic, countries and stakeholders have taken varying precautions to combat and contain the spread of the virus and dampen its collateral economic damage. This review paper aims to synthesize the impact of the virus on a global, micro to macro scale.

Integrated Biophysical Modeling of the SARS-CoV-2 Spike Protein Binding and Allosteric Interactions with Antibodies

Structural and biochemical studies of the severe acute respiratory syndrome (SARS)-CoV-2 spike glycoproteins and complexes with highly potent antibodies have revealed multiple conformation-dependent epitopes highlighting conformational plasticity of spike proteins and capacity for eliciting specific binding and broad neutralization responses. In this study, we used coevolutionary analysis, molecular simulations, and perturbation-based hierarchical network modeling of the SARS-CoV-2 spike protein complexes with a panel of antibodies targeting distinct epitopes to explore molecular mechanisms underlying binding-induced modulation of dynamics and allosteric signaling in the spike proteins. Through coevolutionary analysis of the SARS-CoV-2 spike proteins, we identified highly coevolving hotspots and functional clusters that enable a functional cross-talk between distant allosteric regions in the SARS-CoV-2 spike complexes with antibodies. Coarse-grained and all-atom molecular dynamics simulations combined with mutational sensitivity mapping and perturbation-based profiling of the SARS-CoV-2 receptor-binding domain (RBD) complexes with CR3022 and CB6 antibodies enabled a detailed validation of the proposed approach and an extensive quantitative comparison with the experimental structural and deep mutagenesis scanning data. By combining in silico mutational scanning, perturbation-based modeling, and network analysis of the SARS-CoV-2 spike trimer complexes with H014, S309, S2M11, and S2E12 antibodies, we demonstrated that antibodies can incur specific and functionally relevant changes by modulating allosteric propensities and collective dynamics of the SARS-CoV-2 spike proteins. The results provide a novel insight into regulatory mechanisms of SARS-CoV-2 S proteins showing that antibody-escaping mutations can preferentially target structurally adaptable energy hotspots and allosteric effector centers that control functional movements and allosteric communication in the complexes.