Microscopic mechanisms of cooperative communications within single nanocatalysts.

Catalysis is a method of accelerating chemical reactions that is critically important for fundamental research as well as for industrial applications. It has been recently discovered that catalytic reactions on metal nanoparticles exhibit cooperative effects. The mechanism of these observations, however, remains not well understood. In this work, we present a theoretical investigation on possible microscopic origin of cooperative communications in nanocatalysts. In our approach, the main role is played by positively charged holes on metal surfaces. A corresponding discrete-state stochastic model for the dynamics of holes is developed and explicitly solved. It is shown that the observed spatial correlation lengths are given by the average distances migrated by the holes before they disappear, while the temporal memory is determined by their lifetimes. Our theoretical approach is able to explain the universality of cooperative communications as well as the effect of external electric fields. Theoretical predictions are in agreement with experimental observations. The proposed theoretical framework quantitatively clarifies some important aspects of the microscopic mechanisms of heterogeneous catalysis.

Chemoinformatics Insights on Molecular Jackhammers and Cancer Cells.

One of the most challenging tasks in modern medicine is to find novel efficient cancer therapeutic methods with minimal side effects. The recent discovery of several classes of organic molecules known as "molecular jackhammers" is a promising development in this direction. It is known that these molecules can directly target and eliminate cancer cells with no impact on healthy tissues. However, the underlying microscopic picture remains poorly understood. We present a study that utilizes theoretical analysis together with experimental measurements to clarify the microscopic aspects of jackhammers' anticancer activities. Our physical-chemical approach combines statistical analysis with chemoinformatics methods to design and optimize molecular jackhammers. By correlating specific physical-chemical properties of these molecules with their abilities to kill cancer cells, several important structural features are identified and discussed. Although our theoretical analysis enhances understanding of the molecular interactions of jackhammers, it also highlights the need for further research to comprehensively elucidate their mechanisms and to develop a robust physical-chemical framework for the rational design of targeted anticancer drugs.

Mesoscopic protein-rich clusters host the nucleation of mutant p53 amyloid fibrils.

The protein p53 is a crucial tumor suppressor, often called "the guardian of the genome"; however, mutations transform p53 into a powerful cancer promoter. The oncogenic capacity of mutant p53 has been ascribed to enhanced propensity to fibrillize and recruit other cancer fighting proteins in the fibrils, yet the pathways of fibril nucleation and growth remain obscure. Here, we combine immunofluorescence three-dimensional confocal microscopy of human breast cancer cells with light scattering and transmission electron microscopy of solutions of the purified protein and molecular simulations to illuminate the mechanisms of phase transformations across multiple length scales, from cellular to molecular. We report that the p53 mutant R248Q (R, arginine; Q, glutamine) forms, both in cancer cells and in solutions, a condensate with unique properties, mesoscopic protein-rich clusters. The clusters dramatically diverge from other protein condensates. The cluster sizes are decoupled from the total cluster population volume and independent of the p53 concentration and the solution concentration at equilibrium with the clusters varies. We demonstrate that the clusters carry out a crucial biological function: they host and facilitate the nucleation of amyloid fibrils. We demonstrate that the p53 clusters are driven by structural destabilization of the core domain and not by interactions of its extensive unstructured region, in contradistinction to the dense liquids typical of disordered and partially disordered proteins. Two-step nucleation of mutant p53 amyloids suggests means to control fibrillization and the associated pathologies through modifying the cluster characteristics. Our findings exemplify interactions between distinct protein phases that activate complex physicochemical mechanisms operating in biological systems.

Theoretical understanding of evolutionary dynamics on inhomogeneous networks.

Evolution is the main feature of all biological systems that allows populations to change their characteristics over successive generations. A powerful approach to understand evolutionary dynamics is to investigate fixation probabilities and fixation times of novel mutations on networks that mimic biological populations. It is now well established that the structure of such networks can have dramatic effects on evolutionary dynamics. In particular, there are population structures that might amplify the fixation probabilities while simultaneously delaying the fixation events. However, the microscopic origins of such complex evolutionary dynamics remain not well understood. We present here a theoretical investigation of the microscopic mechanisms of mutation fixation processes on inhomogeneous networks. It views evolutionary dynamics as a set of stochastic transitions between discrete states specified by different numbers of mutated cells. By specifically considering star networks, we obtain a comprehensive description of evolutionary dynamics. Our approach allows us to employ physics-inspired free-energy landscape arguments to explain the observed trends in fixation times and fixation probabilities, providing a better microscopic understanding of evolutionary dynamics in complex systems.

Bacteria-Specific Feature Selection for Enhanced Antimicrobial Peptide Activity Predictions Using Machine-Learning Methods.

There are several classes of short peptide molecules, known as antimicrobial peptides (AMPs), which are produced during the immune responses of living organisms against various infections. In recent years, substantial progress has been achieved in applying machine-learning methods to predict the activities of AMPs against bacteria. In most investigated cases, however, the outcome is not bacterium-specific since the specific features of bacteria, such as chemical composition and structure of membranes, are not considered. To overcome this problem, we developed a new computational approach that allowed us to train several supervised machine-learning models using a specific set of data associated with peptides targeting E. coli bacteria. LASSO regression and Support Vector Machine techniques have been utilized to select, among more than 1500 physicochemical descriptors, the most important features that can be used to classify a peptide as antimicrobial or ineffective against E. coli. We then performed the classification of active versus inactive AMPs using the Support Vector classifiers, Logistic Regression, and Random Forest methods. This computational study allows us to make recommendations of how to design more efficient antibacterial drug therapies.

Predicting Antimicrobial Activity for Untested Peptide-Based Drugs Using Collaborative Filtering and Link Prediction.

The increase of bacterial resistance to currently available antibiotics has underlined the urgent need to develop new antibiotic drugs. Antimicrobial peptides (AMPs), alone or in combination with other peptides and/or existing antibiotics, have emerged as promising candidates for this task. However, given that there are thousands of known AMPs and an even larger number can be synthesized, it is impossible to comprehensively test all of them using standard wet lab experimental methods. These observations stimulated an application of machine-learning methods to identify promising AMPs. Currently, machine learning studies combine very different bacteria without considering bacteria-specific features or interactions with AMPs. In addition, the sparsity of current AMP data sets disqualifies the application of traditional machine-learning methods or makes the results unreliable. Here, we present a new approach, featuring neighborhood-based collaborative filtering, to predict with high accuracy a given bacteria's response to untested AMPs based on similarities between bacterial responses. Furthermore, we also developed a complementary bacteria-specific link prediction approach that can be used to visualize networks of AMP-antibiotic combinations, enabling us to propose new combinations that are likely to be effective.

Molecular machines open cell membranes.

Beyond the more common chemical delivery strategies, several physical techniques are used to open the lipid bilayers of cellular membranes. These include using electric and magnetic fields, temperature, ultrasound or light to introduce compounds into cells, to release molecular species from cells or to selectively induce programmed cell death (apoptosis) or uncontrolled cell death (necrosis). More recently, molecular motors and switches that can change their conformation in a controlled manner in response to external stimuli have been used to produce mechanical actions on tissue for biomedical applications. Here we show that molecular machines can drill through cellular bilayers using their molecular-scale actuation, specifically nanomechanical action. Upon physical adsorption of the molecular motors onto lipid bilayers and subsequent activation of the motors using ultraviolet light, holes are drilled in the cell membranes. We designed molecular motors and complementary experimental protocols that use nanomechanical action to induce the diffusion of chemical species out of synthetic vesicles, to enhance the diffusion of traceable molecular machines into and within live cells, to induce necrosis and to introduce chemical species into live cells. We also show that, by using molecular machines that bear short peptide addends, nanomechanical action can selectively target specific cell-surface recognition sites. Beyond the in vitro applications demonstrated here, we expect that molecular machines could also be used in vivo, especially as their design progresses to allow two-photon, near-infrared and radio-frequency activation.

Dynamics of chemical reactions on single nanocatalysts with heterogeneous active sites.

Modern chemical science and industries critically depend on the application of various catalytic methods. However, the underlying molecular mechanisms of these processes still remain not fully understood. Recent experimental advances that produced highly-efficient nanoparticle catalysts allowed researchers to obtain more quantitative descriptions, opening the way to clarify the microscopic picture of catalysis. Stimulated by these developments, we present a minimal theoretical model that investigates the effect of heterogeneity in catalytic processes at the single-particle level. Using a discrete-state stochastic framework that accounts for the most relevant chemical transitions, we explicitly evaluated the dynamics of chemical reactions on single heterogeneous nanocatalysts with different types of active sites. It is found that the degree of stochastic noise in nanoparticle catalytic systems depends on several factors that include the heterogeneity of catalytic efficiencies of active sites and distinctions between chemical mechanisms on different active sites. The proposed theoretical approach provides a single-molecule view of heterogeneous catalysis and also suggests possible quantitative routes to clarify some important molecular details of nanocatalysts.

Dynamics of single-base editing: Theoretical analysis.

Recent experimental advances led to the development of DNA base editors (BEs) with single-nucleotide precision, which is critical for future progress in various scientific and technological fields. The molecular mechanisms of single-base discrimination, however, remain poorly understood. Using a recently developed stochastic approach, we theoretically investigated the dynamics of single-base editing. More specifically, transient and mean times to edit "TC" motifs by cytosine BEs are explicitly evaluated for correct (target) and incorrect (bystander) locations on DNA. In addition, the effect of mutations on the dynamics of the single-base edition is also analyzed. It is found that for most ranges of parameters, it is possible to temporarily separate target and bystander products of base editing, supporting the idea of dynamic selectivity as a method of improving the precision of single-base editing. We conclude that to improve the efficiency of single-base editing, selecting the probability or selecting the time requires different strategies. Physical-chemical arguments to explain the observed dynamic properties are presented. The theoretical analysis clarifies some important aspects of the molecular mechanisms of selective base editing.

Kinetic control of stationary flux ratios for a wide range of biochemical processes.

One of the most intriguing features of biological systems is their ability to regulate the steady-state fluxes of the underlying biochemical reactions; however, the regulatory mechanisms and their physicochemical properties are not fully understood. Fundamentally, flux regulation can be explained with a chemical kinetic formalism describing the transitions between discrete states, with the reaction rates defined by an underlying free energy landscape. Which features of the energy landscape affect the flux distribution? Here we prove that the ratios of the steady-state fluxes of quasi-first-order biochemical processes are invariant to energy perturbations of the discrete states and are only affected by the energy barriers. In other words, the nonequilibrium flux distribution is under kinetic and not thermodynamic control. We illustrate the generality of this result for three biological processes. For the network describing protein folding along competing pathways, the probabilities of proceeding via these pathways are shown to be invariant to the stability of the intermediates or to the presence of additional misfolded states. For the network describing protein synthesis, the error rate and the energy expenditure per peptide bond is proven to be independent of the stability of the intermediate states. For molecular motors such as myosin-V, the ratio of forward to backward steps and the number of adenosine 5'-triphosphate (ATP) molecules hydrolyzed per step is demonstrated to be invariant to energy perturbations of the intermediate states. These findings place important constraints on the ability of mutations and drug perturbations to affect the steady-state flux distribution for a wide class of biological processes.

Assembly of Synaptic Protein-DNA Complexes: Critical Role of Non-Specific Interactions.

The synaptic protein-DNA complexes, formed by specialized proteins that bridge two or more distant sites on DNA, are critically involved in various genetic processes. However, the molecular mechanism by which the protein searches for these sites and how it brings them together is not well understood. Our previous studies directly visualized search pathways used by SfiI, and we identified two pathways, DNA threading and site-bound transfer pathways, specific to the site-search process for synaptic DNA-protein systems. To investigate the molecular mechanism behind these site-search pathways, we assembled complexes of SfiI with various DNA substrates corresponding to different transient states and measured their stability using a single-molecule fluorescence approach. These assemblies corresponded to specific-specific (synaptic), non-specific-non-specific (non-specific), and specific-non-specific (pre-synaptic) SfiI-DNA states. Unexpectedly, an elevated stability in pre-synaptic complexes assembled with specific and non-specific DNA substrates was found. To explain these surprising observations, a theoretical approach that describes the assembly of these complexes and compares the predictions with the experiment was developed. The theory explains this effect by utilizing entropic arguments, according to which, after the partial dissociation, the non-specific DNA template has multiple possibilities of rebinding, effectively increasing the stability. Such difference in the stabilities of SfiI complexes with specific and non-specific DNA explains the utilization of threading and site-bound transfer pathways in the search process of synaptic protein-DNA complexes discovered in the time-lapse AFM experiments.

Optimal pathways control fixation of multiple mutations during cancer initiation.

Cancer starts after initially healthy tissue cells accumulate several specific mutations or other genetic alterations. The dynamics of tumor formation is a very complex phenomenon due to multiple involved biochemical and biophysical processes. It leads to a very large number of possible pathways on the road to final fixation of all mutations that marks the beginning of the cancer, complicating the understanding of microscopic mechanisms of tumor formation. We present a new theoretical framework of analyzing the cancer initiation dynamics by exploring the properties of effective free-energy landscape of the process. It is argued that although there are many possible pathways for the fixation of all mutations in the system, there are only a few dominating pathways on the road to tumor formation. The theoretical approach is explicitly tested in the system with only two mutations using analytical calculations and Monte Carlo computer simulations. Excellent agreement with theoretical predictions is found for a large range of parameters, supporting our hypothesis and allowing us to better understand the mechanisms of cancer initiation. Our theoretical approach clarifies some important aspects of microscopic processes that lead to tumor formation.

DLPacker: Deep learning for prediction of amino acid side chain conformations in proteins.

Prediction of side chain conformations of amino acids in proteins (also termed "packing") is an important and challenging part of protein structure prediction with many interesting applications in protein design. A variety of methods for packing have been developed but more accurate ones are still needed. Machine learning (ML) methods have recently become a powerful tool for solving various problems in diverse areas of science, including structural biology. In this study, we evaluate the potential of deep neural networks (DNNs) for prediction of amino acid side chain conformations. We formulate the problem as image-to-image transformation and train a U-net style DNN to solve the problem. We show that our method outperforms other physics-based methods by a significant margin: reconstruction RMSDs for most amino acids are about 20% smaller compared to SCWRL4 and Rosetta Packer with RMSDs for bulky hydrophobic amino acids Phe, Tyr, and Trp being up to 50% smaller.

The role of spatial structures of tissues in cancer initiation dynamics.

It is widely believed that biological tissues evolved to lower the risks of cancer development. One of the specific ways to minimize the chances of tumor formation comes from proper spatial organization of tissues. However, the microscopic mechanisms of underlying processes remain not fully understood. We present a theoretical investigation on the role of spatial structures in cancer initiation dynamics. In our approach, the dynamics of single mutation fixations are analyzed using analytical calculations and computer simulations by mapping them to Moran processes on graphs with different connectivity that mimic various spatial structures. It is found that while the fixation probability is not affected by modifying the spatial structures of the tissues, the fixation times can change dramatically. The slowest dynamics is observed in 'quasi-one-dimensional' structures, while the fastest dynamics is observed in 'quasi-three-dimensional' structures. Theoretical calculations also suggest that there is a critical value of the degree of graph connectivity, which mimics the spatial dimension of the tissue structure, above which the spatial structure of the tissue has no effect on the mutation fixation dynamics. An effective discrete-state stochastic model of cancer initiation is utilized to explain our theoretical results and predictions. Our theoretical analysis clarifies some important aspects on the role of the tissue spatial structures in the cancer initiation processes.

Cation-π Interactions and Their Role in Assembling Collagen Triple Helices.

Cation-π interactions play a significant role in the stabilization of globular proteins. However, their role in collagen triple helices is less well understood and they have rarely been used in de novo designed collagen mimetic systems. In this study, we analyze the stabilizing and destabilizing effects in pairwise amino acid interactions between cationic and aromatic residues in both axial and lateral sequential relationships. Thermal unfolding experiments demonstrated that only axial pairs are stabilizing, while the lateral pairs are uniformly destabilizing. Molecular dynamics simulations show that pairs with an axial relationship can achieve a near-ideal interaction distance, but pairs in a lateral relationship do not. Arginine-π systems were found to be more stabilizing than lysine-π and histidine-π. Arginine-π interactions were then studied in more chemically diverse ABC-type heterotrimeric helices, where arginine-tyrosine pairs were found to form the best helix. This work helps elucidate the role of cation-π interactions in triple helices and illustrates their utility in designing collagen mimetic peptides.

Theoretical study of active secondary transport: Unexpected differences in molecular mechanisms for antiporters and symporters.

Successful functioning of biological cells relies on efficient translocation of different materials across cellular membranes. An important part of this transportation system is membrane channels that are known as antiporters and symporters. They exploit the energy stored as a trans-membrane gradient of one type of molecules to transport the other types of molecules against their gradients. For symporters, the directions of both fluxes for driving and driven species coincide, while for antiporters, the fluxes move in opposite directions. There are surprising experimental observations that despite differing only by the direction of transport fluxes, the molecular mechanisms of translocation adopted by antiporters and symporters seem to be drastically different. We present chemical-kinetic models to quantitatively investigate this phenomenon. Our theoretical approach allows us to explain why antiporters mostly utilize a single-site transportation when only one molecule of any type might be associated with the channel. At the same time, the transport in symporters requires two molecules of different types to be simultaneously associated with the channel. In addition, we investigate the kinetic constraints and efficiency of symporters and compare them with the same properties of antiporters. Our theoretical analysis clarifies some important physical-chemical features of cellular trans-membrane transport.

Nanorings to Probe Mechanical Stress of Single-Stranded DNA Mediated by the DNA Duplex.

The interplay between the mechanical properties of double-stranded and single-stranded DNA is a phenomenon that contributes to various genetic processes in which both types of DNA structures coexist. Highly stiff DNA duplexes can stretch single-stranded DNA (ssDNA) segments between the duplexes in a topologically constrained domain. To evaluate such an effect, we designed short DNA nanorings in which a DNA duplex with 160 bp is connected by a 30 nt single-stranded DNA segment. The stretching effect of the duplex in such a DNA construct can lead to the elongation of ssDNA, and this effect can be measured directly using atomic force microscopy (AFM) imaging. In AFM images of the nanorings, the ssDNA regions were identified, and the end-to-end distance of ssDNA was measured. The data revealed a stretching of the ssDNA segment with a median end-to-end distance which was 16% higher compared with the control. These data are in line with theoretical estimates of the stretching of ssDNA by the rigid DNA duplex holding the ssDNA segment within the nanoring construct. Time-lapse AFM data revealed substantial dynamics of the DNA rings, allowing for the formation of transient crossed nanoring formations with end-to-end distances as much as 30% larger than those of the longer-lived morphologies. The generated nanorings are an attractive model system for investigation of the effects of mechanical stretching of ssDNA on its biochemical properties, including interaction with proteins.

Temporal order of mutations influences cancer initiation dynamics.

Cancer is a set of genetic diseases that are driven by mutations. It was recently discovered that the temporal order of genetic mutations affects the cancer evolution and even the nature of the decease itself. The mechanistic origin of these observations, however, remain not well understood. Here we present a theoretical model for cancer initiation dynamics that allows us to quantify the impact of the temporal order of mutations. In our approach, the cancer initiation process is viewed as a set of stochastic transitions between discrete states defined by the different numbers of mutated cells. Using a first-passage analysis, probabilities and times before the cancer initiation are explicitly evaluated for two alternative sequences of two mutations. It is found that the probability of cancer initiation is determined only by the first mutation, while the dynamics depends on both mutations. In addition, it is shown that the acquisition of a mutation with higher fitness before mutation with lower fitness increases the probability of the tumor formation but delays the cancer initiation. Theoretical results are explained using effective free-energy landscapes.

Charge-Free, Stabilizing Amide-π Interactions Can Be Used to Control Collagen Triple-Helix Self-Assembly.

There is a noted lack of understood, controllable interactions for directing the organization of collagen triple helices. While the field has had success using charge-pair interactions and cation-π interactions in helix design, these alone are not adequate for achieving the degree of specificity desirable for these supramolecular structures. Furthermore, because of the reliance on electrostatic interactions, designed heterotrimeric systems have been heavily charged, a property undesirable in some applications. Amide-π interactions are a comparatively understudied class of charge-free interactions, which could potentially be harnessed for triple-helix design. Herein, we propose, validate, and utilize pairwise amino acid amide-π interactions in collagen triple-helix design. Glutamine-phenylalanine pairs, when arranged in an axial geometry, are found to exhibit a moderately stabilizing effect, while in the lateral geometry, this pair is destabilizing. Together this allows glutamine-phenylalanine pairs to effectively set the register of triple helices. In contrast, interactions between asparagine and phenylalanine appear to have little effect on triple-helical stability. After deconvoluting the contributions of these amino acids to triple-helix stability, we demonstrate these new glutamine-phenylalanine interactions in the successful design of a heterotrimeric triple helix. The results of all of these analyses are used to update our collagen triple-helix thermal stability prediction algorithm, Scoring function for Collagen Emulating Peptides' Temperature of Transition (SCEPTTr).

Understanding the molecular mechanisms of transcriptional bursting.

In recent years, it has been experimentally established that transcription, a fundamental biological process that involves the synthesis of messenger RNA molecules from DNA templates, does not proceed continuously as was expected. Rather, it exhibits a distinct dynamic behavior of alternating between productive phases when RNA molecules are actively synthesized and inactive phases when there is no RNA production at all. The bimodal transcriptional dynamics is now confirmed to be present in most living systems. This phenomenon is known as transcriptional bursting and it attracts significant amounts of attention from researchers in different fields. However, despite multiple experimental and theoretical investigations, the microscopic origin and biological functions of the transcriptional bursting remain unclear. Here we discuss the recent developments in uncovering the underlying molecular mechanisms of transcriptional bursting and our current understanding of them. Our analysis presents a physicochemical view of the processes that govern transcriptional bursting in living cells.