The thermodynamics and kinetics of protein folding and protein aggregation in vivo are of great importance in numerous scientific areas including fundamental biophysics research, nanotechnology, and ...medicine. However, these processes remain poorly understood in both in vivo and in vitro systems. Here we extend an established model for protein aggregation that is based on the kinetic equations for the moments of the polymer size distribution by introducing macromolecular crowding particles into the model using scaled-particle and transition-state theories. The model predicts that the presence of crowders can either speed up, cause no change to, or slow down the progress of the aggregation compared to crowder-free solutions, in striking agreement with experimental results from nine different amyloid-forming proteins that utilized dextran as the crowder. These different dynamic effects of macromolecular crowding can be understood in terms of the change of excluded volume associated with each reaction step.
Abstract
Nucleic acid interactions under crowded environments are of great importance for biological processes and nanotechnology. However, the kinetics and thermodynamics of nucleic acid ...interactions in a crowded environment remain poorly understood. We use a coarse-grained model of DNA to study the kinetics and thermodynamics of DNA duplex and hairpin formation in crowded environments. We find that crowders can increase the melting temperature of both an 8-mer DNA duplex and a hairpin with a stem of 6-nt depending on the excluded volume fraction of crowders in solution and the crowder size. The crowding induced stability originates from the entropic effect caused by the crowding particles in the system. Additionally, we study the hybridization kinetics of DNA duplex formation and the formation of hairpin stems, finding that the reaction rate kon is increased by the crowding effect, while koff is changed only moderately. The increase in kon mostly comes from increasing the probability of reaching a transition state with one base pair formed. A DNA strand displacement reaction in a crowded environment is also studied with the model and we find that rate of toehold association is increased, with possible applications to speeding up strand displacement cascades in nucleic acid nanotechnology.
The problem of retrosynthetic planning can be framed as a one-player game, in which the chemist (or a computer program) works backward from a molecular target to simpler starting materials through a ...series of choices regarding which reactions to perform. This game is challenging as the combinatorial space of possible choices is astronomical, and the value of each choice remains uncertain until the synthesis plan is completed and its cost evaluated. Here, we address this search problem using deep reinforcement learning to identify policies that make (near) optimal reaction choices during each step of retrosynthetic planning according to a user-defined cost metric. Using a simulated experience, we train a neural network to estimate the expected synthesis cost or value of any given molecule based on a representation of its molecular structure. We show that learned policies based on this value network can outperform a heuristic approach that favors symmetric disconnections when synthesizing unfamiliar molecules from available starting materials using the fewest number of reactions. We discuss how the learned policies described here can be incorporated into existing synthesis planning tools and how they can be adapted to changes in the synthesis cost objective or material availability.
Selection protocols such as SELEX, where molecules are selected over multiple rounds for their ability to bind to a target of interest, are popular methods for obtaining binders for diagnostic and ...therapeutic purposes. We show that Restricted Boltzmann Machines (RBMs), an unsupervised two-layer neural network architecture, can successfully be trained on sequence ensembles from single rounds of SELEX experiments for thrombin aptamers. RBMs assign scores to sequences that can be directly related to their fitnesses estimated through experimental enrichment ratios. Hence, RBMs trained from sequence data at a given round can be used to predict the effects of selection at later rounds. Moreover, the parameters of the trained RBMs are interpretable and identify functional features contributing most to sequence fitness. To exploit the generative capabilities of RBMs, we introduce two different training protocols: one taking into account sequence counts, capable of identifying the few best binders, and another based on unique sequences only, generating more diverse binders. We then use RBMs model to generate novel aptamers with putative disruptive mutations or good binding properties, and validate the generated sequences with gel shift assay experiments. Finally, we compare the RBM’s performance with different supervised learning approaches that include random forests and several deep neural network architectures.
Celotno besedilo
Dostopno za:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
Holographic cloud probes provide unprecedented information on cloud particle density, size and position. Each laser shot captures particles within a large volume, where images can be computationally ...refocused to determine particle size and location. However, processing these holograms with standard methods or machine learning (ML) models requires considerable computational resources, time and occasional human intervention. ML models are trained on simulated holograms obtained from the physical model of the probe since real holograms have no absolute truth labels. Using another processing method to produce labels would be subject to errors that the ML model would subsequently inherit. Models perform well on real holograms only when image corruption is performed on the simulated images during training, thereby mimicking non-ideal conditions in the actual probe. Optimizing image corruption requires a cumbersome manual labeling effort. Here we demonstrate the application of the neural style translation approach to the simulated holograms. With a pre-trained convolutional neural network, the simulated holograms are "stylized" to resemble the real ones obtained from the probe, while at the same time preserving the simulated image "content" (e.g. the particle locations and sizes). With an ML model trained to predict particle locations and shapes on the stylized data sets, we observed comparable performance on both simulated and real holograms, obviating the need to perform manual labeling. The described approach is not specific to holograms and could be applied in other domains for capturing noise and imperfections in observational instruments to make simulated data more like real world observations.
We propose a kinetic model for the self-aggregation by amyloid proteins. By extending several well-known models for protein aggregation, the time evolution of aggregate concentrations containing r ...proteins, denoted c r (t), can be written in terms of generalized Smoluchowski kinetics. With this approach, we take into account all possible aggregation and fragmentation reactions involving clusters of any size. Correspondingly, an aggregate of size x + y could be formed by or break up into two smaller constituent aggregates of sizes x and y. The rates of each aggregation or fragmentation reaction, called kernels, are specified in terms of the aggregate size, and we solve c r (t) for large cluster sizes using numerical techniques. We show that by using Smoluchowski kinetics many pathways to fibrillation are possible and quantities, such as the aggregate length distribution at an arbitrary time, can be calculated. We show that the predicted results of the model are in agreement with the experimental observations.
Abstract
We use the oxDNA coarse-grained model to provide a detailed characterization of the fundamental structural properties of DNA origami, focussing on archetypal 2D and 3D origami. The model ...reproduces well the characteristic pattern of helix bending in a 2D origami, showing that it stems from the intrinsic tendency of anti-parallel four-way junctions to splay apart, a tendency that is enhanced both by less screened electrostatic interactions and by increased thermal motion. We also compare to the structure of a 3D origami whose structure has been determined by cryo-electron microscopy. The oxDNA average structure has a root-mean-square deviation from the experimental structure of 8.4 Å, which is of the order of the experimental resolution. These results illustrate that the oxDNA model is capable of providing detailed and accurate insights into the structure of DNA origami, and has the potential to be used to routinely pre-screen putative origami designs and to investigate the molecular mechanisms that regulate the properties of DNA origami.
Investigation of protein self-assembly processes is important for understanding the growth processes of functional proteins as well as disease-causing amyloids. Inside cells, intrinsic molecular ...fluctuations are so high that they cast doubt on the validity of the deterministic rate-equation approach. Furthermore, the protein environments inside cells are often crowded with other macromolecules, with volume fractions of the crowders as high as 40%. We have developed a stochastic kinetic framework using Gillespie’s algorithm for general systems undergoing particle self-assembly, including particularly protein aggregation at the cellular level. The effects of macromolecular crowding are investigated using models built on scaled-particle and transition-state theories. The stochastic kinetic method can be formulated to provide information on the dominating aggregation mechanisms in a method called reaction frequency (or propensity) analysis. This method reveals that the change of scaling laws related to the lag time can be directly related to the change in the frequencies of reaction mechanisms. Further examination of the time evolution of the fibril mass and length quantities unveils that maximal fluctuations occur in the periods of rapid fibril growth and the fluctuations of both quantities can be sensitive functions of rate constants. The presence of crowders often amplifies the roles of primary and secondary nucleation and causes shifting in the relative importance of elongation, shrinking, fragmentation, and coagulation of linear aggregates. We also show a dual effect of changing volume on the halftime of aggregation for ApoC2 which is reduced in the presence of crowders. A comparison of the results of stochastic simulations with those of rate equations gives us information on the convergence relation between them and how the roles of reaction mechanisms change as the system volume is varied.
We introduce an extended version of oxDNA, a coarse-grained model of deoxyribonucleic acid (DNA) designed to capture the thermodynamic, structural, and mechanical properties of single- and ...double-stranded DNA. By including explicit major and minor grooves and by slightly modifying the coaxial stacking and backbone-backbone interactions, we improve the ability of the model to treat large (kilobase-pair) structures, such as DNA origami, which are sensitive to these geometric features. Further, we extend the model, which was previously parameterised to just one salt concentration (Na(+) = 0.5M), so that it can be used for a range of salt concentrations including those corresponding to physiological conditions. Finally, we use new experimental data to parameterise the oxDNA potential so that consecutive adenine bases stack with a different strength to consecutive thymine bases, a feature which allows a more accurate treatment of systems where the flexibility of single-stranded regions is important. We illustrate the new possibilities opened up by the updated model, oxDNA2, by presenting results from simulations of the structure of large DNA objects and by using the model to investigate some salt-dependent properties of DNA.
We use Monte Carlo simulations and free-energy techniques to show that binary solutions of penta- and hexavalent two-dimensional patchy particles can form thermodynamically stable quasicrystals even ...at very narrow patch widths, provided their patch interactions are chosen in an appropriate way. Such patchy particles can be thought of as a coarse-grained representation of DNA multi-arm 'star' motifs, which can be chosen to bond with one another very specifically by tuning the DNA sequences of the protruding arms. We explore several possible design strategies and conclude that DNA star tiles that are designed to interact with one another in a specific but not overly constrained way could potentially be used to construct soft quasicrystals in experiment. We verify that such star tiles can form stable dodecagonal motifs using oxDNA, a realistic coarse-grained model of DNA.