Distinctive Chiral Conformations Induced to Poly (naphthalene-1,4-diyl) by Helix-sense-selective Polymerization and Circularly Polarized Light Irradiation.

Optically active poly(naphthalene-1,4-diyl) was prepared through helix-sense-selective polymerization of the corresponding monomers and also through circularly polarized light (CPL) irradiation, resulting in distinctive circular dichroism (CD) spectral patterns. Chirality of the helix-sense-selective polymerization -based polymer is ascribed to preferred-handed helicity while that of the CPL-based polymer to a non-helical, chiral conformation ('biased-dihedral conformation') with preferred-handedness which was stable only in the solid state. The helix of the helix-sense-selective polymerization-based polymer gradually racemized in tetrahydrofuran while it was stabilized by aggregate formation in a hexane-dichloromethane solution. Both helix-sense-selective polymerization- and CPL-based polymers exhibited efficient circularly polarized luminescence.

Simulation workflows to predict the circular dichroism and circularly polarized luminescence of chiral materials.

Chiral materials are attracting considerable interest in various fields in view of their unique properties and optical activity. Indeed, the peculiar features of chiral materials to absorb and emit circularly polarized light enable their use in an extensive range of applications. Motivated by the interest in boosting the development of chiral materials characterized by enhanced chiroptical properties such as circular dichroism (CD) and circular polarized luminescence (CPL), we herein illustrate in this tutorial how theoretical simulations can be used for the predictions and interpretations of chiroptical data and for the identification of chiral geometries. We are focusing on computational frameworks that can be used to investigate the theoretical aspects of chiral materials' photophysical and conformational characteristics. We will then illustrate ab initio methods based on density functional theory (DFT) and its time-dependent extension (TD-DFT) to simulate CD and CPL signals, and we will exemplify a variety of enhanced sampling techniques useful for an adequate sampling of the configurational space for chiral systems.

Amplified Chirality Transfer to Aromatic Molecules through Non-specific Inclusion by Amorphous, Hyperbranched Poly(fluorenevinylene) Derivatives.

Optically active, hyperbranched, poly(fluorene-2,4,7-triylethene-1,2-diyl) [poly(fluorenevinylene)] derivatives bearing a neomenthyl group and a pentyl group at the 9-position of the fluorene backbone at various ratios acted as a chirality donor (host polymers) efficiently included naphthalene, anthracene, pyrene, 9-phenylanthracene, and 9,10-diphenyanthracene as a chirality acceptor (guest molecules) in their interior space in film as well as in solution, with the guest molecules exhibiting intense circular dichroism through chirality transfer with chirality amplification. The efficiency of the chirality transfer was much higher with higher-molar-mass polymers than lower-molar-mass ones as well as with hyperbranched polymers compared to the analogous linear ones. The hyperbranched polymers include the small molecules in their complex structure without any specific interactions at various stoichiometries. The included molecules may have ordered intermolecular arrangement that may be somewhat similar to those of liquid crystals. Naphthalene, anthracene, and pyrene included in the polymer exhibited efficient circularly polarized luminescence, where the chirality was remarkably amplified in excited states, and anthracene exhibited especially high anisotropies in the emission on the order of 10-2 .

Photo racemization of 2,2'-dihydroxy-1,1'-binaphthyl derivatives.

Photo racemization of 2,2'-dihydroxy-1,1'-binaphthyl (BINOL) and its monomethyl ether, monobutyl ether, and dimethyl ether was studied by means of circularly dichroism spectra, chiral HPLC, and theoretical calculations of rotation energy barriers. Racemization was fastest for BINOL and about one seventh as fast for the monomethyl and monobutyl ethers while it was too slow to be detected for the dimethyl ether under the present conditions.

Non-uniform Self-folding of Helical Poly(fluorenevinylene) Derivatives in the Solid State Leading to Amplified Circular Dichroism and Circularly Polarized Light Emission.

An unprecedented non-uniform self-folding of artificial polymer chains composed of turn moieties and stretched segments is presented through the design of a set of optically active poly(fluorene-2,7-diylethene-1,2-diyl) (poly(fluorenevinylene)) derivatives bearing a neomenthyl group and a pentyl group attached at the 9-position of fluorene backbone at various ratios. The folded structure is formed and stabilized through inter-chain interactions in the solid state, leading to remarkably enhanced chiroptical properties (chirality amplification) in terms of circular dichroism (CD) and circularly polarized light (CPL) emission. This phenomenon is rationalized by experimental and theoretical CD and CPL spectral analyses. The polymer arrangements in the solid state were further assessed through transmission electron microscopic observations combined with enhanced sampling molecular dynamics simulations in the solid state revealing the thin film organizations.

Multi-replica biased sampling for photoswitchable π-conjugated polymers.

In recent years, π-conjugated polymers are attracting considerable interest in view of their light-dependent torsional reorganization around the π-conjugated backbone, which determines peculiar light-emitting properties. Motivated by the interest in designing conjugated polymers with tunable photoswitchable pathways, we devised a computational framework to enhance the sampling of the torsional conformational space and, at the same time, estimate ground- to excited-state free-energy differences. This scheme is based on a combination of Hamiltonian Replica Exchange Method (REM), parallel bias metadynamics, and free-energy perturbation theory. In our scheme, each REM samples an intermediate unphysical state between the ground and the first two excited states, which are characterized by time-dependent density functional theory simulations at the B3LYP/6-31G* level of theory. We applied the method to a 5-mer of 9,9-dioctylfluorene and found that upon irradiation, this system can undergo a dihedral inversion from -155° to 155°, crossing a barrier that decreases from 0.1 eV in the ground state (S0) to 0.05 eV and 0.04 eV in the first (S1) and second (S2) excited states. Furthermore, S1 and even more S2 were predicted to stabilize coplanar dihedrals, with a local free-energy minimum located at ±44°. The presence of a free-energy barrier of 0.08 eV for the S1 state and 0.12 eV for the S2 state can trap this conformation in a basin far from the global free-energy minimum located at 155°. The simulation results were compared with the experimental emission spectrum, showing a quantitative agreement with the predictions provided by our framework.

Unraveling the internal conversion process within the Q-bands of a chlorophyll-like-system through surface-hopping molecular dynamics simulations.

The non-radiative relaxation process within the Q-bands of chlorophylls represents a crucial preliminary step during the photosynthetic mechanism. Despite several experimental and theoretical efforts performed in order to clarify the complex dynamics characterizing this stage, a complete understanding of this mechanism is still far to be reached. In this study, non-adiabatic excited-state molecular dynamic simulations have been performed to model the non-radiative process within the Q-bands for a model system of chlorophylls. This system has been considered in the gas phase and then, to have a more representative picture of the environment, with implicit and mixed implicit-explicit solvation models. In the first part of this analysis, absorption spectra have been simulated for each model in order to guide the setup for the non-adiabatic excited-state molecular dynamic simulations. Then, non-adiabatic excited-state molecular dynamic simulations have been performed on a large set of independent trajectories and the population of the Qx and Qy states has been computed as the average of all the trajectories, estimating the rate constant for the process. Finally, with the aim of investigating the possible role played by the solvent in the Qx-Qy crossing mechanism, an essential dynamic analysis has been performed on the generated data, allowing one to find the most important motions during the simulated dynamics.

Role of specific solute-solvent interactions on the photophysical properties of distyryl substituted BODIPY derivatives.

Computational spectroscopy is now a valuable tool to better understand experimental spectroscopic data. Recent advancements in computational methods allow the study of electronic excited states of medium-large molecular systems with greater accuracy, simulating spectra directly comparable with experiments. Besides the system of interest, proper modelling of the environment effects is essential to reach a sufficient level of accuracy. This is especially important for solvents with high polarity or those able to establish specific interactions with the solute. One of the most straightforward and fastest ways to deal with solvent effects is the use of implicit models, like the polarizable continuum model (PCM). However, these models cannot describe specific solute-solvent interactions. In this case, one possible solution is the adoption of a mixed model that includes the solute and few solvent molecules within a PCM cavity. With the aim of investigating the role played by specific solute-solvent interactions, the behaviour of styryl substituted BODIPY systems in methanol has been studied at a theoretical level, considering both implicit and mixed implicit-explicit models. In the first part of our analysis, vibrationally resolved electronic spectra for absorption and emission processes have been simulated, using several sets of coordinates to represent normal modes, and compared with experiments. Then, to verify if specific solute-solvent interactions play a crucial role during the ultrafast intraband relaxation processes, 2D electronic spectra were also simulated, which could provide valuable information on the complex ultrafast dynamics.

The role of the halogen bond in iodothyronine deiodinase: Dependence on chalcogen substitution in naphthyl-based mimetics.

The effects on the activity of thyroxine (T4) due to the chalcogen replacement in a series of peri-substituted naphthalenes mimicking the catalytic function of deiodinase enzymes are computationally examined using density functional theory. In particular, T4 inner-ring deiodination pathways assisted by naphthyl-based models bearing two tellurols and a tellurol-thiol pair in peri-position are explored and compared with the analogous energy profiles for the naphthalene mimic having two selenols. The presence of a halogen bond (XB) in the intermediate formed in the first step and involved in the rate-determining step of the reaction is assumed to facilitate the process increasing the rate of the reaction. The rate-determining step calculated energy barrier heights allow rationalizing the experimentally observed superior catalytic activity of tellurium containing mimics. Charge displacement analysis is used to ascertain the presence and the role of the electron density charge transfer occurring in the rate-determining step of the reaction, suggesting the incipient formation or presence of a XB interaction. © 2019 Wiley Periodicals, Inc.

On the simulation of vibrationally resolved electronic spectra of medium-size molecules: the case of styryl substituted BODIPYs.

BODIPY dyes are used in a variety of applications because of their peculiar spectroscopic and photo-physical properties that vary depending on the stereochemistry of the functional groups attached to the boron-dipyrromethene core structure. In this work, we have applied several computational methods, adapted for semi-rigid molecules based on the Franck-Condon principle, for the study of the optical properties of BODIPY systems and for the understanding of the influence of functional groups on their spectroscopic features. We have analyzed the electronic spectra of two styryl substituted BODIPY molecules of technological interest, properly taking into account the vibronic contribution. For comparison with recently recorded experimental data in methanol, the vibrationally resolved electronic spectra of these systems were computed using both Time-Independent (TI) and Time-Dependent (TD) formalisms. The first step toward the analysis of optical properties of the styryl modified BODIPYs was a benchmark of several density functionals, to select the most appropriate one. We have found that all benchmarked functionals provide good results in terms of band shape but some of them show strong discrepancies in terms of band position. Beyond the issue of the electronic structure calculation method, different levels of sophistication can be adopted for the calculation of vibronic transitions. In this study, the effect of mode couplings and the influence of the Herzberg-Teller terms on the theoretical spectra has been investigated. It has been found that all levels of theory considered give reproducible results for the investigated systems: band positions and shapes are similar at all levels and little improvements have been found in terms of band shape with the inclusion of Herzberg-Teller effect. Inclusion of temperature effects proved to be challenging due to the important impact of large amplitude motions. Better agreement can be achieved by adopting a suitable set of coordinates coupled with a reduced-dimensionality scheme.

Mechanism of thyroxine deiodination by naphthyl-based iodothyronine deiodinase mimics and the halogen bonding role: a DFT investigation.

This paper deals with a systematic density functional theory (DFT) study aiming to unravel the mechanism of the thyroxine (T4) conversion into 3,3',5-triiodothyronine (rT3) by using different bio-inspired naphthyl-based models, which are able to reproduce the catalytic functions of the type-3 deiodinase ID-3. Such naphthalenes, having two selenols, two thiols, and a selenol-thiol pair in peri positions, which were previously synthesized and tested in their deiodinase activity, are able to remove iodine selectively from the inner ring of T4 to produce rT3. Calculations were performed including also an imidazole ring that, mimicking the role of the His residue, plays an essential role deprotonating the selenol/thiol moiety. For all the used complexes, the calculated potential energy surfaces show that the reaction proceeds via an intermediate, characterized by the presence of a X-I-C (X=Se, S) halogen bond, whose transformation into a subsequent intermediate in which the C-I bond is definitively cleaved and the incipient X-I bond is formed represents the rate-determining step of the whole process. The calculated trend in the barrier heights of the corresponding transition states allows us to rationalize the experimentally observed superior deiodinase activity of the naphthyl-based compound with two selenol groups. The role of the peri interactions between chalcogen atoms appears to be less prominent in determining the deiodination activity.

Theoretical study of silver-ion-mediated base pairs: the case of C-Ag-C and C-Ag-A systems.

Silver-mediated base pairs applied to DNA represent a new biomacromolecular nanomaterial useful for generating nanodevices as ion sensors. Reported herein is a full quantum chemical study devoted to give further knowledge on the electronic and energetic properties of C-Ag-C and mixed C-Ag-A mismatched base pairs. The B3LYP functional in conjunction with the dispersion effects (D3) has been applied. Single-point calculations have been also performed by using the M06-L functional. The investigation of their behavior has been extended to the duplex DNA modeled by the (dC-Ag-dC)2 and (dC-Ag-dA)2 more complex systems. The solvent effect has been taken into account by the conductor-like screening model, COSMO. In the case of mixed C-Ag-A and (dC-Ag-dA)2 systems, both the Watson-Crick and Hoogsteen arrangements have been taken into account. Results show that for (dC-Ag-dA)2 systems, the binding energies are almost double that of the corresponding values of C-Ag-A ones.

Unraveling Copper Exchange in the Atox1-Cu(I)-Mnk1 Heterodimer: A Simulation Approach.

Copper, an essential metal for various cellular processes, requires tight regulation to prevent cytotoxicity. Intracellular pathways crucial for maintaining optimal copper levels involve soluble and membrane transporters, namely, metallochaperones and P-type ATPases, respectively. In this study, we used a simulation workflow based on free-energy perturbation (FEP) theory and parallel bias metadynamics (PBMetaD) to predict the Cu(I) exchange mechanism between the human Cu(I) chaperone, Atox1, and one of its two physiological partners, ATP7A. ATP7A, also known as the Menkes disease protein, is a transmembrane protein and one of the main copper-transporting ATPases. It pumps copper into the trans-Golgi network for the maturation of cuproenzymes and is also essential for the efflux of excess copper across the plasma membrane. In this analysis, we utilized the nuclear magnetic resonance (NMR) structure of the Cu(I)-mediated complex between Atox1 and the first soluble domain of the Menkes protein (Mnk1) as a starting point. Independent free-energy simulations were conducted to investigate the dissociation of both Atox1 and Mnk1. The calculations revealed that the two dissociations require free energy values of 6.3 and 6.2 kcal/mol, respectively, following a stepwise dissociation mechanism.

BDNF- and VEGF-Responsive Stimulus to an NGF Mimic Cyclic Peptide with Copper Ionophore Capability and Ctr1/CCS-Driven Signaling.

Neurotrophins are a family of growth factors that play a key role in the development and regulation of the functioning of the central nervous system. Their use as drugs is made difficult by their poor stability, cellular permeability, and side effects. Continuing our effort to use peptides that mimic the neurotrophic growth factor (NGF), the family model protein, and specifically the N-terminus of the protein, here we report on the spectroscopic characterization and resistance to hydrolysis of the 14-membered cyclic peptide reproducing the N-terminus sequence (SSSHPIFHRGEFSV (c-NGF(1-14)). Far-UV CD spectra and a computational study show that this peptide has a rigid conformation and left-handed chirality typical of polyproline II that favors its interaction with the D5 domain of the NGF receptor TrkA. c-NGF(1-14) is able to bind Cu2+ with good affinity; the resulting complexes have been characterized by potentiometric and spectroscopic measurements. Experiments on PC12 cells show that c-NGF(1-14) acts as an ionophore, influencing the degree and the localization of both the membrane transporter (Ctr1) and the copper intracellular transporter (CCS). c-NGF(1-14) induces PC12 differentiation, mimics the protein in TrkA phosphorylation, and activates the kinase cascade, inducing Erk1/2 phosphorylation. c-NGF(1-14) biological activities are enhanced when the peptide interacts with Cu2+ even with the submicromolar quantities present in the culture media as demonstrated by ICP-OES measurements. Finally, c-NGF(1-14) and Cu2+ concur to activate the cAMP response element-binding protein CREB that, in turn, induces the brain-derived neurotrophic factor (BDNF) and the vascular endothelial growth factor (VEGF) release.

Multinuclear Metal-Binding Ability of the N-Terminal Region of Human Copper Transporter Ctr1: Dependence Upon pH and Metal Oxidation State.

The 14mer peptide corresponding to the N-terminal region of human copper transporter Ctr1 was used to investigate the intricate mechanism of metal binding to this plasma membrane permease responsible for copper import in eukaryotic cells. The peptide contains a high-affinity ATCUN Cu(II)/Ni(II)-selective motif, a methionine-only MxMxxM Cu(I)/Ag(I)-selective motif and a double histidine HH(M) motif, which can bind both Cu(II) and Cu(I)/Ag(I) ions. Using a combination of NMR spectroscopy and electrospray mass spectrometry, clear evidence was gained that the Ctr1 peptide, at neutral pH, can bind one or two metal ions in the same or different oxidation states. Addition of ascorbate to a neutral solution containing Ctr11-14 and Cu(II) in 1:1 ratio does not cause an appreciable reduction of Cu(II) to Cu(I), which is indicative of a tight binding of Cu(II) to the ATCUN motif. However, by lowering the pH to 3.5, the Cu(II) ion detaches from the peptide and becomes susceptible to reduction to Cu(I) by ascorbate. It is noteworthy that at low pH, unlike Cu(II), Cu(I) stably binds to methionines of the peptide. This redox reaction could take place in the lumen of acidic organelles after Ctr1 internalization. Unlike Ctr11-14-Cu(II), bimetallic Ctr11-14-2Cu(II) is susceptible to partial reduction by ascorbate at neutral pH, which is indicative of a lower binding affinity of the second Cu(II) ion. The reduced copper remains bound to the peptide, most likely to the HH(M) motif. By lowering the pH to 3.5, Cu(I) shifts from HH(M) to methionine-only coordination, an indication that only the pH-insensitive methionine motif is competent for metal binding at low pH. The easy interconversion of monovalent cations between different coordination modes was supported by DFT calculations.

Neurotrophic Activity and Its Modulation by Zinc Ion of a Dimeric Peptide Mimicking the Brain-Derived Neurotrophic Factor N-Terminal Region.

Brain-derived neurotrophic factor (BDNF) is a neurotrophin (NT) essential for neuronal development and synaptic plasticity. Dysregulation of BDNF signaling is implicated in different neurological disorders. The direct NT administration as therapeutics has revealed to be challenging. This has prompted the design of peptides mimicking different regions of the BDNF structure. Although loops 2 and 4 have been thoroughly investigated, less is known regarding the BDNF N-terminal region, which is involved in the selective recognition of the TrkB receptor. Herein, a dimeric form of the linear peptide encompassing the 1-12 residues of the BDNF N-terminal (d-bdnf) was synthesized. It demonstrated to act as an agonist promoting specific phosphorylation of TrkB and downstream ERK and AKT effectors. The ability to promote TrkB dimerization was investigated by advanced fluorescence microscopy and molecular dynamics (MD) simulations, finding activation modes shared with BDNF. Furthermore, d-bdnf was able to sustain neurite outgrowth and increase the expression of differentiation (NEFM, LAMC1) and polarization markers (MAP2, MAPT) demonstrating its neurotrophic activity. As TrkB activity is affected by zinc ions in the synaptic cleft, we first verified the ability of d-bdnf to coordinate zinc and then the effect of such complexation on its activity. The d-bdnf neurotrophic activity was reduced by zinc complexation, demonstrating the role of the latter in tuning the activity of the new peptido-mimetic. Taken together our data uncover the neurotrophic properties of a novel BDNF mimetic peptide and pave the way for future studies to understand the pharmacological basis of d-bdnf action and develop novel BDNF-based therapeutic strategies.

Two-Dimensional Electronic Spectroscopy Reveals Dynamics and Mechanisms of Solvent-Driven Inertial Relaxation in Polar BODIPY Dyes.

In this work, we demonstrate the use of two-dimensional electronic spectroscopy (2DES) to study the mechanism and time scale of the femtosecond Stokes shift dynamics in molecules characterized by intramolecular charge transfer, such as distyryl-functionalized boron dipyrromethene (BODIPY) molecules. The obtained results demonstrate that 2DES allows clear and direct visualization of the phenomenon. The analysis of the 2D data in terms of 2D frequency-frequency decay associated maps provides indeed not only the time scale of the relaxation process but also the starting and the final point of the energy flow and the associated reorganization energy, identified by looking at the coordinates of a negative signature below the diagonal. The sensitivity of the 2DES technique to vibrational coherence dynamics also allowed the identification of a possible relaxation mechanism involving specific interaction between a vibrational mode of the dye and the solvent.

A DFT investigation of a bulky biomimetic model catalyzing the 5'-outer ring deiodination of thyroxine.

This paper illustrates the outcomes of a density functional theory investigation aimed at unraveling mechanistic aspects of the 5'-outer ring deiodination process of thyroxine (T4) assisted by the sterically protected organoselenol compound BpqSeH. BpqSeH, which was previously synthesized and tested for its deiodinase activity, is able to afford the active hormone 3,5,3'-tetraiodothyronine (T3) by selective outer-ring deiodination of T4, and to protect the SeH moiety inside the nano-sized molecular cavity from further reactivity, allowing its isolation and characterization. Calculations were also performed including an imidazole ring that, mimicking a His residue in the active site of the original enzyme, plays an crucial role in deprotonating the selenol moiety. Both the suggested enol/keto tautomerization and the previously proven formation of an intermediate whose main characteristic is the presence of a Se⋯I⋯C halogen bond, were examined along the pathway leading to 5'-outer ring deiodination. The calculated potential energy surface showed that neither the pathway encompassing enol/keto tautomerism nor the formation of a halogen bond paving the way to C-I bond breaking and chalcogen-I bond forming is viable. The exergonic formation of the final selenenyl iodide product confirms the stabilization effect of the molecular cavity. Graphical Abstract Computed free energy profile describing the 5'-outer deiodination of thyroxine assisted by the steric hindered organoselenol BpqSH compound. The molecular electrostatic potential map reoported for the INT1 intermediate shows the non-covalent Se-I interaction, due to the attraction between charges of opposite sign, that weakens the C-I bond and prepares the formation of the new Se-I bond.