Multiscale Experimental and Theoretical Investigations of Spin Crossover FeII Complexes: Examples of [Fe(phen)2(NCS)2] and [Fe(PM-BiA)2(NCS)2]

For spin crossover (SCO) complexes, computation results are reported and confirmed with experiments at multiscale levels of the isolated molecule and extended solid on the one hand and theory on the other hand. The SCO phenomenon which characterizes organometallics based on divalent iron in an octahedral FeN₆-like environment with high spin (HS) and low spin (LS) states involves the LS/HS switching at the cost of small energies provided by temperature, pressure or light, the latter connected with Light-Induced Excited Spin-State Trapping (LIESST) process. Characteristic infra red (IR) and Raman vibration frequencies are computed within density functional theory (DFT) framework. In [Fe(phen)₂(NCS)₂] a connection of selected frequencies is established with an ultra-fast light-induced LS → HS photoswitching mechanism. In the extended solid, density of state DOS and electron localization function (ELF) are established for both LS and HS forms, leading to characterizion of the compound as an insulator in both spin states with larger gaps for LS configuration, while keeping molecular features in the solid. In [Fe(PM-BiA)₂(NCS)₂], by combining DFT and classical molecular dynamics, the properties and the domains of existence of the different phases are obtained by expressing the potential energy surfaces in a short range potential for Fe–N interactions. Applying such Fe–N potentials inserted in a classical force field and carrying out molecular dynamics (MD) in so-called “semi-classical MD” calculations, lead to the relative energies of HS/LS configurations of the crystal and to the assessment of the experimental (P, T) phase diagram.     

Well-ordered polymer nano-fibers with self-cleaning property by disturbing crystallization process

Bionic self-cleaning surfaces with well-ordered polymer nano-fibers are firstly fabricated by disturbing crystallization during one-step coating-curing process. Orderly thin (100 nm) and long (5–10 μm) polymer nano-fibers with a certain direction are fabricated by external macroscopic force (F_blow) interference introduced by H₂ gas flow, leading to superior superhydrophobicity with a water contact angle (WCA) of 170° and a water sliding angle (WSA) of 0-1°. In contrast, nano-wires and nano-bridges (1–8 μm in length/10-80 nm in width) are generated by “spinning/stretching” under internal microscopic force (F_T) interference due to significant temperature difference in the non-uniform cooling medium. The findings provide a novel theoretical basis for controllable polymer “bionic lotus” surface and will further promote practical application in many engineering fields such as drag-reduction and anti-icing.     

Conformational Dynamics of the Soluble and Membrane-Bound Forms of Interleukin-1 Receptor Type-1: Insights into Linker Flexibility and Domain Orientation

Interleukin-1 receptor type 1 (IL-1R1) is a key player in inflammation and immune responses. This receptor regulates IL-1 activity in two forms: as a membrane-bound form and as a soluble ectodomain. The details and differences between the conformational dynamics of the membrane-bound and the soluble IL-1R1 ectodomains (ECDs) remain largely elusive. Here, we study and compare the structural dynamics of the soluble and membrane-bound IL-1R1-ECDs using molecular dynamics (MD) simulations, focusing on the flexible interdomain linker of the ECD, as well as the spatial rearrangements between the Ig-like domains of the ECD. To explore the membrane-bound conformations, a full-length IL-1R1 structural model was developed and subjected to classical equilibrium MD. Comparative analysis of multiple MD trajectories of the soluble and the membrane-bound IL-1R1-ECDs reveals that (i) as somewhat expected, the extent of the visited “open-to-closed” transitional states differs significantly between the soluble and membrane-bound forms; (ii) the soluble form presents open-closed transitions, sampling a wider rotational motion between the Ig-like domains of the ECD, visiting closed and “twisted” conformations in higher extent, whereas the membrane-bound form is characterized by more conformationally restricted states; (iii) interestingly, the backbone dihedral angles of residues Glu202, Glu203 and Asn204, located in the flexible linker, display the highest variations during the transition between discrete conformational states detected in IL-1R1, thus appearing to work as the “central wheel of a clock’s movement”. The simulations and analyses presented in this contribution offer a deeper insight into the structure and dynamics of IL-1R1, which may be explored in a drug discovery setting.     

Reorienting Mechanism of Harderoheme in Coproheme Decarboxylase—A Computational Study

Coproheme decarboxylase (ChdC) is an important enzyme in the coproporphyrin-dependent pathway (CPD) of Gram-positive bacteria that decarboxylates coproheme on two propionates at position 2 and position 4 sequentially to generate heme b by using H₂O₂ as an oxidant. This work focused on the ChdC from Geobacillus stearothermophilus (GsChdC) to elucidate the mechanism of its sequential two-step decarboxylation of coproheme. The models of GsChdC in a complex with substrate and reaction intermediate were built to investigate the reorienting mechanism of harderoheme. Targeted molecular dynamics simulations on these models validated that harderoheme is able to rotate in the active site of GsChdC with a 19.06-kcal·mol⁻¹ energy barrier after the first step of decarboxylation to bring the propionate at position 4 in proximity of Tyr145 to continue the second decarboxylation step. The harderoheme rotation mechanism is confirmed to be much easier than the release–rebinding mechanism. In the active site of GsChdC, Trp157 and Trp198 comprise a “gate” construction to regulate the clockwise rotation of the harderoheme. Lys149 plays a critical role in the rotation mechanism, which not only keeps the Trp157–Trp198 “gate” from being closed but also guides the propionate at position 4 through the gap between Trp157 and Trp198 through a salt bridge interaction.     

Structural Evolution of Delta (B.1.617.2) and Omicron (BA.1) Spike Glycoproteins

The vast amount of epidemiologic and genomic data that were gathered as a global response to the COVID-19 pandemic that was caused by SARS-CoV-2 offer a unique opportunity to shed light on the structural evolution of coronaviruses and in particular on the spike (S) glycoprotein, which mediates virus entry into the host cell by binding to the human ACE2 receptor. Herein, we carry out an investigation into the dynamic properties of the S glycoprotein, focusing on the much more transmissible Delta and Omicron variants. Notwithstanding the great number of mutations that have accumulated, particularly in the Omicron S glycoprotein, our data clearly showed the conservation of some structural and dynamic elements, such as the global motion of the receptor binding domain (RBD). However, our studies also revealed structural and dynamic alterations that were concentrated in the aa 627–635 region, on a small region of the receptor binding motif (aa 483–485), and the so-called “fusion-peptide proximal region”. In particular, these last two S regions are known to be involved in the human receptor ACE2 recognition and membrane fusion. Our structural evidence, therefore, is likely involved in the observed different transmissibility of these S mutants. Finally, we highlighted the role of glycans in the increased RBD flexibility of the monomer in the up conformation of Omicron.     

Design and Engineering of “Green” Nanoemulsions for Enhanced Topical Delivery of Bakuchiol Achieved in a Sustainable Manner: A Novel Eco-Friendly Approach to Bioretinol

In the present work, we establish novel “environmentally-friendly” oil-in-water nanoemulsions to enhance the transdermal delivery of bakuchiol, the so-called “bioretinol” obtained from powdered Psoralea corylifolia seeds via a sustainable process, i.e., using a supercritical fluid extraction approach with pure carbon dioxide (SC-CO₂). According to Green Chemistry principles, five novel formulations were stabilized by “green” hybrid ionic surfactants such as coco-betaine—surfactin molecules obtained from coconut and fermented rapeseed meal. Preliminary optimization studies involving three dispersion stability tests, i.e., centrifugation, heating, and cooling cycles, indicated the most promising candidates for further physicochemical analysis. Finally, nanoemulsion colloidal characterization provided by scattering (dynamic and electrophoretic light scattering as well as backscattering), microscopic (transmission electron and confocal laser scanning microscopy), and spectroscopic (UV–Vis spectroscopy) methods revealed the most stable nanocarrier for transdermal biological investigation. In vitro, topical experiments provided on human skin cell line HaCaT keratinocytes and normal dermal NHDF fibroblasts indicated high cell viability upon treatment of the tested formulation with a final 0.02–0.2 mg/mL bakuchiol concentration. This excellent biocompatibility was confirmed by ex vivo and in vivo tests on animal and human skin tissue. The improved permeability and antiaging potential of the bakuchiol-encapsulated rich extract were observed, indicating that the obtained ecological nanoemulsions are competitive with commercial retinol formulations.     

Various Approaches to Studying the Phase Transition in an Octamethylcyclotetrasiloxane Crystal: From X-ray Structural Analysis to Metadynamics

The structure, thermodynamic parameters, and the character of thermal motion in octamethylcyclotetrasiloxane (D4) were investigated using the combination of experimental (single-crystal X-ray diffraction, thermochemistry) and theoretical (density functional theory calculations, ab initio molecular dynamics and metadynamics) methods. Single crystals of D4 were grown in a glass capillary in situ and the structures of high- (238–270 K) and low-temperature (100–230 K) phases were studied in detail. In the temperature interval 230–238 K, a phase transition with rather low enthalpy (−1.04(7) kcal/mol) was detected. It was found that phase transition is accompanied by change of conformation of cyclosiloxane moiety from boat-saddle (cradle) to chair. According to PBE0/6-311G(d,p) calculation of isolated D4, such conformation changes are characterized by a low barrier (0.07 kcal/mol). The character of molecular thermal motion and the path of phase transition were established with combination of periodic DFT calculations, including molecular dynamics and metadynamics. The effect of crystal field led to an increase in the calculated phase transition barrier (4.27 kcal/mol from low- to high-temperature phase and 3.20 kcal/mol in opposite direction).     

Secondary Metabolic Profile as a Tool for Distinction and Characterization of Cultivars of Black Pepper (Piper nigrum L.) Cultivated in Pará State,Brazil

This study evaluated the chemical compositions of the leaves and fruits of eight black pepper cultivars cultivated in Pará State (Amazon, Brazil). Hydrodistillation and gas chromatography–mass spectrometry were employed to extract and analyze the volatile compounds, respectively. Sesquiterpene hydrocarbons were predominant (58.5–90.9%) in the cultivars “Cingapura”, “Equador”, “Guajarina”, “Iaçará”, and “Kottanadan”, and “Bragantina”, “Clonada”, and “Uthirankota” displayed oxygenated sesquiterpenoids (50.6–75.0%). The multivariate statistical analysis applied using volatile composition grouped the samples into four groups: γ-Elemene, curzerene, and δ-elemene (“Equador”/“Guajarina”, I); δ-elemene (“Iaçará”/“Kottanadan”/“Cingapura”, II); elemol (“Clonada”/“Uthirankota”, III) and α-muurolol, bicyclogermacrene, and cubebol (“Bragantina”, IV). The major compounds in all fruit samples were monoterpene hydrocarbons such as α-pinene, β-pinene, and limonene. Among the cultivar leaves, phenolics content (44.75–140.53 mg GAE·g⁻¹ FW), the enzymatic activity of phenylalanine-ammonia lyase (20.19–57.22 µU·mL⁻¹), and carotenoids (0.21–2.31 µg·mL⁻¹) displayed significant variations. Due to black pepper’s susceptibility to Fusarium infection, a molecular docking analysis was carried out on Fusarium protein targets using each cultivar’s volatile components. F. oxysporum endoglucanase was identified as the preferential protein target of the compounds. These results can be used to identify chemical markers related to the susceptibility degree of black pepper cultivars to plant diseases prevalent in Pará State.     

Altered Local Interactions and Long-Range Communications in UK Variant (B.1.1.7) Spike Glycoprotein

The COVID-19 pandemic is caused by SARS-CoV-2. Currently, most of the research efforts towards the development of vaccines and antibodies against SARS-CoV-2 were mainly focused on the spike (S) protein, which mediates virus entry into the host cell by binding to ACE2. As the virus SARS-CoV-2 continues to spread globally, variants have emerged, characterized by multiple mutations of the S glycoprotein. Herein, we employed microsecond-long molecular dynamics simulations to study the impact of the mutations of the S glycoprotein in SARS-CoV-2 Variant of Concern 202012/01 (B.1.1.7), termed the “UK variant”, in comparison with the wild type, with the aim to decipher the structural basis of the reported increased infectivity and virulence. The simulations provided insights on the different dynamics of UK and wild-type S glycoprotein, regarding in particular the Receptor Binding Domain (RBD). In addition, we investigated the role of glycans in modulating the conformational transitions of the RBD. The overall results showed that the UK mutant experiences higher flexibility in the RBD with respect to wild type; this behavior might be correlated with the increased transmission reported for this variant. Our work also adds useful structural information on antigenic “hotspots” and epitopes targeted by neutralizing antibodies.     

Dynamics of Endothelial Engagement and Filopodia Formation in Complex 3D Microscaffolds

The understanding of endothelium–extracellular matrix interactions during the initiation of new blood vessels is of great medical importance; however, the mechanobiological principles governing endothelial protrusive behaviours in 3D microtopographies remain imperfectly understood. In blood capillaries submitted to angiogenic factors (such as vascular endothelial growth factor, VEGF), endothelial cells can transiently transdifferentiate in filopodia-rich cells, named tip cells, from which angiogenesis processes are locally initiated. This protrusive state based on filopodia dynamics contrasts with the lamellipodia-based endothelial cell migration on 2D substrates. Using two-photon polymerization, we generated 3D microstructures triggering endothelial phenotypes evocative of tip cell behaviour. Hexagonal lattices on pillars (“open”), but not “closed” hexagonal lattices, induced engagement from the endothelial monolayer with the generation of numerous filopodia. The development of image analysis tools for filopodia tracking allowed to probe the influence of the microtopography (pore size, regular vs. elongated structures, role of the pillars) on orientations, engagement and filopodia dynamics, and to identify MLCK (myosin light-chain kinase) as a key player for filopodia-based protrusive mode. Importantly, these events occurred independently of VEGF treatment, suggesting that the observed phenotype was induced through microtopography. These microstructures are proposed as a model research tool for understanding endothelial cell behaviour in 3D fibrillary networks.     

Systematic Approach to Find the Global Minimum of Relaxation Dispersion Data for Protein-Induced B–Z Transition of DNA

Carr–Purcell–Meiboom–Gill (CPMG) relaxation dispersion spectroscopy is commonly used for quantifying conformational changes of protein in μs-to-ms timescale transitions. To elucidate the dynamics and mechanism of protein binding, parameters implementing CPMG relaxation dispersion results must be appropriately determined. Building an analytical model for multi-state transitions is particularly complex. In this study, we developed a new global search algorithm that incorporates a random search approach combined with a field-dependent global parameterization method. The robust inter-dependence of the parameters carrying out the global search for individual residues (GSIR) or the global search for total residues (GSTR) provides information on the global minimum of the conformational transition process of the Zα domain of human ADAR1 (hZα_ADAR1)–DNA complex. The global search results indicated that a α-helical segment of hZα_ADAR1 provided the main contribution to the three-state conformational changes of a hZα_ADAR1—DNA complex with a slow B–Z exchange process. The two global exchange rate constants, k_ex and k_ZB, were found to be 844 and 9.8 s⁻¹, respectively, in agreement with two regimes of residue-dependent chemical shift differences—the “dominant oscillatory regime” and “semi-oscillatory regime”. We anticipate that our global search approach will lead to the development of quantification methods for conformational changes not only in Z-DNA binding protein (ZBP) binding interactions but also in various protein binding processes.     

Shape Transformations and Self-Assembly of Hairy Particles under Confinement

Molecular dynamics simulations are used to investigate the behavior of polymer-tethered nanoparticles between two inert or attractive walls. The confinement in pores creates new possibilities for controlling the shape transformation of individual hairy particles and their self-organization. We introduce a minimalistic model of the system; only chain-wall interactions are assumed to be attractive, while the others are softly repulsive. We show how the shape of isolated particles can be controlled by changing the wall separation and the strength of the interaction with the surfaces. For attractive walls, we found two types of structures, “bridges” and “mounds”. The first structures are similar to flanged spools in which the chains are connected with both walls and form bridges between them. We observed various bridges, symmetrical and asymmetrical spools, hourglasses, and pillars. The bridge-like structures can be “nano-oscillators” in which the cores jump from one wall to the other. We also study the self-assembly of a dense fluid of hairy particles in slit-like pores and analyze how the system morphology depends on interactions with the surfaces and the wall separation. The hairy particles form layers parallel to the walls. Different ordered structures, resembling two-dimensional crystalline lattices, are reported. We demonstrate that hairy particles are a versatile soft component forming a variety of structures in the slits.     

Analytical model of the network topology and rigidity of calcium aluminosilicate glasses

Topological constraint theory (TCT) has enabled the prediction of various properties of oxide glasses as a function of their composition and structure. However, the robust application of TCT relies on accurate knowledge of the network structure and topology. Here, based on classical molecular dynamics simulations, we derive a fully analytical model describing the topology of the calcium aluminosilicate [(CaO)_x(Al₂O₃)_y(SiO₂)_1−x−y, CAS] ternary system. This model yields the state of rigidity (flexible, isostatic, or stressed-rigid) of CAS systems as a function of composition and temperature. These results reveal the existence of correlations between network topology and glass-forming ability. This study suggests that glass-forming ability is encoded in the network topology of the liquid state rather than that of the glassy state.     

Density Scaling of Translational and Rotational Molecular Dynamics in a Simple Ellipsoidal Model near the Glass Transition

In this paper, we show that a simple anisotropic model of supercooled liquid properly reflects some density scaling properties observed for experimental data, contrary to many previous results obtained from isotropic models. We employ a well-known Gay–Berne model earlier parametrized to achieve a supercooling and glass transition at zero pressure to find the point of glass transition and explore volumetric and dynamic properties in the supercooled liquid state at elevated pressure. We focus on dynamic scaling properties of the anisotropic model of supercooled liquid to gain a better insight into the grounds for the density scaling idea that bears hallmarks of universality, as follows from plenty of experimental data collected near the glass transition for different dynamic quantities. As a result, the most appropriate values of the scaling exponent γ are established as invariants for a given anisotropy aspect ratio to successfully scale both the translational and rotational relaxation times considered as single variable functions of density^γ/temperature. These scaling exponent values are determined based on the density scaling criterion and differ from those obtained in other ways, such as the virial–potential energy correlation and the equation of state derived from the effective short-range intermolecular potential, which is qualitatively in accordance with the results yielded from experimental data analyses. Our findings strongly suggest that there is a deep need to employ anisotropic models in the study of glass transition and supercooled liquids instead of the isotropic ones very commonly exploited in molecular dynamics simulations of supercooled liquids over the last decades.     

Atomistic Simulations of Functionalized Nano-Materials for Biosensors Applications

Nanoscale biosensors, a highly promising technique in clinical analysis, can provide sensitive yet label-free detection of biomolecules. The spatial and chemical specificity of the surface coverage, the proper immobilization of the bioreceptor as well as the underlying interfacial phenomena are crucial elements for optimizing the performance of a biosensor. Due to experimental limitations at the microscopic level, integrated cross-disciplinary approaches that combine in silico design with experimental measurements have the potential to present a powerful new paradigm that tackles the issue of developing novel biosensors. In some cases, computational studies can be seen as alternative approaches to assess the microscopic working mechanisms of biosensors. Nonetheless, the complex architecture of a biosensor, associated with the collective contribution from “substrate–receptor–analyte” conjugate in a solvent, often requires extensive atomistic simulations and systems of prohibitive size which need to be addressed. In silico studies of functionalized surfaces also require ad hoc force field parameterization, as existing force fields for biomolecules are usually unable to correctly describe the biomolecule/surface interface. Thus, the computational studies in this field are limited to date. In this review, we aim to introduce fundamental principles that govern the absorption of biomolecules onto functionalized nanomaterials and to report state-of-the-art computational strategies to rationally design nanoscale biosensors. A detailed account of available in silico strategies used to drive and/or optimize the synthesis of functionalized nanomaterials for biosensing will be presented. The insights will not only stimulate the field to rationally design functionalized nanomaterials with improved biosensing performance but also foster research on the required functionalization to improve biomolecule–surface complex formation as a whole.     

A Deep Dive into VDAC1 Conformational Diversity Using All-Atom Simulations Provides New Insights into the Structural Origin of the Closed States

The voltage-dependent anion channel 1 (VDAC1) is a crucial mitochondrial transporter that controls the flow of ions and respiratory metabolites entering or exiting mitochondria. As a voltage-gated channel, VDAC1 can switch between a high-conducting “open” state and a low-conducting “closed” state emerging at high transmembrane (TM) potentials. Although cell homeostasis depends on channel gating to regulate the transport of ions and metabolites, structural hallmarks characterizing the closed states remain unknown. Here, we performed microsecond accelerated molecular dynamics to highlight a vast region of VDAC1 conformational landscape accessible at typical voltages known to promote closure. Conformers exhibiting durable subconducting properties inherent to closed states were identified. In all cases, the low conductance was due to the particular positioning of an unfolded part of the N-terminus, which obstructed the channel pore. While the N-terminal tail was found to be sensitive to voltage orientation, our models suggest that stable low-conducting states of VDAC1 predominantly take place from disordered events and do not result from the displacement of a voltage sensor or a significant change in the pore. In addition, our results were consistent with conductance jumps observed experimentally and corroborated a recent study describing entropy as a key factor for VDAC gating.     

Contrasting Local and Macroscopic Effects of Collagen Hydroxylation

Collagen is heavily hydroxylated. Experiments show that proline hydroxylation is important to triple helix (monomer) stability, fibril assembly, and interaction of fibrils with other molecules. Nevertheless, experiments also show that even without hydroxylation, type I collagen does assemble into its native D-banded fibrillar structure. This raises two questions. Firstly, even though hydroxylation removal marginally affects macroscopic structure, how does such an extensive chemical change, which is expected to substantially reduce hydrogen bonding capacity, affect local structure? Secondly, how does such a chemical perturbation, which is expected to substantially decrease electrostatic attraction between monomers, affect collagen’s mechanical properties? To address these issues, we conduct a benchmarked molecular dynamics study of rat type I fibrils in the presence and absence of hydroxylation. Our simulations reproduce the experimental observation that hydroxylation removal has a minimal effect on collagen’s D-band length. We also find that the gap-overlap ratio, monomer width and monomer length are minimally affected. Surprisingly, we find that de-hydroxylation also has a minor effect on the fibril’s Young’s modulus, and elastic stress build up is also accompanied by tightening of triple-helix windings. In terms of local structure, de-hydroxylation does result in a substantial drop (23%) in inter-monomer hydrogen bonding. However, at the same time, the local structures and inter-monomer hydrogen bonding networks of non-hydroxylated amino acids are also affected. It seems that it is this intrinsic plasticity in inter-monomer interactions that preclude fibrils from undergoing any large changes in macroscopic properties. Nevertheless, changes in local structure can be expected to directly impact collagen’s interaction with extra-cellular matrix proteins. In general, this study highlights a key challenge in tissue engineering and medicine related to mapping collagen chemistry to macroscopic properties but suggests a path forward to address it using molecular dynamics simulations.     

Molecular dynamics study on the viscosity of glass-forming systems near and below the glass transition temperature

Temperature-dependent viscosity is critical to decipher two profound questions in condensed matter physics, namely the glass transition and the relaxation of amorphous solids. However, direct measurement of viscosity over a large temperature range is extremely difficult. Here, using classical molecular dynamics (MD) simulations, we report a novel method to calculate the equilibrium viscosity of supercooled liquid both above and below the glass transition temperature (T_g) and to estimate the nonequilibrium viscosity of glass down to room temperature. Based on the shoving model, we derived an analytical formula showing that the shear viscosity in logarithmic scale changes linearly with the shear-induced variation in shear modulus or potential energy of the glass-forming system. The shear viscosity as a function of steady-state potential energy of liquid under different shear strain rates can be directly calculated in MD simulations; together with its equilibrium potential energy, one can extrapolate the zero-strain-rate equilibrium viscosity. We verified the proposed model by reliably calculating equilibrium viscosity near T_g of four glass-forming systems (Kob–Andersen system, silica, Cu_45.5Zr_45.5Al₉, and silicon) with different fragilities. Furthermore, our model can estimate the nonequilibrium viscosity of glass below T_g; the upper-bound nonequilibrium viscosity of amorphous silica and silicon at room temperature are calculated to be ~10³² and 10²⁵ Pa·s, respectively.