The COVID-19 pandemic continues to pose significant challenges to global public health, necessitating the development of effective vaccination strategies to mitigate disease transmission. In Thailand, the COVID-19 epidemic has undergone multiple waves, prompting the implementation of various control measures, including vaccination campaigns. Understanding the dynamics of disease transmission and the impact of vaccination strategies is crucial for guiding public health interventions and optimizing epidemic control efforts. In this study, we developed a comprehensive mathematical model, termed SSvIH1CH2RD, to elucidate the dynamics of the COVID-19 epidemic in Thailand. The model incorporates key epidemiological parameters, vaccination rates, and disease progression stages to assess the effectiveness of different vaccination strategies in curbing disease transmission. Parameter estimation and model fitting were conducted using real-world data from COVID-19 patients in Thailand, enabling the simulation of epidemic scenarios and the exploration of optimal vaccination rates. Our results showed that optimizing vaccination strategies, particularly by administering approximately 119,625 doses per day, can significantly reduce the basic reproduction number (R0) below 1, thereby accelerating epidemic control. Simulation results demonstrated that the optimal vaccination rate led to a substantial decrease in the number of infections, with the epidemic projected to be completely eradicated from the population by June 19, 2022. These findings underscore the importance of targeted vaccination efforts and proactive public health interventions in mitigating the spread of COVID-19 and minimizing the burden on healthcare systems. Our study provides valuable insights into the optimization of vaccination strategies for epidemic control, offering guidance for policymakers and healthcare authorities in Thailand and beyond. By leveraging mathematical modeling techniques and real-world data, stakeholders can develop evidence-based strategies to combat the COVID-19 pandemic and safeguard public health.
Citation: Pannathon Kreabkhontho, Watchara Teparos, Thitiya Theparod. Potential for eliminating COVID-19 in Thailand through third-dose vaccination: A modeling approach[J]. Mathematical Biosciences and Engineering, 2024, 21(8): 6807-6828. doi: 10.3934/mbe.2024298
[1] | Mostean Bahreinipour, Hajar Zarei, Fariba Dashtestani, Jamal Rashidiani, Khadijeh Eskandari, Seyed Ali Moussavi Zarandi, Susan Kabudanian Ardestani, Hiroshi Watabe . Radioprotective effect of nanoceria and magnetic flower-like iron oxide microparticles on gamma radiation-induced damage in BSA protein. AIMS Biophysics, 2021, 8(2): 124-142. doi: 10.3934/biophy.2021010 |
[2] | Thuy Hien Nguyen, Catherine C. Moore, Preston B. Moore, Zhiwei Liu . Molecular dynamics study of homo-oligomeric ion channels: Structures of the surrounding lipids and dynamics of water movement. AIMS Biophysics, 2018, 5(1): 50-76. doi: 10.3934/biophy.2018.1.50 |
[3] | Stephanie H. DeLuca, Samuel L. DeLuca, Andrew Leaver-Fay, Jens Meiler . RosettaTMH: a method for membrane protein structure elucidation combining EPR distance restraints with assembly of transmembrane helices. AIMS Biophysics, 2016, 3(1): 1-26. doi: 10.3934/biophy.2016.1.1 |
[4] | Michael B. Sherman, Juan Trujillo, Benjamin E. Bammes, Liang Jin, Matthias W. Stumpf, Scott C. Weaver . Decontamination of digital image sensors and assessment of electron microscope performance in a BSL-3 containment. AIMS Biophysics, 2015, 2(2): 153-162. doi: 10.3934/biophy.2015.2.153 |
[5] | Panisak Boonamnaj, Pornthep Sompornpisut, Piyarat Nimmanpipug, R.B. Pandey . Thermal denaturation of a coronavirus envelope (CoVE) protein by a coarse-grained Monte Carlo simulation. AIMS Biophysics, 2022, 9(4): 330-340. doi: 10.3934/biophy.2022027 |
[6] | Riyaz A. Mir, Jeff Lovelace, Nicholas P. Schafer, Peter D. Simone, Admir Kellezi, Carol Kolar, Gaelle Spagnol, Paul L. Sorgen, Hamid Band, Vimla Band, Gloria E. O. Borgstahl . Biophysical characterization and modeling of human Ecdysoneless (ECD) protein supports a scaffolding function. AIMS Biophysics, 2016, 3(1): 195-210. doi: 10.3934/biophy.2016.1.195 |
[7] | Ashwani Kumar Vashishtha, William H. Konigsberg . Effect of different divalent cations on the kinetics and fidelity of DNA polymerases. AIMS Biophysics, 2018, 5(4): 272-289. doi: 10.3934/biophy.2018.4.272 |
[8] | Vilas K. Mahajan, Sandip P. Patil, Shirish H. Sonawane, Gunvant H. Sonawane . Ultrasonic, photocatalytic and sonophotocatalytic degradation of Basic Red-2 by using Nb2O5 nano catalyst. AIMS Biophysics, 2016, 3(3): 415-430. doi: 10.3934/biophy.2016.3.415 |
[9] | Ganesh Prasad Tiwari, Santosh Adhikari, Hari Prasad Lamichhane, Dinesh Kumar Chaudhary . Natural bond orbital analysis of dication magnesium complexes [Mg(H2O)6]2+ and [[Mg(H2O)6](H2O)n]2+; n=1-4. AIMS Biophysics, 2023, 10(1): 121-131. doi: 10.3934/biophy.2023009 |
[10] | Roman Marks, Piotr H. Pawłowski . Rotational-electric principles of RNA/DNA and viability. AIMS Biophysics, 2023, 10(3): 385-400. doi: 10.3934/biophy.2023023 |
The COVID-19 pandemic continues to pose significant challenges to global public health, necessitating the development of effective vaccination strategies to mitigate disease transmission. In Thailand, the COVID-19 epidemic has undergone multiple waves, prompting the implementation of various control measures, including vaccination campaigns. Understanding the dynamics of disease transmission and the impact of vaccination strategies is crucial for guiding public health interventions and optimizing epidemic control efforts. In this study, we developed a comprehensive mathematical model, termed SSvIH1CH2RD, to elucidate the dynamics of the COVID-19 epidemic in Thailand. The model incorporates key epidemiological parameters, vaccination rates, and disease progression stages to assess the effectiveness of different vaccination strategies in curbing disease transmission. Parameter estimation and model fitting were conducted using real-world data from COVID-19 patients in Thailand, enabling the simulation of epidemic scenarios and the exploration of optimal vaccination rates. Our results showed that optimizing vaccination strategies, particularly by administering approximately 119,625 doses per day, can significantly reduce the basic reproduction number (R0) below 1, thereby accelerating epidemic control. Simulation results demonstrated that the optimal vaccination rate led to a substantial decrease in the number of infections, with the epidemic projected to be completely eradicated from the population by June 19, 2022. These findings underscore the importance of targeted vaccination efforts and proactive public health interventions in mitigating the spread of COVID-19 and minimizing the burden on healthcare systems. Our study provides valuable insights into the optimization of vaccination strategies for epidemic control, offering guidance for policymakers and healthcare authorities in Thailand and beyond. By leveraging mathematical modeling techniques and real-world data, stakeholders can develop evidence-based strategies to combat the COVID-19 pandemic and safeguard public health.
Abbreviations: IDP: Intrinsically disordered protein; IDR: Intrinsically disordered region; YPRA: Yeast proteinase A; TFE: 2,2,2-trifluoroethanol; CD: Circular dichroism; SDSL: Site-directed spin labeling; CW-EPR: Continuous wave electron paramagnetic resonance; IAP: 3-(2-Iodoacetamido)-PROXYL; MTSL: (1-oxyl-2,2,5,5-tetramethyl-Δ3-pyrroline-3-methyl) methanethiosulfonate; MSL: 4-maleimido-TEMPO; DTT: Dithiothreitol
Proteins or protein regions of 50 or more residues, which do not have stable secondary or tertiary structure under physiological conditions, are characterized as intrinsically disordered proteins (IDPs) or intrinsically disordered regions (IDRs) [1,2,3,4,5,6,7]. Although the structure-function paradigm suggests that 3D structure is necessary for function, the discovered importance of intrinsic disorder within proteins or protein regions of higher eukaryotic systems shows how function could arise from various unstructured states [3,4,8,9]. IDPs have important roles in cellular signal transduction, translation, and transcription, influencing the study of disordered proteins, focusing on both the functions provided by IDPs, and the conformational changes associated with target binding [4,7,10].
IA3, an IDP found in Sacchoromyeces cerevisiae, is composed of 68 amino acid residues and acts as an inhibitor of yeast proteinase A (YPRA). Previous studies have shown IA3 to adopt an α-helical conformation when bound in the active site of YPRA and when exposed to the secondary structure stabilizer 2,2,2-trifluoroethanol (TFE). These observed helical tendencies of IA3 differ between the N-terminal residues (2-34) and C-terminal residues (35-68) [11,12,13,14,15]. X-ray crystallographic models of IA3 bound to YPRA indicate α-helical structure in the N-terminal region (residues 2-34). In contrast, the C-terminal region (residues 35-68) has unresolved electron density, suggesting disorder (Figure 1) [11]. Biophysical methods of circular dichroism (CD) spectroscopy, 2D H1 N15 nuclear magnetic resonance spectroscopy, molecular dynamics simulations, laser temperature-jump fluorescence spectroscopy, fluorescence resonance energy transfer spectroscopy, and site directed spin labeling (SDSL), provide data supporting a two-state transition for IA3, where the C-terminus undergoes a helical transition but is less pronounced than the N-terminus [13,14,15,16,17,18].
SDSL, in combination with continuous wave electron paramagnetic resonance (CW-EPR) spectroscopy is a biophysical tool for probing structure, dynamics, and conformational changes in macromolecules [19,20,21,22,23]. For proteins, typically a CW-EPR active chemical group, such as a nitroxide spin label, is incorporated at desired locations by chemical modification of a substituted cysteine. The introduction of an active chemical group provides a reactive side group for covalent bond spin labeling. Recent progress in spin labeling methodology has also demonstrated the ability to incorporate non-natural amino acids for subsequent chemical modification, as well as novel schemes for specifically labeling tyrosine residues; a potentially very useful strategy for spin labeling IDPs [24,25,26,27].
The resultant CW-EPR spectral line shape reflects motional averaging of the nitroxide spin label [20,21,28]. Typically, the effects are described by three main types of motion, occurring in the 0.1-50 ns time scale [29,30]. The corresponding correlation times are referred to as the following: τR, the overall tumbling of the protein, τI, the movement of the spin label about the bonds connecting it to the protein, and τB, the motion of the spin labeled protein backbone. SDSL CW-EPR has been applied to study structure-to-unstructured transitions in proteins as well as to IDP systems [5,17,22,31,32]. For structured proteins ≥ 15 kDa, line shape parameters such as the second moment and scaled mobility readily reveal conformational changes without the need to increase solution viscosity [29,33]. Ordered to disordered transitions in structured proteins are readily characterized by these “traditional” line shape parameter analyses [34,35]. However, in IDPs, the degree of disorder and the dynamics of the system often times results in “isotropic-like” spectra that fall into the fast motional averaging limit. In these cases, for relatively fast motional averaging, the ratio of the intensities of the low field (h(+1)) and central field (h(0)) transitions provide an alternative line shape analysis parameter, h(1)/h(0), that is useful to describe conformational changes [5,17,22,36].
Here, extending upon previous work, we performed a spin-labeled cysteine scanning profile of the N- and C-termini of IA3, where a series of fifteen singly spin-labeled IA3 constructs consisting of ten and five cysteine substitutions in the N- and C-terminus were generated, respectively (Figure 1) [17]. A comparison is provided between the termini on a site by site basis revealing variations in the degree of transition of the two termini, and sensitivity of the helical content of the N-terminus to the spin-labeled cysteine substitution and select amino acid substitutions. As these results show, the degree of site-specific helical propensity of the N-terminus is modulated by the introduction of a chemically modified cysteine residue. By further probing site V8 of IA3, the specific effects on secondary structure due to site-specific amino acid substitution and spin labeling are analyzed. The size of an introduced residue or chemical modification has marked effects on the helical propensity within the N-terminal region of IA3, whether found on the buried (concave) or solvent exposed (convex) side of the helix when bound to YPRA. By comparing SDSL and CD results, we are able to determine sites within IA3 that can tolerate spin-label incorporation and retain WT helical tendencies; thus, we label those our “wild-type mimics”.
BL21(DE3) pLysS cells were purchased from Invitrogen (Carlsbad, CA). E. coli codon-optimized DNA for the IA3 gene from S. cerevisiae, and DNA primers used for site-directed mutagenesis were purchased from DNA2.0 (Newark, CA). DNA taq polymerase and Dpn1 were purchased from New England Biolabs (Ipswich, MA). A 5-mL chelating column was purchased from GE Healthcare (Wauwatosa, WI). A HiPrep 26/10 desalting column was purchased from Amersham (Pittsburgh, PA). The 16.5% Tris-Tricene gels were purchased from Bio-Rad (Hercules, CA). 3-(2-Iodoacetamido)-PROXYL (IAP), (1-oxyl-2,2,5,5-tetramethyl-Δ3-pyrroline-3-methyl) methanethiosulfonate (MTSL), and 3-maleimido-TEMPO (MSL) spin label was purchased from Sigma Aldrich (St. Louis, MO). The 0.60 i.d. × 0.84 o.d. capillary tubes were purchased from Fiber Optic Center (New Bradford, MA). Unless otherwise stated, all other reagents and products were purchased from Fisher Scientific (Pittsburgh, PA) and used as received.
The optimized IA3 gene was cloned into the pET-22b+ vector containing a C-terminal 6× His-tag, adding the sequence LEHHHHHH to the C-terminus as described previously [17]. Site directed mutagenesis via polymerase chain reaction was used to introduce non-native cysteines by designing mutagenic primers. The substituted sequence in each resultant plasmid was confirmed via Sanger DNA sequencing. WT and cysteine substituted constructs were purified as described previously [17]. Following affinity chromatography, residual nickel from the column was removed by addition of EDTA to a concentration of 100 mM and cysteine reduction was ensured by addition of dithiothreitol (DTT) to a concentration of 0.1 mM. Protein purity was assessed using a 16.5% tris-tricene SDS-PAGE gel.
WT and cysteine substituted proteins were buffer exchanged into 50 mM sodium phosphate, 300 mM sodium chloride, pH 7.4 using a desalting column. This step also removed residual nickel, EDTA, and DTT that were added as described above. After desalting, IAP, MTSL, or MSL nitroxide spin label, dissolved in ethanol, was added in excess (5-10× molar ratio) and allowed to react with cysteine substituted constructs for 10-12 h at room temperature, in the dark. Spin labeling schemes are represented in Figure 2A. Excess spin label was removed by repeating the desalting buffer exchange process as described above.
For each spin labeled construct, nine separate samples were prepared with increasing 5% (v/v) increments of 2,2,2-trifluoroethanol (TFE) that spanned 0-40%. Each individual sample was prepared at a final volume of 500 µL: 300 µL of 50 mM sodium phosphate, 300 mM sodium chloride, pH 7.4 buffer, (200−x) µL of TFE, where x ranged from 0-200 µL, and x µL of appropriately prepared solution of 5× stock buffer diluted with water to keep all samples at equal ionic strength and pH. These larger 500 µL sample volumes were prepared to help mitigate pipetting errors of small volumes. Samples of ∼10 µL were drawn into 0.60 i.d. × 0.84 o.d. capillary tubes, which were then flame sealed on both ends ensuring minimal to no TFE evaporation.
For each spin labeled construct, samples were prepared with varying percentages of increasing TFE concentrations. For IAP labeled cysteine scanning and the V8C construct with varying spin label attachments, three samples were prepared with increasing 15% (v/v) increments of TFE from 0-30%. For amino acid scanning at the 8th position of the IA3 sequence and cysteine or alanine substitution at the 11th position, six or nine samples were prepared with increasing 7% (v/v) increments of TFE from 0-35% or increasing 5% (v/v) increments of TFE from 0-40%, respectively. Samples were prepared with 80-120 μL of protein sample in 50 mM sodium phosphate, 300 mM sodium chloride, pH 7.4 to reach a final concentration of 10 μM protein in 1000 μL. The final volume was adjusted using 10 mM sodium phosphate, pH 7.0 buffer or TFE to reach appropriate % (v/v) TFE. For the V8C un-labeled construct, 0.1 mM DTT was present to prevent disulfide bonding of free cysteine residues.
Either a Bruker ER200 spectrometer with an ER023 M signal channel, an ER032 M field control unit or a Bruker E500 with a loop gap resonator (Medical Advances, Milwaukee, WI) was used to collect CW X-band EPR spectra. For all experiments the temperature was kept at a constant 27 ± 0.2 °C using a nitrogen gas passed through a copper coil submersed into a water bath, while monitored with a temperature probe and thermometer from OMEGA Engineering Inc (Norwalk, CT). All spectra are reported as an average of four scans, collected as 100 G sweep width, 0.6 modulation amplitude, 70 ms conversion time, 100 kHz modulation amplitude, and 2 mW incident microwave power. Spectra were normalized using labVIEW software allowing for baseline correction and double integral area normalization. A representative nitroxide EPR spectrum with the typical three-line transition pattern is shown and labeled in Figure 2B. Spectral line shapes were analyzed by calculating the ratio of the intensities of the low field and center field transitions, h(+1)/h(0) as a function of TFE percentage. Although the low field and center field transitions individually are not as sensitive to motion as the high field transition (h(−1)), this ratio has been previously shown to change considerably as the unstructured-to-structured transition occurs and is not sensitive to errors in baseline correction or double integration [17,36]. We have used this analysis previously to monitor the unstructured to helical transition in IA3 by % TFE. Plots of h(+1)/h(0) values versus the % TFE are sigmoidal in shape and are well fit by a two-state Boltzmann transition, given by
y=(A1−A2)1+e(x−x0)/dx+A2 | (1) |
A1 is the initial value of the curve, A2 is the final value of the curve and x0 is the midpoint of the curve—which corresponds well to the midpoint of folding transitions observed by CD and NMR [14]. Data were fit using Origin 8.5 software.
CD measurements were collected using an AVIV model 202 CD spectrometer set at 27 °C. To collect the spectra, 400 μL of each sample were loaded into a 1mm path length quartz cuvette, cleaned between runs with nanopure H2O and ethanol. Typically, four scans of each buffer blank and sample were collected between wavelengths of 200-250 nm with 1 nm increments and a four second averaging time. Each set of CD scans for each construct was averaged, and the averaged buffer baseline was subtracted from each individual construct averaged scan. Units were converted to Δε (molar circular dichroism, M−1 cm−1) from the instruments output unit of Θ (ellipticity, mdeg) using
Δε (M−1cm−1)=Θ (mdeg)xCxL100xM | (2) |
Θ is the output ellipticity value from the CD measurement, C is the concentration of sample in g/L, L is the path length in cm, and M is the average molecular weight in g/mol.
It is well known that the labeling of a naturally occurring or substituted cysteine may impact the native function of a protein. To understand if spin-label and cysteine substitution within IA3 impacted folding or folding propensity, cysteine scanning was performed across the N- and C-termini of IA3. Sites K7, V8, S9, E10, I11, F12, Q13, S14, S15, and S27 were chosen within the N-terminus and span two helical turns containing both solvent exposed and buried residues in the α-helix when bound to YPRA. Sites Y57, N58, K59, L60, and K61, were selected in the C-terminus (Figure 1).
CW X-band EPR spectra were collected for the fifteen IAP spin labeled-IA3 variants (referred to within as P1-IA3) as a function of TFE concentration that ranged from 0-40% (v/v) TFE. Figure 3 shows representative CW-EPR spectra for 0%, 15%, and 30% (v/v) TFE for each P1-IA3 variant. Spectra are plotted with normalized double integral area to allow for easy visualization of changes in mobility. As mobility (defined as both the rate of motion and order of motion) decreases, the line shape broadens. As a result, when plotted with normalized area, the intensity appears to be less. A decrease in mobility upon increasing TFE percentage is observed for nearly all sites and indicates a conformational change upon increasing TFE concentration (earlier work has ruled out the impact of increased solution viscosity) [17]. Inspection of the line shapes shows that the change in mobility differs among these 15 P1-IA3 variants, indicating that some spin-labeled cysteine substitutions affect the conformation of each state, altering the TFE-induced helical propensity of IA3.
A more quantitative way to analyze the TFE-induced conformational changes for each variant is to plot h(+1)/h(0) values as a function of TFE % (Figure 4). For isotropic motion, the h(+1)/h(0) value is expected to be 1. As the mobility decreases as the protein becomes helical, lower values of h(+1)/h(0) are expected. In the unfolded state (absence of TFE), all P1-IA3 variants have similar CW-EPR spectra indicative of similar high nitroxide mobility (indicated by values of h(+1)/h(0) ∼ 0.9). As a function of TFE concentration, the shape of each h(+1)/h(0) plot can be well fit by
Several conclusions can be drawn from the data in Figure 4. Firstly, the sigmoidal behavior of the h(+1)/h(0) parameter for N-terminal P1-IA3 variants (gray shade, solid diamonds) varies both in sharpness and extent of transition, which is defined by the limiting h(+1)/h(0) value. In contrast C-terminal P1-IA3 variants (no shade, open circles) exhibit similar behavior, both in the sharpness and extent of transition. Secondly, there are site specific variations in the behavior of the h(+1)/h(0) parameter in the N-terminus that follow a helical trend that maps to the surface of IA3 in the YRPA bound site. Specifically, sites V8P1, I11P1, and S15P1 which are located on the concave buried face of the peptide when bound to YPRA, have small changes (<0.05 over the 0-40% (v/v) TFE range) in their h(+1)/h(0) values. In contrast, sites S9P1, Q13P1 and S27P1, which reside on the convex solvent exposed face when bound to YPRA, pose large changes in their h(+1)/h(0) values (>0.3 over the 0-40% TFE range). Figure 5 shows the location of these residues in the IA3: YPRA complex crystal structure. These results indicate a sensitivity of the N-terminal helicity to amino acid substitution.
To further explore these observations, CD spectroscopy was performed on each of the N-terminal P1-IA3 variants with 0%, 15%, and 30% (v/v) TFE. Data are shown in Figure 6, which also overlays results from WT IA3 (gray solid line) in 30% (v/v) TFE for comparison. For wild type IA3, CD data reveal a predominantly random coil conformation for 0% (v/v) TFE, whereas with 30% (v/v) TFE the spectra reflect a mostly helical conformation with features at 208 nm and 222 nm indicative of helix formation [37,38]. For all N-terminal P1-IA3 variants, CD data confirm that in the absence of TFE, all constructs adopt a random coil conformation and that the addition of TFE induces a helical conformation. As was observed with the CW-EPR data, the degree of the helical transition, as reflected in the values of Δε222nm in CD spectra, varies when compared to WT for the P1 variants. It can be seen that I11P1 and S14P1 possess CD spectra that very closely match that of WT (i.e., represent SL wild-type mimics). V8P1 is found to be less helical than WT, whereas all other variants have greater helicity than WT. To directly compare the SDSL CW-EPR characterized unstructured-to-structured transition within the N-terminus to the degree of α-helicity from CD spectroscopy, we plotted the values of Δh(+1)/h(0), or degree of transition (h(+1)/h(0) (40% (v/v) TFE) − h(+1)/h(0) (0% (v/v) TFE)) against the Δε222nm (30% (v/v) TFE) (
Because the P1 substitution at site V8 revealed less helical character than WT, we proceeded to generate a series of substitutions at this site to probe how amino acid size and type impacted helical propensity. The following constructs were chosen: V8A and V8I (smaller and larger than V; respectively), V8P (typical helical disrupter), V8D (charged), and V8C (our cysteine mutation without spin-label). CD spectroscopy was performed on each of the V8 variants with 0-40% or 0-35% (v/v) TFE (
The result that V8C contained WT helicity was surprising, which prompted us to examine if and how spin-label structure influenced helicity when incorporated at this site. CD spectroscopy was performed on the V8C spin-labeled variants from 0-40% or 0-30% (v/v) TFE (
Nevertheless, results appear to indicate a dependency of structural bulkiness for disruption or enhancement of helical propensity. Indeed, the structure of IA3 bound to YPRA shows a slight concave bend on the buried interface, making substitution to larger residues disrupt helical behavior and conversely, a smaller residue would enhance helicity. To test this idea, we made two additional constructs, V8P1-I11A and V8A-I11P1, which would place the smaller alanine residue next to the P1 spin-label. CD spectroscopy was performed on each additional construct from 0-40% or 0-35% (v/v) TFE (
It was determined through SDSL CW-EPR that the extent of helical order within the N-terminal region of IA3 is easily perturbed by amino acid substitution and IAP labeling. The C-terminal region did not have similar susceptibility, where changes in labeled amino acid substitution did not reveal drastic changes in the unstructured-to-structured transition, and as we showed previously, retained WT behavior [17]. CD results confirmed that the helical order within the N-terminus varied due to the location of the altered residue in relation to the buried or solvent-exposed, or the concave or convex face of the protein, respectively. The concave face of IA3, is negatively affected by cysteine mutation and labeling, whereas the convex surface increases in helical order. It is proposed that the change in helical propensity on either side of the protein is due to spatial availability when incorporating a cysteine residue and nitroxide probe. Residue V8, located on the concave side of the N-terminal region showed the greatest negative perturbation due to cysteine substitution, and spin labeling. Mutations to the V8 residue, including substitutions in residue, charge, and size did not drastically affect the helical propensity, but rather the addition of the IAP nitroxide probe had the greatest effect on helical formation.
These experiments demonstrate the care that must be taken when spin-labeling IDPs. Proper checks-and-balances are needed to ensure that a labeled site retains WT helical propensity (in addition to any other functional assays). Overall, the results identify key sites for spin-label incorporation that mimic WT helical tendencies that can be used in future studies of IA3 conformational sampling studies.
This work was supported by NSF MCB 1329467 & MCB 1715384 and NIH 1S10-RR031603-01 to GEF. The authors would like to thank Stephen J. Hagen and Robert McKenna for access to the Aviv CD spectrometer, and Christian Altenbach for Labview software to analyze the EPR data.
All authors declare no conflict of interest in this paper.
[1] |
M. Worobey, J. I. Levy, L. M. Serrano, A. Crits-Christoph, J. E. Pekar, S. A. Goldstein, et al., The Huanan Seafood Wholesale Market in Wuhan was the early epicenter of the COVID-19 pandemic, Science, 377 (2022), 951–959. https://doi.org/10.1126/science.abp8715 doi: 10.1126/science.abp8715
![]() |
[2] | World Health Organization, WHO Director-General's statement on IHR Emergency Committee on Novel Coronavirus (2019-nCoV), 2020. Available from: https://www.who.int/director-general/speeches/detail/who-director-general-s-statement-on-ihr-emergency-committee-on-novel-coronavirus-(2019-ncov). |
[3] |
P. Alam, I. Y. H. Chu, H. J. Larson, L. Lin, Social consequences of mass quarantine during epidemics: A systematic review with implications for the COVID-19 response, J. Travel Med., 27 (2020), taaa192. https://doi.org/10.1093/jtm/taaa192 doi: 10.1093/jtm/taaa192
![]() |
[4] |
N. G. Davies, A. J. Kucharski, R. M. Eggo, A. Gimma, W. J. Edmunds, T. Jombart, et al., Effects of non-pharmaceutical interventions on COVID-19 cases, deaths, and demand for hospital services in the UK: A modelling study, Lancet Public Health, 5 (2020), e375–e385. https://doi.org/10.1016/S2468-2667(20)30133-X doi: 10.1016/S2468-2667(20)30133-X
![]() |
[5] |
H. Lau, V. Khosrawipour, P. Kocbach, A. Mikolajczyk, J. Schubert, J. Bania, et al., The positive impact of lockdown in Wuhan on containing the COVID-19 outbreak in China, J. Travel Med., 27 (2020), taaa037. https://doi.org/10.1093/jtm/taaa037 doi: 10.1093/jtm/taaa037
![]() |
[6] | World Health Organization, World Health Organization Coronavirus (COVID-19) Dashboard, 2022. |
[7] | SAT-MOPH, DDC COVID-19 Interactive Dashboard 2022, 2022. Available from: https://ddc.moph.go.th/covid19-dashboard/?dashboard=main (accessed on 16 May 2023). |
[8] |
P. Jarumaneeroj, P. O. Dusadeerungsikul, T. Chotivanich, T. Nopsopon, K. Pongpirul, An epidemiology-based model for the operational allocation of COVID-19 vaccines: A case study of Thailand, Comput. Ind. Eng., 167 (2022), 108031. https://doi.org/10.1016/j.cie.2022.108031 doi: 10.1016/j.cie.2022.108031
![]() |
[9] |
N. Nittayasoot, R. Suphanchaimat, P. Thammawijaya, C. Jiraphongsa, T. Siraprapasiri, K. Ploddi, et al., Real-world effectiveness of COVID-19 vaccines against severe outcomes during the period of omicron predominance in Thailand: a test-negative nationwide case–control study, Vaccines, 10 (2022), 2123. https://doi.org/10.3390/vaccines10122123 doi: 10.3390/vaccines10122123
![]() |
[10] | World Health Organization Country Office for Thailand, Situation of coronavirus disease 2019 in Thailand: 7 February 2024, 2024. Available from: https://www.who.int/thailand/news/detail/07-02-2024-update-on-covid-19-in-thailand--7-february-2024-THA (accessed on 4 March 2022). |
[11] |
P. Ngamchaliew, N. Kaewkuea, N. Nonthasorn, Changes in preventive behaviour after COVID-19 vaccination in Thailand: A cross-sectional study, BMC Public Health, 22 (2022), 2039. https://doi.org/10.1186/s12889-022-14494-x doi: 10.1186/s12889-022-14494-x
![]() |
[12] | World Health Organization Country Office for Thailand, Progress and Challenges in COVID-19 Vaccine Roll-Out, 2021. Available from: https://cdn.who.int/media/docs/default-source/searo/ivd/itag-2021/day2/05-thailand---progress-and-challenges-in-covid-19-vaccine-roll-out.pdf?sfvrsn=2412ff89_5 (accessed on 4 March 2023). |
[13] | Reuters, Thailand seeks 12 mln Sinovac shots for mix-and-match vax strategy, 2021. Available from: https://www.reuters.com/world/asia-pacific/thailand-reports-daily-record-239-new-coronavirus-deaths-2021-08-17/ (accessed on 4 March 2024). |
[14] | Department of disease control, Guidelines for people who need COVID-19 Epidemic Situation in 2021 in Thailand, 2021. Available from: https://ddc.moph.go.th/vaccine-covid19/getFiles/11/1628849610213.pdf (accessed on 15 May 2023). |
[15] | P. Thepgumpanat, P. Wongcha-um, In First, Thailand to Mix Sinovac, AstraZeneca Vaccine Doses, Reuters, 2021. Available from: https://www.reuters.com/world/asia-pacific/thailand-starts-tighter-coronavirus-lockdown-around-capital-2021-07-12/ (accessed on 4 March 2024). |
[16] | Hfocus, Update! Vaccination with 3 doses of vaccine "Children 12–17 years" and 4 doses stimulate immunity, 2021. Available from: https://www.hfocus.org/content/2022/03/24758 (accessed on 4 March 2024). |
[17] | Thaipbs, The Recommended Third Dose and Fourth Dose of COVID-19 Vaccine, 2021. Available from: https://www.thaipbs.or.th/news/content/311590 (accessed on 4 March 2024). |
[18] | Hfocus, Investigative doctor reveals the facts! Effectiveness of each brand of COVID vaccine, 2021. Available from: https://www.hfocus.org/content/2021/05/21582 (accessed on 4 March 2024). |
[19] | C. Zimmer, J. Corum, S. L. Wee, M. Kristoffersen, Coronavirus Vaccine Tracker, 2022. Available from: https://www.nytimes.com/interactive/2020/science/coronavirus-vaccine-tracker.html (accessed on 4 March 2024). |
[20] | K. Katella, Comparing the COVID-19 Vaccines: How Are They Different, YaleMedicine, 2023. Available from: https://www.yalemedicine.org/news/covid-19-vaccine-comparison (accessed on 4 March 2024). |
[21] | Centers for Disease Control and Prevention, Understanding How COVID-19 Vaccines Work, 2023. Available from: https://archive.cdc.gov/www_cdc_gov/coronavirus/2019-ncov/vaccines/different-vaccines/how-they-work.html (accessed on 4 March 2024). |
[22] | World Health Organization, What is COVID-19 vaccine efficacy, 2021. Available from: https://www.afro.who.int/news/what-covid-19-vaccine-efficacy (accessed on 4 March 2024). |
[23] | F. Brauer, C. Castillo-Chavez, Z. Feng, Endemic disease models, in Mathematical Models in Epidemiology, Springer, (2019), 63–116. https://doi.org/10.1007/978-1-4939-9828-9_3 |
[24] |
A. M. Salman, I. Ahmed, M. H. Mohd, M. S. Jamiluddin, M. A. Dheyab, Scenario analysis of COVID-19 transmission dynamics in Malaysia with the possibility of reinfection and limited medical resources scenarios, Comput. Biol. Med., 133 (2021), 104372. https://doi.org/10.1016/j.compbiomed.2021.104372 doi: 10.1016/j.compbiomed.2021.104372
![]() |
[25] |
K. Rajagopal, N. Hasanzadeh, F. Parastesh, I. I. Hamarash, S. Jafari, I. Hussain, A fractional-order model for the novel coronavirus (COVID-19) outbreak, Nonlinear Dyn., 101 (2020), 711–718. https://doi.org/10.1007/s11071-020-05757-6 doi: 10.1007/s11071-020-05757-6
![]() |
[26] |
Q. M. Al-Mdallal, Mathematical modeling and simulation of SEIR model for COVID-19 outbreak: A case study of Trivandrum, Front. Appl. Math. Stat., 9 (2023), 1124897. https://doi.org/10.3389/fams.2023.1124897 doi: 10.3389/fams.2023.1124897
![]() |
[27] |
S. Nana-Kyere, F. A. Boateng, P. Jonathan, A. Donkor, G. K. Hoggar, B. D. Titus, et al., Global analysis and optimal control model of COVID-19, Comput. Math. Methods Med., 2022 (2022), 9491847. https://doi.org/10.1155/2022/9491847 doi: 10.1155/2022/9491847
![]() |
[28] |
S. Boulaaras, R. Ramalingam, A. J. Gnanaprakasam, SEIR model for COVID-19: Stability of the standard coronavirus factor and control mechanism, Eur. Phys. J. Spec. Top., 232 (2023), 2485–2495. https://doi.org/10.1140/epjs/s11734-023-00915-4 doi: 10.1140/epjs/s11734-023-00915-4
![]() |
[29] |
V. E. Papageorgiou, G. Vasiliadis, G. Tsaklidis, Analyzing the asymptotic behavior of an extended SEIR model with vaccination for COVID-19, Mathematics, 12 (2023), 55. https://doi.org/10.3390/math12010055 doi: 10.3390/math12010055
![]() |
[30] |
I. Rahimi, F. Chen, A. H. Gandomi, A review on COVID-19 forecasting models, Neural Comput. Appl., 35 (2023), 23671–23681. https://doi.org/10.1007/s00521-020-05626-8 doi: 10.1007/s00521-020-05626-8
![]() |
[31] |
S. Jitsinchayakul, U. W. Humphries, A. Khan, The SQEIRP mathematical model for the COVID-19 epidemic in Thailand, Axioms, 12 (2023), 75. https://doi.org/10.3390/axioms12010075 doi: 10.3390/axioms12010075
![]() |
[32] |
K. Intawong, S. Chariyalertsak, K. Chalom, T. Wonghirundecha, W. Kowatcharakul, A. Thongprachum, et al., Effectiveness of heterologous third and fourth dose COVID-19 vaccine schedules for SARS-CoV-2 infection during delta and omicron predominance in Thailand: A test-negative, case-control study, Lancet Reg. Health - Southeast Asia, 10 (2023), 100121. https://doi.org/10.1016/j.lansea.2022.100121 doi: 10.1016/j.lansea.2022.100121
![]() |
[33] |
N. E. Aikawa, L. V. Kupa, C. A. Silva, C. G. Saad, S. G. Pasoto, E. F. Yuki, et al., Strong response after fourth dose of mRNA COVID-19 vaccine in autoimmune rheumatic diseases patients with poor response to inactivated vaccine, Rheumatology, 62 (2023), 480–485. https://doi.org/10.1093/rheumatology/keac301 doi: 10.1093/rheumatology/keac301
![]() |
[34] |
S. Y. Tartof, J. M. Slezak, L. Puzniak, V. Hong, T. B. Frankland, F. Xie, et al., Analysis of mRNA COVID-19 vaccine uptake among immunocompromised individuals in a large US health system, JAMA Netw. Open, 6 (2023), e2251833–e2251833. https://doi.org/10.1001/jamanetworkopen.2022.51833 doi: 10.1001/jamanetworkopen.2022.51833
![]() |
[35] |
K. Safa, C. N. Kotton, COVID-19 vaccines and solid organ transplantation: More doses, more protection, Transplantation, 107 (2023), 21–22. https://doi.org/10.1097/TP.0000000000004387 doi: 10.1097/TP.0000000000004387
![]() |
[36] |
K. Bardosh, A. Krug, E. Jamrozik, T. Lemmens, S. Keshavjee, V. Prasad, et al., COVID-19 vaccine boosters for young adults: A risk benefit assessment and ethical analysis of mandate policies at universities, J. Med. Ethics, 50 (2024), 126–138. https://doi.org/10.1136/jme-2022-108852 doi: 10.1136/jme-2022-108852
![]() |
[37] |
S. Ramot, O. Tal, Attitudes of healthcare workers in Israel towards the fourth dose of COVID-19 vaccine, Vaccines, 11 (2023), 385. https://doi.org/10.3390/vaccines11020385 doi: 10.3390/vaccines11020385
![]() |
[38] |
M. H. DarAssi, T. A. Shatnawi, M. A. Safi, Mathematical analysis of a MERS-Cov coronavirus model, Demonstr. Math., 55 (2022), 265–276. https://doi.org/10.1515/dema-2022-0022 doi: 10.1515/dema-2022-0022
![]() |
[39] |
T. Theparod, P. Kreabkhontho, W. Teparos, Booster dose vaccination and dynamics of COVID-19 pandemic in the fifth wave: An efficient and simple mathematical model for disease progression, Vaccines, 11 (2023), 589. https://doi.org/10.3390/vaccines11030589 doi: 10.3390/vaccines11030589
![]() |
[40] |
T. Patalon, Y. Saciuk, A. Peretz, G. Perez, Y. Lurie, Y. Maor, et al., Waning effectiveness of the third dose of the BNT162b2 mRNA COVID-19 vaccine, Nat. Commun., 13 (2022), 3203. https://doi.org/10.1038/s41467-022-30884-6 doi: 10.1038/s41467-022-30884-6
![]() |
[41] | Centers for Disease Control and Prevention, COVID-19 Vaccine Effectiveness Update, 2022. Available from: https://covid.cdc.gov/covid-data-tracker/#vaccine-effectiveness (accessed on 23 July 2023). |
[42] | World Health Organization Country Office for Thailand, COVID-19 Situation, Thailand 27 July 2022, 2022. Available from: https://cdn.who.int/media/docs/default-source/searo/thailand/2022_07_27_tha-sitrep-243-covid-19.pdf?sfvrsn=23bab9e2_1 (accessed on 23 July 2023). |
[43] | BBC NEWS Thailand, Sinovac-Sinopharm: COVID-19 evaluating the two Chinese vaccinations that Thailand uses and their respective efficacies, 2021. Available from: https://www.bbc.com/thai/thailand-57525497 (accessed on 15 March 2024). |
[44] | Ministry of Higher Education Science Research and Innovation, 100 million doses, COVID-19 vaccination statistics in Thailand, 2021. Available from: https://www.mhesi.go.th/index.php/all-media/infographic/5379-6412221Covid_19.html (accessed on 14 March 2024). |
[45] |
K. Dietz, The estimation of the basic reproduction number for infectious diseases, Stat. Methods Med. Res., 2 (1993), 23–41. https://doi.org/10.1177/096228029300200103 doi: 10.1177/096228029300200103
![]() |
[46] | The Bureau of Registration Administration, Official population statistics from the civil registration (monthly), 2022. Available from: https://stat.bora.dopa.go.th/stat/statnew/statMONTH/statmonth/#/displayData (accessed on 15 May 2023). |
[47] | Department of Health Ministry of Public Health, Report on progress of COVID-19 vaccination services, 2022. Available from: https://public.tableau.com/views/SATCOVIDDashboard/1-dash-tiles?:showVizHome=no (accessed on 14 March 2024). |
[48] | UK Health Security Agency, SARS-CoV-2 variants of concern and variants under investigation in England, 2022. Available from: https://assets.publishing.service.gov.uk/government/uploads/system/uploads/attachment_data/file/1050236/technical-briefing-34-14-january-2022.pdf (accessed on 19 May 2023). |
[49] |
Z. Zi, Sensitivity analysis approaches applied to systems biology models, IET Syst. Biol., 5 (2011), 336–346. https://doi.org/10.1049/iet-syb.2011.0015 doi: 10.1049/iet-syb.2011.0015
![]() |
1. | Mohit Jaiswal, Trang T. Tran, Qingjiang Li, Xin Yan, Mingwei Zhou, Krishnendu Kundu, Gail E. Fanucci, Zhongwu Guo, A metabolically engineered spin-labeling approach for studying glycans on cells, 2020, 11, 2041-6520, 12522, 10.1039/D0SC03874A | |
2. | Katie M. Dunleavy, Collin Oi, Tianyan Li, Andrew Secunda, Afnan M. Jaufer, Yinlu Zhu, Lee Friedman, Alexander Kim, Gail E. Fanucci, Hydrogen Bonding Compensation on the Convex Solvent-Exposed Helical Face of IA3, an Intrinsically Disordered Protein, 2023, 0006-2960, 10.1021/acs.biochem.3c00169 | |
3. | Katie M. Dunleavy, Tianyan Li, Eugene Milshteyn, Afnan M. Jaufer, Shamon A. Walker, Gail E. Fanucci, Charge Distribution Patterns of IA3 Impact Conformational Expansion and Hydration Diffusivity of the Disordered Ensemble, 2023, 1520-6106, 10.1021/acs.jpcb.3c06170 |