Homology Models and Molecular Dynamics Simulations of Main Proteinase from Coronavirus Associated with Severe Acute Respiratory Syndrome (SARS)

In this study, two structural models (denoted as M^proST and M^proSH) of the main proteinase (M^pro) from the novel coronavirus associated with severe acute respiratory syndrome (SARS‐CoV) were constructed based on the crystallographic structures of M^pro from transmissible gastroenteritis coronavirus (TGEV) (M^proT) and human coronavirus HcoV‐229E (M^proH), respectively. Various 200 ps molecular dynamics simulations were subsequently performed to investigate the dynamics behaviors of several structural features. Both M^proST and M^proSH exhibit similar folds as their respective template proteins. These structural models reveal three distinct functional domains as well as an intervening loop connecting domains II and III as found in both template proteins. In addition, domain III of these structures exhibits the least secondary structural conservation. A catalytic cleft containing the substrate binding subsites S1 and the S2 between domains I and II are also observed in these structural models. Although these structures share many common features, the most significant difference occurs at the S2 subsite, where the amino acid residues lining up this subsite are least conserved. It may be a critical challenge for designing anti‐SARS drugs by simply screening the known database of proteinase inhibitors.

Molecular dynamics simulations of human cystatin C and its L68Q varient to investigate the domain swapping mechanism

Human cystatin C variant (L68Q), one of the amyloidgenic proteins, has been shown to form dimeric structure spontaneously via domain swapping and easily cause amyloid deposits in the brains of patients suffering from Alzheimer's disease or hereditary cystatin C amyloid angiopathy. The monomeric L68Q and wild-type (wt) HCCs share similar structural feature consisting of a core with a five-stranded anti-parallel beta-sheet (beta-region) wrapped around a central helix. In this study, various molecular dynamics simulations were conducted to investigate the conformational fluctuations of the monomeric L68Q and wt HCCs at various combinations of temperature (300 and 500K) and pH (2 and 7) to gain insights into the domain swapping mechanism. The results show that elevated temperature accelerates the disruption of the hydrophobic core and acidic condition promotes the destruction of three salt bridges between beta2 and beta3 in both HCCs. The results also indicate that the interior hydrophobic core of the L68Q variant is relatively unstable, leading to domain swapping more readily comparing to wt HCC under conditions favoring this process. However, these two monomeric HCCs adopt the same mechanism of domain swapping as follows: (i) first, the interior hydrophobic core is disrupted; (ii) subsequently, the central helix departs from the beta-region; (iii) then, the beta2-L1-beta3 hairpin structure unfolds following the so-called "zip-up" mechanism; and (iv) finally, the open form HCC is generated.

Molecular dynamics simulations to investigate the domain swapping mechanism of human cystatin C

Human cystatin C (HCC), one of the amyloidgenic proteins, has been proved to form a dimeric structure via a domain swapping process and then cause amyloid deposits in the brains of patients suffering from Alzheimer's disease. HCC monomer consists of a core with a five-stranded antiparallel beta-sheet (beta region) wrapped around a central helix. The connectivity of these secondary structures is: (N)-beta1-alpha-beta2-L1-beta3-AS-beta4-L2-beta5-(C). In this study, various molecular dynamics simulations were conducted to investigate the conformational changes of the monomeric HCC at different temperatures (300 and 500 K) and pH levels (2, 4, and 7) to gain insight into the domain swapping mechanism. The results show that high temperature (500 K) and low pH (pH 2) will trigger the domain swapping process of HCC. We further proposed that the domain swapping mechanism of HCC follows four steps: (1) the alpha-helix moves away from the beta region; (2) the contacts between beta2 and beta3-AS disappear; (3) the beta2-L1-beta3 hairpin unfolds following the so-called "zip-up" mechanism; and finally (4) the HCC dimer is formed. Our study shows that high temperature can accelerate the unfolding of HCC and the departure of the alpha-helix from the beta-region, especially at low pH value. This is attributed to the fact that that low pH results in the protonation of the side chains of Asp, Glu, and His residues, which further disrupts the following four salt-bridge interactions stabilizing the alpha-beta interface of the native structure: Asp15-Arg53 (beta1-beta2), Glu21/20-Lys54 (helix-beta2), Asp40-Arg70 (helix-AS), and His43-Asp81 (beta2-AS).

Molecular dynamics simulations of various coronavirus main proteinases

In this study, two homology models (denoted as MproST and MproSH) of main proteinase (Mpro) from the novel coronavirus associated with severe acute respiratory syndrome (SARS-CoV) were constructed based on the crystal structures of Mpro from transmissible gastroenteritis coronavirus (TGEV) (MproT) and human coronavirus HcoV-229E (MproH), respectively. Both MproST and MproSH exhibit similar folds as their respective template proteins. These homology models reveal three distinct functional domains as well as an intervening loop connecting domains II and III as found in both template proteins. A catalytic cleft containing the substrate binding sites S1 and S2 between domains I and II are also observed. S2 undergoes more significant structural fluctuation than S1 during the 400 ps molecular dynamics simulations because it is located at the open mouth of the catalytic cleft, while S1 is situated in the very bottom of this cleft. The thermal unfolding of these proteins begins at domain III, where the structure is least conserved among these proteins. Mpro may still maintain its proteolytic activity while it is partially unfolded. The electrostatic interaction between Arg40 and Asp186 plays an important role in maintaining the structural integrity of both S1 and S2.