ABSTRACT
DNA-encoded libraries (DELs) provide the means to make and screen millions of diverse compounds against a target of interest in a single experiment. However, despite producing large volumes of binding data at a relatively low cost, the DEL selection process is susceptible to noise, necessitating computational follow-up to increase signal-to-noise ratios. In this work, we present a set of informatics tools to employ data from prior DEL screen(s) to gain information about which building blocks are most likely to be productive when designing new DELs for the same target. We demonstrate that similar building blocks have similar probabilities of forming compounds that bind. We then build a model from the inference that the combined behavior of individual building blocks is predictive of whether an overall compound binds. We illustrate our approach on a set of three-cycle OpenDEL libraries screened against soluble epoxide hydrolase (sEH) and report performance of more than an order of magnitude greater than random guessing on a holdout set, demonstrating that our model can serve as a baseline for comparison against other machine learning models on DEL data. Lastly, we provide a discussion on how we believe this informatics workflow could be applied to benefit researchers in their specific DEL campaigns.
Subject(s)
Drug Discovery , Small Molecule Libraries , Small Molecule Libraries/chemistry , DNA/chemistry , Machine LearningABSTRACT
Enzymatic catalysis is essential to cell survival. In many instances, enzymes that participate in reaction cascades have been shown to assemble into metabolons in response to the presence of the substrate for the first enzyme. However, what triggers metabolon formation has remained an open question. Through a combination of theory and experiments, we show that enzymes in a cascade can assemble via chemotaxis. We apply microfluidic and fluorescent spectroscopy techniques to study the coordinated movement of the first four enzymes of the glycolysis cascade: hexokinase, phosphoglucose isomerase, phosphofructokinase and aldolase. We show that each enzyme independently follows its own specific substrate gradient, which in turn is produced by the preceding enzymatic reaction. Furthermore, we find that the chemotactic assembly of enzymes occurs even under cytosolic crowding conditions.
Subject(s)
Fructose-Bisphosphate Aldolase/metabolism , Glucose-6-Phosphate Isomerase/metabolism , Hexokinase/metabolism , Phosphofructokinases/metabolism , Biocatalysis , Chemotaxis , Fructose-Bisphosphate Aldolase/chemistry , Glucose-6-Phosphate Isomerase/chemistry , Glycolysis , Hexokinase/chemistry , Molecular Structure , Phosphofructokinases/chemistry , Substrate SpecificityABSTRACT
Motor proteins such as myosin and kinesin play a major role in cellular cargo transport, muscle contraction, cell division, and engineered nanodevices. Quantifying the collective behavior of coupled motors is critical to our understanding of these systems. An excellent model system is the gliding motility assay, where hundreds of surface-adhered motors propel one cytoskeletal filament such as an actin filament or a microtubule. The filament motion can be observed using fluorescence microscopy, revealing fluctuations in gliding velocity. These velocity fluctuations have been previously quantified by a motional diffusion coefficient, which Sekimoto and Tawada explained as arising from the addition and removal of motors from the linear array of motors propelling the filament as it advances, assuming that different motors are not equally efficient in their force generation. A computational model of kinesin head diffusion and binding to the microtubule allowed us to quantify the heterogeneity of motor efficiency arising from the combination of anharmonic tail stiffness and varying attachment geometries assuming random motor locations on the surface and an absence of coordination between motors. Knowledge of the heterogeneity allows the calculation of the proportionality constant between the motional diffusion coefficient and the motor density. The calculated value (0.3) is within a standard error of our measurements of the motional diffusion coefficient on surfaces with varying motor densities calibrated by landing rate experiments. This allowed us to quantify the loss in efficiency of coupled molecular motors arising from heterogeneity in the attachment geometry.