Search

Discovered in 2016, the enzyme PETase, secreted by bacterial Ideonella Sakaiensis 201-F6, has an excellent hydrolytic activity toward poly(ethylene terephthalate) (PET) at room temperature, while it decreases at higher temperatures due to the low thermostability. Many variants have been engineered to overcome this limitation, which hinders industrial application. In this work, we systematically compare PETase wild-type (WT) and four mutants (DuraPETase, ThermoPETase, FastPETase, and HotPETase) using standard molecular dynamics (MD) simulations and unbinding free energy calculations. In particular, we analyze the enzymes' structural characteristics and binding to a tetrameric PET chain (PET4) under two temperature conditions: T 1─300 K and T 2─350 K. Our results indicate that (i) PET4 forms stable complexes with the five enzymes at room temperature (∼300 K) and (ii) most of the interactions are localized close to the active site of the protein, where the W185 and Y87 residues interact with the aromatic rings of the substrate. Specifically, (iii) the W185 side-chain explores different conformations in each variant (a phenomenon known in the literature as "W185 wobbling"). This suggests that the binding pocket retains structural plasticity and flexibility among the variants, facilitating substrate recognition and localization events at moderate temperatures. Moreover, (iv) PET4 establishes aromatic interactions with the catalytic H237 residue, stabilizing the catalytic triad composed of residues S160-H237-D206, and helping the system achieve an effective configuration for the hydrolysis reaction. Conversely, (v) the binding affinity decreases at a higher temperature (∼350 K), retaining moderate interactions only for HotPETase. Finally, (vi) MD simulations of complexes formed with poly(ethylene-2,5-furan dicarboxylate) (PEF) show no persistent interactions, suggesting that these enzymes are not yet optimized for binding this alternative semiaromatic plastic polymer. Our study offers valuable insights into the structural stability of these enzymes and the molecular determinants driving PET binding onto their surfaces, sheds light on the mechanistic steps that precede the onset of hydrolysis, and provides a foundation for future enzyme optimization.

Acetonitrile, a polar molecule that cannot form hydrogen bonds on its own, interacts with solvent molecules mainly through the lone pair of its nitrogen atom and the π electrons of its CN triple bond [Correction added on 17 July 2024, after first online publication: Acetole has been changed to Acetonitrile in the preceeding sentence.]. Interestingly, acetonitrile exhibits an unexpected strengthening of the triple bond's force constant in an aqueous environment, leading to an upshift (blueshift) in the corresponding stretching vibration: this effect contrasts with the usual consequence of hydrogen bonding on the vibrational frequencies of the acceptor groups, that is, frequency redshift. This investigation elucidates this phenomenon using Raman spectroscopy to examine the behavior of acetonitrile in organic solvent, water, and silver ion aqueous solutions, where an even more pronounced upshift is observed. Raman spectroscopy is particularly well suited for analyzing aqueous solutions due to the minimal scattering effect of water molecules across most of the vibrational spectrum. Computational approaches, both static and dynamical, based on Density Functional Theory and hybrid functionals, are employed here to interpret these findings, and accurately reproduce the vibrational frequencies of acetonitrile in different environments. Our calculations also allow an explanation for this unique behavior in terms of electric charge displacements. On the other hand, the study of the interaction of acetonitrile with water molecules and metal ions is relevant for the use of this molecule as a solvent in both chemical and pharmaceutical applications., (© 2024 Wiley Periodicals LLC.)

This investigation involved an ab initio and Density Functional Theory (DFT) analysis of the hydrolysis mechanism and energetics in a borate network. The focus was on understanding how water molecules interact with and disrupt the borate network, an area where the experimental data are scarce and unreliable. The modeled system consisted of two boron atoms, bridging oxygen atoms, and varying numbers of water molecules. This setup allows for an exploration of hydrolysis under different environmental conditions, including the presence of OH - or H + ions to simulate basic or acidic environments, respectively. Our investigation utilized both ab initio calculations at the MP2 and CCSD(T) levels and DFT with a range of exchange-correlation functionals. The findings indicate that the borate network is significantly more susceptible to hydrolysis in a basic environment, with respect to an acidic or to a neutral pH setting. The inclusion of explicit water molecules in the calculations can significantly affect the results, depending on the nature of the transition state. In fact, some transition states exhibited closed-ring configurations involving water and the boron-oxygen-boron network; in these cases, there were indeed more water molecules corresponding to lower energy barriers for the reaction, suggesting a crucial role of water in stabilizing the transition states. This study provides valuable insights into the hydrolysis process of borate networks, offering a detailed comparison between different computational approaches. The results demonstrate that the functionals B3LYP, PBE0, and wB97Xd closely approximated the reference MP2 and CCSD(T) calculated reaction pathways, both qualitatively in terms of the mechanism, and quantitatively in terms of the differences in the reaction barriers within the 0.1-0.2 eV interval for the most plausible reaction pathways. In addition, CAM-B3LYP also yielded acceptable results in all cases except for the most complicated pathway. These findings are useful for guiding further computational studies, including those employing machine learning approaches, and experimental investigations requiring accurate reference data for hydrolysis reactions in borate networks.

Metadynamics simulations driven by using two X-ray diffraction peaks identified three alternative crystallization pathways of the lithium disilicate crystal from the melt. The most favorable one passes through the formation of disordered layered structures undergoing internal ordering in a second step. The second pathway involves the formation of phase-separated structures composed of nuclei of β-cristobalite crystals surrounded by lithium-rich phases in which metasilicate chains are formed. The conversion of these structures to the stable lithium disilicate crystal involves an intermediate structure whose silicate layers are connected by silicate rings with the energy barrier of 2.5 kJ/mol per formula unit (f.u.). The third pathway is highly unlikely because of the huge energy barrier involved (20 kJ/mol per f.u.). This path also involves the passage through a phase-separated structure of an indefinite silica region surrounded mainly by amorphous lithium oxide.

Density functional theory (DFT) calculations allow us to reproduce the SERS (surface-enhanced Raman scattering) spectra of molecules adsorbed on nanostructured metal surfaces and extract the most information this spectroscopy is potentially able to provide. The latter point mainly concerns the anchoring mechanism and the bond strength between molecule and metal as well as the structural and electronic modifications of the adsorbed molecule. These findings are of fundamental importance for the application of this spectroscopic technique. This review presents and discusses some SERS-DFT studies carried out in Italy as a collaboration between the universities of Modena and Reggio-Emilia and of Florence, giving an overview of the information that we can extract with a combination of experimental SERS spectra and DFT modeling. In addition, a selection of the most recent studies and advancements on the DFT approach to SERS spectroscopy is reported with commentary.

Cancer is a leading cause of death worldwide, its genesis and progression are caused by homeostatic errors, and reactive oxygen species play a major role in promoting aberrant cancer homeostasis. In this scenario, curcumin could be an interesting candidate due to its versatile antioxidant, anti-inflammatory, anti-tumor, anti-HIV, and anti-infection properties. Nonetheless, the major problem related to its use is its poor oral bioavailability, which can be overcome by encapsulating it into small particles, such as hydrogel beads containing mesoporous silica. In this work, various systems have been synthesized: starting from mesoporous silica glasses (MGs), cerium-containing MGs have been produced; then, these systems have been loaded with 4 to 6% of curcumin. Finally, various MGs at different compositions have been included in alginate beads. In vitro studies showed that these hybrid materials enable the stabilization and effective delivery of curcumin and that a synergic effect can be achieved if Ce 3+ /Ce 4+ and curcumin are both part of the beads. From swelling tests, it is possible to confirm a controlled curcumin release compartmentalized into the gastrointestinal tract. For all beads obtained, a curcumin release sufficient to achieve the antioxidant threshold has been reached, and a synergic effect of cerium and curcumin is observed. Moreover, from catalase mimetic activity tests, we confirm the well-known catalytic activity of the couple Ce 3+ /Ce 4+ . In addition, an extremely good radical scavenging effect of curcumin has been demonstrated. In conclusion, these systems, able to promote an enzymatic-like activity, can be used as drug delivery systems for curcumin-targeted dosing., Competing Interests: The authors declare no conflict of interest

The surface-enhanced Raman scattering (SERS) spectra of piperidine adsorbed on silver/chloride colloids were studied by a combined density functional theory (DFT)/time dependent DFT (TD-DFT) approach. The mechanism of chemical enhancement on the Raman signals is due to at least two contributions: the first comes from the changes in the molecular force constants and the dynamic polarizabilities of the normal modes, when the molecule is chemisorbed. DFT calculations satisfactorily reproduce the SERS spectra of piperidine adsorbed on silver, showing that the species formed on the silver particle is a complex formed by a deprotonated piperidine linked to a silver cation. A second contribution to the SERS chemical enhancement is due to a resonance Raman effect occurring when the wavelength of the Raman excitation falls within the electronic excitation band of the molecule/metal complex. Actually, the SERS spectra of piperidine show a significant dependence on the wavelength of the laser excitation, with a marked enhancement in the green-light region. TD-DFT calculations on the most-probable complex explain this behavior, because a strong excitation band of the complex is calculated in the green spectral region. This pinpoints that a resonance between the exciting radiation and the absorption band of this complex is responsible for this enhancement effect.

Metadynamics (MetaD) is a useful technique to study rare events such as crystallization. It has been only recently applied to study nucleation and crystallization in glass-forming liquids such as silicates, but the optimal set of parameters to drive crystallization and obtain converged free energy surfaces is still unexplored. In this work, we systematically investigated the effects of the simulation conditions to efficiently study the thermodynamics and mechanism of crystallization in highly viscous systems. As a prototype system, we used fused silica, which easily crystallizes to β-cristobalite through MetaD simulations, owing to its simple microstructure. We investigated the influence of the height, width, and bias factor used to define the biasing Gaussian potential, as well as the effects of the temperature and system size on the results. Among these parameters, the bias factor and temperature seem to be most effective in sampling the free energy landscape of melt to crystal transition and reaching convergence more quickly. We also demonstrate that the temperature rescaling from T > T m is a reliable approach to recover free energy surfaces below T m , provided that the temperature gap is below 600 K and the configurational space has been properly sampled. Finally, albeit a complete crystallization is hard to achieve with large simulation boxes, these can be reliably and effectively exploited to study the first stages of nucleation.

Protein interactions with engineered gold nanoparticles (AuNPs) and the consequent formation of the protein corona are very relevant and poorly understood biological phenomena. The nanoparticle coverage affects protein binding modalities, and the adsorbed protein sites influence interactions with other macromolecules and cells. Here, we studied four common blood proteins, i.e., hemoglobin, serum albumin, α1-antiproteinase, and complement C3, interacting with AuNPs covered by hydrophobic 11-mercapto-1-undecanesulfonate (MUS). We use Molecular Dynamics and the Martini coarse-grained model to gain quantitative insight into the kinetics of the interaction, the physico-chemical characteristics of the binding site, and the nanoparticle adsorption capacity. Results show that proteins bind to MUS-capped AuNPs through strong hydrophobic interactions and that they adapt to the AuNP surfaces to maximize the contact surface, but no dramatic change in the secondary structure of the proteins is observed. We suggest a new method to calculate the maximum adsorption capacity of capped AuNPs based on the effective surface covered by each protein, which better represents the realistic behavior of these systems.

Translocator protein 18 kDa [TSPO or peripheral-type benzodiazepine receptor (PBR)] was identified in the search of binding sites for benzodiazepine anxiolytic drugs in peripheral regions. In these areas, binding sites for TSPO ligands were recognized in steroid-producing tissues. TSPO plays an important role in many cellular functions, and its coding sequence is highly conserved across species. TSPO is located predominantly on the membrane of mitochondria and is overexpressed in several solid cancers. TSPO basal expression in the CNS is low, but it becomes high in neurodegenerative conditions. Thus, TSPO constitutes not only as an outstanding drug target but also as a valuable marker for the diagnosis of a number of diseases. The aim of the present article is to show the lesson we have learned from our activity in TSPO medicinal chemistry and in approaching the targeted delivery to mitochondria by means of TSPO ligands.

Ceria (CeO 2 ) is a well-known catalytic oxide with many environmental, energy production, and industrial applications, most of them involving water as a reactant, byproduct, solvent, or simple spectator. In this work, we parameterized a Ce/O/H ReaxFF for the study of ceria and ceria/water interfaces. The parameters were fitted to an ab initio training set obtained at the DFT/PBE0 level, including the structures, cohesive energies, and elastic properties of the crystalline phases Ce, CeO 2 , and Ce 2 O 3 ; the O-defective structures and energies of vacancy formation on CeO 2 bulk and CeO 2 (111) surface, as well as the absorption and reaction energies of H 2 and H 2 O molecules on CeO 2 (111). The new potential reproduced reasonably well all the fitted properties as well as the relative stabilities of the different ceria surfaces, the oxygen vacancies formation, and the energies and structures of associative and dissociative water molecules on them. Molecular dynamics simulations of the liquid water on the CeO 2 (111) and CeO 2 (100) surfaces were carried out to study the coverage and the mechanism of water dissociation. After equilibration, on average, 35% of surface sites of CeO 2 (111) are hydroxylated whereas 15% of them are saturated with molecular water associatively adsorbed. As for the CeO 2 (100) surface, we observed that water preferentially dissociates covering 90% of the available surface sites in excellent agreement with recent experimental findings.

The bacterium strain Ideonella sakaiensis 201-F6 is able to hydrolyze low-crystallinity PET films at 30 °C due to two enzymes named PETase and MHETase. Since its discovery, many efforts have been dedicated to elucidating the structure and features of those two enzymes, and various authors have highlighted the necessity to optimize both the substrate binding site and the global structure in order to enhance the stability and catalytic activity of these PET biocatalysts so as to make them more suitable for industrial applications. In this review, the strategies adopted by different research groups to investigate the structure and functionality of both PETase and MHETase in depth are described, emphasizing the advantages provided by the use of computational methods to complement and drive experiments. Subsequently, the modifications implemented with protein engineering are discussed. The versatility of the enzymes secreted by I. sakaiensis enables the prediction that they will find several applications in the disposal of PET debris, encouraging a prioritization of efforts in this prolific research field., (© 2021 Wiley-VCH GmbH.)

Amyloid-β aggregation is one of the principal causes of amyloidogenic diseases that lead to the loss of neuronal cells and to cognitive impairments. The use of gold nanoparticles treating amyloidogenic diseases is a promising approach, because the chemistry of the gold surface can be tuned in order to have a specific binding, obtaining effective tools to control the aggregation. In this paper, we show, by means of Replica Exchange Solute Tempering Molecular Simulations, how electrostatic interactions drive the absorption of Amyloid-β monomers onto citrates-capped gold nanoparticles. Importantly, upon binding, amyloid monomers show a reduced propensity in forming β-sheets secondary structures that are characteristics of mature amyloid fibrils.

Unraveling detailed mechanism of crystal nucleation from amorphous materials is challenging for both experimental and theoretical approaches. In this study, we have examined two methods to understand the initial stage of crystal precipitation from lithium disilicate glasses using molecular dynamics simulations. One of the methods is a modified exploring method to find structurally similar crystalline clusters in the glass models, enabling us to find three different embryos, such as Li 2 Si 2 O 5 (LS 2 ), Li 2 SiO 3 (LS) and Li 3 PO 4 (LP), in the 33Li 2 O·66SiO 2 ·1P 2 O 5 glass (LS 2 P1), in which P 2 O 5 is added as a nucleating agent. Interestingly, LS 2 and LP crystals were found inside the LS 2 P1 glass while LS crystal appeared on the glass surface, which agrees with experimental observations. The other method is free energy calculation using a subnano-scale spherical crystal embedded in the glass model. This method, which we called Free-Energy Seeding Method (FESM), allows us to evaluate free energy change as a function of crystal radius and to identify critical size of the crystal precipitation. The free energy profiles for LS and LS 2 crystal nuclei in the LS 2 glass models possess maximum energy at a critical radius as expected by classical nucleation theory. Furthermore, the critical radius and the energy barrier height agree well with recent experimental investigation, proving the applicability of this method to design glass-ceramics by atomistic modeling.

Amyloids-β (Aβ) fibrils are involved in several neurodegenerative diseases. In this study, atomistic molecular dynamics simulations have been used to investigate how monolayer-protected gold nanoparticles interact with Aβ(1-40) and Aβ(1-42) fibrils. Our results show that small gold nanoparticles bind with the external side of amyloid-β fibrils that is involved in the fibrillation process. The binding affinity, studied for both kinds of fibrils as a function of the monolayer composition and the nanoparticle diameter, is modulated by hydrophobic interactions and ligand monolayer conformation. Our findings thus show that monolayer-protected nanoparticles are good candidates to prevent fibril aggregation and secondary nucleation or to deliver drugs to specific fibril regions.

In this study, we compared the effects of two well-known natural compounds on the early step of the fibrillation process of amyloid-β (1-40), responsible for the formation of plaques in the brains of patients affected by Alzheimer's disease (AD). The use of extensive replica exchange simulations up to the µs scale allowed us to characterize the inhibition activity of (-)-epigallocatechin-3-gallate (EGCG) and curcumin (CUR) on unfolded amyloid fibrils. A reduced number of β-strands, characteristic of amyloid fibrils, and an increased distance between the amino acids that are responsible for the intra- and interprotein aggregations are observed. The central core region of the amyloid-β (Aβ(1-40)) fibril is found to have a high affinity to EGCG and CUR due to the presence of hydrophobic residues. Lastly, the free binding energy computed using the Poisson Boltzmann Surface Ares suggests that EGCG is more likely to bind to unfolded Aβ(1-40) fibrils and that this molecule can be a good candidate to develop new and more effective congeners to treat AD.

Increasing evidence suggests that amyloid polymorphism gives rise to different strains of amyloids with distinct toxicities and pathology-spreading properties. Validating this hypothesis is challenging due to a lack of tools and methods that allow for the direct characterization of amyloid polymorphism in hydrated and complex biological samples. Here, we report on the development of 11-mercapto-1-undecanesulfonate-coated gold nanoparticles (NPs) that efficiently label the edges of synthetic, recombinant, and native amyloid fibrils derived from different amyloidogenic proteins. We demonstrate that these NPs represent powerful tools for assessing amyloid morphological polymorphism, using cryogenic transmission electron microscopy (cryo-EM). The NPs allowed for the visualization of morphological features that are not directly observed using standard imaging techniques, including transmission electron microscopy with use of the negative stain or cryo-EM imaging. The use of these NPs to label native paired helical filaments (PHFs) from the postmortem brain of a patient with Alzheimer's disease, as well as amyloid fibrils extracted from the heart tissue of a patient suffering from systemic amyloid light-chain amyloidosis, revealed a high degree of homogeneity across the fibrils derived from human tissue in comparison with fibrils aggregated in vitro. These findings are consistent with, and strongly support, the emerging view that the physiologic milieu is a key determinant of amyloid fibril strains. Together, these advances should not only facilitate the profiling and characterization of amyloids for structural studies by cryo-EM, but also pave the way to elucidate the structural basis of amyloid strains and toxicity, and possibly the correlation between the pathological and clinical heterogeneity of amyloid diseases., Competing Interests: The authors declare no competing interest., (Copyright © 2020 the Author(s). Published by PNAS.)

The comprehension of the nonlinear effects provided by mixed alkali effect (MAE) in oxide glasses is useful to optimize glass compositions to achieve specific properties that depend on the mobility of ions, such as the chemical durability, glass transition temperature, viscosity and ionic conductivity. Although molecular dynamics (MD) simulations have already been applied to investigate the MAE on silicates, less effort has been devoted to study such phenomenon in mixed alkali aluminosilicate glasses where alkali cations can act both as modifiers, forming non-bridging oxygens and percolation channels, and as charge compensator of the AlO 4 - units present in the network. Moreover, the ionic conductivity has not been computed yet; thus, the accuracy of the atomistic simulations in reproducing the MAE on the property is still open to question. In this work, we have validated five major interatomic potentials for the classical MD simulations by modelling the structure, density, glass transition temperature and ionic conductivity for three aluminosilicate glasses, (25 - x)Na 2 O - x(K 2 O) - 10(Al 2 O 3 ) - 65(SiO 2 ) (x = 0, 12.5, 25). It was observed that only the core-shell (CS) polarizable force field well reproduces the experimentally measured MAE on T g and the ionic conductivity as well as the higher conductivity of single sodium aluminosilicate glass at low temperature and the higher conductivity of single potassium aluminosilicate glass at high temperature. The MAE is related to the suppression of jump events of the alkaline ions between dissimilar sites in the percolation channels consisting of both sodium and potassium ions as in the case of alkaline silicates. The superior reproducibility of the CS potential is originated from the larger and the flexible ring structures due to the smaller Si-O-Si inter-tetrahedra angle, creating appropriate percolation channels for ion conductivity. We also report detailed assessments for using the potential models including the CS potential for investigating MAE on aluminosilicates.

We have demonstrated that designed ankyrin repeat protein (DARPin) _9-29, which specifically targets human epidermal growth factor receptor 2 (HER2), binds tightly to gold mini nanorods (GNRs). Molecular dynamic simulations showed that a single layer of DARPin_9-29 molecules is formed on the surface of the nanorod and that conjugation with the nanorod does not involve the protein's domain responsible for specific binding to HER2. The nanorod-DARPin (DARPin-GNR) conjugate is specifically bound (in nanomolar concentrations) to human breast adenocarcinoma SK-BR-3 cells overexpressing HER2. Illumination by near-infrared light (850 nm) led to almost complete eradication of the conjugate-treated SK-BR-3 cells; the viability of epithelial human breast cancer cells expressing normal amounts of the receptor was much less affected by the illumination. The results reported here pave the way toward application of DARPin-GNR conjugates in phototherapy of cancer.

A Surface-Enhanced Raman Scattering (SERS) spectrum of 4-cyanopyridine (4CNPy) was recorded on silver plasmonic nanoparticles and analyzed by using Density Functional Theory (DFT) calculations. Two simple molecular models of the metal-4CNPy surface complex with a single silver cation or with a neutral dimer (Ag + -4CNPy, Ag 2 -4CNPy), linked through the two possible interacting sites of 4CNPy (aromatic nitrogen, N, and nitrile group, CN), were considered. The calculated vibrational wavenumbers and intensities of the adsorbate and the isolated species are compared with the experimental Raman and SERS results. The analysis of the DFT predictions and the experimental data indicates that 4CNPy adsorbs preferentially on neutral/charged active sites of the silver nanoparticles through the nitrogen atom of the aromatic ring with a perpendicular orientation.

A multiscale molecular dynamics simulation study has been carried out in order to provide in-depth information on the adsorption of hemoglobin, myoglobin, and trypsin over citrate-capped AuNPs of 15 nm diameter. In particular, determinants for single proteins adsorption and simultaneous adsorption of the three types of proteins considered have been studied by Coarse-Grained and Meso-Scale molecular simulations, respectively. The results, discussed in the light of the controversial experimental data reported in the current experimental literature, have provided a detailed description of the (i) recognition process, (ii) number of proteins involved in the early stages of corona formation, (iii) protein competition for AuNP adsorption, (iv) interaction modalities between AuNP and protein binding sites, and (v) protein structural preservation and alteration.

CeO 2 based materials are very attractive as catalytic components for industrial processes and environmentally friendly technologies; therefore, a reliable and computationally affordable theoretical description of the main properties of ceria is needed. In particular, the description of the interconversion between the Ce(IV) and Ce(III) oxidation states, on which lies the main chemical features of the cerium oxide, results in quite a challenge at the Density Functional Theory level. Here, we tested several density functional approximations, spanning from GGA to hybrid (Global, Meta-Global, and Range Separated Corrected) functionals, on the structural, vibrational, electronic, and thermochemical properties of bulk CeO 2 and Ce 2 O 3 . GGA and Meta-GGA xc best predict the thermochemical data, while the discrepancies increase with the introduction of the exact exchange in hybrid functionals. Overall, the Short Range Corrected and Global Hybrid functionals with a percentage of Exact Exchange between 16 and 25 give the best description of the crystal properties. Then, a group of the best performing functionals has been tested on the formation energy of an oxygen vacancy at the (111) CeO 2 surface. In general, increasing the amount of exact exchange in the hybrid functionals leads to a better description of the localized Ce 4 f states, while the energy of formation of the O vacancy decreases, worsening compared to the experiment.

Several neurodegenerative diseases, like Alzheimer's (AD), are characterized by amyloid fibrillar deposition of misfolded proteins, and this feature can be exploited for both diagnosis and therapy design. In this paper, structural modifications of curcumin scaffold were examined in order to improve its bioavailability and stability in physiological conditions, as well as its ability to interfere with β-amyloid fibrils and aggregates. The acid-base behaviour of curcumin derivatives, their pharmacokinetic stability in physiological conditions, and in vitro ability to interfere with Aβ fibrils at different incubation time were investigated. The mechanisms governing these phenomena have been studied at atomic level by means of molecular docking and dynamic simulations. Finally, biological activity of selected curcuminoids has been investigated in vitro to evaluate their safety and efficiency in oxidative stress protection on hippocampal HT-22 mouse cells. Two aromatic rings, π-conjugated structure and H-donor/acceptor substituents on the aromatic rings showed to be the sine qua non structural features to provide interaction and disaggregation activity even at very low incubation time (2h). Computational simulations proved that upon binding the ligands modify the conformational dynamics and/or interact with the amyloidogenic region of the protofibril facilitating disaggregation. Significantly, in vitro results on hippocampal cells pointed out protection against glutamate toxicity and safety when administered at low concentrations (1 μM). On the overall, in view of its higher stability in physiological conditions with respect to curcumin, of his rapid binding to fibrillar aggregates and strong depolymerizing activity, phtalimmide derivative K2F21 appeared a good candidate for both AD diagnostic and therapeutic purposes., (Copyright © 2018 Elsevier Ltd. All rights reserved.)

One of the principal hallmarks of Alzheimer's disease (AD) is related to the aggregation of amyloid-β fibrils in an insoluble form in the brain, also known as amyloidosis. Therefore, a prominent therapeutic strategy against AD consists of either blocking the amyloid aggregation and/or destroying the already formed aggregates. Natural products have shown significant therapeutic potential as amyloid inhibitors from in vitro studies as well as in vivo animal tests. In this study, the interaction of five natural biophenols (curcumin, dopamine, (-)-epigallocatechin-3-gallate, quercetin, and rosmarinic acid) with amyloid-β(1⁻40) fibrils has been studied through computational simulations. The results allowed the identification and characterization of the different binding modalities of each compounds and their consequences on fibril dynamics and aggregation. It emerges that the lateral aggregation of the fibrils is strongly influenced by the intercalation of the ligands, which modulates the double-layered structure stability.

In this work, the adsorption of cytochrome C (CytC) on wet {100}, {111}, {110}, and {120} silver surfaces has been investigated by computational simulations. The effect of polyvinylpyrrolidone (PVP) coating has also been studied. The main results obtained can be summarized as follow: (a) CytC strongly interacts with wet bare high index facets, while the adsorption over the {100} surface is disfavored due to the strong water structuring at the surface; (b) a nonselective protein adsorption mechanism is highlighted; (c) the native structure of CytC is well preserved during adsorption; (d) the heme group of CytC is never found to interact directly with the surface; (e) the interactions with the PVP-capped {100} surface is weak and specific. These results can be exploited to better control biological responses at engineered nanosurface, allowing the development of improved diagnostic tools.

We demonstrate that the designed ankyrin repeat protein (DARPin)_9-29, which specifically targets human epidermal growth factor receptor 2 (HER 2), binds tightly to gold nanoparticles (GNPs). Binding of the protein strongly increases the colloidal stability of the particles. The results of experimental analysis and molecular dynamics simulations show that approximately 35 DARPin_9-29 molecules are bound to the surface of a 5 nm GNP and that the binding does not involve the receptor-binding domain of the protein. The confocal fluorescent microscopy studies show that the DARPin-coated GNP conjugate specifically interacts with the surface of human cancer cells overexpressing epidermal growth factor receptor 2 (HER2) and enters the cells by endocytosis. The high stability under physiological conditions and high affinity to the receptors overexpressed by cancer cells make conjugates of plasmonic gold nanostructures with DARPin molecules promising candidates for cancer therapy.

Strategies for protein detection via surface-enhanced Raman spectroscopy (SERS) currently exploit the formation of randomly generated hot spots at the interfaces of metal colloidal nanoparticles, which are clustered together by intrusive chemical or physical processes in the presence of the target biomolecule. We propose a different approach based on selective and quantitative gathering of protein molecules at regular hot spots generated on the corners of individual silver nanocubes in aqueous medium at physiological pH. Here, the protein, while keeping its native configuration, experiences an intense local E-field, which boosts SERS efficiency and detection sensitivity. Uncontrolled signal fluctuations caused by variable molecular adsorption to different particle areas or inside clustered nanoparticles are circumvented. Advanced electron microscopy analyses and computational simulations outline a strategy relying on a site-selective mechanism with superior Raman signal enhancement, which offers the perspective of highly controlled and reproducible routine SERS detection of proteins.

Two new fluorophenylindenone derivatives were designed as potential p38α MAPK modulators by preserving the key interactions of the vicinal pyridine/fluorophenyl pharmacophore with the enzyme protein. Interestingly, these two fluorophenylindenone isomers showed divergent activities, with compound 6 behaving as an inhibitor and 5 as a putative activator. These results were rationalized by docking studies and molecular dynamics simulations in terms of stabilization of DFG loop, by compound 5 in a conformation more accessible to phosphorylation., (Copyright © 2016 Elsevier Ltd. All rights reserved.)

⁴³Ca MQMAS NMR spectra of three silica-based glasses in which Ca²⁺ ions play different structural roles have been collected and processed in order to extract the underlying NMR parameter distributions. The NMR parameters have been interpreted with the help of molecular dynamics simulations and DFT-GIPAW calculations. This synergetic experimental-computational approach has allowed us to investigate the Ca environment, to estimate Ca coordination numbers from MD-derived models, and to push further the discussion about ⁴³Ca NMR sensitivity to the first and second coordination spheres: ⁴³Ca δiso and Ca-O distance can be successfully correlated as a function of Ca coordination number., (Copyright © 2015 Elsevier Inc. All rights reserved.)

An extensive benchmark of exchange-correlation functionals on the structure of the X-ray resolved phosphine and thiolate-protected Ag14-based nanocluster, named XMC1, is reported. Calculations were performed both on simplified model systems, with the complexity of the ligands greatly reduced, and on the complete XMC1 particle. Most of the density functionals that yielded good relaxed structures on analogous calculations on gold nanoclusters (viz. those employing the generalized gradient approximation) significantly deform the structure of XMC1. On the contrary, some of the exchange-correlation functionals including part of the exact Hartree-Fock exchange (hybrid functionals) reproduce the experimental geometry with minimal errors. In particular, the widely adopted B3LYP yields fairly accurate structures for XMC1, whereas it is outperformed by many other functionals (both hybrids and generalized gradient corrected) in similar calculations on analogous gold-based systems. Time-dependent density functional calculations have been employed to recover the experimental UV-vis spectrum. The present investigation shows that to correctly reproduce the optical feature of XMC1 the ligands cannot be omitted, because they interact with the metal core at energies much closer to the optical gap than in the case of gold-based nanoclusters of similar sizes. Due to this fact, a functional that accurately describes charge-transfer electronic transitions (such as the long-range corrected CAM-B3LYP) has to be adopted.

The ability of Ce-containing bioactive glasses to inhibit oxidative stress in terms of reduction of hydrogen peroxide, by mimicking the catalase enzyme activity is demonstrated here for the first time. The antioxidant properties of three bioactive glasses containing an increasing amount of CeO2 have been evaluated by following the degradation of hydrogen peroxide with time after immersion in H2O2 aqueous solutions with different concentration. XPS and UV-vis measurements allowed us to determine the Ce(3+)/Ce(4+) ratio in the bulk and on the glass surface, and to correlate it with the ability of the samples to show catalase mimetic activity. Interestingly, we have found that the bioactive glass with composition 23.2Na2O-25.7CaO-43.4SiO2-2.4P2O5-5.3CeO2 immersed in 0.1 M H2O2 aqueous solution is able to degrade 90% of it in 1 week. The reduction in bioactivity of the glasses with increasing CeO2 content is here rationalized in terms of a lower amount of phosphate groups available for the hydroxyapatite layer formation, after binding with cerium ions. In fact, classical molecular dynamics simulations revealed that the addition of CeO2 leads to the formation of cerium phosphate rich regions. The formation of an insoluble CePO4 crystalline phase is also observed by XRD analysis after thermal treatment of the glass samples.

A computational protocol that combines periodic and QM/QM' calculations has been applied to investigate the structural (geometrical and electronic) and photophysical absorption properties of the salicylidene aniline (SA) thermochromic molecular crystal. The protocol consists of three different steps, namely (i) the description of the molecular crystal using a periodic approach taking into account dispersion interactions, (ii) the identification of reliable finite models (clusters), and (iii) the calculation of vertical transition energies including environmental effects through the use of an electronic embedding model (QM/QM' ONIOM approach). The encouraging results obtained in this work for the β polymorph of SA, both in terms of accuracy and computational cost, open the way to the simulation and the prediction of the photophysical behavior of other molecular crystals, especially those much less well characterized experimentally.

We present here a detailed time-dependent density-functional theory investigation aimed at systematically dissecting the electronic spectra of two thiolate and phosphine protected undecagold nanoclusters. Calculations performed on the experimental structures of Au11(PPh3)7Cl3 and Au11(PPh3)7(SPyr)3 show that ligands have negligible contributions in the visible region. Metal → ligand charge transfer transitions appear at energies well above the visible threshold, while transitions with some small ligand → metal and ligand → ligand character occur sporadically at even higher energies. Thus, the conjugation effect between the π-electrons of the ligand and electrons of gold, recently hypothesized to interpret the spectra of phosphine and thiolate-protected nanoclusters, is not confirmed by the results of this study.

Multivalency is widely used in nature in specific recognition processes. This paper describes an approach to multivalency in the pentameric 5-HT3 receptor, a ligand-gated ion channel, which constitutes an example of intrinsically multivalent biological receptors. Owing to the picomolar Ki value, TETRA-L represents an outstanding multivalent ligand for the neurotransmitter receptor.

The structural and (13)C/(1)H NMR parameters of the four crystal forms (1α, 1·H2O, 1β, and 1γ) of the solid wheel-and-axle (WAA) metal-organic compound [(p-cymene)Ru(κN-INA)Cl2] have been studied by means of periodic DFT calculations. The quality of the results obtained strongly depends on a correct description of long-range interactions; thus, in the geometry refinement protocol used, the pure DFT functionals need to be coupled with a dispersion-correction term (B3LYP-D2, B3LYP-D*). The solid-state (13)C/(1)H NMR δ(iso) parameters and (13)C MAS NMR spectra, calculated by means of the PBE-GIPAW method, agree well with the experimental data for the four crystal forms (mean absolute deviations of the (13)C and (1)H δ(iso) data values lie in the ranges 1.3-2.9 and 0.3-1.0 ppm, respectively). In this context, some revisions in the experimental assignment of the (13)C/(1)H NMR δ(iso) parameters of the 1·H2O, 1β, and 1γ crystal forms can be suggested. The mismatch in the assignment seems to be due to the rotation of the -COOH moiety, which occurs at the 1α-1·H2O transition and was not considered in the experiments. Finally, the results obtained suggest the presence of two COOH···Cl hydrogen bonds of comparable strength established by the two molecules in the asymmetric unit of the 1γ polymorph, in partial disagreement with previous findings.

A small series of serotonin 5-HT4 receptor ligands has been designed from flexible 2-methoxyquinoline compounds 7a,b by applying the conformational constraint approach. Ligands 7a,b and the corresponding conformationally constrained analogues 8a-g were synthesized and their interactions with the 5-HT4 receptor were examined by measuring both binding affinity and the ability to promote or inhibit receptor-G protein coupling. Ester derivative 7a and conformationally constrained compound 8b were demonstrated to be the most interesting compounds showing a nanomolar 5-HT4R affinity similar to that shown by reference ligands cisapride (1) and RS-23,597-190 (4). The result was rationalized by docking studies in term of high similarity in the binding modalities of flexible 7a and conformationally constrained 8b. The intrinsic efficacy of some selected ligands was determined by evaluating the receptor-G protein coupling and the results obtained demonstrated that the nature and the position of substituents play a critical role in the interaction of these ligands with their receptor., (Copyright © 2014 Elsevier Masson SAS. All rights reserved.)

Palonosetron (Aloxi) is a potent second generation 5-HT(3) receptor antagonist whose mechanism of action is not yet fully understood. Palonosetron acts at the 5-HT(3) receptor binding site but recent computational studies indicated other possible sites of action in the extracellular domain. To test this hypothesis we mutated a series of residues in the 5-HT3A receptor subunit (Tyr(73), Phe(130), Ser(163), and Asp(165)) and in the 5-HT3B receptor subunit (His(73), Phe(130), Glu(170), and Tyr(143)) that were previously predicted by in silico docking studies to interact with palonosetron. Homomeric (5-HT(3)A) and heteromeric (5-HT(3)AB) receptors were then expressed in HEK293 cells to determine the potency of palonosetron using both fluorimetric and radioligand methods to test function and ligand binding, respectively. The data show that the substitutions have little or no effect on palonosetron inhibition of 5-HT-evoked responses or binding. In contrast, substitutions in the orthosteric binding site abolish palonosetron binding. Overall, the data support a binding site for palonosetron at the classic orthosteric binding pocket between two 5-HT3A receptor subunits but not at allosteric sites previously identified by in silico modelling and docking., (Copyright © 2013 The Authors. Published by Elsevier Ltd.. All rights reserved.)

5-hydroxytryptamine type-3 receptor (5-HT₃), an important target of many neuroactive drugs, is a cation selective transmembrane pentamer whose functional stoichiometries and subunit arrangements are still debated, due to the extreme complexity of the system. The three dimensional structure of the 5-HT₃R subunits has not been solved so far. Moreover, most of the available structural and functional data is related to the extracellular ligand-binding domain, whereas the transmembrane and the intracellular receptor domains are far less characterised, although they are crucial for receptor function. Here, for the first time, 3D homology models of the transmembrane and the intracellular receptor domains of all the known human 5-HT₃ subunits have been built and assembled into homopentameric (5-HT(3A)R, 5-HT(3B)R, 5-HT(3C)R, 5-HT(3D)R and 5-HT(3E)R) and heteropentameric receptors (5-HT(3AB), 5-HT(3AC), 5-HT(3AD) and 5-HT(3AE)), on the basis of the known three-dimensional structures of the nicotinic-acetylcholine receptor and of the ligand gated channel from Erwinia chrysanthemi. The comparative analyses of sequences, modelled structures, and computed electrostatic properties of the single subunits and of the assembled pentamers shed new light both on the stoichiometric composition and on the physicochemical requirements of the functional receptors. In particular, it emerges that a favourable environment for the crossing of the pore at the transmembrane and intracellular C terminus domain levels by Ca²⁺ ions is granted by the maximum presence of two B subunits in the 5-HT₃ pentamer.

The structural properties of phosphosilicate glasses based on the 45S5 Bioglass doped with gallium and aluminum (46.2 SiO2·24.3Na2O·26.9CaO·2.6P2O5·1.0X2O3, X = Ga or Al) are investigated by means of classical molecular dynamics simulations. Structural features of the two compositions are compared with those of the original 45S5 Bioglass in order to relate them to the different known bioactivities of these materials. Differences in the coordination environments of Ga and Al, network connectivity, and ion aggregation reveal a microscopic model of these glasses which supports the interpretation of the experimental data and provides new insight into the different biological behaviors of Ga- and Al-containing phosphosilicate glasses. Although Ga is found predominantly in a 4-fold coordination environment, small amounts of 5- and 6-fold coordinated atoms have been detected depending on the interatomic potential model employed. This suggests its possible intermediate role in phosphosilicate glasses. On the contrary, Al plays a network former role and leads to glasses with a more polymerized structure. Interestingly, the results show an increased propensity for aggregation of the Ca(2+) and PO4(3-) ions in the Al-containing phosphosilicate glasses with respect to the Ga-containing ones. This leads to insoluble calcium-phosphate-rich regions not detected in the bioactive glasses.

A new practical route to chaetomellic acid A (ACA), based on the copper catalysed radical cyclization (RC) of (Z)-3-(2,2-dichloropropanoyl)-2-pentadecylidene-1,3-thiazinane, is described. Remarkably, the process entailed: (i) a one-pot preparation of the intermediate N-α-perchloroacyl-2-(Z)-alkyliden-1,3-thiazinanes starting from N-(3-hydroxypropyl)palmitamide, (ii) a two step smooth transformation of the RC products into ACA and (iii) only one intermediate chromatographic purification step. The method offers a versatile approach to the preparation of ACA analogues, through the synthesis of an intermediate maleic anhydride with a vinylic group at the end of the aliphatic tail, a function that can be transformed through a thiol-ene coupling. Serendipitously, the disodium salt of 2-(9-(butylthio)nonyl)-3-methylmaleic acid, that we prepared as a representative sulfurated ACA analogue, was a more competent FTase inhibitor than ACA. This behaviour was analysed by a molecular docking study., (Copyright © 2012 Elsevier Ltd. All rights reserved.)

The functional serotonin type-3 receptor (5-HT(3)-R), which is the target of many neuroactive drugs, is known to be a homopentamer made of five identical subunits A (5-HT(3A)-R) or a binary heteropentamer made of subunits A and B (5-HT(3A/B)-R) with a still debated arrangement and stoichiometry. This complex picture has been recently further complicated by the discovery of additional 5-HT(3)-R subunits, C, D, and E, which, similarly to the B subunit, are apparently able to form functional receptors only if co-expressed with subunit A. Being the binding site for both serotonin and antagonists (i.e. drugs) located at the extracellular interface between two adjacent subunits, the large variability of the 5-HT(3)-R composition becomes a crucial issue, since it can originate many different interfaces providing non-equivalent ligand binding sites and complicating the pharmacological modulation. Here, the different 5-HT(3)-R interfaces are analysed, on the bases of the structural conformations of previously built 3D homology models and of the known subunit sequences, by addressing their physicochemical characterization. The results confirm the presence of an aromatic cluster located in the core of the A-A interface as a key determinant for having an interface both stable and functional. This is used as a discriminant to make hypotheses about the capability of all the other possible interfaces constituted by the known 5-HT(3)-R sequences A, B, C, D, and E to build active receptors.

The p38 mitogen-activated protein kinase is activated by environmental stress and cytokines and plays a role in transcriptional regulation and inflammatory responses. Factors influencing the activity and selectivity of the p38α mitogen-activated protein kinase inhibitors have been investigated in this paper by inspecting the binding orientation and the possible residue-inhibitor interactions in the binding site. The binding pattern of a set of 45 different inhibitors against p38α mitogen-activated protein kinase was studied through Molecular Dynamic Simulations of the protein-inhibitor complexes. Further, Partial Least Squares regression was used to develop a Quantitative Structure Activity Relationship model to predict the binding affinities of ligands. The selected model successfully predicted the test set with a Root Mean Square Error of Prediction of 1.36. The regression coefficients and the Variable Importance in Projection plots highlighted the residue-inhibitor interactions which exhibited the largest absolute effect on the ligand binding, such as the van der Waals interaction with LYS50, ILE81, ASP165; electrostatic interactions with SER29, LEU164; hydrogen bonds with MET106; and total energy interaction with SER29 and LEU83., (© 2012 John Wiley & Sons A/S.)

The results obtained by means of Molecular Dynamics simulations and Multiway Explorative Data Analysis on ADAMTS4, ADAMTS5 and MMP13 complexed with Marimastat and two cis-1(S)2(R)-amino-2-indanol ligands suggest that determinant characteristics for ligand binding and selectivity among the three enzymes are to be found in the different protein conformation flexibility. Moreover, the role of the TS-domain in the inhibitor binding to ADAMTS enzymes has been investigated for the first time in this work. The results obtained suggest that the influence of the TS-domain on the S1' loop fluctuations of ADAMTS4 and ADAMTS5 could be exploited for the design of therapeutics for chronic osteoarthritis diseases., (Copyright © 2012 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.)

The functional serotonin 5-HT type-3 (5-HT(3)) receptor, the target of many neuroactive drugs, is known to be a pseudo-symmetric pentamer made either of five identical subunits A (homomeric 5-HT(3A)-R) or of subunits A and B (heteromeric 5-HT(3A/B)-R) in a still debated arrangement. The serotonin binding site is located in the extracellular region, at the interface between two monomers, called the principal and the complementary subunits. The results of molecular dynamics simulations and computational alanine scanning mutagenesis studies applied here to the homomeric human 5-HT(3A)-R disclose an aromatic "hot" cluster in the centre of the interface formed by residues W178 (principal subunit), Y68, Y83, W85 and Y148 (complementary subunit). Moreover, investigation of the coupling of agonist/antagonist binding to channel activation/inactivation points out the presence of two putative functional pathways at the subunit interface: W116-H180-L179-W178-E124-F125 (principal subunit) and Y136-Y138-Y148-W85-(P150) (complementary subunit), where W178 and Y148 appear to be critical residues for the binding/activation mechanism. Finally, direct comparison of the main features shown by the AA interface in the human 5-HT(3A)-R with those of the BB interface in the homopentameric human 5-HT(3B)-R provides interesting clues about the possible reasons that cause the 5-HT(3B)-R not to be functional.

The serotonin 5-HT3 receptor is a ligand-gated ion channel, which by virtue of its pentameric architecture, can be considered to be an intriguing example of intrinsically multivalent biological receptors. This paper describes a general design approach to the study of multivalency in this multimeric ion channel. Bivalent ligands for 5-HT3 receptor have been designed by linking an arylpiperazine moiety to probes showing different functional features. Both homobivalent and heterobivalent ligands have shown 5-HT3 receptor affinity in the nanomolar range, providing evidence for the viability of our design approach. Moreover, the high affinity shown by homobivalent ligands suggests that bivalency is a promising approach in 5-HT3 receptor modulation and provides the rational basis for applying the concepts of multivalency to the study of 5-HT3 receptor function.