Dynamics of water around the complex structures formed between the KH domains of far upstream element binding protein and single-stranded DNA molecules

Three-Dimensional Modeling of CpG DNA Binding with Matrix Lumican Shows Leucine-Rich Repeat Motif Involvement as in TLR9-CpG DNA Interactions

Lumican is an extracellular matrix proteoglycan known to regulate toll-like receptor (TLR) signaling in innate immune cells. In experimental settings, lumican suppresses TLR9 signaling by binding to and sequestering its synthetic ligand, CpG-DNA, in non-signal permissive endosomes. However, the molecular details of lumican interactions with CpG-DNA are obscure. Here, the 3-D structure of the 22 base-long CpG-DNA (CpG ODN_2395) bound to lumican or TLR9 were modeled using homology modeling and docking methods. Some of the TLR9-CpG ODN_2395 features predicted by our model are consistent with the previously reported TLR9-CpG DNA crystal structure, substantiating our current analysis. Our modeling indicated a smaller buried surface area for lumican-CpG ODN_2395 (1803 Å2) compared to that of TLR9-CpG ODN_2395 (2094 Å2), implying a potentially lower binding strength for lumican and CpG-DNA than TLR9 and CpG-DNA. The docking analysis identified 32 amino acids in lumican LRR1-11 interacting with CpG ODN_2395, primarily through hydrogen bonding, salt-bridges, and hydrophobic interactions. Our study provides molecular insights into lumican and CpG-DNA interactions that may lead to molecular targets for modulating TLR9-mediated inflammation and autoimmunity.

Base pair compositional variability influences DNA structural stability and tunes hydration thermodynamics and dynamics

DNA deformability and differential hydration are crucial determinants of biological processes ranging from genetic material packaging to gene expression; their associative details, however, remain inadequately understood. Herein, we report investigations of the dynamic and thermodynamic responses of the local hydration of a variety of base pair sequences. Leveraging in silico sampling and our in-house analyses, we first report the local conformational propensity of sequences that are either predisposed toward the canonical A- or B-conformations or are restrained to potential transitory pathways. It is observed that the transition from the unrestrained A-form to the B-form leads to lengthwise structural deformation. The insertion of intermittent -(CG)- base pairs in otherwise homogeneous -(AT)- sequences bears dynamical consequences for the vicinal hydration layer. Calculation of the excess (pair) entropy suggests substantially higher values of hydration water surrounding A conformations over the B- conformations. Applying the Rosenfeld approximation, we project that the diffusivity of water molecules proximal to canonical B conformation is least for the minor groove of the canonical B-conformation. We determine that structure, composition, and conformation specific groove dimension together influence the local hydration characteristics and, therefore, are expected to be important determinants of biological processes.

Structure, stability and specificity of the binding of ssDNA and ssRNA with proteins

Recognition of single-stranded DNA (ssDNA) or single-stranded RNA (ssRNA) is important for many fundamental cellular functions. A variety of single-stranded DNA-binding proteins (ssDBPs) and single-stranded RNA-binding proteins (ssRBPs) have evolved that bind ssDNA and ssRNA, respectively, with varying degree of affinities and specificities to form complexes. Structural studies of these complexes provide key insights into their recognition mechanism. However, computational modeling of the specific recognition process and to predict the structure of the complex is challenging, primarily due to the heterogeneity of their binding energy landscape and the greater flexibility of ssDNA or ssRNA compared with double-stranded nucleic acids. Consequently, considerably fewer computational studies have explored interactions between proteins and single-stranded nucleic acids compared with protein interactions with double-stranded nucleic acids. Here, we report a newly developed energy-based coarse-grained model to predict the structure of ssDNA–ssDBP and ssRNA–ssRBP complexes and to assess their sequence-specific interactions and stabilities. We tuned two factors that can modulate specific recognition: base–aromatic stacking strength and the flexibility of the single-stranded nucleic acid. The model was successfully applied to predict the binding conformations of 12 distinct ssDBP and ssRBP structures with their cognate ssDNA and ssRNA partners having various sequences. Estimated binding energies agreed well with the corresponding experimental binding affinities. Bound conformations from the simulation showed a funnel-shaped binding energy distribution where the native-like conformations corresponded to the energy minima. The various ssDNA–protein and ssRNA–protein complexes differed in the balance of electrostatic and aromatic energies. The lower affinity of the ssRNA–ssRBP complexes compared with the ssDNA–ssDBP complexes stems from lower flexibility of ssRNA compared to ssDNA, which results in higher rate constants for the dissociation of the complex (k_off) for complexes involving the former.

Quantifying bimolecular self-assembly is pivotal to understanding cellular function. In recent years, a large progress has been made in understanding the structure and biophysics of protein-protein interactions. Particularly, various computational tools are available for predicting these structures and to estimate their stability and the driving forces of their formation. The understating of the interactions between proteins and nucleic acids, however, is still limited, presumably due to the involvement of non-specific interactions as well as the high conformational plasticity that may demand an induced-fit mechanism. In particular, the interactions between proteins and single-stranded nucleic acids (i.e., single-stranded DNA and RNA) is very challenging due to their high flexibility. Furthermore, the interface between proteins and single-stranded nucleic acids is often chemically more heterogeneous than the interface between proteins and double-stranded DNA. In this study, we developed a coarse-grained computational model to predict the structure of complexes between proteins and single-stranded nucleic acids. The model was applied to estimate binding affinities and the estimated binding energies agreed well with the corresponding experimental binding affinities. The kinetics of association as well as the specificity of the complexes between proteins and ssDNA are different than those with ssRNA, mostly due to differences in their conformational flexibility.

Dynamical crossover of water confined within the amphiphilic nanocores of aggregated amyloid β peptides

It is believed that the self-assembly of amyloid beta (Aβ) peptides in the brain is the cause of Alzheimer's disease. Atomistic molecular dynamics simulations of aqueous solutions of Aβ protofilaments of different sizes at room temperature have been carried out to explore the dynamic properties of water confined within the core and at the exterior surface of the protofilaments. Attempts have been made to understand how the non-uniform distortion of the protofilaments associated with their structural crossover influences the diffusivity and the hydrogen bonding environment of the confined water molecules. In contrast to the homogeneous solvent dynamical environment at the exterior surface, the calculations revealed heterogeneously restricted motions of water confined within the distorted cores of the protofilaments. Importantly, it is demonstrated that the structural crossover of the aggregates observed for the decamer is associated with a dynamical transition of water confined within its core. A direct one-to-one correlation between the heterogeneously restricted core water motions and the kinetics of the breaking and formation of hydrogen bonds quantitatively demonstrated that a modified hydrogen bond arrangement within the cores of higher order Aβ protofilaments is the origin behind the crossover in core water mobility. A fraction of the water molecules forming short-lived water-water hydrogen bonds within the core of the crossover protofilament decamer are believed to diffuse away easily from the core and thus play a crucial role in further growth of the protofilament by facilitating the binding of new peptide monomers.

Dynamic cell responses in Thermoanaerobacterium sp. under hyperosmotic stress

As a nongenetic engineering technique, adaptive evolution is an effective and easy-to-operate approach to strain improvement. In this work, a commercial Thermoanaerobacterium aotearoense SCUT27/Δldh-G58 was successfully isolated via sequential batch fermentation with step-increased carbon concentrations. Mutants were isolated under selective high osmotic pressures for 58 passages. The evolved isolate rapidly catabolized sugars at high concentrations and subsequently produced ethanol with good yield. A 1.6-fold improvement of ethanol production was achieved in a medium containing 120 g/L of carbon substrate using the evolved strain, compared to the start strain. The analysis of transcriptome and intracellular solute pools suggested that the adaptive evolution altered the synthesis of some compatible solutes and activated the DNA repair system in the two Thermoanaerobacterium sp. evolved strains. Overall, the results indicated the potential of adaptive evolution as a simple and effective tool for the modification and optimization of industrial microorganisms.

Thermodynamics of complex structures formed between single-stranded DNA oligomers and the KH domains of the far upstream element binding protein

The noncovalent interaction between protein and DNA is responsible for regulating the genetic activities in living organisms. The most critical issue in this problem is to understand the underlying driving force for the formation and stability of the complex. To address this issue, we have performed atomistic molecular dynamics simulations of two DNA binding K homology (KH) domains (KH3 and KH4) of the far upstream element binding protein (FBP) complexed with two single-stranded DNA (ss-DNA) oligomers in aqueous media. Attempts have been made to calculate the individual components of the net entropy change for the complexation process by adopting suitable statistical mechanical approaches. Our calculations reveal that translational, rotational, and configurational entropy changes of the protein and the DNA components have unfavourable contributions for this protein-DNA association process and such entropy lost is compensated by the entropy gained due to the release of hydration layer water molecules. The free energy change corresponding to the association process has also been calculated using the Free Energy Perturbation (FEP) method. The free energy gain associated with the KH4-DNA complex formation has been found to be noticeably higher than that involving the formation of the KH3-DNA complex.

Characterizing Structural Stability of Amyloid Motif Fibrils Mediated by Water Molecules

In biological systems, structural confinements of amyloid fibrils can be mediated by the role of water molecules. However, the underlying effect of the dynamic behavior of water molecules on structural stabilities of amyloid fibrils is still unclear. By performing molecular dynamics simulations, we investigate the dynamic features and the effect of interior water molecules on conformations and mechanical characteristics of various amyloid fibrils. We find that a specific mechanism induced by the dynamic properties of interior water molecules can affect diffusion of water molecules inside amyloid fibrils, inducing their different structural stabilities. The conformation of amyloid fibrils induced by interior water molecules show the fibrils' different mechanical features. We elucidate the role of confined and movable interior water molecules in structural stabilities of various amyloid fibrils. Our results offer insights not only in further understanding of mechanical features of amyloids as mediated by water molecules, but also in the fine-tuning of the functional abilities of amyloid fibrils for applications.

Effects of protein-DNA complex formation on the intermolecular vibrational density of states of interfacial water

Single-stranded DNAs (ss-DNAs) are formed as intermediates during DNA metabolic processes. ss-DNA binding (SSB) proteins specifically bind to the single-stranded segments of the DNA and protect it from being degraded. We have performed room temperature molecular dynamics simulations of the aqueous solution of two DNA binding K homology (KH) domains (KH3 and KH4) of the far upstream element (FUSE) binding protein (FBP) complexed with two ss-DNA oligomers. Efforts have been made to explore the influence of complex formation on low-frequency vibrational density of states of the surface water molecules. It is revealed that increased back scattering of water confined around the complexed structures leads to significant blue shifts of the band corresponding to the O···O···O bending or restricted transverse motions of water, the effect being more for the bridged water molecules. Importantly, it is demonstrated that the formation of such complexed structures of a similar type may often influence the transverse and longitudinal degrees of freedom of the surrounding water molecules in a nonuniform manner.