Search

Dengue virus (DENV) has been expanding its range to temperate areas that are not usually affected, where the spread of vectors has been facilitated by global trade and climate change. In Europe, there have been many cases of DENV imported from other regions in the past few years, leading to local outbreaks of DENV among people without travel history. Here we describe the epidemiological and molecular investigations of three transmission events locally acquired DENV infections caused by serotypes 1, 2 and 3, respectively, in the Latium Region from August to November 2023. Next-generation or Sanger sequencing was used to obtain the whole genomes, or the complete E-gene of the viruses, respectively. The structure of the DENV-1 and DENV-3 sequences was analysed to identify amino acid changes that were not found in the closest related sequences. The major cluster was supported by DENV-1 (originated in South America), with 42 autochthonous infections almost occurring in the eastern area of Rome, probably due to a single introduction followed by local sustained transmission. Seven DENV-1 subclusters have been identified by mutational and phylogenetic analysis. Structural analysis indicated changes whose meaning can be explained by the adaptation of the virus to human hosts and vectors and their interactions with antibodies and cell receptors.

Mutations in the SARS-CoV-2 Spike glycoprotein can affect monoclonal antibody efficacy. Recent findings report the occurrence of resistant mutations in immunocompromised patients after tixagevimab/cilgavimab treatment. More recently, the Food and Drug Agency revoked the authorization for tixagevimab/cilgavimab, while this monoclonal antibody cocktail is currently recommended by the European Medical Agency. We retrospectively reviewed 22 immunocompetent patients at high risk for disease progression who received intramuscular tixagevimab/cilgavimab as early COVID-19 treatment and presented a prolonged high viral load. Complete SARS-CoV-2 genome sequences were obtained for a deep investigation of mutation frequencies in Spike protein before and during treatment. At seven days, only one patient showed evidence of treatment-emergent cilgavimab resistance. Quasispecies analysis revealed two different deletions on the Spike protein (S:del138-144 or S:del141-145) in combination with the resistance S:K444N mutation. The structural and dynamic impact of the two quasispecies was characterized by using molecular dynamics simulations, showing the conservation of the principal functional movements in the mutated systems and their capabilities to alter the structure and dynamics of the RBD, responsible for the interaction with the ACE2 human receptor. Our study underlines the importance of prompting an early virological investigation to prevent drug resistance or clinical failures in immunocompetent patients.

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike (S) protein binds angiotensin-converting enzyme 2 as its primary infection mechanism. Interactions between S and endogenous proteins occur after infection but are not well understood. We profiled binding of S against >9000 human proteins and found an interaction between S and human estrogen receptor α (ERα). Using bioinformatics, supercomputing, and experimental assays, we identified a highly conserved and functional nuclear receptor coregulator (NRC) LXD-like motif on the S2 subunit. In cultured cells, S DNA transfection increased ERα cytoplasmic accumulation, and S treatment induced ER-dependent biological effects. Non-invasive imaging in SARS-CoV-2-infected hamsters localized lung pathology with increased ERα lung levels. Postmortem lung experiments from infected hamsters and humans confirmed an increase in cytoplasmic ERα and its colocalization with S in alveolar macrophages. These findings describe the discovery of a S-ERα interaction, imply a role for S as an NRC, and advance knowledge of SARS-CoV-2 biology and coronavirus disease 2019 pathology.

Vesicle biogenesis, trafficking and signaling via Endoplasmic reticulum-Golgi network support essential developmental processes and their disruption lead to neurodevelopmental disorders and neurodegeneration. We report that de novo missense variants in ARF3, encoding a small GTPase regulating Golgi dynamics, cause a developmental disease in humans impairing nervous system and skeletal formation. Microcephaly-associated ARF3 variants affect residues within the guanine nucleotide binding pocket and variably perturb protein stability and GTP/GDP binding. Functional analysis demonstrates variably disruptive consequences of ARF3 variants on Golgi morphology, vesicles assembly and trafficking. Disease modeling in zebrafish validates further the dominant behavior of the mutants and their differential impact on brain and body plan formation, recapitulating the variable disease expression. In-depth in vivo analyses traces back impaired neural precursors' proliferation and planar cell polarity-dependent cell movements as the earliest detectable effects. Our findings document a key role of ARF3 in Golgi function and demonstrate its pleiotropic impact on development., (© 2022. The Author(s).)

A large number of SARS-CoV-2 mutations in a short period of time has driven scientific research related to vaccines, new drugs, and antibodies to combat the new variants of the virus. Herein, we present a web portal containing the structural information, the tridimensional coordinates, and the molecular dynamics trajectories of the SARS-CoV-2 spike protein and its main variants. The Spike Mutants website can serve as a rapid online tool for investigating the impact of novel mutations on virus fitness. Taking into account the high variability of SARS-CoV-2, this application can help the scientific community when prioritizing molecules for experimental assays, thus, accelerating the identification of promising drug candidates for COVID-19 treatment. Below we describe the main features of the platform and illustrate the possible applications for speeding up the drug discovery process and hypothesize new effective strategies to overcome the recurrent mutations in SARS-CoV-2 genome.

The vast amount of epidemiologic and genomic data that were gathered as a global response to the COVID-19 pandemic that was caused by SARS-CoV-2 offer a unique opportunity to shed light on the structural evolution of coronaviruses and in particular on the spike (S) glycoprotein, which mediates virus entry into the host cell by binding to the human ACE2 receptor. Herein, we carry out an investigation into the dynamic properties of the S glycoprotein, focusing on the much more transmissible Delta and Omicron variants. Notwithstanding the great number of mutations that have accumulated, particularly in the Omicron S glycoprotein, our data clearly showed the conservation of some structural and dynamic elements, such as the global motion of the receptor binding domain (RBD). However, our studies also revealed structural and dynamic alterations that were concentrated in the aa 627-635 region, on a small region of the receptor binding motif (aa 483-485), and the so-called "fusion-peptide proximal region". In particular, these last two S regions are known to be involved in the human receptor ACE2 recognition and membrane fusion. Our structural evidence, therefore, is likely involved in the observed different transmissibility of these S mutants. Finally, we highlighted the role of glycans in the increased RBD flexibility of the monomer in the up conformation of Omicron.

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike (S) protein binds angiotensin-converting enzyme 2 (ACE2) at the cell surface, which constitutes the primary mechanism driving SARS-CoV-2 infection. Molecular interactions between the transduced S and endogenous proteins likely occur post-infection, but such interactions are not well understood. We used an unbiased primary screen to profile the binding of full-length S against >9,000 human proteins and found significant S-host protein interactions, including one between S and human estrogen receptor alpha (ERα). After confirming this interaction in a secondary assay, we used bioinformatics, supercomputing, and experimental assays to identify a highly conserved and functional nuclear receptor coregulator (NRC) LXD-like motif on the S2 subunit and an S-ERα binding mode. In cultured cells, S DNA transfection increased ERα cytoplasmic accumulation, and S treatment induced ER-dependent biological effects and ACE2 expression. Noninvasive multimodal PET/CT imaging in SARS-CoV-2-infected hamsters using [ 18 F]fluoroestradiol (FES) localized lung pathology with increased ERα lung levels. Postmortem experiments in lung tissues from SARS-CoV-2-infected hamsters and humans confirmed an increase in cytoplasmic ERα expression and its colocalization with S protein in alveolar macrophages. These findings describe the discovery and characterization of a novel S-ERα interaction, imply a role for S as an NRC, and are poised to advance knowledge of SARS-CoV-2 biology, COVID-19 pathology, and mechanisms of sex differences in the pathology of infectious disease.

The COVID-19 pandemic is caused by SARS-CoV-2. Currently, most of the research efforts towards the development of vaccines and antibodies against SARS-CoV-2 were mainly focused on the spike (S) protein, which mediates virus entry into the host cell by binding to ACE2. As the virus SARS-CoV-2 continues to spread globally, variants have emerged, characterized by multiple mutations of the S glycoprotein. Herein, we employed microsecond-long molecular dynamics simulations to study the impact of the mutations of the S glycoprotein in SARS-CoV-2 Variant of Concern 202012/01 (B.1.1.7), termed the "UK variant", in comparison with the wild type, with the aim to decipher the structural basis of the reported increased infectivity and virulence. The simulations provided insights on the different dynamics of UK and wild-type S glycoprotein, regarding in particular the Receptor Binding Domain (RBD). In addition, we investigated the role of glycans in modulating the conformational transitions of the RBD. The overall results showed that the UK mutant experiences higher flexibility in the RBD with respect to wild type; this behavior might be correlated with the increased transmission reported for this variant. Our work also adds useful structural information on antigenic "hotspots" and epitopes targeted by neutralizing antibodies.

Essential to understanding life, the biomolecular phenomena have been an important subject in science, therefore a necessary path to be covered to make progress in human knowledge. To fully comprehend these processes, the non-covalent interactions are the key. In this review, we discuss how specific protein-ligand interactions can be efficiently described by low computational cost methods, such as Molecular Mechanics (MM). We have taken as example the case of the halogen bonds (XB). Albeit generally weaker than the hydrogen bonds (HB), the XBs play a key role to drug design, enhancing the affinity and selectivity toward the biological target. Along with the attraction between two electronegative atoms in XBs explained by the σ-hole model, important orbital interactions, as well as relief of Pauli repulsion take place. Nonetheless, such electronic effects can be only well-described by accurate quantum chemical methods that have strong limitations dealing with supramolecular systems due to their high computational cost. To go beyond the poor description of XBs by MM methods, reparametrizing the force-fields equations can be a way to keep the balance between accuracy and computational cost. Thus, we have shown the steps to be considered when parametrizing force-fields to achieve reliable results of complex non-covalent interactions at MM level for In Silico drug design methods., (Copyright © 2019 Santos, Prandi and Ramalho.)

The brain has a unique biological complexity and is responsible for important functions in the human body, such as the command of cognitive and motor functions. Disruptive disorders that affect this organ, e.g. neurodegenerative diseases (NDDs), can lead to permanent damage, impairing the patients' quality of life and even causing death. In spite of their clinical diversity, these NDDs share common characteristics, such as the accumulation of specific proteins in the cells, the compromise of the metal ion homeostasis in the brain, among others. Despite considerable advances in understanding the mechanisms of these diseases and advances in the development of treatments, these disorders remain uncured. Considering the diversity of mechanisms that act in NDDs, a wide range of compounds have been developed to act by different means. Thus, promising compounds with contrasting properties, such as chelating agents and metal-based drugs have been proposed to act on different molecular targets as well as to contribute to the same goal, which is the treatment of NDDs. This review seeks to discuss the different roles and recent developments of metal-based drugs, such as metal complexes and metal chelating agents as a proposal for the treatment of NDDs., Competing Interests: The authors confirm no conflicts of interest in this article’s contents.

The serine/threonine protein phosphatase type 5 (PP5) is a promising target for designing new antitumor drugs. This enzyme is a member of the PPP phosphatases gene family, which catalyzes a dephosphorylation reaction: a regulatory process in the signal transduction pathway that controls various biological processes. The aim of this work is to study and compare the inhibition of PP5 by ten cantharidin-like inhibitors in order to bring about contributions relevant to the better comprehension of their inhibitory activity. In this theoretical investigation, we used molecular dynamics techniques to understand the role of key interactions that occur in the protein active site; QM calculations were employed to study the interaction mode of these inhibitors in the enzyme. In addition, atoms in molecules (AIM) calculations were carried out to characterize the chemical bonds among the atoms involved and investigate the orbital interactions with their respective energy values. The obtained results suggest that the Arg275, Asn303, His304, His352, Arg400, His427, Glu428, Val429, Tyr451, and Phe446 residues favorably contribute to the interactions between inhibitors and PP5. However, the Asp271 and Asp244 amino acid residues do not favor such interactions for some inhibitors. Through the QM calculations, we can suggest that the reactional energy of the coordination mechanism of these inhibitors in the PP5 active site is quite important and is responsible for the inhibitory activity. The AIM technique employed in this work was essential to get a better comprehension of the transition states acquired from the mechanism simulation. This work offers insights of how cantharidin-like inhibitors interact with human PP5, potentially allowing the design of more specific and even less cytotoxic drugs for cancer treatments. Graphical Abstract Interactions of cantharidin-like inhibitors with human protein phosphatase-5 in a Mg 2+ system.

Photosynthetic antenna proteins can be thought of as "programmed solvents", which bind pigments at specific mutual orientations, thus tuning the overall energetic landscape and ensuring highly efficient light-harvesting. While positioning of chlorophyll cofactors is well understood and rationalized by the principle of an "energy funnel", the carotenoids still pose many open questions. Particularly, their short excited state lifetime (<25 ps) renders them potential energy sinks able to compete with the reaction centers and drastically undermine light-harvesting efficiency. Exploration of the orientational phase-space revealed that the placement of central carotenoids minimizes their interaction with the nearest chlorophylls in the plant antenna complexes LHCII, CP26, CP29 and LHCI. At the same time we show that this interaction is highly sensitive to structural perturbations, which has a profound effect on the overall lifetime of the complex. This links the protein dynamics to the light-harvesting regulation in plants by the carotenoids.

Classical molecular dynamics is more and more often coupled to quantum mechanical based techniques as a statistical tool to sample configurations of molecular systems embedded in complex environments. Nonetheless, the classical potentials describing the molecular systems are seldom parametrized to reproduce electronic processes, such as electronic excitations, which are instead very sensitive to the underlining description of the molecular structure. Here, we analyze the challenging case of the peridinin molecule, a natural apocarotenoid responsible for the light-harvesting process in the PCP antenna protein of dinoflagellates. Ground-state structural and vibrational properties, as well as electronic transitions of the pigment are studied by means of quantum-mechanical static and dynamic calculations. Thereafter, classical molecular dynamics simulations are performed with a number of different force-fields, ranging from a popular, general purpose one to refined potentials of increasing level of complexity. From the comparison of classical results with their quantum mechanical counterparts, it appears that, while very poor results are obtained from standard transferrable force-fields, specifically tuned potentials are able to correctly characterize most of the structural and vibrational features of the pigment. Nonetheless, only an advanced parametrization technique is able to give a semiquantitative description of the coupling between vibrations and electronic excitations, thus suggesting that the use of classical MD in combination of QM calculations for the study of photoinduced processes, albeit possible, should be considered with care.

Carotenoids are important actors both in light-harvesting (LH) and in photoprotection functions of photosynthetic pigment-protein complexes. A deep theoretical investigation of this multiple role is still missing owing to the difficulty of describing the delicate interplay between electronic and nuclear degrees of freedom. A possible strategy is to combine accurate quantum mechanical (QM) methods with classical molecular dynamics. To do this, however, accurate force-fields (FF) are necessary. This article presents a new FF for the different carotenoids present in LH complexes of plants. The results show that all the important structural properties described by the new FF are in very good agreement with QM reference values. This increased accuracy in the simulation of the structural fluctuations is also reflected in the description of excited states. Both the energy order and the different nature of the lowest singlet states are preserved during the dynamics when the new FF is used, whereas an unphysical mixing is found when a standard FF is used., (© 2016 Wiley Periodicals, Inc.)

Photosynthetic organisms employ several photoprotection strategies to avoid damage due to the excess energy in high light conditions. Among these, quenching of triplet chlorophylls by neighboring carotenoids (Cars) is fundamental in preventing the formation of singlet oxygen. Cars are able to accept the triplets from chlorophylls by triplet energy transfer (TET). We have here studied TET rates in CP29, a minor light-harvesting complex (LHC) of the Photosystem II in plants. A fully atomistic strategy combining classical molecular dynamics of the LHC in its natural environment with a hybrid time-dependent density functional theory/polarizable MM description of the TET is used. We find that the structural fluctuations of the pigment-protein complex can largely enhance the transfer rates with respect to those predicted using the crystal structure, reducing the triplet quenching times in the subnanosecond scale. These findings add a new perspective for the interpretation of the photoprotection function and its relation with structural motions of the LHC.

Light-harvesting pigment-protein complexes (PPC) represent the fundamental units through which the photosynthetic organisms absorb sunlight and funnel the energy to the reaction centre for carrying out the primary energy conversion reactions of photosynthesis. Here we apply a multiscale computational strategy to a specific PPC present in the photosystem II of plants and algae (CP29) to investigate in what detail should the environment effects due to protein and membrane/solvent be included for an accurate description of optical spectra. We find that a refinement of the crystal structure is needed before any meaningful quantum chemical calculations of pigment transition energies can be performed. For this purpose we apply classical molecular dynamics simulations of the PPC within its natural environment and we perform ab initio computations of the exciton Hamiltonian of the complex, including the environment either implicitly by the polarizable continuum model (PCM) or explicitly using the polarizable QM/MM methodology (MMPol). However, PCM essentially leads to an unspecific redshift of all transition energies, and MMPol is able to reveal site-specific changes in the optical properties of the pigments. Based on the latter and the excitonic couplings obtained within a polarizable QM/MM methodology, optical spectra are calculated, which are in good qualitative agreement with experimental data. A weakness of the approach is however found in the overestimation of the fluctuations of the excitonic parameters of the pigments along the MD trajectory. An explanation for such a finding in terms of the limits of the force fields commonly used for protein cofactors is presented and discussed.