Next Article in Journal
CRISPR-Based Genome Editing Tools: Insights into Technological Breakthroughs and Future Challenges
Next Article in Special Issue
Mutational Asymmetries in the SARS-CoV-2 Genome May Lead to Increased Hydrophobicity of Virus Proteins
Previous Article in Journal
LTBP4 in Health and Disease
Previous Article in Special Issue
Protective Role of a TMPRSS2 Variant on Severe COVID-19 Outcome in Young Males and Elderly Women
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genome-Scale Metabolic Model of Infection with SARS-CoV-2 Mutants Confirms Guanylate Kinase as Robust Potential Antiviral Target

1
Department of Computer Science, University of Tübingen, 72076 Tübingen, Germany
2
Cluster of Excellence ‘Controlling Microbes to Fight Infections’, University of Tübingen, 72076 Tübingen, Germany
3
Computational Systems Biology of Infections and Antimicrobial-Resistant Pathogens, Institute for Bioinformatics and Medical Informatics (IBMI), University of Tübingen, 72076 Tübingen, Germany
4
Bernhard Nocht Institute for Tropical Medicine, Virus Immunology, 20359 Hamburg, Germany
5
German Center for Infection Research (DZIF), Partner Site Tübingen, 72076 Tübingen, Germany
*
Author to whom correspondence should be addressed.
Genes 2021, 12(6), 796; https://doi.org/10.3390/genes12060796
Submission received: 29 April 2021 / Revised: 19 May 2021 / Accepted: 21 May 2021 / Published: 24 May 2021
(This article belongs to the Special Issue COVID-19 and Molecular Genetics)

Abstract

:
The current SARS-CoV-2 pandemic is still threatening humankind. Despite first successes in vaccine development and approval, no antiviral treatment is available for COVID-19 patients. The success is further tarnished by the emergence and spreading of mutation variants of SARS-CoV-2, for which some vaccines have lower efficacy. This highlights the urgent need for antiviral therapies even more. This article describes how the genome-scale metabolic model (GEM) of the host-virus interaction of human alveolar macrophages and SARS-CoV-2 was refined by incorporating the latest information about the virus’s structural proteins and the mutant variants B.1.1.7, B.1.351, B.1.28, B.1.427/B.1.429, and B.1.617. We confirmed the initially identified guanylate kinase as a potential antiviral target with this refined model and identified further potential targets from the purine and pyrimidine metabolism. The model was further extended by incorporating the virus’ lipid requirements. This opened new perspectives for potential antiviral targets in the altered lipid metabolism. Especially the phosphatidylcholine biosynthesis seems to play a pivotal role in viral replication. The guanylate kinase is even a robust target in all investigated mutation variants currently spreading worldwide. These new insights can guide laboratory experiments for the validation of identified potential antiviral targets. Only the combination of vaccines and antiviral therapies will effectively defeat this ongoing pandemic.

1. Introduction

Since its emergence in December 2019 [1], individual cases of Severe Acute Respiratory Syndrome (SARS) coronavirus (CoV) type 2 (SARS-CoV-2) infections have evolved into an uncontrolled pandemic. As a result, more than 2.8 million people have lost their lives to or with Coronavirus Disease 19 (COVID-19) by March 2021. COVID-19 symptoms range from pneumonia to severe lung, heart, liver, kidney, neurological or gastrointestinal dysfunction [2]. While great efforts have been employed to provide effective SARS-CoV-2 vaccines [3,4], their success is overshadowed by the emergence of viral escape mutants and the shortcomings in developing targeted antiviral treatments. A meta-analysis by [5] demonstrates that in non-severe cases of COVID-19, there is little to no evidence for effective use of ribavirin, hydroxychloroquine, umifenovir, lopinavir/ritonavir, or interferon [5]. Even the putative effectiveness of remdesivir is questionable [6,7].
While antiviral medication development was less fruitful, as of March 2021, there are 13 vaccines for SARS-CoV-2 in use, most of them targeting the spike (S) protein [3,8]. Albeit the successes in vaccine development, reports of mutations are increasing. Some of these mutations are even bypassing the immunity provided by several vaccine candidates. Five mutation variants have prevailed, disseminate rapidly, and are classified as variants of concern or variants of interest: (i) B.1.1.7, first detected in the United Kingdom; (ii) P.1 (also called B.1.1.28), first detected in Japan and Brazil; (iii) B.1.351, first detected in South Africa; (iv) B.1.427/B.1.429, first detected in the US [9,10]; and (v) B.1.617, first detected in India [11]. While the consequences of some of these mutations for vaccine efficacy have been reported, the metabolic implications of them remain unclear.
SARS-CoV-2 is a member of the Betacoronavirus genus within the Nidovirales order [4,12]. The virus has a 27 kb to 32 kb positive sense, single-stranded RNA genome encoding 26 proteins, including the four structural proteins spike (S), envelope (E), membrane (M) and nucleocapsid (N) [4,13]. The S trimers [14,15] scan the host cells surface for the viral entry receptor angiotensin converting enzyme 2 (ACE-2) and therefore initiate the entry process [1,4]. The structural proteins E and M facilitate viral transport, assembly, budding, and release of SARS-CoV-2 virions from infected host cells [1,4]. While N is expressed within the host cytoplasm, the other structural proteins S, E, and M are translated within the endoplasmic reticulum–Golgi intermediate compartment (ERGIC) of the host cell [2,4]. SARS-CoV-2 N supports replicating the viral genome in the cytoplasm and encloses novel viral RNA to form viral ribonucleoprotein complexes (vRNPs) [2]. During the viral replication process’s final steps, these cytoplasmic vRNPs are assembled with S, E, and M proteins within the ERGIC [2,4]. The mature virions bud at the ERGIC membrane, forming vesicles which are subsequently released from the host cell via exocytosis [2,4].
Viral lipid envelopes protect the vRNPs and facilitate the particles’ entry into host cells [16]. They are usually acquired via budding from the plasma membrane or other cellular organelles [16,17]. Viruses specifically modify host membrane structures, the composition, and the whole host lipid metabolism to favor viral replication [16,18,19]. Many viruses exploit spatiotemporally enriched microdomains or rafts containing different lipid species [19]. To this end, cholesterol, for instance, increases host membrane fluidity for efficient viral entry, replication, and budding, while phosphatidylserine supports viral entry [19]. Altogether, various modifications in viral egress areas determine the differing composition of viral envelopes, thereby influencing their stability and infectivity [19]. As SARS-CoV-2 buds from the ERGIC [2,4], its envelope lipid bilayer resembles this host organelle’s composition [2,4]. The viral membrane formation mostly requires cholesterol and phospholipids, while sphingomyelin and cardiolipin are presumably less abundant [20,21].
In our previous work, we have generated an integrated human-virus metabolic model, which combines flux balance analysis (FBA) and flux variability analysis (FVA) to model the metabolic changes within SARS-CoV-2 infected human alveolar macrophages [22]. The GEM is based on the already published and well-developed human alveolar macrophage model iAB-AMØ-1410 by by Bordbar et al. [23]. Disabling viral replication in human alveolar macrophages might be an early way of intervention and prevention of the virus’s further spread. The model was employed to predict putative antiviral targets such as guanylate kinase 1 (GK1) or the availability of l-isoleucine and l-lysine [22]. Some of these potential targets may be directly targeted by small molecules or antivirals [24,25,26]. Increasing knowledge of SARS-CoV-2 facilitates the model’s improvement by incorporating recent findings of the copy number of the structural proteins [22,27]. The stoichiometric coefficients of the metabolic requirements for amino acids and nucleotides and energy requirements can be refined to predict the viral replication capacity better. Additionally, the lipid requirements were now accounted for in the viral biomass objective function (VBOF). This study presents an updated version of the integrated alveolar macrophage SARS-CoV-2 GEM and the consequences of prominent mutations for predicted metabolic targets.

2. Materials and Methods

2.1. Correcting the Copy Number of Structural Proteins

In the previous version of the VBOF, the copy number of structural proteins was not yet known. We conducted extensive literature research to identify the precise copy number of each structural protein individually. The search was mainly focused on SARS-CoV-2 directly. However, if no information was found for the novel coronavirus, we also searched for information on closely related coronaviruses.
With the identified copy numbers (see Table 2), the stoichiometric coefficients of the nucleotides, amino acids, and energy requirements were re-calculated, as Renz et al. state [22]. However, instead of using a general copy number for all structural proteins, as Aller et al. describe [27], the individual copy numbers of the respective structural proteins were used.
After the VBOF was updated with the corrected stoichiometric coefficients, the knock-out and host-derived enforcement analyses were repeated, as Renz et al. describe [22]. The knock-out experiments were performed by subsequently knocking out each reaction and evaluating its effect on the host’s maintenance and viral replication capacity (VBOF). For the host-derived enforcement analyses, the FVA was used to determine flux ranges that allow for 100% maintenance of the host, while decreasing the viral growth by at least 20%. The adapted host-derived enforcement algorithm was used, as Renz et al. describe [22].

2.2. Testing the Targets’ Robustness against for Several Mutations

The Global Initiative on Sharing All Influenza Data (GISAID) database has a collection of more than 1.5 million viral sequences of SARS-CoV-2 (May 2021). We set the following filters for the sequences: (i) variant (VUI202012/01 GRY (B.1.1.7) for variant B.1.1.7; GH/501Y.v2 (B.1.351) for variant B.1.351; GR/501Y.V3 (P.1) for variant B.1.1.28; GH/452R.V1 (B.1.429+B.1.427) for variants B.1.429 and B.1.427; and G/452R.V3 (B.1.617+) for variant B.1.617) and (ii) location (Europe/United Kingdom for variant B.1.1.7; Africa for variant B.1.351; South America for variant B.1.1.28; North America/USA for variants B.1.429 and B.1.427; and Asia/India for variant B.1.617). We randomly downloaded ten sequences from each mutation variant with the filters set as described. In addition to the sequences, we downloaded the mutation information given in the metadata. All tested mutations are listed in the Supplementary Table S1. With this information, the stoichiometric coefficients for the VBOF were calculated for every downloaded mutation. As the calculation of the nulceotides’ stoichiometric coefficients requires the nulceotide sequence, the downloaded sequences were used directly for this step. For the calculation of the amino acids’ coefficients, we used the annotated protein sequence of the SARS-CoV-2 reference sequence (NCBI accession: NC_045512.2) and the mutation information extracted from the metadata files. An algorithm adapted the amino acids from the protein sequence in accordance with the defined mutations, including substitutions, deletions, and introductions of stop codons. With the calculation of the energy requirements and pyrophosphate liberation, all stoichiometric coefficients for the VBOF were available and could be compared. For the first comparison, the mean and standard deviation of all mutations was calculated for each coefficient. These mean values were compared to the wildtype stoichiometric coefficients by calculation the difference. In subsequent analysis, the mean was calculated for the five mutation variants and was then compared to the wildtype. Again, the difference between the coefficients was calculated and visualized. With all generated VBOFs, the reaction knock-out experiments were repeated, as described in the previous section.

2.3. Lipids as Part of the Viral Biomass Objective Function

Literature research was conducted to identify potential fatty acids that occur in the capsid of SARS-CoV-2. As no lipidomics data of SARS-CoV-2 existed at the time of writing, we focused on the five identified lipids phosphatidylcholine, phosphatidylethanolamine, phosphatidylinositol, phosphatidylserine, and cholesterol. The influence of the individual lipids’ inclusion into the VBOF on the objective value when optimizing for the VBOF was evaluated. An overview of the overall procedure for testing the lipids’ influence is given in Figure 1.
As no data were available for the amount of the respective lipids in one virion, we varied the stoichiometric coefficients between 0 and 0.5. The stoichiometric coefficients of the lipids within the macrophage’s biomass maintenance function varied from 0.00102 for phosphatidylserine to 0.0315 for phosphatidylcholine (see also Table 1).
With the variation of the stoichiometric coefficients between 0 and 0.5, we covered the 14 to 490-fold increase of the stoichiometric coefficients, depending on their initial value. In the next step, all lipids were added simultaneously to the VBOF. We evaluated the VBOF’s objective value using both the lipids’ stoichiometric coefficients from the macrophage’s maintenance function and their ten-fold value.
To evaluate the effect of the lipids’ inclusion on the potential antiviral targets, we again used the stoichiometric coefficients of the macrophage’s maintenance function and a multiplication coefficient, ranging from 0 to 10 as the actual coefficient of the lipids is unknown. We conducted the knock-out experiments as Renz et al. describe [22] for each tested coefficient by knocking out each reaction individually and analyzing its effect on both the viral growth and the host’s maintenance function. While varying the multiplication coefficient, two additional reactions occurred, whose knock-out decreased the viral growth rate.
To investigate, which lipid influences the knock-out experiments most, we again analyzed the lipids individually. As done for the effect on the VBOF’s objective value, we first varied the stoichiometric coefficients between 0 and 0.5. Subsequently, we used a multiplication coefficient ranging from 0 to 10, which was multiplied with the coefficient of the macrophage’s maintenance function (see Table 1).

3. Results

3.1. Correcting the Copy Number of Structural Proteins

The single-stranded RNA genome of SARS-CoV-2 has 26 proteins [13], including four structural proteins. These four structural proteins need to be produced by the host in higher amounts than the non-structural proteins. However, the actual number of copies of each structural protein was unknown when the novel coronavirus arose, and the first studies were conducted at the beginning of the year 2020.
After extensive literature research, we collected the latest information about the copy number of the structural proteins of SARS-CoV-2. [15] identify on average 40 copies of the trimeric spike (S) protein on the surface of SARS-CoV-2, resulting in 120 copies of the S protein. [14] estimate the number of S trimers per virion to be 48, resulting in a similar copy number range as [15]. Since [15] use in situ structural analysis and [14] use mathematical estimations, we chose to use a copy number of 120 S proteins for further analysis (see Table 2). The number of the envelope (E) proteins is approximated to 20 copies [29] based on analyses of the OC43 human coronavirus (hCOV) [30] and the transmissible gastroenteritis virus (TGEV) [31]. Exactly like SARS-CoV-2, both viruses belong to the family of Coronaviridae, and hCOV also belongs to the same genus Betacoronavirus as SARS-CoV-2. Currently, no numbers for the E protein are available for SARS-CoV-2. For that reason, the number is approximated from related coronaviruses. The nucleocapsid (N) packs the viral RNA in so-called vRNPs. [14] observe 38 vRNPs per SARS-CoV-2 virion [14]. Approximately 12 copies of the N protein are located in one vRNP in SARS-CoV-2 [32,33]. Multiplying those two numbers results in 456 copies of the N protein. The amount of membrane proteins is not yet determined for SARS-CoV-2. [13] provide key numbers about SARS-CoV-2, including the copy numbers of the S, M, N, and E protein. However, all copy numbers are derived from SARS-CoV-1 or TGEV. We found precise numbers for the copy number of N proteins in SARS-CoV-2, and [34] determine the estimated ratios of M to N proteins ranging from 3M:1N to 1M:1N with 730 to 2200 N proteins per virion [34]. With this information at hand, we estimated the copy number of M proteins to 1000 by doubling the number of N proteins and rounding them up. The ratio of 2M:1N was chosen based on the article of [13], where the number of N proteins is stated as 1000 copies for SARS-CoV-1 and the number of M proteins as 2000. All used copy numbers are listed in Table 2.
With the updated copy numbers, the stoichiometric coefficients of the nucleotides, amino acids, and energy requirements were re-calculated for the viral biomass objective function (VBOF) of SARS-CoV-2. The subsequent analyses for identifying potential antiviral targets consisted of knock-out and host-derived enforcement experiments, as Renz et al. describe [22]. The guanylate kinase 1 (GK1) remains a promising antiviral target after the adaptions of the copy number of structural proteins based on the knock-out experiments.
The results of the host-derived enforcement analyses were dependent on the Copy number of structural proteins [22]. As we now identified more precise copy numbers, we can also determine the host-derived enforcement analysis results more precisely. In total, 21 reactions were identified, whose inhibition decreases the viral replication capacity by at least 20% without harming the host’s maintenance (100%). These reactions, their inhibition range, and the reduction of the VBOF are visualized in Figure 2. Reactions could be inhibited between 72% and 89%. As seen in the knock-out experiments, the guanylate kinase 1 (GK1) is the only reaction where a complete inhibition (100%) is possible. The ribose-5-phosphate isomerase (RPI) and phosphoribosylpyrophosphate synthetase (PRPPS) are part of the pentose phosphate pathway. glutamine phosphoribosyldiphosphate amidotransferase (GLUPRT), phosphoribosylglycinamide synthase (PRAGSr), phosphoribosylglycinamide formyltransferase (GARFT), phosphoribosylformylglycinamidine synthase (PRFGS), phosphoribosylaminoimidazole synthase (PRAIS), Phosphoribosylaminoimidazole carboxylase (AIRCr), phosphoribosylaminoimidazolesuccinocarboxamide synthase (PRASCS), phosphoribosylaminoimidazolecarboxamide formyltransferase (AICART), and inosine monophosphate (IMP) cyclohydrolase (IMPC) are involved in the purines’ biosynthetic pathway, more precisely in the biosynthesis of IMP [35]. Reactions associated with the purine adenosine monophospate (AMP) biosynthesis were also identified as potential targets, namely adenylosuccinate synthase (ADSS), and adenylosuccinate lyase 1 and 2 (ADSL1, ADSL2) [35].
Besides the reactions associated with the purine metabolism, the host-derived enforcement analysis also reported reactions from the pyrimidine biosynthesis, such as the carbamoyl-phosphate synthase (CBPS), aspartate carbamoyltransferase (ASPCTr), dihydroorotase (DHORTS), dihydoorotic acid dehydrogenase (DHORD9), orotate phosphoribosyltransferase (ORPT), and orotidine-5’-phosphate decarboxylase (OMPDC) [36].

3.2. Testing the Targets’ Robustness for Several Mutations

3.2.1. Analysis of Mutant-Specific Variations in the Viral Biomass

Novel mutations of SARS-CoV-2 emerge on a daily basis. Five mutation variants have prevailed, disseminate rapidly, and are classified as variants of concern or variants of interest: (i) B.1.1.7, (ii) P.1 (also called B.1.1.28), (iii) B.1.351, (iv) B.1.427/B.1.429 [9,10], and (v) B.1.617 [11]. The GISAID was launched in 2008 to promote the international sharing of virus data [37,38]. When the novel coronavirus emerged, GISAID was expanded by a database for sharing sequenced viral genomes of SARS-CoV-2 globally. At the time of writing, more than 1.5 million viral sequences of SARS-CoV-2 are collected in the database. To investigate the mutations’ effect on the previously identified potential antiviral targets, sequences of each mutation variant were downloaded from GISAID and analyzed. The stoichiometric coefficients of each variant were calculated as Renz et al. describe [22]: For the calculation of the nucleotides’ stoichiometric coefficients, the downloaded RNA sequence was used. The amino acids’ stoichiometric coefficients were calculated using the provided information about the identified mutations and the reference (wildtype) protein sequence of the first sequenced SARS-CoV-2. With this information, the abundance of the different amino acids in the different proteins was adapted for each mutation variant. The nucleotide and amino acid counts were subsequently used to calculate the pyrophosphate liberation and the adenosine triphosphates (ATPs) requirements. For each downloaded mutation variant, an individualized VBOF was created with the calculated stoichiometric coefficients.
To assess the mutations’ effect on the VBOF’s stoichiometric coefficients, we first calculated the mean and standard deviation from all stoichiometric coefficients for all mutations and compared them to the wildtype (WT) coefficients. The mean stoichiometric coefficients of the mutations are very similar to the wildtpye’s stoichiometric coefficents. The largest difference is observed for the amino acid l-aspartate: The stoichiometric coefficient for l-aspartate is decreased by on average 0.005 in the mutations compared to the wildtype. Figure 3 visualizes the comparison of the mutations’ mean stoichiometric coefficients with the wildtype coefficients.
Since we analyzed five distinct mutation variants, the differences in the stoichiometric coefficients were examined further based on these variants. The mean for each stoichiometric coefficient was calculated variant-wise. With this mean, the deviation from the wildtype coefficient was calculated and visualized as a heat-map in Figure 4. This analysis gives further insight into the properties of the individual mutations.
One can observe a pattern for the stoichiometric coefficients of adenosine diphosphate (ADP) and ATP: While the mutation variants B.1.1.7 and B.1.1.28 have decreased stoichiometric coefficients (−0.01) compared to the wildtype, the variants B.1.351 and B.1.427/429 have increased stoichiometric coefficients (0.019 to 0.021). This pattern is most apparent for ADP and ATP, but can also be observed for other stoichiometric coefficients, such as for diphosphate (PPi), l-lysine, l-threonine, or l-valine. To further investigate this pattern, we examined the calculation for the stoichiometric coefficients. Each coefficient is set in relation to the total viral molar mass (Mv), which is the sum of the total molar mass of all nucleotides (Gi) and amino acids (Gj). The mutation variants B.1.1.7 and B.1.1.28 have a higher total viral molar mass compared to the mutation variants B.1.351 and B.1.427/429. This increased total viral molar mass is based on an increased molar mass of both nucleotides (Gi) and amino acids (Gj). As the stoichiometric coefficients for ADP and ATP larger than the other coefficients, this pattern is more apparent.
However, this pattern does not emerge in all stoichiometric coefficients. There are deviations for, e.g., l-serine. Only the mutation variant B.1.1.7 shows a decreased stoichiometric coefficient compared to the wildtype. We analyzed the documented mutations for this variant and identified two mutations in structural proteins, Spike S982A and N S235F, which only occur in this variant. In both cases, the amino acid l-serine is substituted by another amino acid. As both mutations occur in structural proteins with copy numbers of 120 and 456, respectively, their influence on the amount of amino acid and, thus, the stoichiometric coefficient, is noticeable. Compared to the other mutation variants, variant B.1.1.28 has the highest increase in the stoichiometric coefficient for l-serine. This could be explained by two mutations specific for this variant in the structural spike protein: Spike P26S and Spike R190S. In both cases, other amino acids are replaced by l-serine. As explained for the mutation variant B.1.1.7, the spike protein has a copy number of 120. Changes in these structural proteins can be measurable and influence the stoichiometric coefficient stronger than mutations in non-structural proteins.
The mutation variant B.1.617 does not fit in this pattern. As the variants B.1.315 and B.1.427/429, its stoichiometric coefficients for ADP and ATP are increased, but not as much. Variant B.1.617 has a similar total viral molar mass as B.1.1.7. However, the amount of the nucleotides adenine and uridine is more similar to the variants B.1.315 and B.1.427/429. Having similarities with both pattern groups, variant B.1.617 does fit in neither of the groups. Variant B.1.617 needs less l-isoleucine compared to the wildtype.

3.2.2. Analysis of the Effects of Single Gene Deletions

After highlighting the differences in the stoichiometric coefficients for the different mutation variants, we tested the robustness of our previously identified potential antiviral targets [22]. To do so, we repeated the single-gene-deletion experiments for every mutation variant. Our analysis revealed that in all mutation variants, the guanylate kinase 1 (GK1) is a robust potential antiviral target.

3.3. Lipids as Part of the Viral Biomass Objective Function

The transmembrane domain of the envelope (E) protein is located in lipid bilayers mimicking the ERGIC membrane [21]. Ref. [20] described this ERGIC membrane [20] in 1994. The four phospholipids, phosphatidylcholine, phosphatidylethanolamine, phosphatidylinositol, and phosphatidylserine, were observed in the ERGIC while sphingomyelin and cardiolipin were not present [20]. Ref. [21] use an ERGIC-mimetic consisting of the four described phospholipids and cholesterol to investigate the E-protein’s transmembrane domain [21]. The five lipids are also participating in the macrophage’s maintenance function. Thus, their role and influence on the VBOF and antiviral targets were examined.
As the actual amount of lipids in the SARS-CoV-2 virion is not yet determined, we evaluated varying stoichiometric coefficients. In the first experiments, the individual lipids’ effect on the VBOF’s objective value was analyzed. The objective coefficients from the macrophage’s maintenance function varied between 0.001 for phosphatidylserine and 0.031 for phosphatidylcholine. Therefore, we first varied all lipids’ coefficients between 0 and 0.5 and subsequently used a multiplication coefficient between 0 and 10 to multiply the macrophage’s coefficients. Despite an up to 490-fold increase of the stoichiometric coefficient (for phosphatidylserine) compared to its initial value in the macrophage’s maintenance function, the VBOF’s objective value remained at 0.01886 mmol/(gDW · h). This was also the case when all five lipids were added to the VBOF  simultaneously.
Knock-out experiments were conducted to identify additional potential antiviral targets. All lipids were included in the VBOF, and the coefficients were varied using a multiplication coefficient. At the five-fold increase of the initial stoichiometric coefficients, two novel reactions emerged as new potential antiviral targets: the methionine synthase (METS) and the 5,10-methylenetetrahydrofolate reductase (FADH2) (MTHFR). To identify, which lipids are responsible for the emergence of the novel antiviral target, we repeated the described analysis for every lipid individually, once using absolute stoichiometric coefficients ranging from 0 to 0.5 and once using the above-described multiplication coefficient ranging between 0 and 10. By this approach, we identified phosphatidylcholine to be the responsible lipid for the additional antiviral targets. When increasing the initial macrophage’s stoichiometric coefficient of phosphatidylcholine by at least 4.76, the two enzymes emerge as potential antiviral targets. At a five-fold increase of phosphatidylcholine and the knock-out of either the methionine synthase or the 5,10-methylenetetrahydrofolate reductase (FADH2), the viral growth can be inhibited by approximately 1.5%. With increasing amounts of phosphatidylcholine in the VBOF, the knock-out influence of the two reactions on the viral growth increases, as seen in Figure 5: at an eleven-fold increase of phosphatidylcholine, the viral growth rate is decreased by approximately 50%. A twenty-fold increase of phosphatidylcholine inhibits the viral growth even to 30% of its initial growth rate.
It needs to be highlighted that the guanylate kinase 1 (GK1) was a potential antiviral target during all conducted in silico experiments evaluating the lipids’ effect on potential targets.

4. Discussion

This study presents an updated viral biomass objective function (VBOF) for the novel coronavirus SARS-CoV-2 based on the latest information of its structural proteins. This VBOF was integrated into an already validated model of human alveolar macrophages [23].
The tissue tropism of SARS-CoV-2 comprises most cell types expressing the entry receptor ACE-2, mainly including cell types of the lung, liver, stomach, ileum, kidney, and colon [39,40]. Although SARS-CoV-2 enters the host via the airways, the expression of ACE-2 is comparably low, highlighting the role of possible co-receptors [40]. Nonetheless, human alveolar type 2 cells robustly express ACE-2, while alveolar macrophages possibly express low levels of the entry receptor [40]. It is known that different coronaviruses infect macrophages, such as the human coronavirus strain 229E [41], the Middle East Respiratory Syndrome (MERS) coronavirus [42], and the SARS coronavirus [43]. Also, the novel coronavirus SARS-CoV-2 is reported to infect alveolar macrophages [44]. However, other in vitro studies suggest that challenging alveolar macrophages with SARS-CoV-2 does not lead to a productive infection [45]. However, even without productive infection, alveolar macrophages could serve as Trojan horses, which enable viral anchoring within pulmonary parenchyma [39]. Ref. [45] demonstrate that the tissue-resident alveolar macrophages play a crucial role in SARS-CoV-2 immune evasion [44,45] and are hypothesized to support viral pathogenesis [39]. Disabling viral replication in human alveolar macrophages might be an early way of intervention and prevention of the virus’s further spread.
We corrected the copy number of structural proteins and the stoichiometric coefficients in the viral biomass objective function (VBOF). The amount of the spike (S) and nucleocapsid (N) proteins were derived from studies on SARS-CoV-2 [14,15,32]. The copy number of the envelope (E) protein is derived from the human coronavirus and the transmissible gastroenteritis virus [29]. Numbers for SARS-CoV-2 are currently not available. Same accounts for the copy number of membrane (M) proteins, where information is only available for SARS-CoV-1 [34]. Especially for the M proteins, a range of potential copy numbers exists, as the ratio of M and N proteins ranges from 3M:1N to 1M:1N [34]. With the N protein’s copy number of 456, the M protein’s copy number ranges from 456 to 1368. As soon as additional information on the copy numbers of the E and M protein is available for SARS-CoV-2, the stoichiometric coefficients can be refined further.
However, the current refinement still confirmed the guanylate kinase 1 (GK1) as a potential antiviral target. Even for the investigated mutations, the guanylate kinase seems to be a robust target in human alveolar macrophages to interrupt SARS-CoV-2 replication. Ref. [46] conduct a similar study with the human reconstruction RECON2.2 [47] containing a lung biomass objective function and a viral biomass objective function [46]. They also report the guanylate kinase as a potential target for antiviral therapies [46]. In our previous study, we suggested potential drugs that could be repurposed to fight this SARS-CoV-2 pandemic. Amongst these drugs were cidofovir, brincidofovir, and favipiravir [22]. A virtual screening method identified cidofovir as a potentially effective therapeutic against SARS-CoV-2 [48]. A molecular docking study suggests the repurposing of brincidofovir against SARS-CoV-2 [49]. For favipiravir, several clinical trials are listed in the ClinicalTrials database hosted by the U.S. National Library of Medicine [50], running in several countries, including Italy (NCT04336904), Turkey (NCT04474457), and the United States (NCT04358549). However, these therapeutics are only analogs and do not directly inhibit the guanylate kinase. No direct inhibitor of the guanylate kinase is tested for its antiviral effect on SARS-CoV-2 infections at the time of writing. As the guanylate kinase is a robust target for all currently occurring mutation variants, further investigations could be of high interest to fight this pandemic.
Besides the guanylate kinase, additional potential antiviral targets were identified using the host-derived enforcement analysis. These antiviral targets are located in the pentose phosphate pathway, the purine, and the pyrimidine metabolism. It is shown that the pentose phosphate pathway is remarkably deregulated during SARS-CoV-2 replication, which shows potential implications for antiviral therapies [51]. The purine biosynthesis pathway is enhanced upon SARS-CoV-2 infection to support the de novo synthesis of purines [52]. First in vitro experiments show that the FDA-approved inhibitor of purine biosynthesis methotrexate potently inhibits viral replication [53,54], protein synthesis, and release [53]. The pyrimidine metabolism is also reported as a potential antiviral target, especially the dihydroorotate dehydrogenase. Its inhibition by, for example, brequinar or leflunomide is already demonstrated to have antiviral activity against other viruses [55,56,57], such as rotavirus [58] and Ebola virus [59]. The dihydroorotate dehydrogenase inhibitor PTC299 is shown to arrest SARS-CoV-2 replication in vitro [60]. The dihydroorotate dehydrogenase inhibitors S312 and S416 are validated to have high antiviral efficacy in vivo [61]. To conclude, our identified antiviral targets are currently under discussion in the scientific community, and for some, the influence and relevance for viral replication are confirmed.
Analyses of the documented mutations revealed that virus variant B.1.617 needs less l-isoleucine compared to the wildtype because of a mutation in the membrane protein, M I82T, where l-isoleucine is substituted by l-threonine. As the membrane protein has a copy number of 1000, its replacement could influence the stoichiometric coefficient of the replaced amino acid. Same accounts for the mutation N D402H in the nucleocapsid protein, where l-aspartate is replaced by l-histidine, which might explain the decreased stoichiometric coefficient for l-aspartate. Changes in these structural proteins can be measurable and influence the stoichiometric coefficient stronger than mutations in non-structural proteins.
Alongside the mutation variants that could complicate the fight against SARS-CoV-2 with vaccines, the S protein’s glycosylation could impact antibodies’ ability to bind to a pathogenic S glycoprotein by shielding its surface [62,63]. Currently, this glycosylation process is not reflected in the VBOF or the model. As soon as more information about the glycosylation is available that can be used to determine a range or precise stoichiometric coefficients, the glycosylation of the spike protein can be incorporated into the model simulations.
The inclusion of lipids in the VBOF opens new perspectives for potential antiviral targets. It is shown that virus infections can dramatically impact on lipid metabolism [64,65,66,67]. Upon rhinovirus infection multiple lipid pathways are altered, and changes in phospholipids, lysophospholipids, fatty acids, and inositol phospholipids are observed [66]. For the human coronavirus 229E (hCoV-229E), the host cell lipid response upon infection was comprehensively characterized. Glycerophospholipids and fatty acids were significantly elevated. Lysophosphatidylcholine, which is hydrolyzed from phosphatidylcholine, was significantly elevated and accounted for approximately 60% of all identified lipids with significant elevation [68]. Our study also highlighted phosphatidylcholine as an essential lipid upon SARS-CoV-2 infection, confirming the findings from [68] for hCoV-229E. As metabolic alterations harbor potential antiviral targets, regulating or targeting the lipid metabolism is suggested and discussed [64,66,68]. We identified two novel potential antiviral targets connected with lipid metabolism: the methionine synthase and the 5,10-methylenetetrahydrofolate reductase (FADH2). S-adenosyl-l-methionine is a pivotal methyl donor in the synthesis of phosphatidylcholine [69,70]. Thus, the synthesis of l-methionine by the 5,10-methylenetetrahydrofolate reductase (FADH2) and methionine synthase seem to be an antiviral target to disrupt the synthesis of phosphatidylcholine. These novel insights could guide further laboratory experiments for investigating and validating the lipid’s role in SARS-CoV-2 infections.
This study confirmed the guanylate kinase 1 (GK1) as a robust antiviral target against SARS-CoV-2 and its arising mutation variants. With the refined copy numbers of structural proteins, the list of further potential antiviral targets was improved, and some targets are already under discussion or even under validation. The inclusion of the lipids into the VBOF opened new perspectives for additional metabolic targets to fight against this pandemic.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/genes12060796/s1, Table S1: List of tested variants. The tested variants are listed together with their GISAID accession number and the observed mutations. The stoichiometric coefficients for the VBOF’s compounds in the different virus mutations are listed.

Author Contributions

Conceptualization, A.R. and A.D.; methodology, A.R. and L.W.; investigation, A.R. and L.W.; software, A.R.; visualization, A.R.; supervision, A.D.; funding acquisition, A.D.; writing—original draft preparation, A.R. and L.W.; writing—review and editing, A.R., L.W. and A.D; All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by the Federal Ministry of Education and Research (BMBF) and the Baden-Württemberg Ministry of Science as part of the Excellence Strategy of the German Federal and State Governments, by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) under Germany’s Excellence Strategy—EXC 2124—390838134, and supported by the German Center for Infection Research (DZIF, doi:10.13039/100009139) within the Deutsche Zentren der Gesundheitsforschung (BMBF-DZG, German Centers for Health Research of the Federal Ministry of Education and Research), grant No 8020708703. The authors acknowledge support by the Open Access Publishing Fund of the University of Tübingen (https://uni-tuebingen.de/en/58988, accessed on 21 May 2021).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The genome-scale metabolic model of the human alveolar macrophage infected with SARS-CoV-2 is available in the BioModels Database [71] as an SBML Level 3 Version 1 file [72,73,74] with the flux balance constraints (fbc) extension package [75] within a COMBINE Archive OMEX file [76] under the accession number MODEL2003020001. A supplementary table in Microsoft Excel format is available as Supplementary Materials along with this article.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
ACE-2angiotensin converting enzyme 2
ADPadenosine diphosphate
ADSL1adenylosuccinate lyase 1
ADSL2adenylosuccinate lyase 2
ADSSadenylosuccinate synthase
AICARTphosphoribosylaminoimidazolecarboxamide formyltransferase
AIRCrPhosphoribosylaminoimidazole carboxylase
AMPadenosine monophospate
ASPCTraspartate carbamoyltransferase
ATPadenosine triphosphate
CBPScarbamoyl-phosphate synthase
COMBINEComputational Modeling in Biology Network
COVID-19Coronavirus Disease 19
CspCopy number of structural proteins
DHORD9dihydoorotic acid dehydrogenase
DHORTSdihydroorotase
Eenvelope
ERGICendoplasmic reticulum–Golgi intermediate compartment
FBAflux balance analysis
fbcflux balance constraints
FVAflux variability analysis
GARFTphosphoribosylglycinamide formyltransferase
GEMgenome-scale metabolic model
GISAIDGlobal Initiative on Sharing All Influenza Data
GK1guanylate kinase 1
GLUPRTglutamine phosphoribosyldiphosphate amidotransferase
hCOVhuman coronavirus
hCoV-229Ehuman coronavirus 229E
IDidentifier
IMPinosine monophosphate
IMPCIMP cyclohydrolase
Mmembrane
MERSMiddle East Respiratory Syndrome
Mvviral molar mass
Nnucleocapsid
OMEXOpen Modeling EXchange format
OMPDCorotidine-5’-phosphate decarboxylase
ORPTorotate phosphoribosyltransferase
PPidiphosphate
PRAGSrphosphoribosylglycinamide synthase
PRAISphosphoribosylaminoimidazole synthase
PRASCSphosphoribosylaminoimidazolesuccinocarboxamide synthase
PRFGSphosphoribosylformylglycinamidine synthase
PRPPSphosphoribosylpyrophosphate synthetase
RPIribose-5-phosphate isomerase
Sspike
SARSSevere Acute Respiratory Syndrome
SARS-CoV-2Severe Acute Respiratory Syndrome coronavirus type 2
SBMLSystems Biology Markup Language
TGEVtransmissible gastroenteritis virus
VBOFviral biomass objective function
vRNPviral ribonucleoprotein complexes
WTwildtype

Appendix A

Table A1. Reactions from the host-derived enforcement experiments. The reaction identifiers listed in Figure 2 are BiGG identifiers [28]. In this table, the BiGG reaction identifiers are given, together with the reaction name and the subsystem, they occur in.
Table A1. Reactions from the host-derived enforcement experiments. The reaction identifiers listed in Figure 2 are BiGG identifiers [28]. In this table, the BiGG reaction identifiers are given, together with the reaction name and the subsystem, they occur in.
Reaction-IDReaction NameSubsystem
ADSL1adenylosuccinate lyase 1Purine metabolism
ADSL2adenylosuccinate lyase 2Purine metabolism
ADSSadenylosuccinate synthasePurine metabolism
AICARTphosphoribosylaminoimidazolecarboxamide formyltransferasePurine metabolism
AIRCrPhosphoribosylaminoimidazole carboxylasePurine metabolism
ASPCTraspartate carbamoyltransferasePyrimidine metabolism
CBPScarbamoyl-phosphate synthasePyrimidine metabolism
DHORD9dihydoorotic acid dehydrogenasePyrimidine metabolism
DHORTSdihydroorotasePyrimidine metabolism
GARFTphosphoribosylglycinamide formyltransferasePurine metabolism
GK1guanylate kinase 1Purine metabolism
GLUPRTglutamine phosphoribosyldiphosphate amidotransferasePurine metabolism
IMPCIMP cyclohydrolasePurine metabolism
OMPDCorotidine-5’-phosphate decarboxylasePyrimidine metabolism
ORPTorotate phosphoribosyltransferasePyrimidine metabolism
PRAGSrphosphoribosylglycinamide synthasePurine metabolism
PRAISphosphoribosylaminoimidazole synthasePurine metabolism
PRASCSphosphoribosylaminoimidazolesuccinocarboxamide synthasePurine metabolism
PRFGSphosphoribosylformylglycinamidine synthasePurine metabolism
PRPPSphosphoribosylpyrophosphate synthetasePentose phosphate pathway
RPIribose-5-phosphate isomerasePentose phosphate pathway

References

  1. Zhou, P.; Yang, X.L.; Wang, X.G.; Hu, B.; Zhang, L.; Zhang, W.; Si, H.R.; Zhu, Y.; Li, B.; Huang, C.L.; et al. A pneumonia outbreak associated with a new coronavirus of probable bat origin. Nature 2020, 579, 270–273. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Mirtaleb, M.S.; Mirtaleb, A.H.; Nosrati, H.; Heshmatnia, J.; Falak, R.; Zolfaghari Emameh, R. Potential therapeutic agents to COVID-19: An update review on antiviral therapy, immunotherapy, and cell therapy. Biomed. Pharmacother. 2021, 138, 111518. [Google Scholar] [CrossRef] [PubMed]
  3. Li, D.D.; Li, Q.H. SARS-CoV-2: Vaccines in the pandemic era. Mil. Med. Res. 2021, 8, 1. [Google Scholar] [CrossRef] [PubMed]
  4. Zhao, J.; Zhao, S.; Ou, J.; Zhang, J.; Lan, W.; Guan, W.; Wu, X.; Yan, Y.; Zhao, W.; Wu, J.; et al. COVID-19: Coronavirus Vaccine Development Updates. Front. Immunol. 2020, 11, 602256. [Google Scholar] [CrossRef] [PubMed]
  5. Liu, W.; Zhou, P.; Chen, K.; Ye, Z.; Liu, F.; Li, X.; He, N.; Wu, Z.; Zhang, Q.; Gong, X.; et al. Efficacy and safety of antiviral treatment for COVID-19 from evidence in studies of SARS-CoV-2 and other acute viral infections: A systematic review and meta-analysis. CMAJ 2020, 192, E734–E744. [Google Scholar] [CrossRef]
  6. Valle, C.; Martin, B.; Touret, F.; Shannon, A.; Canard, B.; Guillemot, J.C.; Coutard, B.; Decroly, E. Drugs against SARS-CoV-2: What do we know about their mode of action? Rev. Med. Virol. 2020, 30, 1–10. [Google Scholar] [CrossRef]
  7. Grein, J.; Ohmagari, N.; Shin, D.; Diaz, G.; Asperges, E.; Castagna, A.; Feldt, T.; Green, G.; Green, M.L.; Lescure, F.X.; et al. Compassionate Use of Remdesivir for Patients with Severe Covid-19. N. Engl. J. Med. 2020, 382, 2327–2336. [Google Scholar] [CrossRef]
  8. Tregoning, J.S.; Brown, E.S.; Cheeseman, H.M.; Flight, K.E.; Higham, S.L.; Lemm, N.M.; Pierce, B.F.; Stirling, D.C.; Wang, Z.; Pollock, K.M. Vaccines for COVID-19. Clin. Exp. Immunol. 2020, 202, 162–192. [Google Scholar] [CrossRef]
  9. European Centre for Disease Prevention and Control. Risk Related to Spread of New SARS-CoV-2 Variants of Concern in the EU/EEA, First Update; Technical Report; European Centre for Disease Prevention and Control: Stockholm, Sweden, 2021.
  10. Centers for Disease Control and Prevention (US). SARS-CoV-2 Variants of Concern; U.S. Department of Health and Human Services, Centers for Disease Control and Prevention: Atlanta, GA, USA, 2021.
  11. European Centre for Disease Prevention and Control. Emergence of SARS-CoV-2 B.1.617 Variants in India and Situation in the EU/EEA—11 May 2021; Technical Report; European Centre for Disease Prevention and Control: Stockholm, Sweden, 2021.
  12. Cui, J.; Li, F.; Shi, Z.L. Origin and evolution of pathogenic coronaviruses. Nat. Rev. Microbiol. 2019, 17, 181–192. [Google Scholar] [CrossRef] [Green Version]
  13. Bar-On, Y.M.; Flamholz, A.; Phillips, R.; Milo, R. Sars-cov-2 (Covid-19) by the numbers. eLife 2020, 9, e57309. [Google Scholar] [CrossRef]
  14. Klein, S.; Cortese, M.; Winter, S.L.; Wachsmuth-Melm, M.; Neufeldt, C.J.; Cerikan, B.; Stanifer, M.L.; Boulant, S.; Bartenschlager, R.; Chlanda, P. SARS-CoV-2 structure and replication characterized by Situ Cryo-Electron Tomogr. Nat. Commun. 2020, 11, 5885. [Google Scholar] [CrossRef] [PubMed]
  15. Turoňová, B.; Sikora, M.; Schürmann, C.; Hagen, W.J.; Welsch, S.; Blanc, F.E.; von Bülow, S.; Gecht, M.; Bagola, K.; Hörner, C.; et al. Mediat. Three Hinges. Science 2020, 370, 203–208. [Google Scholar] [CrossRef] [PubMed]
  16. Lenard, J. Viral membranes. In Encyclopedia of Virology; Springer: Boston, MA, USA, 2008; pp. 308–314. [Google Scholar] [CrossRef]
  17. de Armas-Rillo, L.; Valera, M.S.; Marrero-Hernández, S.; Valenzuela-Fernández, A. Membrane dynamics associated with viral infection. Rev. Med. Virol. 2016, 26, 146–160. [Google Scholar] [CrossRef] [PubMed]
  18. Miller, S.; Krijnse-Locker, J. Modification of intracellular membrane structures for virus replication. Nat. Rev. Microbiol. 2008, 6, 363–374. [Google Scholar] [CrossRef] [PubMed]
  19. Ketter, E.; Randall, G. Virus Impact on Lipids and Membranes. Annu. Rev. Virol. 2019, 6, 319–340. [Google Scholar] [CrossRef]
  20. Schweizer, A.; Clausen, H.; Van Meer, G.; Hauri, H.P. Localization of O-glycan initiation, sphingomyelin synthesis, and glucosylceramide synthesis in Vero cells with respect to the endoplasmic reticulum-Golgi intermediate compartment. J. Biol. Chem. 1994, 269, 4035–4041. [Google Scholar] [CrossRef]
  21. Mandala, V.S.; McKay, M.J.; Shcherbakov, A.A.; Dregni, A.J.; Kolocouris, A.; Hong, M. Structure and drug binding of the SARS-CoV-2 envelope protein transmembrane domain in lipid bilayers. Nat. Struct. Mol. Biol. 2020, 27, 1202–1208. [Google Scholar] [CrossRef]
  22. Renz, A.; Widerspick, L.; Dräger, A. FBA reveals guanylate kinase as a potential target for antiviral therapies against SARS-CoV-2. Bioinformatics 2020, 36, i813–i821. [Google Scholar] [CrossRef]
  23. Bordbar, A.; Lewis, N.E.; Schellenberger, J.; Palsson, B.; Jamshidi, N. Insight into human alveolar macrophage and M. tuberculosis interactions via metabolic reconstructions. Mol. Syst. Biol. 2010, 6, 422. [Google Scholar] [CrossRef]
  24. Hible, G.; Daalova, P.; Gilles, A.M.; Cherfils, J. Crystal structures of GMP kinase in complex with ganciclovir monophosphate and Ap5G. Biochimie 2006, 88, 1157–1164. [Google Scholar] [CrossRef]
  25. Navé, J.F.; Eschbach, A.; Halazy, S. 9-(Phosphonoalkyl)guanine derivatives as substrates or inhibitors of guanylate kinase. Arch. Biochem. Biophys. 1992, 295, 253–257. [Google Scholar] [CrossRef]
  26. Navé, J.F.; Taylor, D.; Tyms, S.; Kenny, M.; Eggenspiller, A.; Eschbach, A.; Dulworth, J.; Brennan, T.; Piriou, F.; Halazy, S. Synthesis, antiviral activity and enzymatic phosphorylation of 9-phosphonopentenyl derivatives of guanine. Antivir. Res. 1995, 27, 301–316. [Google Scholar] [CrossRef]
  27. Aller, S.; Scott, A.; Sarkar-Tyson, M.; Soyer, O.S. Integrated human-virus metabolic stoichiometric modelling predicts host-based antiviral targets against Chikungunya, Dengue and Zika viruses. J. R. Soc. Interface 2018, 15, 20180125. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Norsigian, C.J.; Pusarla, N.; McConn, J.L.; Yurkovich, J.T.; Dräger, A.; Palsson, B.O.; King, Z. BiGG Models 2020: Multi-strain genome-scale models and expansion across the phylogenetic tree. Nucleic Acids Res. 2019, 48, D402–D406. [Google Scholar] [CrossRef] [PubMed]
  29. Tilocca, B.; Soggiu, A.; Sanguinetti, M.; Babini, G.; De Maio, F.; Britti, D.; Zecconi, A.; Bonizzi, L.; Urbani, A.; Roncada, P. Immunoinformatic analysis of the SARS-CoV-2 envelope protein as a strategy to assess cross-protection against COVID-19. Microbes Infect. 2020, 22, 182–187. [Google Scholar] [CrossRef]
  30. Stodola, J.K.; Dubois, G.; Le Coupanec, A.; Desforges, M.; Talbot, P.J. The OC43 human coronavirus envelope protein is critical for infectious virus production and propagation in neuronal cells and is a determinant of neurovirulence and CNS pathology. Virology 2018, 515, 134–149. [Google Scholar] [CrossRef]
  31. Godet, M.; L’Haridon, R.; Vautherot, J.F.; Laude, H. TGEV corona virus ORF4 encodes a membrane protein that is incorporated into virions. Virology 1992, 188, 666–675. [Google Scholar] [CrossRef]
  32. Yao, H.; Song, Y.; Chen, Y.; Wu, N.; Xu, J.; Sun, C.; Zhang, J.; Weng, T.; Zhang, Z.; Wu, Z.; et al. Molecular Architecture of the SARS-CoV-2 Virus. Cell 2020, 183, 730–738.e13. [Google Scholar] [CrossRef]
  33. Lu, S.; Ye, Q.; Singh, D.; Cao, Y.; Diedrich, J.K.; Yates, J.R.; Villa, E.; Cleveland, D.W.; Corbett, K.D. The SARS-CoV-2 nucleocapsid phosphoprotein forms mutually exclusive condensates with RNA and the membrane-associated M protein. Nat. Commun. 2021, 12, 1–15. [Google Scholar] [CrossRef]
  34. Neuman, B.W.; Kiss, G.; Kunding, A.H.; Bhella, D.; Baksh, M.F.; Connelly, S.; Droese, B.; Klaus, J.P.; Makino, S.; Sawicki, S.G.; et al. A structural analysis of M protein in coronavirus assembly and morphology. J. Struct. Biol. 2011, 174, 11–22. [Google Scholar] [CrossRef]
  35. Pedley, A.M.; Benkovic, S.J. A New View into the Regulation of Purine Metabolism: The Purinosome. Trends Biochem. Sci. 2017, 42, 141–154. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Löffler, M.; Fairbanks, L.D.; Zameitat, E.; Marinaki, A.M.; Simmonds, H.A. Pyrimidine pathways in health and disease. Trends Mol. Med. 2005, 11, 430–437. [Google Scholar] [CrossRef] [PubMed]
  37. Elbe, S.; Buckland-Merrett, G. Data, disease and diplomacy: GISAID’s innovative contribution to global health. Glob. Challenges 2017, 1, 33–46. [Google Scholar] [CrossRef] [Green Version]
  38. Shu, Y.; McCauley, J. GISAID: Global initiative on sharing all influenza data—From vision to reality. Eurosurveillance 2017, 22, 30494. [Google Scholar] [CrossRef] [Green Version]
  39. Abassi, Z.; Knaney, Y.; Karram, T.; Heyman, S.N. The Lung Macrophage in SARS-CoV-2 Infection: A Friend or a Foe? Front. Immunol. 2020, 11, 1312. [Google Scholar] [CrossRef]
  40. Qi, F.; Qian, S.; Zhang, S.; Zhang, Z. Single cell RNA sequencing of 13 human tissues identify cell types and receptors of human coronaviruses. Biochem. Biophys. Res. Commun. 2020, 526, 135–140. [Google Scholar] [CrossRef] [PubMed]
  41. Joel Funk, C.; Wang, J.; Ito, Y.; Travanty, E.A.; Voelker, D.R.; Holmes, K.V.; Mason, R.J. Infection of human alveolar macrophages by human coronavirus strain 229E. J. Gen. Virol. 2012, 93, 494–503. [Google Scholar] [CrossRef] [PubMed]
  42. Al-Qahtani, A.A.; Lyroni, K.; Aznaourova, M.; Tseliou, M.; Al-Anazi, M.R.; Al-Ahdal, M.N.; Alkahtani, S.; Sourvinos, G.; Tsatsanis, C. Middle east respiratory syndrome corona virus spike glycoprotein suppresses macrophage responses via DPP4-mediated induction of IRAK-M and PPARγ. Oncotarget 2017, 8, 9053–9066. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Gu, J.; Gong, E.; Zhang, B.; Zheng, J.; Gao, Z.; Zhong, Y.; Zou, W.; Zhan, J.; Wang, S.; Xie, Z.; et al. Multiple organ infection and the pathogenesis of SARS. J. Exp. Med. 2005, 202, 415–424. [Google Scholar] [CrossRef]
  44. Chu, H.; Chan, J.F.W.; Wang, Y.; Yuen, T.T.T.; Chai, Y.; Hou, Y.; Shuai, H.; Yang, D.; Hu, B.; Huang, X.; et al. Comparative Replication and Immune Activation Profiles of SARS-CoV-2 and SARS-CoV in Human Lungs: An Ex Vivo Study with Implications for the Pathogenesis of COVID-19. Clin. Infect. Dis. 2020, 71, 1400–1409. [Google Scholar] [CrossRef] [Green Version]
  45. Dalskov, L.; Møhlenberg, M.; Thyrsted, J.; Blay-Cadanet, J.; Poulsen, E.T.; Folkersen, B.H.; Skaarup, S.H.; Olagnier, D.; Reinert, L.; Enghild, J.J.; et al. SARS-CoV-2 evades immune detection in alveolar macrophages. EMBO Rep. 2020, 21, e51252. [Google Scholar] [CrossRef] [PubMed]
  46. Delattre, H.; Sasidharan, K.; Soyer, O.S. Inhibiting the reproduction of SARS-CoV-2 through perturbations in human lung cell metabolic network. Life Sci. Alliance 2021, 4, e202000869. [Google Scholar] [CrossRef] [PubMed]
  47. Swainston, N.; Smallbone, K.; Hefzi, H.; Dobson, P.D.; Brewer, J.; Hanscho, M.; Zielinski, D.C.; Ang, K.S.; Gardiner, N.J.; Gutierrez, J.M.; et al. Recon 2.2: From reconstruction to model of human metabolism. Metabolomics 2016, 12, 109. [Google Scholar] [CrossRef] [PubMed]
  48. Rahman, M.R.; Banik, A.; Chowdhury, I.M.; Sajib, E.H.; Sarkar, S. Identification of potential antivirals against SARS-CoV-2 using virtual screening method. Inform. Med. Unlocked 2021, 23, 100531. [Google Scholar] [CrossRef] [PubMed]
  49. Hussien, M.A.; Abdelaziz, A.E. Molecular docking suggests repurposing of brincidofovir as a potential drug targeting SARS-CoV-2 ACE2 receptor and main protease. Netw. Model. Anal. Health Inform. Bioinform. 2020, 9, 1–18. [Google Scholar] [CrossRef]
  50. National Library of Medicine (U.S.). Available online: https://clinicaltrials.gov/ (accessed on 29 May 2012).
  51. Bojkova, D.; Costa, R.; Bechtel, M.; Ciesek, S.; Michaelis, M.; Cinatl, J. Targeting pentose phosphate pathway for SARS-CoV-2 therapy. bioRxiv 2020. [Google Scholar] [CrossRef]
  52. Zhang, Y.; Guo, R.; Kim, S.H.; Shah, H.; Zhang, S.; Liang, J.H.; Fang, Y.; Gentili, M.; Leary, C.N.; Elledge, S.J.; et al. SARS-CoV-2 hijacks folate and one-carbon metabolism for viral replication. Nat. Commun. 2021, 12, 1676. [Google Scholar] [CrossRef]
  53. Caruso, A.; Caccuri, F.; Bugatti, A.; Zani, A.; Vanoni, M.; Bonfanti, P.; Cazzaniga, M.E.; Perno, C.F.; Messa, C.; Alberghina, L. Methotrexate inhibits SARS-CoV-2 virus replication “in vitro”. J. Med. Virol. 2021, 93, 1780–1785. [Google Scholar] [CrossRef]
  54. Stegmann, K.M.; Dickmanns, A.; Gerber, S.; Nikolova, V.; Klemke, L.; Manzini, V.; Schlösser, D.; Bierwirth, C.; Freund, J.; Sitte, M.; et al. The folate antagonist methotrexate diminishes replication of the coronavirus SARS-CoV-2 and enhances the antiviral efficacy of remdesivir in cell culture models. bioRxiv 2020. [Google Scholar] [CrossRef]
  55. Coelho, A.R.; Oliveira, P.J. Dihydroorotate dehydrogenase inhibitors in SARS-CoV-2 infection. Eur. J. Clin. Investig. 2020, 50, e13366. [Google Scholar] [CrossRef]
  56. Hoffmann, H.H.; Kunz, A.; Simon, V.A.; Palese, P.; Shaw, M.L. Broad-spectrum antiviral that interferes with De Novo Pyrimidine Biosynthesis. Proc. Natl. Acad. Sci. USA 2011, 108, 5777–5782. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Cheung, N.N.; Lai, K.K.; Dai, J.; Kok, K.H.; Chen, H.; Chan, K.H.; Yuen, K.Y.; Tsun Kao, R.Y. Broad-spectrum inhibition of common respiratory RNA viruses by a pyrimidine synthesis inhibitor with involvement of the host antiviral response. J. Gen. Virol. 2017, 98, 946–954. [Google Scholar] [CrossRef] [PubMed]
  58. Chen, S.; Ding, S.; Yin, Y.; Xu, L.; Li, P.; Peppelenbosch, M.P.; Pan, Q.; Wang, W. Suppression of pyrimidine biosynthesis by targeting DHODH enzyme robustly inhibits rotavirus replication. Antivir. Res. 2019, 167, 35–44. [Google Scholar] [CrossRef] [PubMed]
  59. Luthra, P.; Naidoo, J.; Pietzsch, C.A.; De, S.; Khadka, S.; Anantpadma, M.; Williams, C.G.; Edwards, M.R.; Davey, R.A.; Bukreyev, A.; et al. Inhibiting pyrimidine biosynthesis impairs Ebola virus replication through depletion of nucleoside pools and activation of innate immune responses. Antivir. Res. 2018, 158, 288–302. [Google Scholar] [CrossRef] [PubMed]
  60. Luban, J.; Sattler, R.A.; Mühlberger, E.; Graci, J.D.; Cao, L.; Weetall, M.; Trotta, C.; Colacino, J.M.; Bavari, S.; Strambio-De-Castillia, C.; et al. The DHODH inhibitor PTC299 arrests SARS-CoV-2 replication and suppresses induction of inflammatory cytokines. Virus Res. 2021, 292, 198246. [Google Scholar] [CrossRef]
  61. Xiong, R.; Zhang, L.; Li, S.; Sun, Y.; Ding, M.; Wang, Y.; Zhao, Y.; Wu, Y.; Shang, W.; Jiang, X.; et al. Novel and potent inhibitors targeting DHODH are broad-spectrum antivirals against RNA viruses including newly-emerged coronavirus SARS-CoV-2. Protein Cell 2020, 11, 723–739. [Google Scholar] [CrossRef] [PubMed]
  62. Grant, O.C.; Montgomery, D.; Ito, K.; Woods, R.J. Analysis of the SARS-CoV-2 spike protein glycan shield reveals implications for immune recognition. Sci. Rep. 2020, 10, 14991. [Google Scholar] [CrossRef]
  63. Watanabe, Y.; Allen, J.D.; Wrapp, D.; McLellan, J.S.; Crispin, M. Site-specific glycan analysis of the SARS-CoV-2 spike. Science 2020, 369, 330–333. [Google Scholar] [CrossRef]
  64. Abu-Farha, M.; Thanaraj, T.A.; Qaddoumi, M.G.; Hashem, A.; Abubaker, J.; Al-Mulla, F. The Role of Lipid Metabolism in COVID-19 Virus Infection and as a Drug Target. Int. J. Mol. Sci. 2020, 21, 3544. [Google Scholar] [CrossRef]
  65. Wu, Q.; Zhou, L.; Sun, X.; Yan, Z.; Hu, C.; Wu, J.; Xu, L.; Li, X.; Liu, H.; Yin, P.; et al. Altered Lipid Metabolism in Recovered SARS Patients Twelve Years after Infection. Sci. Rep. 2017, 7, 9110. [Google Scholar] [CrossRef]
  66. Nguyen, A.; Guedán, A.; Mousnier, A.; Swieboda, D.; Zhang, Q.; Horkai, D.; Le Novere, N.; Solari, R.; Wakelam, M.J. Host lipidome analysis during rhinovirus replication in HBECs identifies potential therapeutic target. J. Lipid Res. 2018, 59, 1671–1684. [Google Scholar] [CrossRef] [Green Version]
  67. Yan, B.; Zou, Z.; Chu, H.; Chan, G.; Tsang, J.O.L.; Lai, P.M.; Yuan, S.; Yip, C.C.Y.; Yin, F.; Kao, R.Y.T.; et al. Lipidomic Profiling Reveals Significant Perturbations of Intracellular Lipid Homeostasis in Enterovirus-Infected Cells. Int. J. Mol. Sci. 2019, 20, 5952. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  68. Yan, B.; Chu, H.; Yang, D.; Sze, K.H.; Lai, P.M.; Yuan, S.; Shuai, H.; Wang, Y.; Kao, R.Y.T.; Chan, J.F.W.; et al. Characterization of the lipidomic profile of human coronavirus-infected cells: Implications for lipid metabolism remodeling upon coronavirus replication. Viruses 2019, 11, 73. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  69. Visioli, F.; Colombo, C.; Monti, S.; Giulidori, P.; Galli, C. S-adenosyl-L-methionine: Role in phosphatidylcholine synthesis and in vitro effects on the ethanol-induced alterations of lipid metabolism. Pharmacol. Res. 1998, 37, 203–206. [Google Scholar] [CrossRef] [PubMed]
  70. Ye, C.; Sutter, B.M.; Wang, Y.; Kuang, Z.; Tu, B.P. A Metabolic Function for Phospholipid and Histone Methylation. Mol. Cell 2017, 66, 180–193.e8. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  71. Malik-Sheriff, R.S.; Glont, M.; Nguyen, T.V.N.; Tiwari, K.; Roberts, M.G.; Xavier, A.; Vu, M.T.; Men, J.; Maire, M.; Kananathan, S.; et al. BioModels—15 years of sharing computational models in life science. Nucleic Acids Res. 2020, 48, D407–D415. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  72. Keating, S.M.; Waltemath, D.; König, M.; Zhang, F.; Dräger, A.; Chaouiya, C.; Bergmann, F.T.; Finney, A.; Gillespie, C.S.; Helikar, T.; et al. SBML Level 3: An extensible format for the exchange and reuse of biological models. Mol. Syst. Biol. 2020, 16, e9110. [Google Scholar] [CrossRef]
  73. Renz, A.; Mostolizadeh, R.; Dräger, A. Clinical Applications of Metabolic Models in SBML Format. In Systems Medicine; Wolkenhauer, O., Ed.; Academic Press: Oxford, UK, 2020; Volume 3, pp. 362–371. [Google Scholar] [CrossRef]
  74. Hucka, M.; Bergmann, F.T.; Dräger, A.; Hoops, S.; Keating, S.M.; Le Novère, N.; Myers, C.J.; Olivier, B.G.; Sahle, S.; Schaff, J.C.; et al. Systems Biology Markup Language (SBML) Level 3 Version 1 Core. J. Integr. Bioinform. 2018, 15, 1. [Google Scholar] [CrossRef] [Green Version]
  75. Olivier, B.G.; Bergmann, F.T. SBML Level 3 Package: Flux Balance Constraints version 2. J. Integr. Bioinform. 2018, 15. [Google Scholar] [CrossRef]
  76. Bergmann, F.T.; Adams, R.; Moodie, S.; Cooper, J.; Glont, M.; Golebiewski, M.; Hucka, M.; Laibe, C.; Miller, A.K.; Nickerson, D.P.; et al. COMBINE archive and OMEX format: One file to share all information to reproduce a modeling project. BMC Bioinform. 2014, 15, 369. [Google Scholar] [CrossRef]
Figure 1. Workflow for the investigation of lipids’ influence on the VBOF. The five lipids phosphatidylcholine (pchol), phosphatidylethanolamine (pe), phosphatidylinositol (pail), phosphatidylserine (ps), and cholesterol (chstol) were added together and individually to the VBOF. The stoichiometric coefficients were either an absolute value identical for all lipids, or the initial stoichiometric coefficient from the macrophage biomass function factorized with a multiplication-coefficient. For all scenarios, the influence of the different VBOFs on the objective value was analyzed. Additionally, potential antiviral targets were examined using reaction knock-outs.
Figure 1. Workflow for the investigation of lipids’ influence on the VBOF. The five lipids phosphatidylcholine (pchol), phosphatidylethanolamine (pe), phosphatidylinositol (pail), phosphatidylserine (ps), and cholesterol (chstol) were added together and individually to the VBOF. The stoichiometric coefficients were either an absolute value identical for all lipids, or the initial stoichiometric coefficient from the macrophage biomass function factorized with a multiplication-coefficient. For all scenarios, the influence of the different VBOFs on the objective value was analyzed. Additionally, potential antiviral targets were examined using reaction knock-outs.
Genes 12 00796 g001
Figure 2. Results of the host-derived enforcement experiments. With the help of the host-derived enforcement, the range and effect of reaction inhibitions on the VBOF can be investigated while keeping the host’s maintenance at 100%. The minimum possible reaction inhibition rate to reduce the viral replication capacity (VBOF) is given in dark blue. The maximum inhibition of the reaction does not harm the host’s maintenance and is indicated in light blue. The reduction of the VBOF is given in comparison to the un-inhibited state. All reaction identifiers are BiGG identifiers [28]. Table A1 lists all reaction identifiers with their corresponding reaction name and the subsystem they occur in.
Figure 2. Results of the host-derived enforcement experiments. With the help of the host-derived enforcement, the range and effect of reaction inhibitions on the VBOF can be investigated while keeping the host’s maintenance at 100%. The minimum possible reaction inhibition rate to reduce the viral replication capacity (VBOF) is given in dark blue. The maximum inhibition of the reaction does not harm the host’s maintenance and is indicated in light blue. The reduction of the VBOF is given in comparison to the un-inhibited state. All reaction identifiers are BiGG identifiers [28]. Table A1 lists all reaction identifiers with their corresponding reaction name and the subsystem they occur in.
Genes 12 00796 g002
Figure 3. Difference of stoichiometric coefficients between wildtype (WT) and all mutations. The stoichiometric coefficients for all metabolites participating in the viral biomass objective function (VBOF) are compared. WT stoichiometric coefficients are indicated in red, the mean stoichiometric coefficients of all mutation variants are indicated in orange, including standard deviations (black). If the difference of the stoichiometric coefficients between WT and mutation variants was more than 0.001, the difference is indicated above the bars. The stoichiometric coefficients for the metabolites ATP, ADP and PPi, are higher compared to the other coefficients. The mutation variants’ mean coefficients show little deviation. Additionally, the differences between the stoichiometric coefficients of WT and mutation variants are very small.
Figure 3. Difference of stoichiometric coefficients between wildtype (WT) and all mutations. The stoichiometric coefficients for all metabolites participating in the viral biomass objective function (VBOF) are compared. WT stoichiometric coefficients are indicated in red, the mean stoichiometric coefficients of all mutation variants are indicated in orange, including standard deviations (black). If the difference of the stoichiometric coefficients between WT and mutation variants was more than 0.001, the difference is indicated above the bars. The stoichiometric coefficients for the metabolites ATP, ADP and PPi, are higher compared to the other coefficients. The mutation variants’ mean coefficients show little deviation. Additionally, the differences between the stoichiometric coefficients of WT and mutation variants are very small.
Genes 12 00796 g003
Figure 4. Difference of stoichiometric coefficients between wildtype (WT) and the individual mutations. The deviation between WT and the mean of the individual mutation variants was calculated. Higher stoichiometric coefficients in the mutation compared to the WT are indicated in blue, while lower stoichiometric coefficients are indicated in red. Based on similar sequence length for the mutation variants B.1.1.7 and B.1.1.28 and resulting similar total viral molar masses, a pattern emerges, which is most apparent for the stoichiometric coefficients of ATP and ADP. This pattern, however, is not present for all stoichiometric coefficients. The coefficient for l-serine, for example, is only decreased in the mutation variant B.1.1.7 based on two mutations in two structural proteins. Overall, the deviations from the WT are very small.
Figure 4. Difference of stoichiometric coefficients between wildtype (WT) and the individual mutations. The deviation between WT and the mean of the individual mutation variants was calculated. Higher stoichiometric coefficients in the mutation compared to the WT are indicated in blue, while lower stoichiometric coefficients are indicated in red. Based on similar sequence length for the mutation variants B.1.1.7 and B.1.1.28 and resulting similar total viral molar masses, a pattern emerges, which is most apparent for the stoichiometric coefficients of ATP and ADP. This pattern, however, is not present for all stoichiometric coefficients. The coefficient for l-serine, for example, is only decreased in the mutation variant B.1.1.7 based on two mutations in two structural proteins. Overall, the deviations from the WT are very small.
Genes 12 00796 g004
Figure 5. Influence of stoichiometric coefficient on reduction of VBOF during METS knock-out. With increasing factorization of phosphatidylcholine’s stoichiometric coefficient, the objective value of the VBOF’s optimization decreases during the knock-out of the methionine synthase (METS) reaction. The hosts growth maintenance stays at 100%. At an eleven-fold increase of the initial stoichiometric coefficient extracted from the host’s maintenance function results in a 50% decrease of the viral growth rate.
Figure 5. Influence of stoichiometric coefficient on reduction of VBOF during METS knock-out. With increasing factorization of phosphatidylcholine’s stoichiometric coefficient, the objective value of the VBOF’s optimization decreases during the knock-out of the methionine synthase (METS) reaction. The hosts growth maintenance stays at 100%. At an eleven-fold increase of the initial stoichiometric coefficient extracted from the host’s maintenance function results in a 50% decrease of the viral growth rate.
Genes 12 00796 g005
Table 1. Stoichiometric coefficients of the five lipids in the macrophage’s maintenance function. The stoichiometric coefficients of the five lipids were extracted from the macrophage’s maintenance function. Additionally, the BiGG identifiers [28] of the lipids are given. These stoichiometric coefficients formed the starting point for evaluating the lipids’ influence on the viral biomass objective function (VBOF).
Table 1. Stoichiometric coefficients of the five lipids in the macrophage’s maintenance function. The stoichiometric coefficients of the five lipids were extracted from the macrophage’s maintenance function. Additionally, the BiGG identifiers [28] of the lipids are given. These stoichiometric coefficients formed the starting point for evaluating the lipids’ influence on the viral biomass objective function (VBOF).
LipidBiGG IDCoefficient
Phosphatidylcholinepchol_hs_c0.03152
Phosphatidylethanolaminepe_hs_c0.02110
Phosphatidylinositolpail_hs_c0.00374
Phosphatidylserineps_hs_c0.00102
Cholesterolchsterol_c0.02093
Table 2. Copy number of structural proteins. The Copy number of structural proteins (Csp) was determined based on extensive literature research. Besides the reference and the copy number of structural proteins, the investigated organism is given as a source.
Table 2. Copy number of structural proteins. The Copy number of structural proteins (Csp) was determined based on extensive literature research. Besides the reference and the copy number of structural proteins, the investigated organism is given as a source.
ProteinNameReferenceSourceCsp
SS[15]SARS-CoV-2120
EE[29]hCOV, TGEV20
NN[14,32]SARS-CoV-2456
MM[34]SARS-CoV-11000
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Renz , A.; Widerspick , L.; Dräger , A. Genome-Scale Metabolic Model of Infection with SARS-CoV-2 Mutants Confirms Guanylate Kinase as Robust Potential Antiviral Target. Genes 2021, 12, 796. https://doi.org/10.3390/genes12060796

AMA Style

Renz  A, Widerspick  L, Dräger  A. Genome-Scale Metabolic Model of Infection with SARS-CoV-2 Mutants Confirms Guanylate Kinase as Robust Potential Antiviral Target. Genes. 2021; 12(6):796. https://doi.org/10.3390/genes12060796

Chicago/Turabian Style

Renz , Alina, Lina Widerspick , and Andreas Dräger . 2021. "Genome-Scale Metabolic Model of Infection with SARS-CoV-2 Mutants Confirms Guanylate Kinase as Robust Potential Antiviral Target" Genes 12, no. 6: 796. https://doi.org/10.3390/genes12060796

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop