Conformational response to ligand binding of TMPRSS2, a protease involved in SARS-CoV-2 infection: Insights through computational modeling.

Thanks to the considerable research which has been undertaken in the last few years to improve our understanding of the biology and mechanism of action of SARS-CoV-2, we know how the virus uses its surface spike protein to infect host cells. The transmembrane prosthesis, serine 2 (TMPRSS2) protein, located on the surface of human cells, recognizes the cleavage site in the spike protein, leading to the release of the fusion peptide and entry of the virus into the host cells. Because of its role, TMPRSS2 has been proposed as a drug target to prevent infection by the virus. In this study, we aim to increase our understanding of TMPRSS2 using long scale microsecond atomistic molecular dynamics simulations, focusing on the conformational changes over time. The comparison between simulations conducted on the protein in the native (apo) and inhibited form (holo), has shown that in the holo form the inhibitor stabilizes the catalytic site and induces rearrangements in the extracellular domain of the protein. In turn, it leads to the formation of a new cavity in the vicinity of the ligand binding pocket that is stable in the microsecond time scale. Given the low specificity of known protease inhibitors, these findings suggest a new potential drug target site that can be used to improve TMPRSS2 specific recognition by newly designed inhibitors.

Computational Study of Helicase from SARS-CoV-2 in RNA-Free and Engaged Form.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the causative agent of the pandemic that broke out in 2020 and continues to be the cause of massive global upheaval. Coronaviruses are positive-strand RNA viruses with a genome of ~30 kb. The genome is replicated and transcribed by RNA-dependent RNA polymerase together with accessory factors. One of the latter is the protein helicase (NSP13), which is essential for viral replication. The recently solved helicase structure revealed a tertiary structure composed of five domains. Here, we investigated NSP13 from a structural point of view, comparing its RNA-free form with the RNA-engaged form by using atomistic molecular dynamics (MD) simulations at the microsecond timescale. Structural analyses revealed conformational changes that provide insights into the contribution of the different domains, identifying the residues responsible for domain-domain interactions in both observed forms. The RNA-free system appears to be more flexible than the RNA-engaged form. This result underlies the stabilizing role of the nucleic acid and the functional core role of these domains.

Computational Studies of SARS-CoV-2 3CLpro: Insights from MD Simulations.

Given the enormous social and health impact of the pandemic triggered by severe acute respiratory syndrome 2 (SARS-CoV-2), the scientific community made a huge effort to provide an immediate response to the challenges posed by Coronavirus disease 2019 (COVID-19). One of the most important proteins of the virus is an enzyme, called 3CLpro or main protease, already identified as an important pharmacological target also in SARS and Middle East respiratory syndrome virus (MERS) viruses. This protein triggers the production of a whole series of enzymes necessary for the virus to carry out its replicating and infectious activities. Therefore, it is crucial to gain a deeper understanding of 3CLpro structure and function in order to effectively target this enzyme. All-atoms molecular dynamics (MD) simulations were performed to examine the different conformational behaviors of the monomeric and dimeric form of SARS-CoV-2 3CLpro apo structure, as revealed by microsecond time scale MD simulations. Our results also shed light on the conformational dynamics of the loop regions at the entry of the catalytic site. Studying, at atomic level, the characteristics of the active site and obtaining information on how the protein can interact with its substrates will allow the design of molecules able to block the enzymatic function crucial for the virus.

Comparing the Montreal Cognitive Assessment and Rowland Universal Dementia Assessment Scale in a multicultural rehabilitation setting.

In Australia it is recommended that all older people undergoing rehabilitation have a cognitive screen. We performed a longitudinal study comparing the correlation of two cognitive screening tools - the Rowland Universal Dementia Assessment Scale (RUDAS) and Montreal Cognitive Assessment (MoCA) with discharge outcomes in a geriatric inpatient setting. The RUDAS cut-off (<23/30) was associated with discharge to a nursing home (sensitivity 52%, specificity 70%). This was also noted with a MoCA cut-off <18/30 (sensitivity 57%, specificity 69%). Furthermore the association between the RUDAS and discharge destination was independent of its association with the Functional Independence Measure (r = 0.116; P = 0.275) and had a shorter administration time. Both RUDAS and MoCA scores could be used as predictors of discharge destination in a multicultural population.

In silico pharmacogenetic approach: The natalizumab case study.

Natalizumab is a humanized monoclonal antibody to α4ß1 integrin and is approved for the treatment of Multiple Sclerosis. In patients there is a great variation in drug response and there is much evidence that genetic contributors play an important role in defining an individual's susceptibility. Natalizumab binds to α4-residues Gln-152, Lys-201, Lys256, and these seem to be essential for its activity. Studies on a range of species in disease model have showed a loss of reactivity when any one of those three residues were different to human. Based on these animal studies, we thought that the single nucleotide polymorphism in the ITGA4 human gene causing a lysine to arginine transversion at amino acid position 256 require further investigations in the context of individual drug susceptibility. So, the aim of our study was to investigate the association between this genetic polymorphism and the resistance to natalizumab. We had applied molecular dynamics simulation to study the possible conformational changes induced by Lys256Arg transversion on the overall structure of integrin and we have analyzed the binding affinities of natalizumab in the non-mutated and mutated structures through HINT score. We found that this SNP does not affect the VLA4-natalizumab interaction. Instead, the binding affinities are slightly higher in the mutated complex than in the wild-type. We reported one of the first work in which MD simulation was applied in the pharmacogenetic context, and this approach is rapid and cost effective, since a population survey is carried out only after the positive prediction of simulation.

ImmunoGrid, an integrative environment for large-scale simulation of the immune system for vaccine discovery, design and optimization.

Vaccine research is a combinatorial science requiring computational analysis of vaccine components, formulations and optimization. We have developed a framework that combines computational tools for the study of immune function and vaccine development. This framework, named ImmunoGrid combines conceptual models of the immune system, models of antigen processing and presentation, system-level models of the immune system, Grid computing, and database technology to facilitate discovery, formulation and optimization of vaccines. ImmunoGrid modules share common conceptual models and ontologies. The ImmunoGrid portal offers access to educational simulators where previously defined cases can be displayed, and to research simulators that allow the development of new, or tuning of existing, computational models. The portal is accessible at .

Fusobacterium necrophorum Bacteremia With Evidence of Cavitary Pulmonary Lesion.

Fusobacterium necrophorum (F. necrophorum) is a gram-negative anaerobic bacterium and a known etiologic agent in Lemierre's syndrome. This rare disease commonly presents with persistent sore throat and dysphagia, which can spread to involve the internal jugular vein. Presented in this report is an interesting case of a patient who presented with a progressively worsening sore throat, dysphagia, and productive cough on admission. Blood cultures were positive for F. necrophorum and computed tomography angiogram (CTA) of the chest detected cavitation in the left lower lobe and a large consolidation within the right lower lobe without evidence of a vascular defect. CT of the neck with IV contrast demonstrated no findings of abnormal vascular structures. This patient was diagnosed with pneumonia secondary to F. necrophorum bacteremia and treated successfully with antibiotics and was discharged home. Clinical suspicion is warranted in patients with worsening symptoms of sore throat and dysphagia, as this rare syndrome may be present.

ImmunoGrid - the virtual human immune system project.

ImmunoGrid is a 3 year project funded by the European Union which began in February 2006 and establishes an infrastructure for the simulation of the immune system that integrates processes at molecular, cellular and organ levels. It is designed for applications that support clinical outcomes such as the design of vaccines, immunotherapies and optimization of immunization protocols. The first phase of the project concentrated on improving and extending current models of the immune system. We are now entering the second phase which will design and implement a human immune system simulator. Since the new models are orders of magnitude more complex than the previous ones, grid technologies will be essential in providing the necessary computer infrastructure. The final phase of the project will validate the simulator with pre-clinical trials using mouse models.

Long-range transcriptome sequencing reveals cancer cell growth regulatory chimeric mRNA.

mRNA chimeras from chromosomal translocations often play a role as transforming oncogenes. However, cancer transcriptomes also contain mRNA chimeras that may play a role in tumor development, which arise as transcriptional or post-transcriptional events. To identify such chimeras, we developed a deterministic screening strategy for long-range sequence analysis. High-throughput, long-read sequencing was then performed on cDNA libraries from major tumor histotypes and corresponding normal tissues. These analyses led to the identification of 378 chimeras, with an unexpectedly high frequency of expression (≈2 x 10(-5) of all mRNA). Functional assays in breast and ovarian cancer cell lines showed that a large fraction of mRNA chimeras regulates cell replication. Strikingly, chimeras were shown to include both positive and negative regulators of cell growth, which functioned as such in a cell-type-specific manner. Replication-controlling chimeras were found to be expressed by most cancers from breast, ovary, colon, uterus, kidney, lung, and stomach, suggesting a widespread role in tumor development.

Simulation of urea-induced protein unfolding: a lesson from bovine ß-lactoglobulin.

To investigate the molecular mechanisms involved in the very initial stages of protein unfolding, we carried out one long (1 µs) simulation of bovine ß-lactoglobulin (BLG) together with three (500 ns) supporting MD runs, in which the unfolding conditions were produced by adding the osmolyte urea to the simulated systems and/or by increasing the thermal energy raising the temperature from 300 to 350 K. BLG was chosen, since it is a well-characterized model protein, for which structural and folding properties have been widely investigated by X-ray and NMR. MD trajectories were analyzed not only in terms of standard progress variables, such as backbone H-bonds, gyration radius width, secondary structure elements, but also through the scrutiny of interactions and dynamical behavior of specific key residues previously pointed out and investigated by NMR and belonging to a well known hydrophobic cluster. MD trajectories simulated in different unfolding conditions suggest that urea destabilizes BLG structure weakening protein::protein hydrophobic interactions and the hydrogen bond network. The early unfolding events, better observed at higher temperature, affect both secondary and tertiary structure of the protein.

ImmunoGrid: towards agent-based simulations of the human immune system at a natural scale.

The ultimate aim of the EU-funded ImmunoGrid project is to develop a natural-scale model of the human immune system-that is, one that reflects both the diversity and the relative proportions of the molecules and cells that comprise it-together with the grid infrastructure necessary to apply this model to specific applications in the field of immunology. These objectives present the ImmunoGrid Consortium with formidable challenges in terms of complexity of the immune system, our partial understanding about how the immune system works, the lack of reliable data and the scale of computational resources required. In this paper, we explain the key challenges and the approaches adopted to overcome them. We also consider wider implications for the present ambitious plans to develop natural-scale, integrated models of the human body that can make contributions to personalized health care, such as the European Virtual Physiological Human initiative. Finally, we ask a key question: How long will it take us to resolve these challenges and when can we expect to have fully functional models that will deliver health-care benefits in the form of personalized care solutions and improved disease prevention?

High-pressure dissociation of crystalline para-diiodobenzene: optical experiments and Car-Parrinello calculations.

We have investigated the high-pressure properties of the molecular crystal para-diiodobenzene, by combining optical absorption, reflectance, and Raman experiments with Car-Parrinello simulations. The optical absorption edge exhibits a large red shift from 4 eV at ambient conditions to about 2 eV near 30 GPa. Reflectance measurements up to 80 GPa indicate a redistribution of oscillator strength toward the near-infrared. The calculations, which describe correctly the two known molecular crystal phases at ambient pressure, predict a nonmolecular metallic phase, stable at high pressure. This high-density phase is characterized by an extended three-dimensional network, in which chemically bound iodine atoms form layers connected by hydrocarbon bridges. Experimentally, Raman spectra of samples recovered after compression show vibrational modes of elemental solid iodine. This result points to a pressure-induced molecular dissociation process which leads to the formation of domains of iodine and disordered carbon.