Modeling solvent effects and convergence of 31P-NMR shielding calculations with COBRAMM.

Solvent effects on 31P-NMR parameters for triphenylphosphine oxide and triphenylphosphine in chloroform have been extensively investigated by testing different solvation models. The solvent is described implicitly, mixed implicitly/explicitly, and using full explicit models. Polarizable continuum model (PCM), molecular dynamic (MD) simulations, and hybrid quantum mechanics/molecular mechanics (QM/MM) calculations are used to disclose the effects of solute/solvent interactions and, more generally, the role of the embedding in NMR simulations. The results show the beneficial effect of carrying out QM/MM optimizations on top of geometries directly extracted from classical MD simulations, used to ensure representative conformational sampling. The nuclear shielding convergence has been tested against a different number of snapshots and with the inclusion of solvent shells into the QM region. An automated MD//QM/MM//GIAO protocol, implemented in the COBRAMM package, is here proposed and tested on trimethyl phosphite showing that our approach boosts the convergence of nuclear shielding satisfactorily. The present work aims to be a stepping-stone to assess proper QM/MM computational strategies in simulating chemical shifts in non-homogeneous systems like supramolecular and biological systems.

Characterizing Counterion-Dependent Aggregation of Rhodamine B by Classical Molecular Dynamics Simulations.

The aggregation in a solution of charged dyes such as Rhodamine B (RB) is significantly affected by the type of counterion, which can determine the self-assembled structure that in turn modulates the optical properties. RB aggregation can be boosted by hydrophobic and bulky fluorinated tetraphenylborate counterions, such as F5TPB, with the formation of nanoparticles whose fluorescence quantum yield (FQY) is affected by the degree of fluorination. Here, we developed a classical force field (FF) based on the standard generalized Amber parameters that allows modeling the self-assembling process of RB/F5TPB systems in water, consistent with experimental evidence. Namely, the classical MD simulations employing the re-parametrized FF reproduce the formation of nanoparticles in the RB/F5TPB system, while in the presence of iodide counterions, only RB dimeric species can be formed. Within the large, self-assembled RB/F5TPB aggregates, the occurrence of an H-type RB-RB dimer can be observed, a species that is expected to quench RB fluorescence, in agreement with the experimental data of FQY. The outcome provides atomistic details on the role of the bulky F5TPB counterion as a spacer, with the developed classical FF representing a step towards reliable modeling of dye aggregation in RB-based materials.

Distinct allosteric pathways in imidazole glycerol phosphate synthase from yeast and bacteria.

Understanding the relationship between protein structures and their function is still an open question that becomes very challenging when allostery plays an important functional role. Allosteric proteins, in fact, exploit different ranges of motions (from sidechain local fluctuations to long-range collective motions) to effectively couple distant binding sites, and of particular interest is whether allosteric proteins of the same families with similar functions and structures also necessarily share the same allosteric mechanisms. Here, we compared the early dynamics initiating the allosteric communication of a prototypical allosteric enzyme from two different organisms, i.e., the imidazole glycerol phosphate synthase (IGPS) enzymes from the thermophilic bacteria and the yeast, working at high and room temperatures, respectively. By combining molecular dynamics simulations and network models derived from graph theory, we found rather distinct early allosteric dynamics in the IGPS from the two organisms, involving significatively different allosteric pathways in terms of both local and collective motions. Given the successful prediction of key allosteric residues in the bacterial IGPS, whose mutation disrupts its allosteric communication, the outcome of this study paves the way for future experimental studies on the yeast IGPS that could foster therapeutic applications by exploiting the control of IGPS enzyme allostery.

Singular Interface Dynamics of the SARS-CoV-2 Delta Variant Explained with Contact Perturbation Analysis.

Emerging SARS-CoV-2 variants raise concerns about our ability to withstand the Covid-19 pandemic, and therefore, understanding mechanistic differences of those variants is crucial. In this study, we investigate disparities between the SARS-CoV-2 wild type and five variants that emerged in late 2020, focusing on the structure and dynamics of the spike protein interface with the human angiotensin-converting enzyme 2 (ACE2) receptor, by using crystallographic structures and extended analysis of microsecond molecular dynamics simulations. Dihedral angle principal component analysis (PCA) showed the strong similarities in the spike receptor binding domain (RBD) dynamics of the Alpha, Beta, Gamma, and Delta variants, in contrast with those of WT and Epsilon. Dynamical perturbation networks and contact PCA identified the peculiar interface dynamics of the Delta variant, which cannot be directly imputable to its specific L452R and T478K mutations since those residues are not in direct contact with the human ACE2 receptor. Our outcome shows that in the Delta variant the L452R and T478K mutations act synergistically on neighboring residues to provoke drastic changes in the spike/ACE2 interface; thus a singular mechanism of action eventually explains why it dominated over preceding variants.

Modelling quenching mechanisms of disordered molecular systems in the presence of molecular aggregates.

Exciton density dynamics recorded in time-resolved spectroscopic measurements is a useful tool to recover information on energy transfer (ET) processes that can occur at different timescales, up to the ultrafast regime. Macroscopic models of exciton density decays, involving both direct Förster-like ET and diffusion mechanisms for exciton-exciton annihilation, are largely used to fit time-resolved experimental data but generally neglect contributions from molecular aggregates that can work as quenching species. In this work, we introduce a macroscopic model that includes contributions from molecular aggregate quenchers in a disordered molecular system. As an exemplifying case, we considered a homogenous distribution of rhodamine B dyes embedded in organic nanoparticles to set the initial parameters of the proposed model. The influence of such model parameters is systematically analysed, showing that the presence of molecular aggregate quenchers can be monitored by evaluating the exciton density long time decays. We showed that the proposed model can be applied to molecular systems with ultrafast decays, and we anticipated that it could be used in future studies for global fitting of experimental data with potential support from first-principles simulations.

Control of Protonated Schiff Base Excited State Decay within Visual Protein Mimics: A Unified Model for Retinal Chromophores.

Artificial biomimetic chromophore-protein complexes inspired by natural visual pigments can feature color tunability across the full visible spectrum. However, control of excited state dynamics of the retinal chromophore, which is of paramount importance for technological applications, is lacking due to its complex and subtle photophysics/photochemistry. Here, ultrafast transient absorption spectroscopy and quantum mechanics/molecular mechanics simulations are combined for the study of highly tunable rhodopsin mimics, as compared to retinal chromophores in solution. Conical intersections and transient fluorescent intermediates are identified with atomistic resolution, providing unambiguous assignment of their ultrafast excited state absorption features. The results point out that the electrostatic environment of the chromophore, modified by protein point mutations, affects its excited state properties allowing control of its photophysics with same power of chemical modifications of the chromophore. The complex nature of such fine control is a fundamental knowledge for the design of bio-mimetic opto-electronic and photonic devices.

Heterometallic rhodium clusters as electron reservoirs: Chemical, electrochemical, and theoretical studies of the centered-icosahedral [Rh12E(CO)27]n- atomically precise carbonyl compounds.

In this paper, we present a comparative study of the redox properties of the icosahedral [Rh12E(CO)27]n- (n = 4 when E = Ge or Sn and n = 3 when E = Sb or Bi) family of clusters through in situ infrared spectroelectrochemistry experiments and density functional theory computational studies. These clusters show shared characteristics in terms of molecular structure, being all E-centered icosahedral species, and electron counting, possessing 170 valence electrons as predicted by the electron-counting rules, based on the cluster-borane analogy, for compounds with such metal geometry. However, in some cases, clusters of similar nuclearity, and beyond, may show multivalence behavior and may be stable with a different electron counting, at least on the time scale of the electrochemical analyses. The experimental results, confirmed by theoretical calculations, showed a remarkable electron-sponge behavior for [Rh12Ge(CO)27]4- (1), [Rh12Sb(CO)27]3- (3), and [Rh12Bi(CO)27]3- (4), with a cluster charge going from -2 to -6 for 1 and 3 and from -2 to -7 for cluster 4, making them examples of molecular electron reservoirs. The [Rh12Sn(CO)27]4- (2) derivative, conversely, presents a limited ability to exist in separable reduced cluster species, at least within the experimental conditions, while in the gas phase it appears to be stable both as a penta- and hexa-anion, therefore showing a similar redox activity as its congeners. As a fallout of those studies, during the preparation of [Rh12Sb(CO)27]3-, we were able to isolate a new species, namely, [Rh11Sb(CO)26]2-, which presents a Sb-centered nido-icosahedral metal structure possessing 158 cluster valence electrons, in perfect agreement with the polyhedral skeletal electron pair theory.

Eigenvector centrality for characterization of protein allosteric pathways.

Determining the principal energy-transfer pathways responsible for allosteric communication in biomolecules remains challenging, partially due to the intrinsic complexity of the systems and the lack of effective characterization methods. In this work, we introduce the eigenvector centrality metric based on mutual information to elucidate allosteric mechanisms that regulate enzymatic activity. Moreover, we propose a strategy to characterize the range of correlations that underlie the allosteric processes. We use the V-type allosteric enzyme imidazole glycerol phosphate synthase (IGPS) to test the proposed methodology. The eigenvector centrality method identifies key amino acid residues of IGPS with high susceptibility to effector binding. The findings are validated by solution NMR measurements yielding important biological insights, including direct experimental evidence for interdomain motion, the central role played by helix h[Formula: see text], and the short-range nature of correlations responsible for the allosteric mechanism. Beyond insights on IGPS allosteric pathways and the nature of residues that could be targeted by therapeutic drugs or site-directed mutagenesis, the reported findings demonstrate the eigenvector centrality analysis as a general cost-effective methodology to gain fundamental understanding of allosteric mechanisms at the molecular level.

Polymerization Isomerism in Co-M (M = Cu, Ag, Au) Carbonyl Clusters: Synthesis, Structures and Computational Investigation.

The reaction of [Co(CO)4]- (1) with M(I) compounds (M = Cu, Ag, Au) was reinvestigated unraveling an unprecedented case of polymerization isomerism. Thus, as previously reported, the trinuclear clusters [M{Co(CO)4}2]- (M = Cu, 2; Ag, 3; Au, 4) were obtained by reacting 1 with M(I) in a 2:1 molar ratio. Their molecular structures were corroborated by single-crystal X-ray diffraction (SC-XRD) on isomorphous [NEt4][M{Co(CO)4}2] salts. [NEt4](3)represented the first structural characterization of 3. More interestingly, changing the crystallization conditions of solutions of 3, the hexanuclear cluster [Ag2{Co(CO)4}4]2- (5) was obtained in the solid state instead of 3. Its molecular structure was determined by SC-XRD as Na2(5)·C4H6O2, [PPN]2(5)·C5H12 (PPN = N(PPh3)2]+), [NBu4]2(5) and [NMe4]2(5) salts. 5 may be viewed as a dimer of 3 and, thus, it represents a rare case of polymerization isomerism (that is, two compounds having the same elemental composition but different molecular weights) in cluster chemistry. The phenomenon was further studied in solution by IR and ESI-MS measurements and theoretically investigated by computational methods. Both experimental evidence and density functional theory (DFT) calculations clearly pointed out that the dimerization process occurs in the solid state only in the case of Ag, whereas Cu and Au related species exist only as monomers.

Key role of the REC lobe during CRISPR-Cas9 activation by 'sensing', 'regulating', and 'locking' the catalytic HNH domain.

Understanding the conformational dynamics of CRISPR (clustered regularly interspaced short palindromic repeat)-Cas9 is of the utmost importance for improving its genome editing capability. Here, molecular dynamics simulations performed using Anton-2 - a specialized supercomputer capturing micro-to-millisecond biophysical events in real time and at atomic-level resolution - reveal the activation process of the endonuclease Cas9 toward DNA cleavage. Over the unbiased simulation, we observe that the spontaneous approach of the catalytic domain HNH to the DNA cleavage site is accompanied by a remarkable structural remodeling of the recognition (REC) lobe, which exerts a key role for DNA cleavage. Specifically, the significant conformational changes and the collective conformational dynamics of the REC lobe indicate a mechanism by which the REC1-3 regions 'sense' nucleic acids, 'regulate' the HNH conformational transition, and ultimately 'lock' the HNH domain at the cleavage site, contributing to its catalytic competence. By integrating additional independent simulations and existing experimental data, we provide a solid validation of the activated HNH conformation, which had been so far poorly characterized, and we deliver a comprehensive understanding of the role of REC1-3 in the activation process. Considering the importance of the REC lobe in the specificity of Cas9, this study poses the basis for fully understanding how the REC components control the cleavage of off-target sequences, laying the foundation for future engineering efforts toward improved genome editing.

Decrypting the Information Exchange Pathways across the Spliceosome Machinery.

Intron splicing of a nascent mRNA transcript by spliceosome (SPL) is a hallmark of gene regulation in eukaryotes. SPL is a majestic molecular machine composed of an entangled network of proteins and RNAs that meticulously promotes intron splicing through the formation of eight intermediate complexes. Cross-communication among the critical distal proteins of the SPL assembly is pivotal for fast and accurate directing of the compositional and conformational readjustments necessary to achieve high splicing fidelity. Here, molecular dynamics (MD) simulations of an 800 000 atom model of SPL C complex from yeast Saccharomyces cerevisiae and community network analysis enabled us to decrypt the complexity of this huge molecular machine, by identifying the key channels of information transfer across long distances separating key protein components. The reported study represents an unprecedented attempt in dissecting cross-communication pathways within one of the most complex machines of eukaryotic cells, supporting the critical role of Clf1 and Cwc2 splicing cofactors and specific domains of the Prp8 protein as signal conveyors for pre-mRNA maturation. Our findings provide fundamental advances into mechanistic aspects of SPL, providing a conceptual basis for controlling the SPL via small-molecule modulators able to tackle splicing-associated diseases by altering/obstructing information-exchange paths.

Sodium diffusion in ionic liquid-based electrolytes for Na-ion batteries: the effect of polarizable force fields.

Understanding the transport of sodium ions in ionic liquids is key to designing novel electrolyte materials for sodium-ion batteries. In this work, we combine molecular dynamics simulation and experiments to study how molecular interactions and local ordering affect relevant physico-chemical properties. Ionic transport and local solvation environments are investigated in electrolytes composed of sodium bis(fluorosulfonyl)imide, (Na[FSI]), in N,N-methylpropylpyrrolidinium bis(fluorosulfonyl)imide, [C3C1pyr][FSI], at different salt concentrations. The electrolyte systems are modelled by means of molecular dynamic simulations using a polarizable force field. We show that including polarization effects explicitly in the molecular simulations is required in order to attain a reliable description of the transport properties of sodium in the [C3C1pyr][FSI] electrolyte. The validation of the computational results upon comparison with experimental data allows us to assess the suitability of polarizable force fields in describing and interpreting the structure and dynamics of the sodium salt-ionic liquid system, which is essential to enable the application of IL-based electrolytes in novel energy-storage technologies.

First-principles characterization of the singlet excited state manifold in DNA/RNA nucleobases.

An extensive theoretical characterization of the singlet excited state manifold of the five canonical DNA/RNA nucleobases (thymine, cytosine, uracil, adenine and guanine) in gas-phase is carried out with time-dependent density functional theory (TD-DFT) and restricted active space second-order perturbation theory (RASPT2) approaches. Both ground state and excited state absorptions are analyzed and compared between these different theoretical approaches, assessing the performance of the hybrid B3LYP and CAM-B3LYP (long-range corrected) functionals with respect to the RASPT2 reference. By comparing the TD-DFT estimates with our reference for high-lying excited states, we are able to narrow down specific energetic windows where TD-DFT may be safely employed to qualitatively reproduce the excited state absorption (ESA) signals registered in non-linear and time-resolved spectroscopy for monitoring photoinduced phenomena. Our results show a qualitative agreement between the RASPT2 reference and the B3LYP computed ESAs of pyrimidines in the near-IR/Visible spectral probing window while for purines the agreement is limited to the near-IR ESAs, with generally larger discrepancies obtained with the CAM-B3LYP functional. This outcome paves the way for appropriate application of cost-effective TD-DFT approaches to simulate linear and non-linear spectroscopies of realistic multichromophoric DNA/RNA systems with biological and nanotechnological relevance.

Modeling multidimensional spectral lineshapes from first principles: application to water-solvated adenine.

In this discussion we present a methodology to describe spectral lineshape from first principles, providing insight into the solvent-solute molecular interactions in terms of static and dynamic disorder and how these shape the signals recorded experimentally in linear and nonlinear optical spectroscopies, including two-dimensional electronic spectroscopy (2DES). Two different strategies for simulating the lineshape are compared: both rely on the same evaluation of the coupling between the electronic states and the intra-molecular vibrations, while they differ in describing the influence exerted by the diverse water configurations attained along a molecular dynamics (MD) simulation. The first method accounts for such water arrangements as first order perturbations on the adenine energies computed for a single reference (gas phase) quantum calculation. The second method requires computation of the manifold of excited states explicitly at each simulation snapshot, employing a hybrid quantum mechanics/molecular mechanics (QM/MM) scheme. Both approaches are applied to a large number of states of the adenine singlet excited manifold (chosen because of its biological role), and compared with available experimental data. They give comparable results but the first approach is two orders of magnitude faster. We show how the various contributions (static/dynamic disorder, intra-/inter-molecular interactions) sum up to build the total broadening observed in experiments.

Ultrafast Carotenoid to Retinal Energy Transfer in Xanthorhodopsin Revealed by the Combination of Transient Absorption and Two-Dimensional Electronic Spectroscopy.

By comparing two-dimensional electronic spectroscopy (2DES) and Pump-Probe (PP) measurements on xanthorhodopsin (XR) and reduced-xanthorhodopsin (RXR) complexes, the ultrafast carotenoid-to-retinal energy transfer pathway is revealed, at very early times, by an excess of signal amplitude at the associated cross-peak and by the carotenoid bleaching reduction due to its ground state recovery. The combination of the measured 2DES and PP spectroscopic data with theoretical modelling allows a clear identification of the main experimental signals and a comprehensive interpretation of their origin and dynamics. The remarkable velocity of the energy transfer, despite the non-negligible energy separation between the two chromophores, and the analysis of the underlying transport mechanism, highlight the role played by the ground state carotenoid vibrations in assisting the process.

Two-dimensional electronic spectroscopy as a tool for tracking molecular conformations in DNA/RNA aggregates.

A computational strategy to simulate two-dimensional electronic spectra (2DES) is introduced, which allows us to analyse ground state dynamics and to sample and measure different conformations attained by flexible molecular systems in solution. An explicit mixed quantum mechanics/molecular mechanics (QM/MM) approach is employed for the evaluation of the necessary electronic excited state energies and transition dipole moments. The method is applied towards a study of the highly flexible water-solvated adenine-adenine monophosphate (ApA), a system featuring two interacting adenine moieties that display various intermolecular arrangements, known to deeply affect their photochemical outcome. Molecular dynamics simulations and cluster analysis have been used to select the molecular conformations, reducing the complexity of the flexible ApA conformational space. By using our sum-over-states (SOS) approach to obtain the 2DES spectra for each of these selected conformations, we can discern spectral changes and relate them to specific nuclear arrangements: close lying π-stacked bases exhibit a splitting of their respective 1La signal traces; T-stacked bases exhibit the appearance of charge transfer states in the low-energy Vis probing window while displaying no 1La splitting, being particularly favoured when promoting amino to 5-ring interactions; unstacked and distant adenine moieties exhibit signals similar to those of the adenine monomer, as is expected for non-interacting nucleobases. 2DES maps reveal the spectral fingerprints associated with specific molecular conformations, and are thus a promising option to enable their quantitative spectroscopic detection beyond standard 1D pump-probe techniques. This is expected to aid the understanding of how nucleobase aggregation controls and modulates the photostability and photo-damage of extended DNA/RNA systems.

The effect of solvent relaxation in the ultrafast time-resolved spectroscopy of solvated benzophenone.

Benzophenone (BP) despite its relatively simple molecular structure is a paradigmatic sensitizer, featuring both photocatalytic and photobiological effects due to its rather complex photophysical properties. In this contribution we report an original theoretical approach to model realistic, ultra-fast spectroscopy data, which requires describing intra- and intermolecular energy and structural relaxation. In particular we explicitly simulate time-resolved pump-probe spectra using a combination of state-of-the art hybrid quantum mechanics/molecular mechanics dynamics to treat relaxation and vibrational effects. The comparison with experimental transient absorption data demonstrates the efficiency and accuracy of our approach. Furthermore the explicit inclusion of the solvent, water for simulation and methanol for experiment, allows us, despite the inherent different behavior of the two, to underline the role played by the H-bonding relaxation in the first hundreds of femtoseconds after optical excitation. Finally we predict for the first time the two-dimensional electronic spectrum (2DES) of BP taking into account the vibrational effects and hence modelling partially symmetric and asymmetric ultrafast broadening.

UV-induced long-lived decays in solvated pyrimidine nucleosides resolved at the MS-CASPT2/MM level.

The most relevant 'dark' electronic excited states in DNA/RNA pyrimidine nucleosides are mapped in water employing hybrid MS-CASPT2/MM optimisations with explicit solvation and including the sugar. Conical intersections (CIs) between initially accessed bright 1ππ* and the lowest energy dark 1nπ* excited states, involving the lone pair localised on the oxygen and/or nitrogen atoms are characterised. They are found in the vicinities of the Franck-Condon (FC) region and are shown to facilitate non-adiabatic population transfer. The excited state population of the 1nOπ* state, localised in the carbonyl moiety on all pyrimidine nucleosides, is predicted to rapidly evolve to its minimum, displaying non-negligible potential energy barriers along its non-radiative decay, and accounting for the ps signal registered in pump-probe experiments as well as for an efficient population of the triplet state. Cytidine displays an additional 1nNπ* state localised in the N3 atom and that leads to its excited state minimum displaying large potential energy barriers in the pathway connecting to the CI with the ground state. Sugar-to-base hydrogen/proton transfer processes are assessed in solution for the first time, displaying a sizable barrier along its decay and thus being competitive with other slow decay channels in the ps and ns timescales. A unified deactivation scheme for the long-lived channels of pyrimidine nucleosides is delivered, where the 1nOπ* state is found to mediate the long-lived decay in the singlet manifold and act as the doorway for triplet population and thus accounting for the recorded phosphorescence and, more generally, for the transient/photoelectron spectral signals registered up to the ns timescale.

Characterization of Protein Tyrosine Phosphatase 1B Inhibition by Chlorogenic Acid and Cichoric Acid.

Protein tyrosine phosphatase 1B (PTP1B) is a known regulator of the insulin and leptin signaling pathways and is an active target for the design of inhibitors for the treatment of type II diabetes and obesity. Recently, cichoric acid (CHA) and chlorogenic acid (CGA) were predicted by docking methods to be allosteric inhibitors that bind distal to the active site. However, using a combination of steady-state inhibition kinetics, solution nuclear magnetic resonance experiments, and molecular dynamics simulations, we show that CHA is a competitive inhibitor that binds in the active site of PTP1B. CGA, while a noncompetitive inhibitor, binds in the second aryl phosphate binding site, rather than the predicted benzfuran binding pocket. The molecular dynamics simulations of the apo enzyme and cysteine-phosphoryl intermediate states with and without bound CGA suggest CGA binding inhibits PTP1B by altering hydrogen bonding patterns at the active site. This study provides a mechanistic understanding of the allosteric inhibition of PTP1B.

Protospacer Adjacent Motif-Induced Allostery Activates CRISPR-Cas9.

CRISPR-Cas9 is a genome editing technology with major impact in life sciences. In this system, the endonuclease Cas9 generates double strand breaks in DNA upon RNA-guided recognition of a complementary DNA sequence, which strictly requires the presence of a protospacer adjacent motif (PAM) next to the target site. Although PAM recognition is essential for cleavage, it is unknown whether and how PAM binding activates Cas9 for DNA cleavage at spatially distant sites. Here, we find evidence of a PAM-induced allosteric mechanism revealed by microsecond molecular dynamics simulations. PAM acts as an allosteric effector and triggers the interdependent conformational dynamics of the Cas9 catalytic domains (HNH and RuvC), responsible for concerted cleavage of the two DNA strands. Targeting such an allosteric mechanism should enable control of CRISPR-Cas9 functionality.