A Portrait of the Chromophore as a Young System-Quantum-Derived Force Field Unraveling Solvent Reorganization upon Optical Excitation of Cyclocurcumin Derivatives.

The study of fast non-equilibrium solvent relaxation in organic chromophores is still challenging for molecular modeling and simulation approaches, and is often overlooked, even in the case of non-adiabatic dynamics simulations. Yet, especially in the case of photoswitches, the interaction with the environment can strongly modulate the photophysical outcomes. To unravel such a delicate interplay, in the present contribution we resorted to a mixed quantum-classical approach, based on quantum mechanically derived force fields. The main task is to rationalize the solvent reorganization pathways in chromophores derived from cyclocurcumin, which are suitable for light-activated chemotherapy to destabilize cellular lipid membranes. The accurate and reliable decryption delivered by the quantum-derived force fields points to important differences in the solvent's reorganization, in terms of both structure and time scale evolution.

Interaction between Yersinia pestis Ail Outer Membrane Protein and the C-Terminal Domain of Human Vitronectin.

Yersinia pestis, the causative agent of plague, is capable of evading the human immune system response by recruiting the plasma circulating vitronectin proteins, which act as a shield and avoid its lysis. Vitronectin recruitment is mediated by its interaction with the bacterial transmembrane protein Ail, protruding from the Y. pestis outer membrane. By using all-atom long-scale molecular dynamic simulations of Ail embedded in a realistic model of the bacterial membrane, we have shown that vitronectin forms a stable complex, mediated by interactions between the disordered moieties of the two proteins. The main amino acids driving the complexation have also been evidenced, thus favoring the possible rational design of specific peptides which, by inhibiting vitronectin recruitment, could act as original antibacterial agents.

Interaction of the Immune System TIM-3 Protein with a Model Cellular Membrane Containing Phosphatidyl-Serine Lipids.

T cell transmembrane, Immunoglobulin, and Mucin (TIM) are important immune system proteins which are especially present in T-cells and regulated the immune system by sensing cell engulfment and apoptotic processes. Their role is exerted by the capacity to detect the presence of phosphatidyl-serine lipid polar head in the outer leaflet of cellular membranes (correlated with apoptosis). In this contribution by using equilibrium and enhanced sampling molecular dynamics simulation we unravel the molecular bases and the thermodynamics of TIM, and in particular TIM-3, interaction with phosphatidyl serine in a lipid bilayer. Since TIM-3 deregulation is an important factor of pro-oncogenic tumor micro-environment understanding its functioning at a molecular level may pave the way to the development of original immunotherapeutic approaches.

Synthesis of Original Cyclic Dinucleotide Analogues Using the Sulfo-click Reaction.

The stimulator of interferon genes (STING) protein plays a crucial role in the activation of the innate immune response. Activation of STING is initiated by cyclic dinucleotides (CDNs) which prompted the community to synthesize structural analogues to enhance their biological properties. We present here the synthesis and biological evaluation of four novel CDN analogues composed of an N-acylsulfonamide linkage. These CDNs were obtained in high overall yields via the sulfo-click reaction as a key step.

Electronic Effects in Phosphino-Iminophosphorane PdII Complexes upon Varying the N Substituent.

Three series of palladium(II) complexes supported by a phosphine-iminophosphorane ligand built upon an ortho-phenylene core were investigated to study the influence of the iminophosphorane N substituent. Cis-dichloride palladium(II) complexes 1 in which the N atom bears an isopropyl (iPr, 1 a), a phenyl (Ph, 1 b), a trimethylsilyl (TMS, 1 c) group or an H atom (1 d) were synthesized in high yield. They were characterized by NMR, IR spectroscopy, HR-mass spectrometry, elemental analysis, and X-ray diffraction. A substantial bond length difference between the Pd-Cl bonds was observed in 1. Complexes 1 a-d were converted into [Pd(LR )Cl(CNt Bu)](OTf)] 2 a-d whose isocyanide is located trans to the iminophosphorane. The corresponding dicationic complexes [Pd(LR )(CNt Bu)2 ](OTf)2 3 a-d were also synthesized, however they exhibited lower stability in solution than 2, the isopropyl derivative 3 a being the most stable of the series. Molecular modeling was performed to rationalize the regioselectivity of the substitution of the single chloride by isocyanide (from 1 to 2) and to study the electronic distribution in the complexes. In particular differences between the TMS and H containing complexes vs. the iPr and Ph ones were found. This suggests that the nature of the N substituent is far from innocent and can help tune the reactivity of iminophosphorane complexes.

Molecular Photoswitches Regulating the Activity of the Human Serotonin Transporter.

Serotonin is an essential mediator regulating diverse neural processes, and its deregulation is related to the development of debilitating neurological diseases. In particular, the human serotonin transporter (hSERT) is fundamental in completing the synaptic neural cycle by allowing reuptake of serotonin. Its inhibition is particularly attractive, especially as a pharmacological target against depressive syndrome. Here, we analyze, by using long-range molecular dynamic simulations, the behavior of a molecular photoswitch whose cis- and trans-isomers inhibit the hSERT differently. In particular, we evidence the structural and molecular basis behind the higher inhibiting capacity of the cis-isomer, which blocks more efficiently the hSERT conformational cycle, leading to serotonin uptake.

Resolving a guanine-quadruplex structure in the SARS-CoV-2 genome through circular dichroism and multiscale molecular modeling.

The genome of SARS-CoV-2 coronavirus is made up of a single-stranded RNA fragment that can assume a specific secondary structure, whose stability can influence the virus's ability to reproduce. Recent studies have identified putative guanine quadruplex sequences in SARS-CoV-2 genome fragments that are involved in coding for both structural and non-structural proteins. In this contribution, we focus on a specific G-rich sequence referred to as RG-2, which codes for the non-structural protein 10 (Nsp10) and assumes a guanine-quadruplex (G4) arrangement. We provide the secondary structure of RG-2 G4 at atomistic resolution by molecular modeling and simulation, validated by the superposition of experimental and calculated electronic circular dichroism spectra. Through both experimental and simulation approaches, we have demonstrated that pyridostatin (PDS), a widely recognized G4 binder, can bind to and stabilize RG-2 G4 more strongly than RG-1, another G4 forming sequence that was previously proposed as a potential target for antiviral drug candidates. Overall, this study highlights RG-2 as a valuable target to inhibit the translation and replication of SARS-CoV-2, paving the way towards original therapeutic approaches against emerging RNA viruses.

Mechanism of the Covalent Inhibition of Human Transmembrane Protease Serine 2 as an Original Antiviral Strategy.

The Transmembrane Protease Serine 2 (TMPRSS2) is a human enzyme which is involved in the maturation and post-translation of different proteins. In addition to being overexpressed in cancer cells, TMPRSS2 plays a further fundamental role in favoring viral infections by allowing the fusion of the virus envelope with the cellular membrane, notably in SARS-CoV-2. In this contribution, we resort to multiscale molecular modeling to unravel the structural and dynamical features of TMPRSS2 and its interaction with a model lipid bilayer. Furthermore, we shed light on the mechanism of action of a potential inhibitor (nafamostat), determining the free-energy profile associated with the inhibition reaction and showing the facile poisoning of the enzyme. Our study, while providing the first atomistically resolved mechanism of TMPRSS2 inhibition, is also fundamental in furnishing a solid framework for further rational design targeting transmembrane proteases in a host-directed antiviral strategy.

Modelling the effects of E/Z photoisomerization of a cyclocurcumin analogue on the properties of cellular lipid membranes.

The use of photosensitive molecules capable of isomerizing under light stimuli, and thus induce perturbation in biological systems, is becoming increasingly popular for potential light-activated chemotherapeutic purposes. We recently show that a cyclocurcumin derivative (CCBu), may be suitable for light-activated chemotherapy and may constitute a valuable alternative to traditional photodynamic therapy, due to its oxygen-independent mechanism of action, which allows the treatment of hypoxic solid tumors. In particular, we have shown that the E/Z photoisomerization of CCBu correlates with strong perturbations of model lipid bilayers. In this work, we perform all-atom classical molecular dynamics for a more complex bilayer, whose composition is, thus, much closer to eukaryotic outer cell membranes. We have evidenced important differences in the interaction pathway between CCBu and the complex lipid bilayer as compared to previous models, concerning both the membrane penetration capacity and the isomerization-induced perturbations. While we confirm that structural perturbations of the lipid membrane are induced by isomerization, we also show how the use of a simplified membrane model can result in an oversimplification of the system and hinder key physical and biological phenomena. Although, CCBu may be considered as a suitable candidate for light-activated chemotherapy, we also underline how the inclusion of bulkier substituents, inducing larger perturbations upon photoisomerization, may enhance its efficiency.

Rationalizing the environment-dependent photophysical behavior of a DNA luminescent probe by classical and non-adiabatic molecular dynamics simulations.

Environment-sensitive fluorescent nucleoside analogs are of utmost importance to investigate the structure of nucleic acids, their intrinsic flexibility, and sequence-specific DNA- and RNA-binding proteins. The latter play indeed a key role in transcription, translation as well as in the regulation of RNA stability, localization and turnover, and many other cellular processes. The sensitivity of the embedded fluorophore to polarity, hydration, and base stacking is clearly dependent on the specific excited-state relaxation mechanism and can be rationalized combining experimental and computational techniques. In this work, we elucidate the mechanisms leading to the population of the triplet state manifold for a versatile nucleobase surrogate, namely the 2-thienyl-3-hydroxychromone in gas phase, owing to non-adiabatic molecular dynamics simulations. Furthermore, we analyze its behavior in the B-DNA environment via classical molecular dynamics simulations, which evidence a rapid extrusion of the adenine facing the 2-thienyl-3-hydroxychromone nucleobase surrogate. Our simulations provide new insights into the dynamics of this family of chromophores, which could give rise to an integrated view and a fine tuning of their photochemistry, and namely the role of excited-state intramolecular proton transfer for the rational design of the next generation of fluorescent nucleoside analogs.

Predicting the Three-Dimensional Structure of the c-KIT Proto-Oncogene Promoter and the Dynamics of Its Strongly Coupled Guanine Quadruplexes.

Guanine quadruplexes (G4s) play essential protective and regulatory roles within cells, influencing gene expression. In several gene-promoter regions, multiple G4-forming sequences are in close proximity and may form three-dimensional arrangements. We analyze the interplay among the three neighboring G4s in the c-KIT proto-oncogene promoter (WK1, WSP, and WK2). We highlight that the three G4s are structurally linked and their cross-talk favors the formation of a parallel structure for WSP. Relying on all-atom molecular dynamic simulations exceeding the µs time scale and using enhanced sampling methods, we provide the first computationally resolved structure of a well-organized G4 cluster in the promoter of a crucial gene involved in cancer development. Our results indicate that neighboring G4s influence their mutual three-dimensional arrangement and provide a powerful tool to predict and interpret complex DNA structures that can ultimately be used as a starting point for drug discovery.

How SARS-CoV-2 Alters the Regulation of Gene Expression in Infected Cells.

Nonstructural accessory proteins in viruses play a key role in hijacking the basic cellular mechanisms, which is essential to promote the virus survival and evasion of the immune system. The immonuglobulin-like open reading frame 8 (ORF8) protein expressed by SARS-CoV-2 accumulates in the nucleus and may influence the regulation of the gene expression in infected cells. In this contribution, by using microsecond time-scale all-atom molecular dynamics simulations, we unravel the structural bases behind the epigenetic action of ORF8. In particular, we highlight how the protein is able to form stable aggregates with DNA through a histone tail-like motif, and how this interaction is influenced by post-translational modifications, such as acetylation and methylation, which are known epigenetic markers in histones. Our work not only clarifies the molecular mechanisms behind the perturbation of the epigenetic regulation caused by the viral infection but also offers an unusual perspective which may foster the development of original antivirals.

Molecular Basis of the pH-Controlled Maturation of the Tick-Borne Encephalitis Flavivirus.

Flaviviruses are enveloped viruses causing high public concerns. Their maturation spans several cellular compartments having different pH. Thus, complex control mechanisms are in place to avoid premature maturation. Here we report the dynamical behavior at neutral and acidic pH of the precursor of the membrane fusion protein E of tick-borne encephalitis, showing the different stabilizations of the E dimer and the role played by the small fusion-assisting protomer (pr). The comprehension, at atomic resolution, of the fine regulation of viral maturation will be fundamental to the development of efficient strategies against emerging viral threats.

Revealing the Molecular Interactions between Human ACE2 and the Receptor Binding Domain of the SARS-CoV-2 Wild-Type, Alpha and Delta Variants.

After a sudden and first spread of the pandemic caused by the novel SARS-CoV-2 (Severe Acute Respiratory Syndrome-Coronavirus 2) wild-type strain, mutants have emerged which have been associated with increased infectivity, inducing surges in the contagions. The first of the so-called variants of concerns, was firstly isolated in the United Kingdom and later renamed Alpha variant. Afterwards, in the middle of 2021, a new variant appeared called Delta. The latter is characterized by the presence of point mutations in the Spike protein of SARS-CoV-2, especially in the Receptor Binding Domain (RBD). When in its active conformation, the RBD can interact with the human receptor Angiotensin-Converting Enzyme 2 (ACE2) to allow the entry of the virions into cells. In this contribution, by using extended all-atom molecular dynamic simulations, complemented with machine learning post-processing, we analyze the changes in the molecular interaction network induced by these different strains in comparison with the wild-type. On one hand, although relevant variations are evidenced, only limited changes in the global stability indicators and in the flexibility profiles have been observed. On the other hand, key differences were obtained by tracking hydrophilic and hydrophobic molecular interactions, concerning both positioning at the ACE2/RBD interface and formation/disruption dynamic behavior.

Understanding the Interactions of Guanine Quadruplexes with Peptides as Novel Strategies for Diagnosis or Tuning Biological Functions.

Guanine quadruplexes (G4s) are nucleic acid structures exhibiting a complex structural behavior and exerting crucial biological functions in both cells and viruses. The specific interactions of peptides with G4s, as well as an understanding of the factors driving the specific recognition are important for the rational design of both therapeutic and diagnostic agents. In this review, we examine the most important studies dealing with the interactions between G4s and peptides, highlighting the strengths and limitations of current analytic approaches. We also show how the combined use of high-level molecular simulation techniques and experimental spectroscopy is the best avenue to design specifically tuned and selective peptides, thus leading to the control of important biological functions.

Phosphonium Ylides vs Iminophosphoranes: The Role of the Coordinating Ylidic Atom in cis-[Phosphine-Ylide Rh(CO)2] Complexes.

The coordinating properties of two families of ylides, namely, phosphonium ylides and iminophosphoranes, differently substituted at the ylidic center (CH2- vs NiPr-), have been investigated in structurally related cationic phosphine-ylide Rh(CO)2 complexes obtained from readily available phosphine-phosphonium salt precursors derived from an ortho-phenylene bridge. However, while the Rh(CO)2 complex bearing the P+-CH2- donor moiety proved to be stable, the P═NiPr donor end appeared to induce lability to one of the CO groups. All of the RhI carbonyl complexes in both ylide series were fully characterized, including through X-ray diffraction analysis. Based on the experimental and calculated infrared (IR) CO stretching frequencies in Rh(CO)2 complexes, we evidenced that the phosphonium ylide ligand is a stronger donor than the iminophosphorane ligand. However, we also found that the difference in the intrinsic electronic properties can be largely compensated by the introduction of an iPr substituent on the N atom of the iminophosphorane, hence pointing to the noninnocent role of the peripheral substituent and opening novel possibilities to tune the properties of metal complexes containing ylide ligands.

A posteriori localization of many-body excited states through simultaneous diagonalization.

In this paper we propose a numerical method to localize many-electron excited states. To characterize the electronic structure of the electronic excited states of a system, quantum chemistry methods typically yield a delocalized description of the excitations. Some a priori localization methods have been developed to provide an intuitive local picture of the excited states. They typically require a good strategy to separate the system of interest from its environment, or a set of a priori localized orbitals, that may reduce their computational accuracy. Here, we introduce an a posteriori method to localize delocalized many-body excited states directly obtained from quantum chemistry calculations. A localization metric for the excited states is defined from their representation as electron-hole pairs, which is encoded in the transition density matrix. This novel a posteriori strategy thus allows to localize excitons within a volume around selected fragments of a complex molecular system without tempering with its quantum chemical treatment. The method is tested on π-stacked oligomers of phenanthrenes and pyrenes. It is found to efficiently localize and separate the excitons according to their character while preserving the information about delocalized many-body states at a low computational cost.

Perturbation of Lipid Bilayers by Biomimetic Photoswitches Based on Cyclocurcumin.

The use of photoswitches which may be activated by suitable electromagnetic radiation is an attractive alternative to conventional photodynamic therapy. Here, we report all-atom molecular dynamics simulation of a biomimetic photoswitch derived from cyclocurcumin and experiencing E/Z photoisomerization. In particular, we show that the two isomers interact persistently with a lipid bilayer modeling a cellular membrane. Furthermore, the interaction with the membrane is strongly dependent on the concentration, and a transition between ordered and disordered arrangements of the photoswitches is observed. We also confirm that the structural parameters of the bilayer are differently affected by the two isomers and hence can be modulated through photoswitching, offering interesting perspectives for future applications.

Biomimetic Photo-Switches Softening Model Lipid Membranes.

We report the synthesis and characterization of a novel photo-switch based on biomimetic cyclocurcumin analogous and interacting with the lipid bilayer, which can be used in the framework of oxygen-independent light-induced therapy. More specifically, by using molecular dynamics simulations and free energy techniques, we show that the inclusion of hydrophobic substituents is needed to allow insertion in the lipid membrane. After having confirmed experimentally that the substituents do not preclude the efficient photoisomerization, we show through UV-vis and dynamic light scattering measurements together with compression isotherms that the chromophore is internalized in both lipid vesicles and monomolecular film, respectively, inducing their fluidification. The irradiation of the chromophore-loaded lipid aggregates modifies their properties due to the different organization of the two diastereoisomers, E and Z. In particular, a competition between a fast structural reorganization and a slower expulsion of the chromophore after isomerization can be observed in the kinetic profiles recorded during E to Z photoisomerization. This report paves the way for future investigations in the optimization of biomimetic photoswitches potentially useful in modern light-induced therapeutic strategies.

Modeling the Enzymatic Mechanism of the SARS-CoV-2 RNA-Dependent RNA Polymerase by DFT/MM-MD: An Unusual Active Site Leading to High Replication Rates.

Viral infection relies on the hijacking of cellular machineries to enforce the reproduction of the infecting virus and its subsequent diffusion. In this context, the replication of the viral genome is a key step performed by specific enzymes, i.e., polymerases. The replication of SARS-CoV-2, the causative agent of the COVID-19 pandemics, is based on the duplication of its RNA genome, an action performed by the viral RNA-dependent RNA polymerase. In this contribution, by using highly demanding DFT/MM-MD computations coupled to 2D-umbrella sampling techniques, we have determined the chemical mechanisms leading to the inclusion of a nucleotide in the nascent viral RNA strand. These results highlight the high efficiency of the polymerase, which lowers the activation free energy to less than 10 kcal/mol. Furthermore, the SARS-CoV-2 polymerase active site is slightly different from those usually found in other similar enzymes, and in particular, it lacks the possibility to enforce a proton shuttle via a nearby histidine. Our simulations show that this absence is partially compensated by lysine whose proton assists the reaction, opening up an alternative, but highly efficient, reactive channel. Our results present the first mechanistic resolution of SARS-CoV-2 genome replication at the DFT/MM-MD level and shed light on its unusual enzymatic reactivity paving the way for the future rational design of antivirals targeting emerging RNA viruses.