Unrevealing the Nitrogen Elusive Chirality of 3-Sulfanyl and 3-Sulfinyl N-Tosyl Isoindolinones by ECD Spectra: An Experimental and Theoretical Investigation.

The nitrogen-hybridization/pyramidalization of two solvated N-tosylisoindolinone derivatives having chiral residues in adjacent (I) or adjacent and distal (II) position has been investigated by a theoretical-computational procedure based on Molecular Dynamics simulations and Quantum-Chemical calculations. After validation of our methodology in providing a reliable repertory of conformations by modeling the electronic circular dichroism (EDC) spectra, the electronic features associated with N-pyramidalization were further characterized through Natural Bond Order (NBO) analysis. Comparing against the N-geometry observed in crystal structures as a reference, our findings reveal that the presence of neighbouring chiral centers induces a more pronounced N-pyramidalization in solution than in the solid state, both in I and II. Furthermore, NBO analysis confirms that the N-lactam mostly retains the sp2 character but exhibits slight configurational distortion (ξI=13°; ξII=21°), which significantly influences the chiroptical activities observed in ECD spectra of I and II. This substantiates the N-lactams as configurationally stable chiral centers.

On the Statistical Regime, Coherence versus Incoherence and Ergodicity of Quantum Vibrational Trajectories in Soft Condensed Molecular Systems.

A theoretical-computational procedure, recently proposed for modelling Vibrational Energy Relaxation (VER) processes of a molecule (Quantum Center, QC) embedded in a complex atomic-molecular system, is extended and applied for analyzing in detail the features of the QC density matrix (DM) temporal evolution. The results, obtained using aqueous azide ion as a case study, show the total lack of coherence in the DM, when the system is prepared to be initially in a pure vibrational eigenstate. This finding is fully in line with the statistical interpretation of the process typically adopted also in the experimental studies where the relaxation processes are all described within the typical schemes of chemical kinetics. Consistently, when the initial vibrational state corresponds to an eigenstate mixture, although initially coherent, the DM relaxes to a fully incoherent condition with a mean lifetime related to the one of the diagonal elements relaxation. These specific DM features turn out to be essentially governed by the thermal equilibrium condition of the atomic-molecular classical coordinates which drive the ensemble of the quantum-trajectories toward the observed statistical regime. Finally, from the analysis of a single long timescale quantum vibrational trajectory it also clearly emerges its ergodic behaviour.

Direct Regioselective Hydro(hetero)arylation/Cyclocondensation Reactions of ß-(2-Aminophenyl)-α,ß-ynones by Means of Transition-Metal Catalysis/Brønsted Acid Synergism: Experimental Results and Computational Insights.

Experimental results and computational insights explain the key role of transition-metal catalysis/Brønsted acid synergism in the achievement of the sequential regioselective direct heteroarylation/cyclocondensation reactions of ß-(2-aminophenyl)-α,ß-ynones with a variety of electron-rich aromatic heterocyclic/arenes to afford quinoline-(hetero)aromatic hybrids. The first approach to the synthesis of 4-(1H-pyrrol-2-yl)quinolines is described. The effectiveness of various transition metals is compared.

Plasmon-enhanced circular dichroism spectroscopy of chiral drug solutions.

We investigate the potential of surface plasmon polaritons at noble metal interfaces for surface-enhanced chiroptical sensing of dilute chiral drug solutions with nl volume. The high quality factor of surface plasmon resonances in both Otto and Kretschmann configurations enables the enhancement of circular dichroism differenatial absorption thanks to the large near-field intensity of such plasmonic excitations. Furthermore, the subwavelength confinement of surface plasmon polaritons is key to attain chiroptical sensitivity to small amounts of drug volumes placed around ≃100 nm by the metal surface. Our calculations focus on reparixin, a pharmaceutical molecule currently used in clinical studies for patients with community-acquired pneumonia, including COVID-19 and acute respiratory distress syndrome. Considering realistic dilute solutions of reparixin dissolved in water with concentration ≤5 mg/ml and nl volume, we find a circular-dichroism differential absorption enhancement factor of the order ≃20 and chirality-induced polarization distortion upon surface plasmon polariton excitation. Our results are relevant for the development of innovative chiroptical sensors capable of measuring the enantiomeric imbalance of chiral drug solutions with nl volume.

Theoretical-Computational Modeling of CD Spectra of Aqueous Monosaccharides by Means of Molecular Dynamics Simulations and Perturbed Matrix Method.

The electronic circular dichroism (ECD) spectra of aqueous d-glucose and d-galactose were modeled using a theoretical-computational approach combining molecular dynamics (MD) simulations and perturbed matrix method (PMM) calculations, hereafter termed MD-PMM. The experimental spectra were reproduced with a satisfactory accuracy, confirming the good performances of MD-PMM in modeling different spectral features in complex atomic-molecular systems, as already reported in previous studies. The underlying strategy of the method was to perform a preliminary long timescale MD simulation of the chromophore followed by the extraction of the relevant conformations through essential dynamics analysis. On this (limited) number of relevant conformations, the ECD spectrum was calculated via the PMM approach. This study showed that MD-PMM was able to reproduce the essential features of the ECD spectrum (i.e., the position, the intensity, and the shape of the bands) of d-glucose and d-galactose while avoiding the rather computationally expensive aspects, which were demonstrated to be important for the final outcome, such as (i) the use of a large number of chromophore conformations; (ii) the inclusion of quantum vibronic coupling; and (iii) the inclusion of explicit solvent molecules interacting with the chromophore atoms within the chromophore itself (e.g., via hydrogen bonds).

Alternative Fast and Slow Primary Charge-Separation Pathways in Photosystem II.

Photosystem-II (PSII) is a multi-subunit protein complex that harvests sunlight to perform oxygenic photosynthesis. Initial light-activated charge separation takes place at a reaction centre consisting of four chlorophylls and two pheophytins. Understanding the processes following light excitation remains elusive due to spectral congestion, the ultrafast nature, and multi-component behaviour of the charge-separation process. Here, using advanced computational multiscale approaches which take into account the large-scale configurational flexibility of the system, we identify two possible primary pathways to radical-pair formation that differ by three orders of magnitude in their kinetics. The fast (short-range) pathway is dominant, but the existence of an alternative slow (long-range) charge-separation pathway hints at the evolution of redundancy that may serve other purposes, adaptive or protective, related to formation of the unique oxidative species that drives water oxidation in PSII.

A computational approach for modeling electronic circular dichroism of solvated chromophores.

The present study consists in a novel computational protocol to model the UV-circular dichroism spectra of solvated species. It makes use of quantum-chemical calculations on a series of conformations of a flexible chromophore or on a series of chromophore/solvent clusters extracted from molecular dynamic simulations. The protocol is described and applied to the aqueous cationic tripeptide GAG+ and to the aqueous neutral decapeptide (GVGVP)2 . The protocol has proven able to: (i) properly consider the conformational motion of solute in the given environment; (ii) give the actual statistical weight of each conformational state; (iii) provide a reliable quantum mechanical method able to reproduce the spectral features. Temperature effects on conformations and spectral properties are properly taken into account. The role of explicit solvent on the conformational analysis and the spectra calculation is discussed. The comparison of the calculated circular dichroism spectra with experimental ones recorded at different temperatures represents a strict validation test of the method.

Correction: Theoretical and experimental study on the O(3P) + 2,5-dimethylfuran reaction in the gas phase.

Correction for 'Theoretical and experimental study on the O(3P) + 2,5-dimethylfuran reaction in the gas phase' by Andrea Giustini et al., Phys. Chem. Chem. Phys., 2021, 23, 19424-19434, DOI: 10.1039/D1CP01724A.

A Simplified Treatment for Efficiently Modeling the Spectral Signal of Vibronic Transitions: Application to Aqueous Indole.

In this paper, we introduce specific approximations to simplify the vibronic treatment in modeling absorption and emission spectra, allowing us to include a huge number of vibronic transitions in the calculations. Implementation of such a simplified vibronic treatment within our general approach for modelling vibronic spectra, based on molecular dynamics simulations and the perturbed matrix method, provided a quantitative reproduction of the absorption and emission spectra of aqueous indole with higher accuracy than the one obtained when using the existing vibronic treatment. Such results, showing the reliability of the approximations employed, indicate that the proposed method can be a very efficient and accurate tool for computational spectroscopy.

Electrochemical-Induced Cascade Reaction of 2-Formyl Benzonitrile with Anilines: Synthesis of N-Aryl Isoindolinones.

An electrochemical initiated tandem reaction of anilines with 2-formyl benzonitrile has been developed. Thus, unprecedented 3-N-aryl substituted isoindolinones have been conveniently achieved by constant current electrolysis in a divided cell using catalytic amount of electricity and supporting electrolyte and a Pt-cathode as working electrode. The origin of the electrochemical activation as well as the mechanism of the subsequent chemical cascade reactions have been investigated by DFT calculations.

Theoretical-Computational Modeling of Gas-State Thermodynamics in Flexible Molecular Systems: Ionic Liquids in the Gas Phase as a Case Study.

A theoretical-computational procedure based on the quasi-Gaussian entropy (QGE) theory and molecular dynamics (MD) simulations is proposed for the calculation of thermodynamic properties for molecular and supra-molecular species in the gas phase. The peculiarity of the methodology reported in this study is its ability to construct an analytical model of all the most relevant thermodynamic properties, even within a wide temperature range, based on a practically automatic sampling of the entire conformational repertoire of highly flexible systems, thereby bypassing the need for an explicit search for all possible conformers/rotamers deemed relevant. In this respect, the reliability of the presented method mainly depends on the quality of the force field used in the MD simulations and on the ability to discriminate in a physically coherent way between semi-classical and quantum degrees of freedom. The method was tested on six model systems (n-butane, n-butane, n-octanol, octadecane, 1-butyl-3-methylimidazolium hexafluorophosphate and 1-butyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide ionic pairs), which, being experimentally characterized and already addressed by other theoretical-computational methods, were considered as particularly suitable to allow us to evaluate the method's accuracy and efficiency, bringing out advantages and possible drawbacks. The results demonstrate that such a physically coherent yet relatively simple method can represent a further valid computational tool that is alternative and complementary to other extremely efficient computational methods, as it is particularly suited for addressing the thermodynamics of gaseous systems with a high conformational complexity over a large range of temperature.

Theoretical and experimental study on the O(3P) + 2,5-dimethylfuran reaction in the gas phase.

In this work we report a joint experimental and computational study on the 2,5-dimethylfuran oxidation reaction in the gas phase initiated by atomic oxygen O(3P). The experiments have been performed by using vacuum-ultraviolet synchrotron radiation at the Advanced Light Source (ALS) of the Lawrence Berkeley National Laboratory (LBNL), at a temperature of 550 K and a pressure of 8 Torr. The experimental data were supported by quantum-chemical calculations along with a kinetic model, also taking into account the possible involvement of different magnetic states, performed in the framework of the RRKM theory. Propyne, acetaldehyde, methylglyoxal, dimethylglyoxal, 3-penten-2-one, 2,5-dimethylfuran-3(2H)-one, and 1,2-diacetyl ethylene have been identified as the main primary products arising under the conditions of the experiment. Our computational model suggests that these species can be formed at the concentration and branching ratio experimentally observed only in the presence of a non-negligible fraction of non-thermalized intermediates.

A Combined Experimental and Computational Study of Halogen and Hydrogen Bonding in Molecular Salts of 5-Bromocytosine.

Although natural or artificial modified pyrimidine nucleobases represent important molecules with valuable properties as constituents of DNA and RNA, no systematic analyses of the structural aspects of bromo derivatives of cytosine have appeared so far in the literature. In view of the biochemical and pharmaceutical relevance of these compounds, six different crystals containing proton-transfer derivatives of 5-bromocytosine are prepared and analyzed in the solid-state by single crystal X-ray diffraction. All six compounds are organic salts, with proton transfer occurring to the Nimino atom of the pyridine ring. Experimental results are then complemented with Hirshfeld surface analysis to quantitively evaluate the contribution of different intermolecular interactions in the crystal packing. Furthermore, theoretical calculations, based on different arrangements of molecules extracted from the crystal structure determinations, are carried out to analyze the formation mechanism of halogen bonds (XBs) in these compounds and provide insights into the nature and strength of the observed interactions. The results show that the supramolecular architectures of the six molecular salts involve extensive classical intermolecular hydrogen bonds. However, in all but one proton-transfer adducts, weak to moderate XBs are revealed by C-Br…O short contacts between the bromine atom in the fifth position, which acts as XB donor (electron acceptor). Moreover, the lone pair electrons of the oxygen atom of adjacent pyrimidine nucleobases and/or counterions or water molecules, which acts as XB acceptor (electron donor).

Experimental and computational evidence on gold-catalyzed regioselective hydration of phthalimido-protected propargylamines: an entry to ß-amino ketones.

The results of our investigations on the Au-catalyzed regioselective hydration reaction of both alkyl- and aryl-substituted N-propargyl phthalimides directed to the selective formation of the corresponding ß-phthalimido ketones are described. Experimental data, in particular the observed regioselectivity, have been qualitatively supported by quantum-chemical calculations carried out on model systems in the framework of Density Functional Theory (DFT) followed by quantum theory of atoms in molecules (QTAIMS). Our results suggest that the electronic features of the initial adduct between the propargyl triple bond and the Au(i) catalyst, in particular the character of the gold-triple bond interaction, are essential for the observed regioselectivity. Other effects, such as the presence of the solvent and the formation of a H-bond between the water molecule and the phthalimido moiety, although apparently irrelevant for the regioselectivity, have proven to be kinetically and catalytically rather important.

Experimental results and computational insight into sequential reactions of ß-(2-aminophenyl)-α,ß-ynones with aryl isocyanates/benzoyl isothiocyanate.

The investigation on tandem addition/cyclization reactions of ß-(2-aminophenyl)-α,ß-ynones with aryl isocyanates/benzoyl isothiocyanate is reported. Experimental results show the suitable conditions to selectively direct the reaction outcome towards the product of 6-exo-dig N-, O-, or S-annulation of the in situ generated alkynyl urea/thiourea intermediate. The reaction of a variety of ß-(2-aminophenyl)-α,ß-ynones with aryl isocyanates/benzoyl isothiocyanate led to the selective formation of quinazoline or benzoxazine/benzothiazine derivatives, respectively. Density functional theory calculations provide a plausible rationale for the reaction outcome.

Cooperative protein-solvent tuning of proton transfer energetics: carbonic anhydrase as a case study.

We investigate the coupling between the proton transfer (PT) energetics and the protein-solvent dynamics using the intra-molecular PT in wild type (wt) human carbonic anhydrase II and its ten-fold faster mutant Y7F/N67Q as a test case. We calculate the energy variation upon PT, and from that we also calculate the PT reaction free energy, making use of a hybrid quantum mechanics/molecular dynamics approach. In agreement with the experimental data, we obtain that the reaction free energy is basically the same in the two systems. Yet, we show that the instantaneous PT energy is on average lower in the mutant possibly contributing to the faster PT rate. Analysis of the contribution to the PT energetics of the solvent and of each protein residue, also not in the vicinity of the active site, provides evidence for electrostatic tuning of the PT energy arising from the combined effect of the solvent and the protein environment. These findings open up a way to the more general task of the rational design of mutants with either enhanced or reduced PT rate.

Kinetic Profile and Molecular Dynamic Studies Show that Y229W Substitution in an NDM-1/L209F Variant Restores the Hydrolytic Activity of the Enzyme toward Penicillins, Cephalosporins, and Carbapenems.

The New Delhi metallo-ß-lactamase-1 (NDM-1) enzyme is the most common metallo-ß-lactamase identified in many Gram-negative bacteria causing severe nosocomial infections. The aim of this study was to focus the attention on non-active-site residues L209 and Y229 of NDM-1 and to investigate their role in the catalytic mechanism. Specifically, the effect of the Y229W substitution in the L209F variant was evaluated by antimicrobial susceptibility testing, kinetic, and molecular dynamic (MD) studies. The Y229W single mutant and L209F-Y229W double mutant were generated by site-directed mutagenesis. The Km , kcat, and kcat/Km kinetic constants, calculated for the two mutants, were compared with those of (wild-type) NDM-1 and the L209F variant. Compared to the L209F single mutant, the L209F-Y229W double mutant showed a remarkable increase in kcat values of 100-, 240-, 250-, and 420-fold for imipenem, meropenem, benzylpenicillin, and cefepime, respectively. In the L209F-Y229W enzyme, we observed a remarkable increase in kcat/Km of 370-, 140-, and 80-fold for cefepime, meropenem, and cefazolin, respectively. The same behavior was noted using the antimicrobial susceptibility test. MD simulations were carried out on both L209F and L209F-Y229W enzymes complexed with benzylpenicillin, focusing attention on the overall mechanical features and on the differences between the two systems. With respect to the L209F variant, the L209F-Y229W double mutant showed mechanical stabilization of loop 10 and the N-terminal region. In addition, Y229W substitution destabilized both the C-terminal region and the region from residues 149 to 154. The epistatic effect of the Y229W mutation jointly with the stabilization of loop 10 led to a better catalytic efficiency of ß-lactams. NDM numbering is used in order to facilitate the comparison with other NDM-1 studies.

Modelling vibrational relaxation in complex molecular systems.

In this paper we show how it is possible to treat the quantum vibrational relaxation of a chromophore, embedded in a complex atomic-molecular environment, via the explicit solution of the time-dependent Schroedinger equation once using a proper separation between quantum and semiclassical degrees of freedom. The rigorous theoretical framework derived, based on first principles and making use of well defined approximations/assumptions, is utilized to construct a general model for the kinetics of the vibrational relaxation as obtained by the direct evaluation of the density matrix for all the relevant quantum state transitions. Application to (deuterated) N-methylacetamide (the typical benchmark used as a model for the amino acids) shows that the obtained theoretical-computational approach captures the essential features of the experimental process, unveiling the basic relaxation mechanism involving several vibrational state transitions.

Effect of the Incorporation of Functionalized Cyclodextrins in the Liposomal Bilayer.

Liposomes loaded with drug­cyclodextrin complexes are widely used as drug delivery systems, especially for species with low aqueous solubility and stability. Investigation of the intimate interactions of macrocycles with liposomes are essential for formulation of efficient and stable drug-in-cyclodextrin-in-liposome carriers. In this work, we reported the preparation of unilamellar vesicles of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) embedded with native ß-cyclodextrin and two synthetic derivatives: heptakis(2,3,6-tri-O-methyl)-ß-cyclodextrin (TMCD) and heptakis(2,3-di-O-acetyl)-ß-cyclodextrin (DACD). We then studied the effect of these macrocycles on the liposomal size, membrane viscosity, and liposomal stability at different temperatures and concentrations. We observed that TMCD and DACD affected vesicle size and the change of size was related to CD concentration. Irrespective of its nature, the macrocycle established interactions with the phospholipidic head groups, preventing cyclodextrins to diffuse into the lipid bilayer, as confirmed by molecular dynamics simulations. Such supramolecular structuring improves liposome stability making these colloid systems promising carriers for biologically active compounds.

P174E Substitution in GES-1 and GES-5 ß-Lactamases Improves Catalytic Efficiency toward Carbapenems.

GES-type ß-lactamases are a group of enzymes that have evolved their hydrolytic activity against carbapenems. In this study, the role of residue 174 inside the Ω-loop of GES-1 and GES-5 was investigated. GES-1P174E and GES-5P174E mutants, selected by site saturation mutagenesis, were purified and kinetically characterized. In comparison with GES-1 and GES-5 wild-type enzymes, GES-1P174E and GES-5P174E mutants exhibited lower kcat and kcat/Km values for cephalosporins and penicillins. Concerning carbapenems, GES-1P174E shared higher kcat values but lower Km values than those calculated for GES-1. The GES-1P174E and GES-5P174E mutants showed high hydrolytic efficiency for imipenem, with kcat/Km values 100- and 660-fold higher, respectively, than those of GES-1. Clavulanic acid and tazobactam are good inhibitors for both GES-1P174E and GES-5P174E Molecular dynamic (MD) simulations carried out for GES-1, GES-5, GES-1P174E, and GES-5P174E complexed with imipenem and meropenem have shown that mutation at position 174 induces a drastic increase of enzyme flexibility, in particular in the Ω-loop. The circular dichroism (CD) spectroscopy spectra of the four enzymes indicate that the P174E substitution in GES-1 and GES-5 does not affect the secondary structural content of the enzymes.