Distinctive RNA expression profiles in blood associated with Alzheimer disease after accounting for white matter hyperintensities.

BACKGROUND: Defining the RNA transcriptome in Alzheimer Disease (AD) will help understand the disease mechanisms and provide biomarkers. Though the AD blood transcriptome has been studied, effects of white matter hyperintensities (WMH) were not considered. This study investigated the AD blood transcriptome and accounted for WMH. METHODS: RNA from whole blood was processed on whole-genome microarrays. RESULTS: A total of 293 probe sets were differentially expressed in AD versus controls, 5 of which were significant for WMH status. The 288 AD-specific probe sets classified subjects with 87.5% sensitivity and 90.5% specificity. They represented 188 genes of which 29 have been reported in prior AD blood and 89 in AD brain studies. Regulated blood genes included MMP9, MME (Neprilysin), TGFß1, CA4, OCLN, ATM, TGM3, IGFR2, NOV, RNF213, BMX, LRRN1, CAMK2G, INSR, CTSD, SORCS1, SORL1, and TANC2. CONCLUSIONS: RNA expression is altered in AD blood irrespective of WMH status. Some genes are shared with AD brain.

Quantitative and Longitudinal Assessment of Systemic Innate Immunity in Health and Disease Using a 2D Gene Model.

Dysregulation of innate immunity is deeply involved in infectious and autoimmune diseases. For a better understanding of pathogenesis and improved management of these diseases, it is of vital importance to implement convenient monitoring of systemic innate immunity. Built upon our previous works on the host transcriptional response to infection in peripheral blood, we proposed a 2D gene model for the simultaneous assessment of two major components of systemic innate immunity, including VirSig as the signature of the host response to viral infection and BacSig as the signature of the host response to bacterial infection. The revelation of dysregulation in innate immunity by this 2D gene model was demonstrated with a wide variety of transcriptome datasets. In acute infection, distinctive patterns of VirSig and BacSig activation were observed in viral and bacterial infection. In comparison, both signatures were restricted to a defined range in the vast majority of healthy adults, regardless of age. In addition, BacSig showed significant elevation during pregnancy and an upward trend during development. In tuberculosis (TB), elevation of BacSig and VirSig was observed in a significant portion of active TB patients, and abnormal BacSig was also associated with a longer treatment course. In cystic fibrosis (CF), abnormal BacSig was observed in a subset of patients, and no overall change in BacSig abnormality was observed after the drug treatment. In systemic sclerosis-associated interstitial lung disease (SSc-ILD), significant elevation of VirSig and BacSig was observed in some patients, and treatment with a drug led to the further deviation of BacSig from the control level. In systemic lupus erythematosus (SLE), positivity for the anti-Ro autoantibody was associated with significant elevation of VirSig in SLE patients, and the additive effect of VirSig/BacSig activation was also observed in SLE patients during pregnancy. Overall, these data demonstrated that the 2D gene model can be used to assess systemic innate immunity in health and disease, with the potential clinical applications including patient stratification, prescription of antibiotics, understanding of pathogenesis, and longitudinal monitoring of treatment response.

Hypoxia and Activation of Neutrophil Degranulation-Related Genes in the Peripheral Blood of COVID-19 Patients.

Severe COVID-19 is characterized by systematic hyper-inflammation and subsequent damage to various organs. Therefore, it is critical to trace this cascade of hyper-inflammation. Blood transcriptome has been routinely utilized in the interrogation of host immune response in COVID-19 and other infectious conditions. In this study, consensus gene dysregulation in the blood was obtained from 13 independent transcriptome studies on COVID-19. Among the up-regulated genes, the most prominent functional categories were neutrophil degranulation and cell cycle, which is clearly different from the classical activation of interferon signaling pathway in seasonal flu. As for the potential upstream causal factors of the atypical gene dysregulation, systemic hypoxia was further examined because it is much more widely reported in COVID-19 than that in seasonal flu. It was found that both physiological and pathological hypoxia can induce activation of neutrophil degranulation-related genes in the blood. Furthermore, COVID-19 patients with different requirement for oxygen intervention showed distinctive levels of gene expression related to neutrophil degranulation in the whole blood, which was validated in isolated neutrophils. Thus, activation of neutrophil degranulation-related genes in the blood of COVID-19 could be partially attributed to hypoxia. Interestingly, similar pattern was also observed in H1N1 infection (the cause of Spanish flu) and several other severe respiratory viral infections. As for the molecular mechanism, both HIF-dependent and HIF-independent pathways have been examined. Since the activation of neutrophil degranulation-related genes is highly correlated with disease severity in COVID-19, early detection of hypoxia and active intervention may prevent further activation of neutrophil degranulation-related genes and other harmful downstream hyper-inflammation. This common mechanism is applicable to current and future pandemic as well as the severe form of common respiratory infection.

A two-gene marker for the two-tiered innate immune response in COVID-19 patients.

For coronavirus disease 2019 (COVID-19), a pandemic disease characterized by strong immune dysregulation in severe patients, convenient and efficient monitoring of the host immune response is critical. Human hosts respond to viral and bacterial infections in different ways, the former is characterized by the activation of interferon stimulated genes (ISGs) such as IFI27, while the latter is characterized by the activation of anti-bacterial associated genes (ABGs) such as S100A12. This two-tiered innate immune response has not been examined in COVID-19. In this study, the activation patterns of this two-tiered innate immune response represented by IFI27 and S100A12 were explored based on 1421 samples from 17 transcriptome datasets derived from the blood of COVID-19 patients and relevant controls. It was found that IFI27 activation occurred in most of the symptomatic patients and displayed no correlation with disease severity, while S100A12 activation was more restricted to patients under severe and critical conditions with a stepwise activation pattern. In addition, most of the S100A12 activation was accompanied by IFI27 activation. Furthermore, the activation of IFI27 was most pronounced within the first week of symptom onset, but generally waned after 2-3 weeks. On the other hand, the activation of S100A12 displayed no apparent correlation with disease duration and could last for several months in certain patients. These features of the two-tiered innate immune response can further our understanding on the disease mechanism of COVID-19 and may have implications to the clinical triage. Development of a convenient two-gene protocol for the routine serial monitoring of this two-tiered immune response will be a valuable addition to the existing laboratory tests.

The protein folding network indicates that the ultrafast folding mutant of villin headpiece subdomain has a deeper folding funnel.

Protein folding is a dynamic process with continuous transitions among different conformations. In this work, the dynamics in the protein folding network of villin headpiece subdomain (HP35) has been investigated based on multiple reversible folding trajectories of HP35 and its ultrafast folding mutant where sub-angstrom folding was achieved. The four folding states were clearly separated on the network, validating the classification of the states. Examination of the eight conformers with different formation of the individual helices revealed high plasticity of the three helices in all the four states. A consistent feature between the wild type and mutant protein is the dominant conformer 111 (all three helices formed) in the folded state and conformers 111 and 011 (helices II and III formed) in the major intermediate state, indicating the critical role of helices II and III in the folding mechanism. When compared to the wild type, the folding landscape of the ultrafast folding mutant exhibited a deeper folding funnel towards the folded state. The very beginning of the folding (0-10 ns) was very similar for both protein variants but it soon diverged and displayed different folding pathways. Although going through the major intermediate state is the dominant pathway for both, it was also observed that some folding went through the minor intermediate state for the mutant. The intriguing difference resulting from the mutation at two residues in helix III has been carefully analyzed and discussed in details.

Trends in template/fragment-free protein structure prediction.

Predicting the structure of a protein from its amino acid sequence is a long-standing unsolved problem in computational biology. Its solution would be of both fundamental and practical importance as the gap between the number of known sequences and the number of experimentally solved structures widens rapidly. Currently, the most successful approaches are based on fragment/template reassembly. Lacking progress in template-free structure prediction calls for novel ideas and approaches. This article reviews trends in the development of physical and specific knowledge-based energy functions as well as sampling techniques for fragment-free structure prediction. Recent physical- and knowledge-based studies demonstrated that it is possible to sample and predict highly accurate protein structures without borrowing native fragments from known protein structures. These emerging approaches with fully flexible sampling have the potential to move the field forward.

A single transcript for the prognosis of disease severity in COVID-19 patients.

With many countries strapped for medical resources due to the COVID-19 pandemic, it is highly desirable to allocate the precious resources to those who need them the most. Several markers have been found to be associated with the disease severity in COVID-19 patients. However, the established markers only display modest prognostic power individually and better markers are urgently needed. The aim of this study is to investigate the potential of S100A12, a prominent marker gene for bacterial infection, in the prognosis of disease severity in COVID-19 patients. To ensure the robustness of the association, a total of 1695 samples from 14 independent transcriptome datasets on sepsis, influenza infection and COVID-19 infection were examined. First, it was demonstrated that S100A12 was a marker for sepsis and severity of sepsis. Then, S100A12 was found to be a marker for severe influenza infection, and there was an upward trend of S100A12 expression as the severity level of influenza infection increased. As for COVID-19 infection, it was found that S100A12 expression was elevated in patients with severe and critical COVID-19 infection. More importantly, S100A12 expression at hospital admission was robustly correlated with future quantitative indexes of disease severity and outcome in COVID-19 patients, superior to established prognostic markers including CRP, PCT, d-dimer, ferritin, LDH and fibrinogen. Thus, S100A12 is a valuable novel prognostic marker for COVID-19 severity and deserves more attention.

Single-cell RNA-Seq revealed profound immune alteration in the peripheral blood of patients with bacterial infection.

OBJECTIVES: Bacterial infection remains one of the greatest threats to human health. However, how human hosts respond to bacterial infection has not been thoroughly understood. Better understanding of this response will improve human health. METHODS: Here, we conducted an investigation on host response to bacterial infection using unperturbed clinical samples and single-cell RNA-Seq (scRNA-Seq) technology. To evaluate immune alteration upon bacterial infection in scRNA-Seq data of peripheral blood mononuclear cells (PBMCs), we developed a barcode analytical framework named PBMCode. RESULTS: Using this PBMCode framework, we revealed profound immune alteration in peripheral blood under bacterial infection, including the emergence of natural killer T (NKT) cell cluster, reduction of B cell population, and considerable changes in T cells and monocytes. In addition, we also observed a large quantity of low-density neutrophils. CONCLUSIONS: Our investigation on single cells provided unprecedented details in the alteration of both cell population and cell state under bacterial infection. These findings may be relevant to clinical decisions. The complexity of host response to bacterial infection revealed by scRNA-Seq deserves further attention in future studies.

A host-based two-gene model for the identification of bacterial infection in general clinical settings.

OBJECTIVES: In this study, we aimed to develop a simple gene model to identify bacterial infection, which can be implemented in general clinical settings. METHODS: We used a clinically availablereal-time quantitative polymerase chain reaction platform to conduct focused gene expression assays on clinical blood samples. Samples were collected from 2 tertiary hospitals. RESULTS: We found that the 8 candidate genes for bacterial infection were significantly dysregulated in bacterial infection and displayed good performance in group classification, whereas the 2 genes for viral infection displayed poor performance. A two-gene model (S100A12 and CD177) displayed 93.0% sensitivity and 93.7% specificity in the modeling stage. In the independent validation stage, 87.8% sensitivity and 96.6% specificity were achieved in one set of case-control groups, and 93.6% sensitivity and 97.1% specificity in another set. CONCLUSIONS: We have validated the signature genes for bacterial infection and developed a two-gene model to identify bacterial infection in general clinical settings.

Folding network of villin headpiece subdomain.

Protein folding is a complex multidimensional process that is difficult to illustrate by the traditional analyses based on one- or two-dimensional profiles. Analyses based on transition networks have become an alternative approach that has the potential to reveal detailed features of protein folding dynamics. However, due to the lack of successful reversible folding of proteins from conventional molecular-dynamics simulations, this approach has rarely been utilized. Here, we analyzed the folding network from several 10 µs conventional molecular-dynamics reversible folding trajectories of villin headpiece subdomain (HP35). The folding network revealed more complexity than the traditional two-dimensional map and demonstrated a variety of conformations in the unfolded state, intermediate states, and the native state. Of note, deep enthalpic traps at the unfolded state were observed on the folding landscape. Furthermore, in contrast to the clear separation of the native state and the primary intermediate state shown on the two-dimensional map, the two states were mingled on the folding network, and prevalent interstate transitions were observed between these two states. A more complete picture of the folding mechanism of HP35 emerged when the traditional and network analyses were considered together.

Improved sampling methods for molecular simulation.

Molecular simulation has broad application in the biological sciences. One of the greatest challenges in molecular simulation is the limited conformational sampling due to slow barrier crossing on the rugged energy landscape of complex biomolecules and to the relatively short simulation time. Many enhanced sampling techniques have been developed over the years to alleviate this problem. Significant progress has been made in the past couple of years, with emerging methods targeting specific aspects of the potential energy surface and new variants of the replica exchange method.

Dual folding pathways of an alpha/beta protein from all-atom ab initio folding simulations.

Successful ab initio folding of proteins with both alpha-helix and beta-sheet requires a delicate balance among a variety of forces in the simulation model, which may explain that the successful folding of any alpha/beta proteins to within experimental error has yet to be reported. Here we demonstrate that it is an achievable goal to fold alpha/beta proteins with a force field emphasizing the balance between the two major secondary structures. Using our newly developed force field, we conducted extensive ab initio folding simulations on an alpha/beta protein full sequence design (FSD) employing both conventional molecular dynamics and replica exchange molecular dynamics in combination with a generalized-Born solvation model. In these simulations, the folding of FSD to the native state with high population (>64.2%) and high fidelity (C(alpha)-Root Mean Square Deviation of 1.29 A for the most sampled conformation when compared to the experimental structure) was achieved. The folding of FSD was found to follow two pathways. In the major pathway, the folding started from the formation of the helix. In the minor pathway, however, folding of the beta-hairpin started first. Further examination revealed that the helix initiated from the C-terminus and propagated toward the N-terminus. The formation of the hydrophobic contacts coincided with the global folding. Therefore the hydrophobic force does not appear to be the driving force of the folding of this protein.

Gene dysregulation in peripheral blood of moyamoya disease and comparison with other vascular disorders.

OBJECTIVE: Moyamoya disease (MMD) is a chronic occlusive cerebrovascular disease with unknown etiology, sharing many similar clinical symptoms with other vascular disorders. This study aimed to investigate gene dysregulation in peripheral blood of MMD and compare it with other vascular disorders. METHODS: Transcriptomic profiles of 12 MMD patients and 8 healthy controls were obtained using RNA sequencing. Differentially expressed genes (DEGs) were identified and several were validated by quantitative real-time PCR in independent samples. Biological pathway enrichment analysis of DEGs and deconvolution of leukocyte subsets in peripheral blood were performed. Expression profiles for other vascular diseases were downloaded from public database and consistent DEGs were calculated. Gene set enrichment analysis (GSEA) was conducted to compare gene dysregulation pattern between MMD and other vascular diseases. RESULTS: A total of 533 DEGs were identified for MMD. Up-regulated genes were mainly involved in extracellular matrix (ECM) organization, whereas down-regulated genes were primarily associated with inflammatory and immune responses. As for cell populations, significantly increased naïve B cells and naïve CD4 cells as well as obviously decreased resting natural killer cells were observed in peripheral blood of MMD patients. GSEA analysis indicated that only up-regulated genes of ischemic stroke and down-regulated genes of coronary artery disease and myocardial infarction were enriched in up-regulated and down-regulated genes of MMD, respectively. CONCLUSION: Dysregulated genes in peripheral blood of MMD mainly played key roles in ECM organization, inflammatory and immune responses. This gene dysregulation pattern was specific compared with other vascular diseases. Besides, naïve B cells, naïve CD4 cells and resting natural killer cells were aberrantly disrupted in peripheral blood of MMD patients. These results will help elucidate the complicated pathogenic mechanism of MMD.

Evaluation of Peripheral Immune Dysregulation in Alzheimer's Disease and Vascular Dementia.

Immune dysregulation has been observed in the brain and blood of patients with Alzheimer's disease (AD). However, a convenient assay to evaluate peripheral immune dysregulation in AD has not been developed, partly due to the inconsistent observations from different studies. We hypothesized that peripheral immune dysregulation may only exist in a subpopulation of AD patients; therefore it may be valuable to identify this subpopulation with a convenient assay. Along this line, we selected 14 candidate genes based on our analysis of microarray data on peripheral blood of AD and other diseases. We used RT-qPCR to examine the expression of these 14 genes in a cohort of 288 subjects, including 74 patients with AD, 64 patients with mild cognitive impairment (MCI), 51 patients with vascular dementia (VaD), and 99 elderly controls with no cognitive dysfunction/impairment. Seven of these 14 genes displayed significant difference in group comparison. Switching from group comparison to individualized evaluation revealed more in-depth information. First, there existed a wide dynamic range for the expression of these immune genes in peripheral blood even within the control group. Second, for the vast majority of the patients (AD, VaD, and MCI patients), the expression of these genes fell within the dynamic range of the control group. Third, a small portion of outliers were observed in the patient groups, more so in the VaD group than that in the AD or MCI groups. This is our first attempt to conduct personalized evaluation of peripheral immune dysregulation in AD and VaD. These findings may be applicable to the identification of peripheral immune dysregulation in AD and VaD patients which may lead to tailored treatment toward those patients.

Two-stage folding of HP-35 from ab initio simulations.

Accurate ab initio simulation of protein folding is a critical step toward elucidation of protein-folding mechanisms. Here, we demonstrate highly accurate folding of the 35 residue villin headpiece subdomain (HP35) by all-atom molecular dynamics simulations using AMBER FF03 and the generalized-Born solvation model. In a set of 20 micros long simulations, the protein folded to the native state in multiple trajectories, with the lowest C(alpha) RMSD being 0.39 A for residues 2-34 (excluding residues 1 and 35). The native state had the highest population among all sampled conformations, and the center of most populated cluster had a C(alpha) RMSD of 1.63 A. Folding of this protein can be described as a two-stage process that followed a well-defined pathway. In the first stage, formation of helices II and III as a folding intermediate constituted the rate-limiting step and was initiated at a folding nucleus around residues Phe17 and Pro21. The folding intermediate further acted as a template that facilitated the folding and docking of helix I in the second stage. Detailed descriptions of the folding kinetics and the roles of key residues are presented.

Protein folding and unfolding by all-atom molecular dynamics simulations.

Computational protein folding can be classified into pathway and sampling approaches. Here, we use the AMBER simulation package as an example to illustrate the protocols for all-atom molecular simulations of protein folding, including system setup, simulation, and analysis. We introduced two traditional pathway approaches: ab inito folding and high-temperature unfolding. The popular replica exchange method was chosen to represent sampling approaches. Our emphasis is placed on the analysis of the simulation trajectories, and some in-depth discussions are provided for commonly encountered problems.

Conformational changes in protein function.

Conformational changes are the hallmarks of protein dynamics and are often intimately related to protein functions. Molecular dynamics (MD) simulation is a powerful tool to study the time-resolved properties of protein structure in atomic details. In this chapter, we discuss the various applications of MD simulation to the study of protein conformational changes, and introduce several selected advanced techniques that may significantly increase the sampling efficiencies, including locally enhanced sampling (LES), and grow-to-fit molecular dynamics (G2FMD).

The fast-folding HP35 double mutant has a substantially reduced primary folding free energy barrier.

The LYS24/29NLE double mutant of villin headpiece subdomain (HP35) is the fastest folding protein known so far with a folding time constant of 0.6 micros. In this work, the folding mechanism of the mutant has been investigated by both conventional and replica exchange molecular dynamics (CMD and REMD) simulations with AMBER FF03 force field and a generalized-Born solvation model. Direct comparison to the ab initio folding of the wild type HP35 enabled a close examination on the mutational effect on the folding process. The mutant folded to the native state, as demonstrated by the 0.50 A C(alpha)-root mean square deviation (RMSD) sampled in both CMD and REMD simulations and the high population of the folded conformation compared with the denatured conformations. Consistent with experiments, the significantly reduced primary folding free energy barrier makes the mutant closer to a downhill folder than the wild type HP35 that directly leads to the faster transition and higher melting temperature. However, unlike the proposed downhill folding which envisages a smooth shift between unfolded and folded states without transition barrier, we observed a well-defined folding transition that was consistent with experiments. Further examination of the secondary structures revealed that the two mutated residues have higher intrinsic helical preference that facilitated the formation of both helix III and the intermediate state which contains the folded segment helix II/III. Other factors contributing to the faster folding include the more favorable electrostatic interactions in the transition state with the removal of the charged NH(3)(+) groups from LYS. In addition, both transition state ensemble and denatured state ensemble are shifted in the mutant.

Folding processes of the B domain of protein A to the native state observed in all-atom ab initio folding simulations.

Reaching the native states of small proteins, a necessary step towards a comprehensive understanding of the folding mechanisms, has remained a tremendous challenge to ab initio protein folding simulations despite the extensive effort. In this work, the folding process of the B domain of protein A (BdpA) has been simulated by both conventional and replica exchange molecular dynamics using AMBER FF03 all-atom force field. Started from an extended chain, a total of 40 conventional (each to 1.0 micros) and two sets of replica exchange (each to 200.0 ns per replica) molecular dynamics simulations were performed with different generalized-Born solvation models and temperature control schemes. The improvements in both the force field and solvent model allowed successful simulations of the folding process to the native state as demonstrated by the 0.80 A C(alpha) root mean square deviation (RMSD) of the best folded structure. The most populated conformation was the native folded structure with a high population. This was a significant improvement over the 2.8 A C(alpha) RMSD of the best nativelike structures from previous ab initio folding studies on BdpA. To the best of our knowledge, our results demonstrate, for the first time, that ab initio simulations can reach the native state of BdpA. Consistent with experimental observations, including Phi-value analyses, formation of helix II/III hairpin was a crucial step that provides a template upon which helix I could form and the folding process could complete. Early formation of helix III was observed which is consistent with the experimental results of higher residual helical content of isolated helix III among the three helices. The calculated temperature-dependent profile and the melting temperature were in close agreement with the experimental results. The simulations further revealed that phenylalanine 31 may play critical to achieve the correct packing of the three helices which is consistent with the experimental observation. In addition to the mechanistic studies, an ab initio structure prediction was also conducted based on both the physical energy and a statistical potential. Based on the lowest physical energy, the predicted structure was 2.0 A C(alpha) RMSD away from the experimentally determined structure.

Molecular dynamics simulations and free energy analyses on the dimer formation of an amyloidogenic heptapeptide from human beta2-microglobulin: implication for the protofibril structure.

Amyloid formation is associated with many neurodegenerative diseases. Recent findings suggest that early oligomeric aggregates could be major sources of toxicity. We present a computational investigation of the first step of amyloid initiation-dimer formation of a seven residue peptide (NHVTLSQ) from human beta2-microglobulin at pH 2.0, which renders +2.0 units charges to each peptide. A total of over 1.2 micros of simulations with explicit solvent and 1.0 micros of simulations with implicit solvent were conducted. Main-chain conformational restraint was applied to facilitate the formation of ordered dimers. An antiparallel beta-sheet with six main-chain hydrogen bonds was dominant in the implicit solvent simulations. In contrast, no stable dimers were observed in the two negative controls, the mouse heptapeptide (KHDSMAE, +3.0 units charges) and the scrambled human heptapeptide (QVLHTSN). Explicit solvent simulations presented a more complex scenario. The wild-type human heptapeptide formed predominantly antiparallel beta-sheets ( approximately 38%) although parallel ones ( approximately 12%) were also observed. Hydrophobic contacts preceded hydrogen bond saturation in the majority of the association events in the explicit solvent simulations, highlighting the important role of hydrophobic interaction in amyloid initiation. The fact that the mouse dimer dissociated immediately after the removal of conformational restraint suggests that the higher conformational entropy barrier, along with the stronger charge repulsion and weaker hydrophobic interaction, contributed to its inability to form amyloid fibril. The closeness of positive charge pairs in the dimers of the scrambled human heptapeptide may prohibit further beta-sheet extension and fibril growth. Combining the results from simulations and free energy analyses, we propose that the building block for this amyloid fibril is an antiparallel dimer with a two-residue register shift and six main-chain hydrogen bonds. A double-layer protofibril structure is also proposed in which two antiparallel beta-sheets face each other and are held together by hydrophobic staples and hydrogen bonds of the polar side-chains.