A computational approach for modeling electronic circular dichroism of solvated chromophores

The Immunosuppressive Properties of Cyclo-[D-Pro-Pro-β3-HoPhe-Phe-] Tetrapeptide Selected from Stereochemical Variants of Cyclo-[Pro-Pro-β3-HoPhe-Phe-] Peptide

The anti-inflammatory, antiviral, and anti-cancer properties, as well as the mechanism of action of cyclo-[Pro-Pro-β3-HoPhe-Phe-] tetrapeptide (denoted as 4B8M), were recently described. The aim of this work was to synthesize and evaluate the immunosuppressive actions of the stereochemical variants of 4B8M by sequential substitution of L-amino acids by D-amino acids (a series of peptides denoted as P01-P07) using parent 4B8M as a reference compound. In addition, diverse available bioinformatics tools using machine learning and artificial intelligence were tested to find the bio-pharmacokinetic and polypharmacological attributes of analyzed stereomers. All peptides were non-toxic to human peripheral blood mononuclear cells (PBMCs) and only cyclo-[D-Pro-Pro-β3-HoPhe-Phe-] peptide (P03) was capable of inhibiting mitogen-induced PBMC proliferation. The peptides inhibited the lipopolysaccharide (LPS)-induced production of tumor necrosis factor-alpha (TNF-α) to various degrees, with P04 (cyclo-[Pro-Pro-D-β3-HoPhe-Phe-]) and P03 being the most potent. For further in vivo studies, P03 was selected because it had the combined properties of inhibiting cell proliferation and TNF-α production. P03 demonstrated a comparable ability to 4B8M in the inhibition of auricle edema and lymph node cell number and in the normalization of a distorted blood cell composition in contact sensitivity to the oxazolone mouse model. In the mouse model of carrageenan-induced inflammation of the air pouch, P03 exhibited a similar inhibition of the cell number in the air pouches as 4B8M, but its inhibitory effects on the percentage of neutrophils and eosinophils in the air pouches and blood, as well as on mastocyte degranulation in the air pouches, were stronger in comparison to 4B8M. Lastly, in a mouse model of dextran sulfate-induced colitis, similar effects to 4B8M regarding thymocyte number restoration and normalization of the blood cell pictures by P03 were observed. In summary, depending on either experimental findings or in silico predictions, P03 demonstrated comparable, or even better, anti-inflammatory and bio-pharmacokinetic properties to 4B8M and may be considered as a potential therapeutic. The possibility of P00 and P03 identification by circular dichroism measurements was tested by quantum-chemical calculations.

A Joint Experimental and Theoretical Study on the Structural and Spectroscopic Properties of the Piv-Pro-d-Ser-NHMe Peptide

In this paper, we investigate the secondary structure of the Piv-Pro-d-Ser-NHMe peptide by means of nuclear magnetic resonance (NMR) and electronic circular dichroism (ECD) experiments, in conjunction with theoretical simulations based on molecular dynamics and time-dependent density functional theory calculations including polarizable embedding to account for solvent effects. The various experimental and theoretical protocols are assessed and validated, and are shown to provide a consistent description of the turn structure adopted by this peptide in solution. In addition, a simple fitting procedure is proposed to make the simulated and experimental ECD almost perfectly match. This full methodology is finally tested on another small peptide, enlightening its efficiency and robustness.

Close-Up Look at Electronic Spectroscopic Signatures of Common Pharmaceuticals in Solution

Simulating electronic properties and spectral signals requires robust computational approaches that need tuning with the system's peculiarities. In this paper, we test implicit and fully atomistic solvation models for the calculation of UV-vis and electronic circular dichroism (ECD) spectra of two pharmaceutically relevant molecules, namely, (2S)-captopril and (S)-naproxen, dissolved in aqueous solution. Room temperature molecular dynamics simulations reveal that these two drugs establish strong contacts with the surrounding solvent molecules via hydrogen bonds. Such specific interactions, which play a major role in the spectral response and are neglected in implicit approaches, are further characterized and quantified with natural bond orbital methods. Our calculations show that simulated spectra, and especially ECD, are in good agreement with experiments solely when conformational and configurational dynamics, mutual polarization, and solute-solvent repulsion effects are considered.

Directing Intrinsic Chirality in Gold Nanoclusters: Preferential Formation of Stable Enantiopure Clusters in High Yield and Experimentally Unveiling the "Super" Chirality of Au144

Chiral gold nanoclusters offer significant potential for exploring chirality at a fundamental level and for exploiting their applications in sensing and catalysis. However, their widespread use is impeded by low yields in synthesis, tedious separation procedures of their enantiomeric forms, and limited thermal stability. In this study, we investigated the direct synthesis of enantiopure chiral nanoclusters using the chiral ligand 2-MeBuSH in the fabrication of Au25, Au38, and Au144 nanoclusters. Notably, this approach leads to the unexpected formation of intrinsically chiral clusters with high yields for chiral Au38 and Au144 nanoclusters. Experimental evaluation of chiral activity by circular dichroism (CD) spectroscopy corroborates previous theoretical calculations, highlighting the stronger CD signal exhibited by Au144 compared to Au38 or Au25. Furthermore, the formation of a single enantiomeric form is experimentally confirmed by comparing it with intrinsically chiral Au38(2-PET)24 (2-PET: 2-phenylethanethiol) and is supported theoretically for both Au38 and Au144. Moreover, the prepared chiral clusters show stability against diastereoisomerization, up to temperatures of 80 °C. Thus, our findings not only demonstrate the selective preparation of enantiopure, intrinsically chiral, and highly stable thiolate-protected Au nanoclusters through careful ligand design but also support the predicted "super" chirality in the Au144 cluster, encompassing hierarchical chirality in ligands, staple configuration, and core structure.

Diverging conformations guide dipeptide self-assembly into crystals or hydrogels

The prediction of dipeptide assembly into crystals or gels is challenging. This work reveals the diverging conformational landscape that guides self-organization towards different outcomes. In silico and experimental data enabled deciphering of the electronic circular dichroism (ECD) spectra of self-assembling dipeptides to reveal folded or extended conformers as key players.

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).

Electronic circular dichroism from real-time propagation in state space

In this paper, we propose to compute the electronic circular dichroism (ECD) spectra of chiral molecules using a real-time propagation of the time-dependent Schrödinger equation (TDSE) in the space of electronic field-free eigenstates, by coupling TDSE with a given treatment of the electronic structure of the target. The time-dependent induced magnetic moment is used to compute the ECD spectrum from an explicit electric perturbation. The full matrix representing the transition magnetic moment in the space of electronic states is generated from that among pairs of molecular orbitals. In the present work, we show the ECD spectra of methyloxirane, of several conformers of L-alanine, and of the Λ-Co(acac)₃ complex, computed from a singly excited ansatz of time-dependent density functional theory eigenstates. The time-domain ECD spectra properly reproduce the frequency-domain ones obtained in the linear-response regime and quantitatively agree with the available experimental data. Moreover, the time-domain approach to ECD allows us to naturally go beyond the ground-state rotationally averaged ECD spectrum, which is the standard outcome of the linear-response theory, e.g., by computing the ECD spectra from electronic excited states.

The Conformational Dynamics of the Ligands Determines the Electronic Circular Dichroism of the Chiral Au38(SC2H4Ph)24 Cluster

Effects of the conformational dynamics of 2-PET protective ligands on the electronic circular dichroism (ECD) of the chiral Au38(SC2H4Ph)24 cluster are investigated. We adopt a computational protocol in which ECD spectra are calculated via the first principle polTDDFT approach on a series of conformations extracted from MD simulations by using Essential Dynamics (ED) analysis, and then properly weighted to predict the final spectrum. We find that the experimental spectral features are well reproduced, whereas significant discrepancies arise when the spectrum is calculated using the experimental X-ray structure. This result unambiguously demonstrates the need to account for the conformational effects in the ECD modeling of chiral nanoclusters. The present procedure proved to be able of capturing the essential conformational features of the dynamic Au38(SC2H4Ph)24 system, opening the possibility to model the ECD of soluble chiral nanoclusters in a realistic way.

Adhesion of Bis-Salphen-Based Coordination Polymers to Graphene: Insights from Free Energy Perturbation Study

Manipulation of nanoscale objects using molecular self-assembly is a potent tool to achieve large scale nanopatterning with small effort. Coordination polymers of bis-salphen compounds based on zinc have demonstrated their ability to align carbon nanotubes into micro-scale networks with an unusual "rings-and-rods" pattern. This paper investigates how the compounds interact with pristine and functionalized graphene using density functional theory calculations and molecular dynamic simulations. Using the free energy perturbation method we will show how the addition of phenyl side groups to the core compound and functionalization of graphene affect the stability, mobility and conformation adopted by a dimer of bis-(Zn)salphen compound adsorbed on graphene surface and what it can reveal about the arrangement of chains of bis-(Zn)salphen polymer around carbon nanotubes during the self-assembly of microscale networks.

Tailoring the interaction between a gold nanocluster and a fluorescent dye by cluster size: creating a toolbox of range-adjustable pH sensors

We present a novel strategy for tailoring the fluorescent azadioxatriangulenium (KU) dye-based pH sensor to the target pH range by regulating the pK _a value of the gold nanoclusters. Based on the correlation between the pK _a and surface curvature of ligand-protected nanoparticles, the pK _a value of the gold nanoclusters was controlled by size. In particular, three different-sized para-mercaptobenzoic acid (p-MBA) protected gold nanoclusters, Au₂₅(p-MBA)₁₈, Au₁₀₂(p-MBA)₄₄, and Au_210-230(p-MBA)_70-80 were used as the regulator for the pH range of the KU response. The negatively charged gold nanoclusters enabled the positively charged KU to bind to the surface, forming a complex and quenching the fluorescence of the KU by the energy transfer process. The fluorescence was restored after adjusting the surface charge of the gold nanocluster by controlling the solution pH. In addition, the KU exhibited a significantly different pH response behaviour for each gold nanocluster. Au_210-230(p-MBA)_70-80 showed a higher pH response range than Au₁₀₂(p-MBA)_44, which was intuitive. However, Au₂₅(p-MBA)₁₈ showed an unexpectedly high pH response behaviour. pK _a titration measurement, molecular dynamics simulations, and essential dynamics analysis showed that small nanoclusters do not follow the scaling between the curvature and the pK _a value. Instead, the behaviour is governed by the distribution and interaction of p-MBA ligands on the nanocluster surface. This work presents an effective design strategy for fabricating a range adjustable pH sensor by understanding the protonation behaviour of the ultrasmall gold nanoclusters in an atomic range.