Multiscale modeling of molecular structure and optical properties of complex supramolecular aggregates.
Ontology highlight
ABSTRACT:
SUBMITTER: Bondarenko AS
PROVIDER: S-EPMC8162738 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Dye aggregates are of interest for excitonic applications, including biomedical imaging, organic photovoltaics, and quantum information systems. Dyes with large transition dipole moments (μ) are necessary to optimize coupling within dye aggregates. Extinction coefficients (ε) can be used to determine the μ of dyes, and so dyes with a large ε (>150,000 M-1cm-1) should be engineered or identified. However, dye properties leading to a large ε are not fully understood, and low-throughput methods of dye screening, such as experimental measurements or density functional theory (DFT) calculations, can be time-consuming. In order to screen large datasets of molecules for desirable properties (i.e., large ε and μ), a computational workflow was established using machine learning (ML), DFT, time-dependent (TD-) DFT, and molecular dynamics (MD). ML models were developed through training and validation on a dataset of 8802 dyes using structural features. A Classifier was developed with an accuracy of 97% and a Regressor was constructed with an R2 of above 0.9, comparing between experiment and ML prediction. Using the Regressor, the ε values of over 18,000 dyes were predicted. The top 100 dyes were further screened using DFT and TD-DFT to identify 15 dyes with a μ relative to a reference dye, pentamethine indocyanine dye Cy5. Two benchmark MD simulations were performed on Cy5 and Cy5.5 dimers, and it was found that MD could accurately capture experimental results. The results of this study exhibit that our computational workflow for identifying dyes with a large μ for excitonic applications is effective and can be used as a tool to develop new dyes for excitonic applications.
| S-EPMC9182218 | biostudies-literature
Project description:The ability of phenylalanine to form fibrillar nanostructures was demonstrated on multiple occasions, and such an oligomerization reaction could be the cause of cytotoxicity in patients with phenylketonuria. These findings were supported by claims that L-phenylalanine (Phe) fibrils have amyloid properties and can be detected using thioflavin T fluorescence assay. However, a part of Phe aggregation studies reported the opposite data, suggesting no amyloid structures to be formed. Due to the contradicting reports, the amyloid nature of Phe aggregates remains uncertain. In this work we tested Phe aggregation under conditions where amyloid formation was previously reported. We show the emergence of Phe aggregates with visible light optical properties that overlap with the spectra of dyes used in amyloid fibril assays, which could lead to false-positive identifications.
| S-EPMC6394350 | biostudies-literature
Project description:Lignin is an aromatic biopolymer found in ubiquitous sources of woody biomass. Designing and optimizing lignin valorization processes requires a fundamental understanding of lignin structures. Experimental characterization techniques, such as 2D-heteronuclear single quantum coherence (HSQC) nuclear magnetic resonance (NMR) spectra, could elucidate the global properties of the polymer molecules. Computer models could extend the resolution of experiments by representing structures at the molecular and atomistic scales. We introduce a graph-based multiscale modeling framework for lignin structure generation and visualization. The framework employs accelerated rejection-free polymerization and hierarchical Metropolis Monte Carlo optimization algorithms. We obtain structure libraries for various lignin feedstocks based on literature and new experimental NMR data for poplar wood, pinewood, and herbaceous lignin. The framework could guide researchers towards feasible lignin structures, efficient space exploration, and future kinetics modeling. Its software implementation in Python, LigninGraphs, is open-source and available on GitHub.
| S-EPMC9261032 | biostudies-literature
Project description:Neisseria gonorrhoeae is the causative agent of one of the most common sexually transmitted diseases, gonorrhea. Over the past two decades there has been an alarming increase of reported gonorrhea cases where the bacteria were resistant to the most commonly used antibiotics thus prompting for alternative antimicrobial treatment strategies. The crucial step in this and many other bacterial infections is the formation of microcolonies, agglomerates consisting of up to several thousands of cells. The attachment and motility of cells on solid substrates as well as the cell-cell interactions are primarily mediated by type IV pili, long polymeric filaments protruding from the surface of cells. While the crucial role of pili in the assembly of microcolonies has been well recognized, the exact mechanisms of how they govern the formation and dynamics of microcolonies are still poorly understood. Here, we present a computational model of individual cells with explicit pili dynamics, force generation and pili-pili interactions. We employ the model to study a wide range of biological processes, such as the motility of individual cells on a surface, the heterogeneous cell motility within the large cell aggregates, and the merging dynamics and the self-assembly of microcolonies. The results of numerical simulations highlight the central role of pili generated forces in the formation of bacterial colonies and are in agreement with the available experimental observations. The model can quantify the behavior of multicellular bacterial colonies on biologically relevant temporal and spatial scales and can be easily adjusted to include the geometry and pili characteristics of various bacterial species. Ultimately, the combination of the microbiological experimental approach with the in silico model of bacterial colonies might provide new qualitative and quantitative insights on the development of bacterial infections and thus pave the way to new antimicrobial treatments.
| S-EPMC8132470 | biostudies-literature
Project description:Chemical synapses exhibit a diverse array of internal mechanisms that affect the dynamics of transmission efficacy. Many of these processes, such as release of neurotransmitter and vesicle recycling, depend strongly on activity-dependent influx and accumulation of Ca2+. To model how each of these processes may affect the processing of information in neural circuits, and how their dysfunction may lead to disease states, requires a computationally efficient modelling framework, capable of generating accurate phenomenology without incurring a heavy computational cost per synapse. Constructing a phenomenologically realistic model requires the precise characterization of the timing and probability of neurotransmitter release. Difficulties arise in that functional forms of instantaneous release rate can be difficult to extract from noisy data without running many thousands of trials, and in biophysical synapses, facilitation of per-vesicle release probability is confounded by depletion. To overcome this, we obtained traces of free Ca2+ concentration in response to various action potential stimulus trains from a molecular MCell model of a hippocampal Schaffer collateral axon. Ca2+ sensors were placed at varying distance from a voltage-dependent calcium channel (VDCC) cluster, and Ca2+ was buffered by calbindin. Then, using the calcium traces to drive deterministic state vector models of synaptotagmin 1 and 7 (Syt-1/7), which respectively mediate synchronous and asynchronous release in excitatory hippocampal synapses, we obtained high-resolution profiles of instantaneous release rate, to which we applied functional fits. Synchronous vesicle release occurred predominantly within half a micron of the source of spike-evoked Ca2+ influx, while asynchronous release occurred more consistently at all distances. Both fast and slow mechanisms exhibited multi-exponential release rate curves, whose magnitudes decayed exponentially with distance from the Ca2+ source. Profile parameters facilitate on different time scales according to a single, general facilitation function. These functional descriptions lay the groundwork for efficient mesoscale modelling of vesicular release dynamics.
| S-EPMC9119629 | biostudies-literature
Project description:Supramolecular copolymers, non-covalent analogues of synthetic copolymers, constitute a new and promising class of polymers. In contrast to their covalent counterparts, the details of their mechanism of formation, as well as the factors determining their composition and length, are still poorly understood. Here, the supramolecular copolymerization between two slightly structurally different benzene-1,3,5-tricarboxamide (BTA) monomers functionalized with either oligodimethylsiloxane (oDMSi) or alkyl side chains is unraveled by combining experimental and theoretical approaches. By applying the "sergeant-and-soldiers" approach using circular dichroism (CD) experiments, we are able to obtain detailed insights into the structure and composition of these supramolecular copolymers. Moreover, we observe an unexpected chiral induction upon mixing two independently CD-silent solutions of the achiral (soldier) and chiral (sergeant) monomers. We find that the subtle differences in the chemical structure of the two monomers impact their homopolymerization mechanism: whereas alkyl-BTAs cooperatively self-assemble, oDMSi-BTAs self-assemble in an isodesmic manner. The effect of these mechanistic differences in the supramolecular copolymerization process is investigated as a function of the composition of the two monomers and explicitly rationalized by mathematical modeling. The results show that, at low fractions of oDMSi-BTA sergeants (<10 mol%), the polymerization process is cooperative and the supramolecular helicity is biased toward the helical preference of the sergeant. However, at higher fractions of oDMSi-BTA sergeant (>25 mol%), the isodesmic assembly of the increasing amounts of sergeant becomes more dominant, and different species start to coexist in the copolymerization process. The analysis of the experimental data with a newly developed theoretical model allows us to quantify the thermodynamic parameters, the distribution of different species, and the compositions and stack lengths of the formed supramolecular copolymers existing at various feed ratios of the two monomers.
| S-EPMC5445503 | biostudies-literature
Project description:Cryogenic transmission electron microscopy (cryo-TEM) studies suggest that TTBC molecules self-assemble in aqueous solution to form single-walled tubes with a diameter of about 35 Å. In order to reveal the arrangement and mutual orientations of the individual molecules in the tube, we combine information from crystal structure data of this dye with a calculation of linear absorbance and linear dichroism spectra and molecular dynamics simulations. We start with wrapping crystal planes in different directions to obtain tubes of suitable diameter. This set of tube models is evaluated by comparing the resulting optical spectra with experimental data. The tubes that can explain the spectra are investigated further by molecular dynamics simulations, including explicit solvent molecules. From the trajectories of the most stable tube models, the short-range ordering of the dye molecules is extracted and the optimization of the structure is iteratively completed. The final structural model is a tube of rings with 6-fold rotational symmetry, where neighboring rings are rotated by 30° and the transition dipole moments of the chromophores form an angle of 74° with respect to the symmetry axis of the tube. This model is in agreement with cryo-TEM images and can explain the optical spectra, consisting of a sharp red-shifted J-band that is polarized parallel to to the symmetry axis of the tube and a broad blue-shifted H-band polarized perpendicular to this axis. The general structure of the homogeneous spectrum of this hybrid HJ-aggregate is described by an analytical model that explains the difference in redistribution of oscillator strength inside the vibrational manifolds of the J- and H-bands and the relative intensities and excitation energies of those bands. In addition to the particular system investigated here, the present methodology can be expected to aid the structure prediction for a wide range of self-assembled dye aggregates.
| S-EPMC5021387 | biostudies-literature
Project description:Adaptive hydrogels tailor-made from silk-elastin-like proteins (SELPs) possess excellent biocompatibility and biodegradability with properties that are tunable and responsive to multiple simultaneous external stimuli. To unravel the molecular mechanisms of their physical response to external stimuli in tandem with experiments, here we predict and measure the variation in structural properties as a function of temperature through coarse-grained (CG) modeling of individual and crosslinked SE8Y and S4E8Y molecules, which have ratios of 1:8 and 4:8 of silk to elastin blocks respectively. Extensive structural reshuffling in single SE8Y molecules led to the increased compactness of the structure, whereas S4E8Y molecules did not experience any significant changes as they already adopted very compact structures at low temperatures. Crosslinking of SE8Y molecules at high concentrations impeded their structural transition at high temperatures that drastically reduced the degree of deswelling through extensive suppression of the structural shuffling and the trapping of the molecules in high potential energy states due to inter-molecular constraints. This integrative experimental and computational understanding of the thermal response in single molecules of SELPs and their crosslinked networks should lead to further improvements in the properties of SELP hydrogels through predictive designs and their wider applications in biomaterials and tissue engineering.
| S-EPMC6241539 | biostudies-literature
Project description:The need for more accessible energy resources makes shale formations increasingly important. Characterization of such low-permeability formations is complicated, due to the presence of multiscale features, and defies conventional methods. High-quality 3D imaging may be an ultimate solution for revealing the complexities of such porous media, but acquiring them is costly and time consuming. High-quality 2D images, on the other hand, are widely available. A novel three-step, multiscale, multiresolution reconstruction method is presented that directly uses 2D images in order to develop 3D models of shales. It uses a high-resolution 2D image representing the small-scale features to reproduce the nanopores and their network, a large scale, low-resolution 2D image to create the larger-scale characteristics, and generates stochastic realizations of the porous formation. The method is used to develop a model for a shale system for which the full 3D image is available and its properties can be computed. The predictions of the reconstructed models are in excellent agreement with the data. The method is, however, quite general and can be used for reconstructing models of other important heterogeneous materials and media. Two biological examples and from materials science are also reconstructed to demonstrate the generality of the method.
| S-EPMC4642334 | biostudies-literature
Project description:Cardiovascular disease is often related to defects of subcellular components in cardiac myocytes, specifically in the dyadic cleft, which include changes in cleft geometry and channel placement. Modeling of these pathological changes requires both spatially resolved cleft as well as whole cell level descriptions. We use a multiscale model to create dyadic structure-function relationships to explore the impact of molecular changes on whole cell electrophysiology and calcium cycling. This multiscale model incorporates stochastic simulation of individual L-type calcium channels and ryanodine receptor channels, spatially detailed concentration dynamics in dyadic clefts, rabbit membrane potential dynamics, and a system of partial differential equations for myoplasmic and lumenal free Ca2+ and Ca2+-binding molecules in the bulk of the cell. We found action potential duration, systolic, and diastolic [Ca2+] to respond most sensitively to changes in L-type calcium channel current. The ryanodine receptor channel cluster structure inside dyadic clefts was found to affect all biomarkers investigated. The shape of clusters observed in experiments by Jayasinghe et al. and channel density within the cluster (characterized by mean occupancy) showed the strongest correlation to the effects on biomarkers.
| S-EPMC6990380 | biostudies-literature