Water-mediated Interactions Research Articles

Comparative Bindings of Glycosaminoglycans with CXCL8 Monomer and Dimer: Insights from Conformational Dynamics and Kinetics of Hydrogen Bonds.

GAGs bind to both the monomeric and dimeric forms of CXCL8, helping to form a concentration gradient of the chemokine that facilitates the recruitment of neutrophils to an injury site and supports other biological functions. In this study, atomistic molecular dynamics simulations were conducted to investigate the binding behavior of two hexameric GAGs sulfated at two different positions, chondroitin sulfate (CS) and heparan sulfate (HS), with the monomer (SIL8) and dimer (DIL8) forms of the CXCL8 protein. The results support that the conformational diversity of CS and HS appeared to be more when binding with monomer SIL8 than dimer DIL8. CS gained more configurational entropy from glycosidic linkage flexibility when bound to SIL8 than DIL8, with a higher energy barrier, whereas HS exhibited a lower energy barrier for configurational entropy when bound to SIL8 and DIL8. The monomer SIL8 exhibited more favorable and preferential binding with GAGs compared to DIL8. Formation of hydrogen bonds with the basic amino acids of SIL8 and GAG was more rigid and required higher activation energy to break than the other identified hydrogen bondings. Water molecules involved in hydrogen bonding with GAGs, excluding those with basic amino acids of DIL8, showed longer lifetimes and slower relaxation compared to SIL8. This suggests that water-mediated interactions also favor binding of DIL8 with GAGs. Despite having more basic amino acids, DIL8 did not display stronger binding than SIL8, indicating the significant role of basic residues in stabilizing the GAG-protein interactions in the monomers. The reason could be that the greater number of nonbasic amino acids in dimeric CXCL8 stabilizes the complex by forming water-mediated hydrogen bonds, reducing the conformational preferences for binding with GAGs. In contrast, the monomeric form of CXCL8 exhibits a higher conformational preference for protein-GAG interactions.

Mechanistic Insights of Amino Acid Binding to Hydroxyapatite: Molecular Dynamics Charts Future Directions in Biomaterial Design.

Extensive efforts have been made to improve the understanding of hard tissue regeneration, essential for advancing medical applications like bone graft materials. However, the mechanisms of bone biomineralization, particularly the regulation of hydroxyapatite growth by proteins/peptides, remain debated. Small biomolecules such as amino acids are ideal for studying these mechanisms due to their simplicity and relevance as protein/peptide building blocks. This study investigates the binding affinity of four amino acids including glycine (Gly), proline (Pro), lysine (Lys), and aspartic acid (Asp) to the hydroxyapatite (HAP) (100) surface through molecular dynamics simulations. Our findings reveal that aspartic acid exhibits the most energetically favorable binding affinity, attributed to its additional carboxylate group (-COO-), which facilitates stronger interactions with Ca2+ ions on the HAP surface compared to other amino acids with single carboxylate groups. This highlights the critical role of specific functional groups in modulating binding strength, emphasizing that the presence of multiple binding sites in amino acids enhances binding stability. Interestingly, the study also uncovers the significance of water-mediated interactions, as the compact water layer above the HAP surface acts as a barrier, complicating direct binding and underscoring the need to consider solvation effects in simulations. Glycine, due to its small size, demonstrates a unique ability to penetrate this tightly bound water monolayer, suggesting that molecular size influences binding dynamics. These simulations offer detailed insights into the atomic-level interactions, providing a deeper understanding of binding affinity and stability. These insights are pertinent for designing peptides or proteins with enhanced interactions with biomaterials, particularly in mimicking natural bone-binding processes.

Effect of Direct and Water-Mediated Interactions on the Identification of Hotspots in Biomolecular Complexes with Multiple Subsystems.

Identification of important residues in biochemical complexes is often a crucial step for many problems in molecular biology and biochemistry. A method is proposed to identify hotspots in biomolecular complexes based on a new metric, derived from networks representing molecular subunits (residues, bridging solvent molecules, ligands etc.) connected by interactions. A singular value decomposition of the weighted adjacency matrix is used to construct a scalar rank for each subunit that reflects its importance in the residue interaction network. This metric is called the singular value centrality. In addition, a new formalism is proposed to account for water-mediated interactions in the ranking of residues. Interactions for a residue network can be provided by various computational methods. In this work interactions are obtained from full quantum-mechanical calculations of protein-protein complexes using the fragment molecular orbital method. The ranking results are shown to be in good agreement with earlier computational and experimental studies. The developed method can be used to gain a deeper insight into the role of subunits in complex biomolecular systems.

Mechanistic Basis for the Translation Inhibition of Cutibacterium acnes by Clindamycin

Inflammation and the Gram-positive anaerobic bacterium Cutibacterium acnes, which is implicated in acne pathogenesis and pilosebaceous unit inflammation, are the main targets of antibiotic-based therapy against acne vulgaris (acne). The most widely used antibiotics in acne therapy are tetracyclines, macrolides, and lincosamides. Unfortunately, C. acnes bacteria over the past several decades have demonstrated increased resistance to these antibiotics, particularly to clindamycin. The precise knowledge of how antibiotics interact with their clinical target is needed to overcome this problem. Toward this goal, we determined the structure of clindamycin in complex with the ribosome of Cutibacterium acnes at 2.53 Å resolution using cryogenic electron microscopy. The galactose sugar moiety of clindamycin interacts with nucleotides of the 23S rRNA directly or through a conserved network of water-mediated interactions. Its propyl pyrrolidinyl group interacts with the 23S rRNA through van der Waals forces. Clindamycin binding to the Cutibacterium acnes ribosome interferes with both: proper orientation of the aminoacyl group of the A-site bound tRNA that is needed for peptide bond formation and with the extension of the nascent peptide. Our data are important for advancing understanding of antibiotic resistance and development of narrow- spectrum antibacterial drugs, which is an urgent need for contemporary antibiotic stewardship.

On the role of water in regulating the mechanics of collagen fibers

The mechanical properties of collagen fiber–based materials have garnered considerable attention owing to their application in various hydration environments. Herein, we conducted a comprehensive investigation on the mediatory role of water in the mechanical properties and structure of collagen fibers. Lower hydration levels are conducive to the mechanical strength, particularly the tensile and tear properties, and reduce the softness of collagen fibers as well as the D-period and the lateral spacing between collagen molecules. Atomic force microscopy measurements indicate that lower hydration levels also enhance the roughness, adhesion force and lateral forces between collagen fibers. Theoretical calculations further reveal that the hydration level regulates weaker components of water-mediated interactions, such as electrostatic and van der Waals forces, which dominate the resistance against mechanical shock. Our findings provide valuable insights into the role of water in regulating the mechanical behaviour of collagen fibers, thereby contributing to the development of novel high-performance collagen fiber materials.

Targeting DNA junction sites by bis-intercalators induces topological changes with potent antitumor effects.

Targeting inter-duplex junctions in catenated DNA with bidirectional bis-intercalators is a potential strategy for enhancing anticancer effects. In this study, we used d(CGTATACG)2, which forms a tetraplex base-pair junction that resembles the DNA-DNA contact structure, as a model target for two alkyl-linked diaminoacridine bis-intercalators, DA4 and DA5. Cross-linking of the junction site by the bis-intercalators induced substantial structural changes in the DNA, transforming it from a B-form helical end-to-end junction to an over-wounded side-by-side inter-duplex conformation with A-DNA characteristics and curvature. These structural perturbations facilitated the angled intercalation of DA4 and DA5 with propeller geometry into two adjacent duplexes. The addition of a single carbon to the DA5 linker caused a bend that aligned its chromophores with CpG sites, enabling continuous stacking and specific water-mediated interactions at the inter-duplex contacts. Furthermore, we have shown that the different topological changes induced by DA4 and DA5 lead to the inhibition of topoisomerase 2 activities, which may account for their antitumor effects. Thus, this study lays the foundations for bis-intercalators targeting biologically relevant DNA-DNA contact structures for anticancer drug development.

Interactions of adsorbing cosolutes with hydrophobic hydration shells.

The analysis of water density fluctuations in the hydration shell of nonpolar solutes provides insights into water-mediated interactions, especially hydrophobic interactions. These fluctuations are sensitive to small perturbations due to changes in thermodynamic conditions, such as temperature and pressure, but also to the presence of cosolutes, such as salts or small organic molecules. Herein, we investigate the effect of two classes of adsorbing cosolutes, using urea and methanol as representatives, on the fluctuations in energy and solvent density within the solvation shell of a model extended hydrophobic solute. We focus on the interactions of the cosolutes with the hydrophobic hydration shell, rather than with the solute itself, which though important remain largely unexplored. We calculate and analyze the interfacial partial molar energy of the cosolute, using a methodology based on the small system method. This approach provides correlated solvent density and energy fluctuations and allows us to decompose them into contributions due to interactions between the different components present in the solvation shell of the solute. The results show that adsorbed urea molecules interact more favorably with water than nonadsorbed urea molecules, which leads to the attenuation of interfacial density fluctuations and thus to the stabilization of the solvation shell. By contrast, the adsorbed methanol molecules interact preferably with other methanol molecules in the solvation shell, leading to a nano-phase segregated structure, which enhances interfacial fluctuations.

Inverse design of pore wall chemistry and topology through active learning of surface group interactions.

Design of next-generation membranes requires a nanoscopic understanding of the effect of biologically inspired heterogeneous surface chemistries and topologies (roughness) on local water and solute behavior. In particular, the rejection of small, neutral solutes, such as boric acid, poses a heretofore unsolved challenge. In prior work, a computational inverse design technique using an evolutionary optimization successfully uncovered new surface design strategies for optimized transport of water over solutes in smooth, model pores consisting of two surface chemistries. However, extending such an approach to more complex (and realistic) scenarios involving many surface chemistries as well as surface roughness is challenging due to the expanded design space. In this work, we develop a new approach that uses active learning to optimize in a reduced feature space of surface group interactions, finding parameters that lead to their assembly into ordered, optimal patterns. This approach rapidly identifies novel surface functionalizations that maximize the difference in water and boric acid transport through the nanopore. Moreover, we find that the roughness of the nanopore wall, independent of its chemistry, can be leveraged to enhance transport selectivity: oscillations in the pore wall diameter optimally inhibit boric acid transport by creating energetic wells from which the solute must escape to transport down the pore. This proof-of-concept demonstrates the potential for active learning strategies, in concert with molecular simulations, to rapidly navigate complex design spaces of aqueous interfaces and is promising as a tool for engineering water-mediated surface interactions for a broad range of applications.

Speciation and Organic Phase Structure in Nitric Acid Extraction with Trioctylamine.

Understanding chemical speciation and intermolecular interactions in multicomponent liquids is essential to understanding their phase and chemical equilibria, which underpin chemical separation processes, including solvent extraction. Here we report on the extraction of nitric acid from its aqueous solutions into organic solutions of trioctylamine (TOA) in toluene, investigated with spectroscopic, X-ray scattering, and computational tools to understand molecular speciation in the organic phase and its relationship with the nanoscale structure of the organic phase. Trends in acid and water extraction clearly show two and three regimes, respectively, indicating different stoichiometric relationships, but speciation of HNO3, water, and amine in these regimes is not apparent. 1H NMR of the organic phase shows that there are at least two distinct acidic protons in the organic phase while ATR-FTIR results show that the organic phase with excess acid extraction is a mixture of trioctylammonium-nitrate ion pairs (TOAH·NO3), and undissociated HNO3 molecules. Comparison with DFT-computed IR spectra show that the chain-like configurations of TOAH·NO3·HNO3·H2O are favored over TOAH·NO3·H2O·HNO3, i.e., direct interaction between the nitrate and HNO3 molecules is more favored compared to a water-mediated interaction. SAXS of the organic phases were modeled as sums of Ornstein-Zernike (O-Z) scattering and a prepeak feature in the higher Q region that corresponds to extractant packing. The extraction of undissociated HNO3 by the ion pairs leads to an increased X-ray scattering contrast in the organic phase without any significant change in the correlation length. These results show that the organic phase nanostructure is more sensitive to the concentration of TOAH·NO3 and is relatively unaffected by excess acid extraction. These findings will enable a molecular understanding of the mechanisms behind metal extraction from acidic media with basic extractants.

Long-range, water-mediated interaction between a moderately active antifreeze protein molecule and the surface of ice.

Using molecular dynamics simulations, we show that a molecule of moderately active antifreeze protein (type III AFP, QAE HPLC-12 isoform) is able to interact with ice in an indirect manner. This interaction occurs between the ice binding site (IBS) of the AFP III molecule and the surface of ice, and it is mediated by liquid water, which separates these surfaces. As a result, the AFP III molecule positions itself at a specific orientation and distance relative to the surface of ice, which enables the effective binding (via hydrogen bonds) of the molecule with the nascent ice surface. Our results show that the final adsorption of the AFP III molecule on the surface of ice is not achieved by chaotic diffusion movements, but it is preceded by a remote, water-mediated interaction between the IBS and the surface of ice. The key factor that determines the existence of this interaction is the ability of water molecules to spontaneously form large, high-volume aggregates that can be anchored to both the IBS of the AFP molecule and the surface of ice. The results presented in this work for AFP III are in full agreement with the ones obtained by us previously for hyperactive CfAFP, which indicates that the mechanism of the remote interaction of these molecules with ice remains unchanged despite significant differences in the molecular structure of their ice binding sites. For that reason, we can expect that also other types of AFPs interact with the ice surface according to an analogous mechanism.

α-Methylacyl-CoA Racemase from Mycobacterium tuberculosis-Detailed Kinetic and Structural Characterization of the Active Site.

α-Methylacyl-CoA racemase in M. tuberculosis (MCR) has an essential role in fatty acid metabolism and cholesterol utilization, contributing to the bacterium's survival and persistence. Understanding the enzymatic activity and structural features of MCR provides insights into its physiological and pathological significance and potential as a therapeutic target. Here, we report high-resolution crystal structures for wild-type MCR in a new crystal form (at 1.65 Å resolution) and for three active-site mutants, H126A, D156A and E241A, at 2.45, 1.64 and 1.85 Å resolutions, respectively. Our analysis of the new wild-type structure revealed a similar dimeric arrangement of MCR molecules to that previously reported and details of the catalytic site. The determination of the structures of these H126A, D156A and E241A mutants, along with their detailed kinetic analysis, has now allowed for a rigorous assessment of their catalytic properties. No significant change outside the enzymatic active site was observed in the three mutants, establishing that the diminution of catalytic activity is mainly attributable to disruption of the catalytic apparatus involving key hydrogen bonding and water-mediated interactions. The wild-type structure, together with detailed mutational and biochemical data, provide a basis for understanding the catalytic properties of this enzyme, which is important for the design of future anti-tuberculosis drug molecules.

Interplay between water-mediated and direct interactions and the organization of peptide condensates

Interplay between water-mediated and direct interactions and the organization of peptide condensates

Phase separation of multicomponent peptide mixtures into dehydrated clusters with hydrophilic cores

Phase separation of biomolecules underlies the formation and regulation of various membraneless condensates in cells. How condensates function reliably while surrounded by heterogeneous and dynamic mixtures of biomolecular components with specific and nonspecific interactions is yet to be understood. Studying multicomponent biomolecular mixtures with designer peptides has recently become an attractive avenue for learning about physicochemical principles governing cellular condensates. In this work, we employed long-timescale atomistic simulations of multicomponent tripeptide mixtures with all residue substitutions to illuminate the nature of direct and water-mediated interactions in a prototypical cellular condensate environment. We find that peptide mixtures form clusters with inverse hydrophobic order. Most multivalent and charged residues are localized in the cluster’s core, with a large fraction of nonaromatic hydrophobic residues remaining on the surface. This inverse hydrophobic order in peptide clusters is partly driven by the expulsion of nonspecifically bound water molecules following peptide cluster growth. The growth of clusters is also accompanied by the formation of increasing numbers of specific water-mediated interactions between polar and charged residues. While the present study focused on the condensation of short peptide motifs, the general findings and analysis techniques should be helpful for future studies on larger peptides and protein condensation.

New model for aspartic acid species in aqueous calcium carbonate growth environments: challenges and perspectives.

The lack of experimental data on the dynamics of aspartic acid species in water for its range of protonation states and the details of their atomic-level interaction with aqueous calcium carbonate species is a driver for accurate force field development. A classical model that is consistent with the few pieces of experimental data available and with first principles calculations has been developed. The complex dynamics of the aspartate anions relevant to biomineralization and calcium carbonate crystal growth has been explored in water, providing a quantitative description of solvation structure and free energies, including conformational free energy profiles and pairing free energies. The model has been used to probe the structure and dynamics of aqueous calcium aspartate homo- and hetero-chiral clusters, confirming their unlikelihood due to weak and water-mediated interactions. This supports the hypothesis that the formation of such clusters, observed while growing vaterite in the presence of acidic chiral amino acids, is favoured by the presence of the crystal surface.

How do water-mediated interactions and osmotic second virial coefficients vary with particle size?

We examine quantitatively the solute-size dependences of the effective interactions between nonpolar solutes in water and in a simple liquid. The potential w(r) of mean force and the osmotic second virial coefficients B are calculated with high accuracy from molecular dynamics simulations. As the solute diameter increases from methane's to C60's with the solute-solute and solute-solvent attractive interaction parameters fixed to those for the methane-methane and methane-water interactions, the first minimum of w(r) lowers from -1.1 to -4.7 in units of the thermal energy kT. Correspondingly, the magnitude of B (<0) increases proportional to σα with some power close to 6 or 7, which reinforces the solute-size dependence of B found earlier for a smaller range of σ [H. Naito, R. Okamoto, T. Sumi and K. Koga, J. Chem. Phys., 2022, 156, 221104]. We also demonstrate that the strength of the attractive interactions between solute and solvent molecules can qualitatively change the characteristics of the effective pair interaction between solute particles, both in water and in a simple liquid. If the solute-solvent attractive force is set to be weaker (stronger) than a threshold, the effective interaction becomes increasingly attractive (repulsive) with increasing solute size.

Molecular Dynamics Investigations of Human DNA-Topoisomerase I Interacting with Novel Dewar Valence Photo-Adducts: Insights into Inhibitory Activity.

Chronic exposure to ultraviolet (UV) radiation is known to induce the formation of DNA photo-adducts, including cyclobutane pyrimidine dimers (CPDs) and Dewar valence derivatives (DVs). While CPDs usually occur at higher frequency than DVs, recent studies have shown that the latter display superior selectivity and significant stability in interaction with the human DNA/topoisomerase 1 complex (TOP1). With the aim to deeply investigate the mechanism of interaction of DVs with TOP1, we report here four all-atom molecular dynamic simulations spanning one microsecond. These simulations are focused on the stability and conformational changes of two DNA/TOP1-DV complexes in solution, the data being compared with the biomimetic thymine dimer counterparts. Results from root-mean-square deviation (RMSD) and root-mean-square fluctuation (RMSF) analyses unequivocally confirmed increased stability of the DNA/TOP1-DV complexes throughout the simulation duration. Detailed interaction analyses, uncovering the presence of salt bridges, hydrogen bonds, water-mediated interactions, and hydrophobic interactions, as well as pinpointing the non-covalent interactions within the complexes, enabled the identification of specific TOP1 residues involved in the interactions over time and suggested a potential TOP1 inhibition mechanism in action.

Dynamics of water-mediated interaction effects on the stability and transmission of Omicron

SARS-Cov-2 Omicron variant and its highly transmissible sublineages amidst news of emerging hybrid variants strengthen the evidence of its ability to rapidly spread and evolve giving rise to unprecedented future waves. Owing to the presence of isolated RBD, monomeric and trimeric Cryo-EM structures of spike protein in complex with ACE2 receptor, comparative analysis of Alpha, Beta, Gamma, Delta, and Omicron assist in a rational assessment of their probability to evolve as new or hybrid variants in future. This study proposes the role of hydration forces in mediating Omicron function and dynamics based on a stronger interplay between protein and solvent with each Covid wave. Mutations of multiple hydrophobic residues into hydrophilic residues underwent concerted interactions with water leading to variations in charge distribution in Delta and Omicron during molecular dynamics simulations. Moreover, comparative analysis of interacting moieties characterized a large number of mutations lying at RBD into constrained, homologous and low-affinity groups referred to as mutational drivers inferring that the probability of future mutations relies on their function. Furthermore, the computational findings reveal a significant difference in angular distances among variants of concern due 3 amino acid insertion (EPE) in Omicron variant that not only facilitates tight domain organization but also seems requisite for characterization of mutational processes. The outcome of this work signifies the possible relation between hydration forces, their impact on conformation and binding affinities, and viral fitness that will significantly aid in understanding dynamics of drug targets for Covid-19 countermeasures. The emerging scenario is that hydration forces and hydrophobic interactions are crucial variables to probe in mutational analysis to explore conformational landscape of macromolecules and reveal the molecular origins of protein behaviors.

Graphene oxide incorporated PSU/SPEKEKK for potential application as proton exchange membranes in fuel cells

This study aims to develop a new PEM for fuel cell applications consisting of PSU/SPEKEKK/SGO membranes with improved physiochemical properties to surmount the low power generation of Nafion. The DFT calculations using NCI plots show that there are network interactions between PSU and SPEKEKK, established by the hydrogen bonding and vdW interactions. SGO was incorporated into the as-prepared PSU/SPEKEKK membrane. A solution casting method was employed to fabricate SGO-loaded composite membranes in the range of 0.1% to 1.5% (wt%). Structural evidence of the SGO and composite membranes was characterized by ATR-FTIR and FESEM, respectively. Adding SGO to PSU/SPEKEKK shows a significant increase in proton conductivity due to the presence of more protonated sites (SO3H) and water-mediated interactions for proton conduction. Of all the compositions considered, PSU/SPEKEKK/SGO-0.6-1.5% loading was found to have the highest proton conductivity, 0.03Scm−1. Although adding SGO to the composite membrane has reduced water absorption, the hydrophilic oxygenated functional groups in SGO, the bound water molecules, and the increased –SO3H groups have successfully improved the proton conductivity. The findings reveal that all composite membranes show excellent oxidative stabilities for 32 h in Fenton's reagent. To summarize, all composite membranes fabricated in this study show acceptable water uptake and swelling ratio, high oxidation, excellent thermal stability, and high proton conductivity.

Elucidating the Competitive Adsorption of H2O and CO2 in CALF-20: New Insights for Enhanced Carbon Capture Metal-Organic Frameworks.

In light of the pressing need for efficient carbon capture solutions, our study investigates the simultaneous adsorption of water (H2O) and carbon dioxide (CO2) as a function of relative humidity in CALF-20, a highly scalable and stable metal-organic framework (MOF). Advanced computer simulations reveal that due to their similar interactions with the framework, H2O and CO2 molecules compete for the same binding sites, occupying similar void regions within the CALF-20 pores. This competition results in distinct thermodynamic and dynamical behaviors of H2O and CO2 molecules, depending on whether one or both guest species are present. Notably, the presence of CO2 molecules forces the H2O molecules to form more connected hydrogen-bond networks within smaller regions, slowing water reorientation dynamics and decreasing water entropy. Conversely, the presence of water speeds up the reorientation of CO2 molecules, decreases the CO2 entropy, and increases the propensity for CO2 to be adsorbed within the framework due to stronger water-mediated interactions. Due to the competition for the same void spaces, both H2O and CO2 molecules exhibit slower diffusion when molecules of the other guest species are present. These findings offer valuable strategies and insights into enhancing the differential affinity of H2O and CO2 for MOFs specifically designed for carbon capture applications.

De novo prediction of explicit water molecule positions by a novel algorithm within the protein design software MUMBO

By mediating interatomic interactions, water molecules play a major role in protein–protein, protein–DNA and protein–ligand interfaces, significantly affecting affinity and specificity. This notwithstanding, explicit water molecules are usually not considered in protein design software because of high computational costs. To challenge this situation, we analyzed the binding characteristics of 60,000 waters from high resolution crystal structures and used the observed parameters to implement the prediction of water molecules in the protein design and side chain-packing software MUMBO. To reduce the complexity of the problem, we incorporated water molecules through the solvation of rotamer pairs instead of relying on solvated rotamer libraries. Our validation demonstrates the potential of our algorithm by achieving recovery rates of 67% for bridging water molecules and up to 86% for fully coordinated waters. The efficacy of our algorithm is highlighted further by the prediction of 3 different proteinligand complexes. Here, 91% of water-mediated interactions between protein and ligand are correctly predicted. These results suggest that the new algorithm could prove highly beneficial for structure-based protein design, particularly for the optimization of ligand-binding pockets or protein–protein interfaces.