Feature Selection Enhances Peptide Binding Predictions for TCR-Specific Interactions.

T-cell receptors (TCRs) play a critical role in the immune response by recognizing specific ligand peptides presented by major histocompatibility complex (MHC) molecules. Accurate prediction of peptide binding to TCRs is essential for advancing immunotherapy, vaccine design, and understanding mechanisms of autoimmune disorders. This study presents a novel theoretical method that explores the impact of feature selection techniques on enhancing the predictive accuracy of peptide binding models tailored for specific TCRs. To evaluate the universality of our approach across different TCR systems, we utilized a dataset that includes peptide libraries tested against three distinct murine TCRs. A broad range of physicochemical properties, including amino acid composition, dipeptide composition, and tripeptide features, were integrated into the machine learning-based feature selection framework to identify key features contributing to binding affinity. Our analysis reveals that leveraging optimized feature subsets not only simplifies the model complexity but also enhances predictive performance, enabling more precise identification of TCR-peptide interactions. The results of our feature selection method are consistent with findings from hybrid approaches that utilize both sequence and structural data as input as well as experimental data. Our theoretical approach highlights the role of feature selection in peptide-TCR interactions, providing a powerful tool for uncovering the molecular mechanisms of the T-cell response and assisting in the design of more advanced targeted therapeutics.

Linear-Decoupling Enables Accurate Speed and Accuracy Predictions for Copolymerization Processes.

Biological processes exhibit remarkable accuracy and speed and can be theoretically explored through various approaches. The Markov-chain copolymerization theory, describing polymer growth kinetics as a Markov chain, provides an exact set of equations to solve for error and speed. Still, due to nonlinearity, these equations are hard to solve. Alternatively, the enzyme-kinetics approach, which formulates a set of linear equations, simplifies the biological processes as transitions between discrete chemical states, but generally, it might not be accurate. Here, we show that the enzyme-kinetic approach can lead to inaccurate fluxes, even for first-order polymerization processes. To address the problem, we propose a simplified linear-decoupling approximation for steady-state probabilities of higher-order copolymer chains under biologically relevant conditions. Our findings demonstrate that the stationary speed and error rate obtained from the linear-decoupling method align closely with exact values from the Markov-chain (nonlinear) approximation. Extending the technique to higher-order processes with proofreading and internal states shows that it works equally well to describe trade-offs between speed and accuracy for DNA replication and transcription elongation. Our work underscores the proposed linear-decoupling approximation's efficacy in addressing the nonlinear behavior of the Markov-chain approach and the enzyme-kinetic approach's limitations, ensuring accurate predictions for high-fidelity biological processes.

How Transcription Factors Binding Stimulates Transcriptional Bursting.

Transcription is a fundamental biological process of transferring genetic information which often occurs in stochastic bursts when periods of intense activity alternate with quiescent phases. Recent experiments identified strong correlations between the association of transcription factors (TFs) to gene promoters on DNA and transcriptional activity. However, the underlying molecular mechanisms of this phenomenon remain not well understood. Here, we present a theoretical framework that allowed us to investigate how binding dynamics of TF influences transcriptional bursting. Our minimal theoretical model incorporates the most relevant physical-chemical features, including TF exchange among multiple binding sites at gene promoters and TF association/dissociation dynamics. Using analytical calculations supported by Monte Carlo computer simulations, it is demonstrated that transcriptional bursting dynamics depends on the strength of TF binding and the number of binding sites. Stronger TF binding affinity prolongs burst duration but reduces variability, while an optimal number of binding sites maximizes transcriptional noise, facilitating cellular adaptation. Our theoretical method explains available experimental observations quantitatively, confirming the model's predictive accuracy. This study provides important insights into molecular mechanisms of gene expression and regulation, offering a new theoretical tool for understanding complex biological processes.

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.

Unraveling the role of physicochemical differences in predicting protein-protein interactions.

The ability to accurately predict protein-protein interactions is critically important for understanding major cellular processes. However, current experimental and computational approaches for identifying them are technically very challenging and still have limited success. We propose a new computational method for predicting protein-protein interactions using only primary sequence information. It utilizes the concept of physicochemical similarity to determine which interactions will most likely occur. In our approach, the physicochemical features of proteins are extracted using bioinformatics tools for different organisms. Then they are utilized in a machine-learning method to identify successful protein-protein interactions via correlation analysis. It was found that the most important property that correlates most with the protein-protein interactions for all studied organisms is dipeptide amino acid composition (the frequency of specific amino acid pairs in a protein sequence). While current approaches often overlook the specificity of protein-protein interactions with different organisms, our method yields context-specific features that determine protein-protein interactions. The analysis is specifically applied to the bacterial two-component system that includes histidine kinase and transcriptional response regulators, as well as to the barnase-barstar complex, demonstrating the method's versatility across different biological systems. Our approach can be applied to predict protein-protein interactions in any biological system, providing an important tool for investigating complex biological processes' mechanisms.

Elucidating Physicochemical Features of Holin Proteins Responsible for Bacterial Cell Lysis.

Bacterial resistance to conventional antibiotics stimulated the development of so-called "phage therapies" that rely on cell lysis, which is a process of destroying bacterial cells due to their infections by bacterial viruses. For λ bacteriophages, it is known that the critical role in this process is played by holin proteins that aggregate in cellular membranes before breaking them apart. While multiple experimental studies probed various aspects of cell lysis, the underlying molecular mechanisms remain not well understood. Here we investigate what physicochemical properties of holin proteins are the most relevant for these processes by employing statistical correlation analysis of cell lysis dynamics for different experimentally observed mutant species. Our findings reveal significant correlations between various physicochemical features and cell lysis dynamics. Notably, we uncover a strong inverse correlation between local hydrophobicity and cell lysis times, underscoring the crucial role of hydrophobic interactions in membrane disruption. Stimulated by these observations, a predictive model capable of explicitly estimating cell lysis times for any holin protein mutants based on their mean hydrophobicity values is developed. Our study not only provides important microscopic insights into cell lysis phenomena but also proposes specific routes to optimize medical and biotechnological applications of bacteriophages.

Heterotrimeric collagen helix with high specificity of assembly results in a rapid rate of folding.

The most abundant natural collagens form heterotrimeric triple helices. Synthetic mimics of collagen heterotrimers have been found to fold slowly, even compared to the already slow rates of homotrimeric helices. These prolonged folding rates are not understood. Here we compare the stabilities, specificities and folding rates of three heterotrimeric collagen mimics designed through a computationally assisted approach. The crystal structure of one ABC-type heterotrimer verified a well-controlled composition and register and elucidated the geometry of pairwise cation-π and axial and lateral salt bridges in the assembly. This collagen heterotrimer folds much faster (hours versus days) than comparable, well-designed systems. Circular dichroism and NMR data suggest the folding is frustrated by unproductive, competing heterotrimer species and these species must unwind before refolding into the thermodynamically favoured assembly. The heterotrimeric collagen folding rate is inhibited by the introduction of preformed competing triple-helical assemblies, which suggests that slow heterotrimer folding kinetics are dominated by the frustration of the energy landscape caused by competing triple helices.

Molecular mechanisms of precise timing in cell lysis.

Many biological systems exhibit precise timing of events, and one of the most known examples is cell lysis, which is a process of breaking bacterial host cells in the virus infection cycle. However, the underlying microscopic picture of precise timing remains not well understood. We present a novel theoretical approach to explain the molecular mechanisms of effectively deterministic dynamics in biological systems. Our hypothesis is based on the idea of stochastic coupling between relevant underlying biophysical and biochemical processes that lead to noise cancellation. To test this hypothesis, we introduced a minimal discrete-state stochastic model to investigate how holin proteins produced by bacteriophages break the inner membranes of gram-negative bacteria. By explicitly solving this model, the dynamic properties of cell lysis are fully evaluated, and theoretical predictions quantitatively agree with available experimental data for both wild-type and holin mutants. It is found that the observed threshold-like behavior is a result of the balance between holin proteins entering the membrane and leaving the membrane during the lysis. Theoretical analysis suggests that the cell lysis achieves precise timing for wild-type species by maximizing the number of holins in the membrane and narrowing their spatial distribution. In contrast, for mutated species, these conditions are not satisfied. Our theoretical approach presents a possible molecular picture of precise dynamic regulation in intrinsically random biological processes.

Insights into Error Control Mechanisms in Biological Processes: Copolymerization and Enzyme-Kinetics Revisited.

The high fidelity observed in biological information processing ranging from replication to translation has stimulated significant research efforts to clarify the underlying microscopic picture. Theoretically, several approaches to analyze the error rates have been proposed. The copolymerization theory describes the addition and removal of monomers at the growing tip of a copolymer, leading to a closed set of nonlinear equations. On the other hand, enzyme-kinetics approaches formulate linear equations of biochemical networks, describing transitions between discrete chemical states. However, it is still unclear whether the error values computed by the two approaches agree. Moreover, there are conflicting interpretations on whether the error is under thermodynamic or kinetic discrimination control. In this work, we examine the error rate in persistent copying biochemical processes by specifically analyzing both theoretical approaches. The initial disagreement of the results between the two theories motivated us to rederive the formula for the error rate in the kinetic model. The error computed with the new method resulted in excellent agreement between both theoretical approaches and with Monte Carlo simulations. Furthermore, our theoretical analysis shows that the kinetic discrimination controls the error, even when the energy difference between adding the right and wrong products is very small. Our theoretical investigation gives important insights into the physical-chemical properties of complex biological processes by providing the quantitative framework to evaluate them.

Physical-Chemical Features Selection Reveals That Differences in Dipeptide Compositions Correlate Most with Protein-Protein Interactions.

The ability to accurately predict protein-protein interactions is critically important for our understanding of major cellular processes. However, current experimental and computational approaches for identifying them are technically very challenging and still have limited success. We propose a new computational method for predicting protein-protein interactions using only primary sequence information. It utilizes a concept of physical-chemical similarity to determine which interactions will most probably occur. In our approach, the physical-chemical features of protein are extracted using bioinformatics tools for different organisms, and then they are utilized in a machine-learning method to identify successful protein-protein interactions via correlation analysis. It is found that the most important property that correlates most with the protein-protein interactions for all studied organisms is dipeptide amino acid compositions. The analysis is specifically applied to the bacterial two-component system that includes histidine kinase and transcriptional response regulators. Our theoretical approach provides a simple and robust method for quantifying the important details of complex mechanisms of biological processes.

Differences in Relevant Physicochemical Properties Correlate with Synergistic Activity of Antimicrobial Peptides.

With the urgent need for new medical approaches due to increased bacterial resistance to antibiotics, antimicrobial peptides (AMPs) have been considered as potential treatments for infections. Experiments indicate that combinations of several types of AMPs might be even more effective at inhibiting bacterial growth with reduced toxicity and a lower likelihood of inducing bacterial resistance. The molecular mechanisms of AMP-AMP synergistic antimicrobial activity, however, remain not well understood. Here, we present a theoretical approach that allows us to relate the physicochemical properties of AMPs and their antimicrobial cooperativity. It utilizes correlation and bioinformatics analysis. A concept of physicochemical similarity is introduced, and it is found that less similar AMPs with respect to certain physicochemical properties lead to greater synergy because of their complementary antibacterial actions. The analysis of correlations between the similarity and the antimicrobial properties allows us to effectively separate synergistic from nonsynergistic AMP pairs. Our theoretical approach can be used for the rational design of more effective AMP combinations for specific bacterial targets, for clarifying the mechanisms of bacterial elimination, and for a better understanding of cooperativity phenomena in biological systems.

Why Are Nucleosome Breathing Dynamics Asymmetric?

In eukaryotic cells, DNA is bound to nucleosomes, but DNA segments occasionally unbind in the process known as nucleosome breathing. Although DNA can unwrap simultaneously from both ends of the nucleosome (symmetric breathing), experiments indicate that DNA prefers to dissociate from only one end (asymmetric breathing). However, the molecular origin of the asymmetry is not understood. We developed a new theoretical approach that gives microscopic explanations of asymmetric breathing. It is based on a stochastic description that leads to a comprehensive evaluation of dynamics by using effective free-energy landscapes. It is shown that asymmetric breathing follows the kinetically preferred pathways. In addition, it is also found that asymmetric breathing leads to a faster target search by transcription factors. Theoretical predictions, supported by computer simulations, agree with experiments. It is proposed that nature utilizes the symmetry of nucleosome breathing to achieve a better dynamic accessibility of chromatin for more efficient genetic regulation.

Brownian Diffusion of Hexagonal Boron Nitride Nanosheets and Graphene in Two Dimensions.

Two-dimensional (2D) nanomaterials have numerous interesting chemical and physical properties that make them desirable building blocks for the manufacture of macroscopic materials. Liquid-phase processing is a common method for forming macroscopic materials from these building blocks including wet-spinning and vacuum filtration. As such, assembling 2D nanomaterials into ordered functional materials requires an understanding of their solution dynamics. Yet, there are few experimental studies investigating the hydrodynamics of disk-like materials. Herein, we report the lateral diffusion of hexagonal boron nitride nanosheets (h-BN and graphene) in aqueous solution when confined in 2-dimensions. This was done by imaging fluorescent surfactant-tagged nanosheets and visualizing them by using fluorescence microscopy. Spectroscopic studies were conducted to characterize the interactions between h-BN and the fluorescent surfactant, and atomic force microscopy (AFM) was conducted to characterize the quality of the dispersion. The diffusion data under different gap sizes and viscosities displayed a good correlation with Kramers' theory. We propose that the yielded activation energies by Kramers' equation express the magnitude of the interaction between fluorescent surfactant tagged h-BN and glass because the energies remain constant with changing viscosity and decrease with increasing confinement size. The diffusion of graphene presented a similar trend with similar activation energy as the h-BN. This relationship suggests that Kramers' theory can also be applied to simulate the diffusion of other 2D nanomaterials.

How to Build Plasmon-Driven Molecular Jackhammers that Disassemble Cell Membranes and Cytoskeletons in Cancer.

Plasmon-driven molecular machines with ultrafast motion at the femtosecond scale are effective for the treatment of cancer and other diseases. It is recently shown that cyanine dyes act as molecular jackhammers (MJH) through vibronic (vibrational and electronic mode coupling) driven activation that causes the molecule to stretch longitudinally and axially through concerted whole molecule vibrations. However, the theoretical and experimental underpinnings of these plasmon-driven motions in molecules are difficult to assess. Here the use of near-infrared (NIR) light-activated plasmons in a broad array of MJH that mechanically disassemble membranes and cytoskeletons in human melanoma A375 cells is described. The characteristics of plasmon-driven molecular mechanical disassembly of supramolecular biological structures are observed and recorded using real-time fluorescence confocal microscopy. Molecular plasmon resonances in MJH are quantified through a new experimental plasmonicity index method. This is done through the measurement of the UV-vis-NIR spectra in various solvents, and quantification of the optical response as a function of the solvent polarity. Structure-activity relationships are used to optimize the synthesis of plasmon-driven MJH, applying them to eradicate human melanoma A375 cells at low lethal concentrations of 75 nm and 80 mW cm-2 of 730 nm NIR-light for 10 min.

Increasing Heterogeneity in Antimicrobial Peptide Combinations Enhances Their Synergistic Activities.

Antimicrobial peptides (AMPs) are short biopolymers produced by living organisms as an immune system defense against infections. They have been considered as potential alternatives to conventional antibiotics. Experiments suggest that combining several types of different AMPs might enhance their antimicrobial activity more effectively than using single-component AMPs. However, a clear understanding of the underlying microscopic mechanisms is still lacking. We present a theoretical investigation of antibacterial cooperativity mechanisms involving several types of AMPs. It is argued that synergy results from intermolecular interactions when the presence of one type of AMP stimulates the association of another type of AMP to bacteria. It is found that increasing the number of different AMPs in the mixtures increases the number of such interactions, making them more efficient in eliminating infections. Our theoretical framework provides valuable insights into the mechanisms of antimicrobial action.

Role of Nucleosome Sliding in the Protein Target Search for Covered DNA Sites.

Associations of transcription factors (TFs) with specific sites on DNA initiate major cellular processes. But DNA in eukaryotic cells is covered by nucleosomes which prevent TFs from binding. However, nucleosome structures on DNA are not static and exhibit breathing and sliding. We develop a theoretical framework to investigate the effect of nucleosome sliding on a protein target search. By analysis of a discrete-state stochastic model of nucleosome sliding, search dynamics are explicitly evaluated. It is found that for long sliding lengths the target search dynamics are faster for normal TFs that cannot enter the nucleosomal DNA. But for more realistic short sliding lengths, the so-called pioneer TFs, which can invade nucleosomal DNA, locate specific sites faster. It is also suggested that nucleosome breathing, which is a faster process, has a stronger effect on protein search dynamics than that of nucleosome sliding. Theoretical arguments to explain these observations are presented.

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.

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.

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.

Nucleosome Breathing Facilitates the Search for Hidden DNA Sites by Pioneer Transcription Factors.

Transfer of genetic information starts with transcription factors (TFs) binding to specific sites on DNA. But in living cells, DNA is mostly covered by nucleosomes. There are proteins, known as pioneer TFs, that can efficiently reach the DNA sites hidden by nucleosomes, although the underlying mechanisms are not understood. Using the recently proposed idea of interaction-compensation mechanism, we develop a stochastic model for the target search on DNA with nucleosome breathing. It is found that nucleosome breathing can significantly accelerate the search by pioneer TFs in comparison to situations without breathing. We argue that this is the result of the interaction-compensation mechanism that allows proteins to enter the inner nucleosome region through the outer DNA segment. It is suggested that nature optimized pioneer TFs to take advantage of nucleosome breathing. The presented theoretical picture provides a possible microscopic explanation for the successful invasion of nucleosome-buried genes.