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.

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.

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.

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.

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.

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.

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.

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.

Dynamics of Single-Base Editing: Theoretical Analysis.

Recent experimental advances led to the development of DNA base editors (BEs) with a single-nucleotide precision that is critical for future progress in various scientific and technological fields. The molecular mechanisms of single-base discrimination, however, remain not well 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 molecular mechanisms of selective base editing.

Synergy among Pausing, Intrinsic Proofreading, and Accessory Proteins Results in Optimal Transcription Speed and Tolerable Accuracy.

Cleavage of dinucleotides after the misincorporational pauses serves as a proofreading mechanism that increases transcriptional elongation accuracy. The accuracy is further improved by accessory proteins such as GreA and TFIIS. However, it is not clear why RNAP pauses and why cleavage-factor-assisted proofreading is necessary despite transcriptional errors in vitro being of the same order as those in downstream translation. Here, we developed a chemical-kinetic model that incorporates most relevant features of transcriptional proofreading and uncovers how the balance between speed and accuracy is achieved. We found that long pauses are essential for high accuracy, whereas cleavage-factor-stimulated proofreading optimizes speed. Moreover, in comparison to the cleavage of a single nucleotide or three nucleotides, RNAP backtracking and dinucleotide cleavage improve both speed and accuracy. Our results thereby show how the molecular mechanism and the kinetic parameters of the transcriptional process were evolutionarily optimized to achieve maximal speed and tolerable accuracy.

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.

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.

Beyond Sequence: Internucleosomal Interactions Dominate Array Assembly.

The organization of the nucleosome array is a critical component of the chromatin assembly into higher order structure as well as its function. Here, we investigated the contributions of the DNA sequence and internucleosomal interactions on the organization of the nucleosomal arrays in compact structures using atomic force microscopy. We assembled nucleosomes on DNA substrates allowing for the formation of tetranucleosomes. We found that nucleosomes are capable of close positioning with no discernible space between them, even in the case of assembled dinucleosomes. This morphology of the array is in contrast with that observed for arrays assembled with repeats of the nucleosome positioning motifs separated by uniform spacers. Simulated assembly of tetranucleosomes by random placement along the substrates revealed that nucleosome array compaction is promoted by the interaction of the nucleosomes. We developed a theoretical model to account for the role of DNA sequence and internucleosomal interactions in the formation of the nucleosome structures. These findings suggest that, in the chromatin assembly, the affinity of the nucleosomes to the DNA sequence and the strengths of the internucleosomal interactions are the two major factors defining the compactness of the chromatin.

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.

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.