Search

Multilayer vesicles with an onion-like architecture can form by self-assembly of organic amphiphiles such as dendrimers, small-molecule surfactants, and block copolymers. Thus far, there are limited reports about multilayer vesicles based on coordination compounds. Herein, we show that nano onions are obtained by aggregation of an amphiphilic Au 3 (pyrazolate) 3 complex in aqueous solution. The nanostructures were characterized by cryogenic and transition electron microscopy, dynamic light scattering, and energy-dispersive X-ray analysis. Control experiments with analogous Ag(I) and Cu(I) complexes revealed the importance of Au(I) for the formation of well-defined nano onions. A structurally related Au(I) complex without solubilizing polyethylene glycol side chains was analyzed by single-crystal X-ray diffraction. In the solid state, columns of offset, π-stacked Au 3 (pyrazolate) 3 complexes are observed, but short intermolecular Au⋯Au contacts were not found. Molecular dynamics simulations provided further insights into the aggregation process in aqueous solution, supporting the formation of nano-onion structures through lateral interactions between stacked complexes.

Even at low temperatures, metal nanoparticles (NPs) possess atomic dynamics that are key for their properties but challenging to elucidate. Recent experimental advances allow obtaining atomic-resolution snapshots of the NPs in realistic regimes, but data acquisition limitations hinder the experimental reconstruction of the atomic dynamics present within them. Molecular simulations have the advantage that these allow directly tracking the motion of atoms over time. However, these typically start from ideal/perfect NP structures and, suffering from sampling limits, provide results that are often dependent on the initial/putative structure and remain purely indicative. Here, by combining state-of-the-art experimental and computational approaches, how it is possible to tackle the limitations of both approaches and resolve the atomistic dynamics present in metal NPs in realistic conditions is demonstrated. Annular dark-field scanning transmission electron microscopy enables the acquisition of ten high-resolution images of an Au NP at intervals of 0.6 s. These are used to reconstruct atomistic 3D models of the real NP used to run ten independent molecular dynamics simulations. Machine learning analyses of the simulation trajectories allow resolving the real-time atomic dynamics present within the NP. This provides a robust combined experimental/computational approach to characterize the structural dynamics of metal NPs in realistic conditions., (© 2024 The Authors. Advanced Science published by Wiley‐VCH GmbH.)

In the early Earth, rudimentary enzymes must have utilized the available light energy source to modulate protometabolic processes. Herein, we report the light-responsive C-C bond manipulation via short peptide-based assemblies bound to the photosensitive molecular cofactor (azo-based photoswitch) where the energy of the light source regulated the binding sites which subsequently modulated the retro-aldolase activity. In the presence of a continual source of high-energy photons, temporal realization of a catalytically more proficient state could be achieved under nonequilibrium conditions. Further, the hydrophobic surface of peptide assemblies facilitated the binding of an orthogonal molecular catalyst that showed augmented activity (promiscuous hydrolytic activity) upon binding. This latent activity was utilized for the in situ generation of light-sensitive cofactor that subsequently modulated the retro-aldolase activity, thus creating a reaction network.

The most advanced structure prediction methods are powerless in exploring the conformational ensemble of disordered peptides and proteins and for this reason the "protein folding problem" remains unsolved. We present a novel methodology that enables the accurate prediction of spectroscopic fingerprints (circular dichroism, infrared, Raman, and Raman optical activity), and by this allows for "tidying up" the conformational ensembles of disordered peptides and disordered regions in proteins. This concept is elaborated for and applied to a dodecapeptide, whose spectroscopic fingerprint is measured and theoretically predicted by means of enhanced-sampling molecular dynamics coupled with quantum mechanical calculations. Following this approach, we demonstrate that peptides lacking a clear propensity for ordered secondary-structure motifs are not randomly, but only conditionally disordered. This means that their conformational landscape, or phase-space, can be well represented by a basis-set of conformers including about 10 to 100 structures. The implications of this finding have profound consequences both for the interpretation of experimental electronic and vibrational spectral features of peptides in solution and for the theoretical prediction of these features using accurate and computationally expensive techniques. The here-derived methods and conclusions are expected to fundamentally impact the rationalization of so-far elusive structure-spectra relationships for disordered peptides and proteins, towards improved and versatile structure prediction methods., Competing Interests: There are no conflicts to declare., (This journal is © The Royal Society of Chemistry.)

It is known that metal nanoparticles (NPs) may be dynamic and atoms may move within them even at fairly low temperatures. Characterizing such complex dynamics is key for understanding NPs' properties in realistic regimes, but detailed information on, e.g., the stability, survival, and interconversion rates of the atomic environments (AEs) populating them are non-trivial to attain. In this study, we decode the intricate atomic dynamics of metal NPs by using a machine learning approach analyzing high-dimensional data obtained from molecular dynamics simulations. Using different-shape gold NPs as a representative example, an AEs' dictionary allows us to label step-by-step the individual atoms in the NPs, identifying the native and non-native AEs and populating them along the MD simulations at various temperatures. By tracking the emergence, annihilation, lifetime, and dynamic interconversion of the AEs, our approach permits estimating a "statistical equivalent identity" for metal NPs, providing a comprehensive picture of the intrinsic atomic dynamics that shape their properties., (© 2023. The Author(s).)

Metals are traditionally considered hard matter. However, it is well known that their atomic lattices may become dynamic and undergo reconfigurations even well below the melting temperature. The innate atomic dynamics of metals is directly related to their bulk and surface properties. Understanding their complex structural dynamics is, thus, important for many applications but is not easy. Here, we report deep-potential molecular dynamics simulations allowing to resolve at an atomic resolution the complex dynamics of various types of copper (Cu) surfaces, used as an example, near the Hüttig (∼1/3 of melting) temperature. The development of deep neural network potential trained on density functional theory calculations provides a dynamically accurate force field that we use to simulate large atomistic models of different Cu surface types. A combination of high-dimensional structural descriptors and unsupervized machine learning allows identifying and tracking all the atomic environments (AEs) emerging in the surfaces at finite temperatures. We can directly observe how AEs that are non-native in a specific (ideal) surface, but that are, instead, typical of other surface types, continuously emerge/disappear in that surface in relevant regimes in dynamic equilibrium with the native ones. Our analyses allow estimating the lifetime of all the AEs populating these Cu surfaces and to reconstruct their dynamic interconversions networks. This reveals the elusive identity of these metal surfaces, which preserve their identity only in part and in part transform into something else under relevant conditions. This also proposes a concept of "statistical identity" for metal surfaces, which is key to understanding their behaviors and properties.

Spatial confinement is widely employed by nature to attain unique efficiency in controlling chemical reactions. Notable examples are enzymes, which selectively bind reactants and exquisitely regulate their conversion into products. In an attempt to mimic natural catalytic systems, supramolecular metal-organic cages capable of encapsulating guests in their cavity and of controlling/accelerating chemical reactions under confinement are attracting increasing interest. However, the complex nature of these systems, where reactants/products continuously exchange in-and-out of the host, makes it often difficult to elucidate the factors controlling the reactivity in dynamic regimes. As a case study, here we focus on a coordination cage that can encapsulate amide guests and enhance their hydrolysis by favoring their mechanical twisting towards reactive molecular configurations under confinement. We designed an advanced multiscale simulation approach that allows us to reconstruct the reactivity in such host-guest systems in dynamic regimes. In this way, we can characterize amide encapsulation/expulsion in/out of the cage cavity (thermodynamics and kinetics), coupling such host-guest dynamic equilibrium with characteristic hydrolysis reaction constants. All computed kinetic/thermodynamic data are then combined, obtaining a statistical estimation of reaction acceleration in the host-guest system that is found in optimal agreement with the available experimental trends. This shows how, to understand the key factors controlling accelerations/variations in the reaction under confinement, it is necessary to take into account all dynamic processes that occur as intimately entangled in such host-guest systems. This also provides us with a flexible computational framework, useful to build structure-dynamics-property relationships for a variety of reactive host-guest systems., Competing Interests: There are no conflicts to declare., (This journal is © The Royal Society of Chemistry.)

This work reports for the first time a quantum mechanical study of the interactions of a model benzodiazepine drug, i.e., nitrazepam, with various models of amorphous silica surfaces, differing in structural and interface properties. The interest in these systems is related to the use of mesoporous silica as carrier in drug delivery. The adopted computational procedure has been chosen to investigate whether silica-drug interactions favor the drug degradation mechanism or not, hindering the beneficial pharmaceutical effect. Computed structural, energetics, and vibrational properties represent a relevant comparison for future experiments. Our simulations demonstrate that adsorption of nitrazepam on amorphous silica is a strongly exothermic process in which a partial proton transfer from the surface to the drug is observed, highlighting a possible catalytic role of silica in the degradation reaction of benzodiazepines.

Melatonin (MT) is a molecule of paramount importance in all living organisms, due to its presence in many biological activities, such as circadian (sleep-wake cycle) and seasonal rhythms (reproduction, fattening, molting, etc.). Unfortunately, it suffers from poor solubility and, to be used as a drug, an appropriate transport vehicle has to be developed, in order to optimize its release in the human tissues. As a possible drug-delivery system, β-cyclodextrin (βCD) represents a promising scaffold which can encapsulate the melatonin, releasing when needed. In this work, we present a computational study supported by experimental IR spectra on inclusion MT/βCD complexes. The aim is to provide a robust, accurate and, at the same time, low-cost methodology to investigate these inclusion complexes both with static and dynamic simulations, in order to study the main actors that drive the interactions of melatonin with β-cyclodextrin and, therefore, to understand its release mechanism.

We undertake steps to overcome four challenges that have hindered the understanding of ZnO/biomolecule interfaces at the atomic scale: parametrization of a classical force field, ZnO surface termination and amino acid protonation state in methanol, and convergence of enhanced sampling molecular dynamics simulations. We predict adsorption free energies for histidine, serine, cysteine, and tryptophan in remarkable agreement with experimental measurements obtained via a novel indicator-displacement assay. Adsorption is driven by direct surface/amino-acid interactions mediated by terminal hydroxyl groups and stabilized by strongly structured methanol solvation shells.

Although successful in the structural determination of ordered biomolecules, the spectroscopic investigation of oligopeptides in solution is hindered by their complex and rapidly changing conformational ensemble. The measured circular dichroism (CD) spectrum of an oligopeptide is an ensemble average over all microstates, severely limiting its interpretation, in contrast to ordered biomolecules. Spectral deconvolution methods to estimate the secondary structure contributions in the ensemble are still mostly based on databases of larger ordered proteins. Here, we establish how the interpretation of CD spectra of oligopeptides can be enhanced by the ability to compute the same observable from a set of atomic coordinates. Focusing on two representative oligopeptides featuring a known propensity toward an α-helical and β-hairpin motif, respectively, we compare and cross-validate the structural information coming from deconvolution of the experimental CD spectra, sequence-based de novo structure prediction, and molecular dynamics simulations based on enhanced sampling methods. We find that small conformational variations can give rise to significant changes in the CD signals. While for the simpler conformational landscape of the α-helical peptide de novo structure prediction can already give reasonable agreement with the experiment, an extended ensemble of conformers needs to be considered for the β-hairpin sequence.

Adsorption of enzymes on solid surfaces may lead to conformational changes that reduce their catalytic conversion activity and are thus detrimental to the efficiency of biotechnology or biosensing applications. This work is a joint theoretical and experimental endeavor in which we identify and quantify the conformational changes that chymotrypsin undergoes when in contact with the surface of amorphous silica nanoparticles. For this purpose, we use circular dichroism spectroscopy, standard molecular dynamics, and advanced-sampling methods. Only the combination of these techniques allowed us to pinpoint a destabilization effect of silica on specific structural motifs of chymotrypsin. They are linked by the possibility of theoretically predicting CD spectra, allowing us to elucidate the source of the experimentally observed spectral changes. We find that chymotrypsin loses part of its helical content upon adsorption, with minor perturbation of its overall tertiary structure, associated with changes in the aromatic interactions. We demonstrate that the C-terminal helical fragment is unfolded as an isolated oligopeptide in pure water, folded as an α-helix as terminus of chymotrypsin in solution, and again partly disordered when the protein is adsorbed on silica. We believe that the joint methodology introduced in this manuscript has a direct general applicability to investigate any biomolecule-inorganic surface system. Methods to theoretically predict circular dichroism spectra from atomistic simulations were compared and improved. The drawbacks of the approaches are discussed; in particular, the limited capability of advanced-sampling MD schemes to explore the conformational phase space of large proteins and the dependency of the predicted ellipticity bands on the choice of calculation parameters.

Li 7 La 3 Zr 2 O 12 (LLZO) and related compounds are considered as promising candidates for future all-solid-state Li-ion battery applications. Still, the processing of those materials into thin membranes with the right stoichiometry and crystal structure is difficult and laborious. The sensitivity of the Li-ion conductive garnets against moisture and the associated Li + /H + cation exchange makes their processing even more difficult. Formulation of suitable polymer/ceramic hybrid solid state electrolytes could be a prosperous way to reach the future large scale production of solid state Li-ion batteries. In fact, solvent mediated and/or slurry based wet-processing of the LLZO, e.g., tape-casting, could result in irreversible Li-ion loss of the pristine material due to Li + /H + cation exchange. The concomitant structural changes and loss in functionality in terms of Li-ion conductivity are the results of the above process. Therefore, in the present work a systematic study on the chemical stability and structural retention of Al-substituted LLZO in different solvents is reported. It was found that Li + /H + exchange in LLZO occurs upon solvent immersion, and its magnitude is dependent on the availability of -OH functional groups of the solvent molecules. As a result, a larger degree of Li + /H + exchange causes higher increase of the lattice parameter of the LLZO, determined by synchrotron diffraction analyses. The expansion of the cubic unit cell was ascertained, when Li + was replaced by H + in the host lattice, by ab initio computational studies. The application of the most common solvent as dispersion medium, i.e., high purity water, causes the most significant Li + /H + exchange and, therefore, structural change, while acetonitrile was proven to be the best suitable solvent for wet postprocessing of LLZO. Finally, computational calculations suggested that the Li + /H + exchange could result in diminished ionic, i.e., mixed Li + -H + , conductivity due to the insertion of protons with lower mobility than that of Li-ions.

Objectives: To estimate the incidence of acute rheumatic fever (ARF) in a metropolitan area of Northern Italy and study how the introduction of the 2015 revised Jones criteria affects the epidemiology in a region with moderate to high incidence of ARF., Study Design: The incidence of ARF in children 5-14 years old living in the Province of Turin was estimated using low-risk criteria in a 10-year period (group A patients). The proportion of patients fulfilling only high-risk (HR) criteria (group B patients) was also calculated both prospectively (from July 2015 through December 2016) and retrospectively (from January 2007 through June 2015)., Results: One hundred thirty-five group A patients were identified for an annual incidence of 3.2-9.6 out of 100 000 children. The use of HR criteria identified an additional 28 patients (group B), resulting in a 20.7% increase in the incidence of ARF. Age, sex annual incidence, and seasonal distribution pattern were comparable between group A and group B patients., Conclusions: HR criteria should be used for the diagnosis ARF in our region. The application of these criteria led to a 20% increase in patients with the diagnosis of ARF. The characteristics of patients fulfilling only HR criteria are similar to the remaining patients, suggesting that these criteria are sensitive and specific., (Copyright © 2018 Elsevier Inc. All rights reserved.)

A density functional theory (PBE functional) investigation is carried out, in which a model of an amorphous silica surface is functionalized by ortho-benzoquinone. Surface functionalization with catechol and quinone-based compounds is relevant in biomedical fields, from prosthetic implants to dentistry, to develop multifunctional coatings with antimicrobial properties. The present study provides atomistic information on the specific interactions between the functionalizing agent and the silanol groups at the silica surface. The distinct configurations of the functional groups, the hydrogen bond pattern, the role of dispersion forces and the simulated IR spectra provide detailed insight into the features of this model surface coating. Ab initio molecular dynamics gives further insights into the mobility of the functionalizing groups. As a final step, we studied the condensation reaction with allylamine, via Schiff base formation, to ground subsequent simulations on condensation with model peptides of antimicrobial activity.

Molecular simulations of proteins have been usually accomplished through empirical or semi-empirical potentials, due to the large size and inherent complexity of these biological systems. On the other hand, a theoretical description of proteins based on quantum-mechanical methods would however provide an unbiased characterization of their electronic properties, possibly offering a link between these and the ultimate biological activity. Yet, such approaches have been historically hindered by the large amount of requested computational power. Here we demonstrate the feasibility of periodic all-electron density functional theory calculations in the description of the crystal of the protein crambin (46 aminoacids), which is determined with exceptional structural accuracy. We have employed the hybrid B3LYP functional, coupled to an empirical description of London interactions (D*) to simulate the crambin crystal with an increasing amount of lattice water molecules in the cell (up to 172H 2 O per cell). The agreement with the experiment is good for both protein geometry and protein-water interactions. The energetics was computed to predict crystal formation energies, protein-water and protein-protein interaction energies. We studied the role of dispersion interactions which are crucial for holding the crambin crystal in place. B3LYP-D* electrostatic potential and dipole moment of crambin as well as the electronic charge flow from crambin to the solvating water molecules (0.0015 e per H 2 O) have also been predicted. These results proved that quantum-mechanical simulations of small proteins, both free and in their crystalline state, are now feasible in a reasonable amount of time, by programs capable of exploiting high performance computing architectures, allowing the study of protein properties not easily amenable through classical force fields.

In silico modeling of acidic (CH2COOH) or basic (CH2NH2) functionalized silica surfaces has been carried out by means of a density functional approach based on a gradient-corrected functional to provide insight into the characterization of experimentally functionalized surfaces via a plasma method. Hydroxylated surfaces of crystalline cristobalite (sporting 4.8 OH/nm(2)) mimic an amorphous silica interface as unsubstituted material. To functionalize the silica surface we transformed the surface Si-OH groups into Si-CH2COOH and Si-CH2NH2 moieties to represent acidic/basic chemical character for the substitution. Structures, energetics, electronic, and vibrational properties were computed and compared as a function of the increasing loading of the functional groups (from 1 to 4 per surface unit cell). Classical molecular dynamics simulations of selected cases have been performed through a Reax-FF reactive force field to assess the mobility of the surface added chains. Both DFT and force field calculations identify the CH2NH2 moderate surface loading (1 group per unit cell) as the most stable functionalization, at variance with the case of the CH2COOH group, where higher loadings are preferred (2 groups per unit cell). The vibrational fingerprints of the surface functionalities, which are the ν(C═O) stretching and δ(NH2) bending modes for acidic/basic cases, have been characterized as a function of substitution percentage in order to guide the assignment of the experimental data. The final results highlighted the different behavior of the two types of functionalization. On the one hand, the frequency associated with the ν(C═O) mode shifts to lower wavenumbers as a function of the H-bond strength between the surface functionalities (both COOH and SiOH groups), and on the other hand, the δ(NH2) frequency shift seems to be caused by a subtle balance between the H-bond donor and acceptor abilities of the NH2 moiety. Both sets of data are in general agreement with experimental measurements on the corresponding silica-functionalized materials and provide finer details for a deeper interpretation of experimental spectra.

Silica-based materials find applications as excipients and, particularly for those of mesoporous nature, as drug delivery agents for pharmaceutical formulations. Their performance can be crucially affected by water moisture, as it can modify the behavior of these formulations, by limiting their shelf life. Here we describe the role of water microsolvation on the features of ibuprofen adsorbed on a model of amorphous silica surface by means of density functional theory (DFT) simulations. Starting from the results of the simulation of ibuprofen in interaction with a dry hydrophobic amorphous silica surface, a limited number of water molecules has been added to study the configurational landscape of the microsolvated system. Structural and energetics properties, as well as the role of dispersive forces, have been investigated. Our simulations have revealed that the silica surface exhibits a higher affinity for water than for ibuprofen, even if several structures coexist at room temperature, with an active competition of ibuprofen and water for the exposed surface silanols. Dispersive interactions play a key role in this system, as pure DFT fails to correctly describe its potential energy surface. Indeed, van der Waals forces are the leading contribution to adsorption, independently of whether the drug is hydrogen-bonded directly to the surface or via water molecules.

Amorphous silica is widely employed in pharmaceutical formulations both as a tableting, anticaking agent and as a drug delivery system, whereas MCM-41 mesoporous silica has been recently proposed as an efficient support for the controlled release of drugs. Notwithstanding the relevance of this topic, the atomistic details about the specific interactions between the surfaces of the above materials and drugs and the energetic of adsorption are almost unknown. In this work, we resort to a computational ab initio approach, based on periodic Density Functional Theory (DFT), to study the adsorption behavior of two popular drugs (aspirin and ibuprofen) on two models of an amorphous silica surface characterized by different hydrophilic/hydrophobic properties due to different SiOH surface groups' density. Particular effort was devoted to understand the role of dispersive (vdW) interactions in the adsorption mechanism and their interplay with H-bond interactions. On the hydrophilic silica surface, the H-bond pattern of the Si-OH groups rearranges to comply with the formation of new H-bond interactions triggered by the adsorbed drug. The interaction energy of ibuprofen with the hydrophilic model of the silica surface is computed to be very close to the sublimation energy of the ibuprofen molecular crystal, accounting for the experimental evidence of ibuprofen crystal amorphization induced by the contact with the mesoporous silica material. For both surface models, dispersion interactions play a crucial role in dictating the features of the drug/silica system, and they become dominant for the hydrophobic surface. It was proved that a competition may exist between directional H-bonds and nonspecific dispersion driven interactions, with important structural and energetic consequences for the adsorption. The results of this work emphasize the inadequacy of plain DFT methods to model adsorption processes involving inorganic surfaces and drugs of moderate size, due to the missing term accounting for London dispersion interactions.

Fully ab initio treatment of complex solid systems needs computational software which is able to efficiently take advantage of the growing power of high performance computing (HPC) architectures. Recent improvements in CRYSTAL, a periodic ab initio code that uses a Gaussian basis set, allows treatment of very large unit cells for crystalline systems on HPC architectures with high parallel efficiency in terms of running time and memory requirements. The latter is a crucial point, due to the trend toward architectures relying on a very high number of cores with associated relatively low memory availability. An exhaustive performance analysis shows that density functional calculations, based on a hybrid functional, of low-symmetry systems containing up to 100,000 atomic orbitals and 8000 atoms are feasible on the most advanced HPC architectures available to European researchers today, using thousands of processors., (Copyright © 2012 Wiley Periodicals, Inc.)

Background: For many years surgery was the cornerstone of treatment for head and neck cancers and radiotherapy was the treatment of choice in adjuvant and advanced inoperable settings. Recently, induction sequential chemotherapy followed by radiotherapy has shown good tolerability and has prolonged the median overall survival. This phase II trial explored the feasibility of the concurrent association with radiotherapy of a full-dose chemotherapy based on an original schedule of docetaxel and cisplatin., Patients and Methods: Twenty-four patients with head and neck squamous cell carcinoma (HNSCC) were enrolled. Taxotere (docetaxel) was administered on day 1, weekly for 6 weeks. The dose was 33 mg/m2 /w. Cisplatin was administered on day 2 at the dose of 70 mg/m2. Radiotherapy delivered was 60 Gy divided in 30 administrations over 6 weeks., Results and Conclusion: This schedule of treatment for HNSCC proved feasible. Appropriate support treatment, however, appears to be necessary for the feasibility of this concurrent chemo-radiotherapy.

Usually head and neck cancer is treated with combined therapy, applying surgery, if possible, and then radiotherapy and chemotherapy in a sequential or concomitant way. Sequential approach seems to be preferred, because of the high toxicity rate of concomitant therapy. Platinum compounds and 5-fluorouracil are the standard drugs, but new drugs are entering therapeutic arena: gemcitabine and taxanes are the most promising ones. The efficacy of these drugs, especially in association with radiotherapy, must be assessed; moreover it is essential to ascertain how to associate these drugs to radiotherapy and to evaluate drug toxicity when combined with the latter. End point of the study here presented is a preliminar assessment of toxicity and feasibility of concurrent radio-chemoterapy with docetaxel and cisplatinum in patients with head and neck cancer. The number of enrolled patients and the relatively short time of follow up do not allow to evaluate treatment efficacy.