Effects of electrostatic interactions on ligand dissociation kinetics

Multiscale modeling reveals the ion-mediated phase separation of nucleosome core particles

Due to the vast length scale inside the cell nucleus, multiscale models are required to understand chromatin folding, structure, and dynamics and how they regulate genomic activities such as DNA transcription, replication, and repair. We study the interactions and structure of condensed phases formed by the universal building block of chromatin, the nucleosome core particle (NCP), using bottom-up multiscale coarse-grained (CG) simulations with a model extracted from all-atom MD simulations. In the presence of the multivalent cations Mg(H2O)62+ or CoHex3+, we analyze the internal structures of the NCP aggregates and the contributions of histone tails and ions to the aggregation patterns. We then derive a "super" coarse-grained (SCG) NCP model to study the macroscopic scale phase separation of NCPs. The SCG simulations show the formation of NCP aggregates with Mg(H2O)62+ concentration-dependent densities and sizes. Variation of the CoHex3+ concentrations results in highly ordered lamellocolumnar and hexagonal columnar phases in agreement with experimental data. The results give detailed insights into nucleosome interactions and for understanding chromatin folding in the cell nucleus.

Potent activity of polymyxin B is associated with long-lived super-stoichiometric accumulation mediated by weak-affinity binding to lipid A

Polymyxins are gram-negative antibiotics that target lipid A, the conserved membrane anchor of lipopolysaccharide in the outer membrane. Despite their clinical importance, the molecular mechanisms underpinning polymyxin activity remain unresolved. Here, we use surface plasmon resonance to kinetically interrogate interactions between polymyxins and lipid A and derive a phenomenological model. Our analyses suggest a lipid A-catalyzed, three-state mechanism for polymyxins: transient binding, membrane insertion, and super-stoichiometric cluster accumulation with a long residence time. Accumulation also occurs for brevicidine, another lipid A-targeting antibacterial molecule. Lipid A modifications that impart polymyxin resistance and a non-bactericidal polymyxin derivative exhibit binding that does not evolve into long-lived species. We propose that transient binding to lipid A permeabilizes the outer membrane and cluster accumulation enables the bactericidal activity of polymyxins. These findings could establish a blueprint for discovery of lipid A-targeting antibiotics and provide a generalizable approach to study interactions with the gram-negative outer membrane.

Dental biomaterials redefined: molecular docking and dynamics-driven dental resin composite optimization

BACKGROUND

Dental resin-based composites are widely recognized for their aesthetic appeal and adhesive properties, which make them integral to modern restorative dentistry. Despite their advantages, adhesion and biomechanical performance challenges persist, necessitating innovative strategies for improvement. This study addressed the challenges associated with adhesion and biomechanical properties in dental resin-based composites by employing molecular docking and dynamics simulation.

METHODS

Molecular docking assesses the binding energies and provides valuable insights into the interactions between monomers, fillers, and coupling agents. This investigation prioritizes SiO2 and TRIS, considering their consistent influence. Molecular dynamics simulations, executed with the Forcite module and COMPASS II force field, extend the analysis to the mechanical properties of dental composite complexes. The simulations encompassed energy minimization, controlled NVT and NPT ensemble simulations, and equilibration stages. Notably, the molecular dynamics simulations spanned a duration of 50 ns.

RESULTS

SiO2 and TRIS consistently emerged as influential components, showcasing their versatility in promoting solid interactions. A correlation matrix underscores the significant roles of van der Waals and desolvation energies in determining the overall binding energy. Molecular dynamics simulations provide in-depth insights into the mechanical properties of dental composite complexes. HEMA-SiO2-TRIS excelled in stiffness, BisGMA-SiO2-TRIS prevailed in terms of flexural strength, and EBPADMA-SiO2-TRIS offered a balanced combination of mechanical properties.

CONCLUSION

These findings provide valuable insights into optimizing dental composites tailored to diverse clinical requirements. While EBPADMA-SiO2-TRIS demonstrates distinct strengths, this study emphasizes the need for further research. Future investigations should validate the computational findings experimentally and assess the material's response to dynamic environmental factors.

Assessment of binding interaction to salmon sperm DNA of two antiviral agents and ecofriendly nanoparticles: comprehensive spectroscopic study

The direct binding of antiviral agents; Daclatasvir and valacyclovir and green synthesized nanoparticles to salmon sperm DNA have been assessed in a comparative study. The nanoparticles were synthesized by the hydrothermal autoclave method and have been fully characterized. The interactive behavior and competitive binding of the analytes to DNA in addition to the thermodynamic properties were deeply investigated by the UV-visible spectroscopy. The binding constants were monitored in the physiological pH conditions to be 1.65 × 10⁶, 4.92 × 10⁵ and 3.12 × 10⁵ for daclatasvir,valacyclovir and quantum dots, respectively. The significant changes in the spectral features of all analytes have proven intercalative binding. The competitive study has confirmed that, daclatasvir, valacyclovir, and the quantum dots have exhibited groove binding. All analytes have shown good entropy and enthalpy values indicating stable interactions. The electrostatic and non-electrostatic kinetic parameters have been determined through studying the binding interactions at different concentrations of KCl solutions. A molecular modelling study has been applied to demonstrate the binding interactions and their mechanisms. The obtained results were complementary and afforded new eras for the therapeutic applications.

Polymorphisms in alternative oxidase genes from ecotypes of Arabidopsis and rice revealed an environment-induced linkage to altitude and rainfall

We investigated SNPs in alternative oxidase (AOX) genes and their connection to ecotype origins (climate, altitude, and rainfall) by using genomic data sets of Arabidopsis and rice populations from 1190 and 90 ecotypes, respectively. Parameters were defined to detect non-synonymous SNPs in the AOX ORF, which revealed amino acid (AA) changes in AOX1c, AOX1d, and AOX2 from Arabidopsis and AOX1c from rice in comparison to AOX references from Columbia-0 and Japonica ecotypes, respectively. Among these AA changes, Arabidopsis AOX1c_A161E&G165R and AOX1c_R242S revealed a link to high rainfall and high altitude, respectively, while all other changes in Arabidopsis and rice AOX was connected to high altitude and rainfall. Comparative 3D modeling showed that all mutant AOX presented structural differences in relation to the respective references. Molecular docking analysis uncovered lower binding affinity values between AOX and the substrate ubiquinol for most of the identified structures compared to their reference, indicating better enzyme-substrate binding affinities. Thus, our in silico data suggest that the majority of the AA changes found in the available ecotypes will confer better enzyme-subtract interactions and thus indicate environment-related, more efficient AOX activity.

Local Ion Densities can Influence Transition Paths of Molecular Binding

Improper reaction coordinates can pose significant problems for path-based binding free energy calculations. Particularly, omission of long timescale motions can lead to over-estimation of the energetic barriers between the bound and unbound states. Many methods exist to construct the optimal reaction coordinate using a pre-defined basis set of features. Although simulations are typically conducted in explicit solvent, the solvent atoms are often excluded by these feature sets—resulting in little being known about their role in reaction coordinates, and ultimately, their role in determining (un)binding rates and free energies. In this work, analysis is done on an extensive set of host-guest unbinding trajectories, working to characterize differences between high and low probability unbinding trajectories with a focus on solvent-based features, including host-ion interactions, guest-ion interactions and location-dependent ion densities. We find that differences in ion densities as well as guest-ion interactions strongly correlate with differences in the probabilities of reactive paths that are used to determine free energies of (un)binding and play a significant role in the unbinding process.

The Role of Ligand Rebinding and Facilitated Dissociation on the Characterization of Dissociation Rates by Surface Plasmon Resonance (SPR) and Benchmarking Performance Metrics

Surface plasmon resonance (SPR) is a real-time kinetic measurement principle that can probe the kinetic interactions between ligands and their binding sites, and lies at the backbone of pharmaceutical, biosensing, and biomolecular research. The extraction of dissociation rates from SPR-response signals often relies on several commonly adopted assumptions, one of which is the exponential decay of the dissociation part of the response signal. However, certain conditions, such as high density of binding sites or high concentration fluctuations near the surface as compared to the bulk, can lead to non-exponential decays via ligand rebinding or facilitated dissociation. Consequently, fitting the data with an exponential function can underestimate or overestimate the measured dissociation rates. Here, we describe a set of alternative fit functions that can take such effects into consideration along with plasmonic sensor design principles with key performance metrics, thereby suggesting methods for error-free high-precision extraction of the dissociation rates.

How do DNA-bound proteins leave their binding sites? The role of facilitated dissociation

Dissociation of a protein from DNA is often assumed to be described by an off rate that is independent of other molecules in solution. Recent experiments and computational analyses have challenged this view by showing that unbinding rates (residence times) of DNA-bound proteins can depend on concentrations of nearby molecules that are competing for binding. This 'facilitated dissociation' (FD) process can occur at the single-binding site level via formation of a ternary complex, and can dominate over 'spontaneous dissociation' at low (submicromolar) concentrations. In the crowded intracellular environment FD introduces new regulatory possibilities at the level of individual biomolecule interactions.

Force-Dependent Facilitated Dissociation Can Generate Protein-DNA Catch Bonds

Cellular structures are continually subjected to forces, which may serve as mechanical signals for cells through their effects on biomolecule interaction kinetics. Typically, molecular complexes interact via "slip bonds," so applied forces accelerate off rates by reducing transition energy barriers. However, biomolecules with multiple dissociation pathways may have considerably more complicated force dependencies. This is the case for DNA-binding proteins that undergo "facilitated dissociation," in which competitor biomolecules from solution enhance molecular dissociation in a concentration-dependent manner. Using simulations and theory, we develop a generic model that shows that proteins undergoing facilitated dissociation can form an alternative type of molecular bond, known as a "catch bond," for which applied forces suppress protein dissociation. This occurs because the binding by protein competitors responsible for the facilitated dissociation pathway can be inhibited by applied forces. Within the model, we explore how the force dependence of dissociation is regulated by intrinsic factors, including molecular sensitivity to force and binding geometry and the extrinsic factor of competitor protein concentration. We find that catch bonds generically emerge when the force dependence of the facilitated unbinding pathway is stronger than that of the spontaneous unbinding pathway. The sharpness of the transition between slip- and catch-bond kinetics depends on the degree to which the protein bends its DNA substrate. This force-dependent kinetics is broadly regulated by the concentration of competitor biomolecules in solution. Thus, the observed catch bond is mechanistically distinct from other known physiological catch bonds because it requires an extrinsic factor-competitor proteins-rather than a specific intrinsic molecular structure. We hypothesize that this mechanism for regulating force-dependent protein dissociation may be used by cells to modulate protein exchange, regulate transcription, and facilitate diffusive search processes.

Crystallizing protein assemblies via free and grafted linkers

Porous protein superlattices have plausible catalytic applications in biotechnology and nanotechnology. They are solid yet open structures with the potential for preserving the activity of enzymes. However, there is still a lack of understanding of the design parameters that are required to arrange proteins in a periodic porous fashion. Here, we introduce a coarse-grained molecular dynamics (MD) simulation approach to study the effects of the lengths and geometries of linkers on the stability of 3D crystalline assemblies of metal ion anchored ferritin proteins. By simulating a system of proteins (eight metal ion anchored sites per protein) and linkers (two free ends per linker), we find that there is a range of optimal linker lengths for crystalline order. The optimal linker length is found to depend on the linker to protein concentration ratio and binding energy. We also examine the case of grafted flexible linkers on the protein surface as an alternative route for constructing highly porous crystalline structures. Our study demonstrates that the length of grafted linkers is a better tunable parameter than the length of free linkers to achieve high porosity protein superlattices. The computational study developed here provides guidelines to assemble biomolecules into crystals with high porosity.

Receptor-Ligand Rebinding Kinetics in Confinement

Rebinding kinetics of molecular ligands plays a key role in the operation of biomachinery, from regulatory networks to protein transcription, and is also a key factor in design of drugs and high-precision biosensors. In this study, we investigate initial release and rebinding of ligands to their binding sites grafted on a planar surface, a situation commonly observed in single-molecule experiments and that occurs in vivo, e.g., during exocytosis. Via scaling arguments and molecular dynamic simulations, we analyze the dependence of nonequilibrium rebinding kinetics on two intrinsic length scales: the average separation distance between the binding sites and the total diffusible volume (i.e., height of the experimental reservoir in which diffusion takes place or average distance between receptor-bearing surfaces). We obtain time-dependent scaling laws for on rates and for the cumulative number of rebinding events. For diffusion-limited binding, the (rebinding) on rate decreases with time via multiple power-law regimes before the terminal steady-state (constant on-rate) regime. At intermediate times, when particle density has not yet become uniform throughout the diffusible volume, the cumulative number of rebindings exhibits a novel, to our knowledge, plateau behavior because of the three-dimensional escape process of ligands from binding sites. The duration of the plateau regime depends on the average separation distance between binding sites. After the three-dimensional diffusive escape process, a one-dimensional diffusive regime describes on rates. In the reaction-limited scenario, ligands with higher affinity to their binding sites (e.g., longer residence times) delay entry to the power-law regimes. Our results will be useful for extracting hidden timescales in experiments such as kinetic rate measurements for ligand-receptor interactions in microchannels, as well as for cell signaling via diffusing molecules.