Heterogeneous and Allosteric Role of Surface Hydration for Protein-Ligand Binding

Investigating the Bromoform Membrane Interactions Using Atomistic Simulations and Machine Learning: Implications for Climate Change Mitigation

Methane emissions from livestock contribute to global warming. Seaweeds used as food additive offer a promising emission mitigation strategy because seaweeds are enriched in bromoform─a methanogenesis inhibitor. Therefore, understanding bromoform storage and production in seaweeds and particularly in a cell-like environment is crucial. As a first step toward this aim, we present an atomistic description of bromoform dynamics, diffusion, and aggregation in the presence of lipid membranes. Using all-atom molecular dynamics simulations with customized CHARMM-formatted bromoform force field files, we investigate the interactions of bromoform and lipid bilayer across various concentrations. Bromoform penetrates membranes and at high concentrations forms aggregates outside the membrane without affecting membrane thickness or lipid tail order. Aggregates outside the membrane influence the membrane curvature. Within the membrane, bromoform preferentially localizes in the membrane hydrophobic core and diffuses the slowest along the membrane normal. Employing general local-atomic descriptors and unsupervised machine learning, we demonstrate the similarity of bromoform local structures between the liquid and aggregated forms.

Solvation free energy in governing equations for DNA hybridization, protein-ligand binding, and protein folding

This work examines the thermodynamics of model biomolecular interactions using a governing equation that accounts for the participation of bulk water in the equilibria. In the first example, the binding affinities of two DNA duplexes, one of nine and one of 10 base pairs in length, are measured and characterized by isothermal titration calorimetry (ITC) as a function of concentration. The results indicate that the change in solvation free energy that accompanies duplex formation (ΔGS) is large and unfavorable. The duplex with the larger number of G:C pairings yields the largest change in solvation free energy, ΔGS = +460 kcal·mol-1per base pair at 25 °C. A van't Hoff analysis of the data is complicated by the varying degree of intramolecular base stacking within each DNA strand as a function of temperature. A modeling study reveals how the solvation free energy alters the output of a typical ITC experiment and leads to a good, though misleading, fit to the classical equilibrium equation. The same thermodynamic framework is applied to a model protein-ligand interaction, the binding of ribonuclease A with the nucleotide inhibitor 3'-UMP, and to a conformational equilibrium, the change in tertiary structure of α-lactalbumin in molar guanidinium chloride solutions. The ribonuclease study yields a value of ΔGS = +160 kcal·mol-1, whereas the folding equilibrium yields ΔGS ≈ 0, an apparent characteristic of hydrophobic interactions. These examples provide insight on the role of solvation energy in binding equilibria and suggest a pivot in the fundamental application of thermodynamics to solution chemistry.

CHARMM at 45: Enhancements in Accessibility, Functionality, and Speed

Since its inception nearly a half century ago, CHARMM has been playing a central role in computational biochemistry and biophysics. Commensurate with the developments in experimental research and advances in computer hardware, the range of methods and applicability of CHARMM have also grown. This review summarizes major developments that occurred after 2009 when the last review of CHARMM was published. They include the following: new faster simulation engines, accessible user interfaces for convenient workflows, and a vast array of simulation and analysis methods that encompass quantum mechanical, atomistic, and coarse-grained levels, as well as extensive coverage of force fields. In addition to providing the current snapshot of the CHARMM development, this review may serve as a starting point for exploring relevant theories and computational methods for tackling contemporary and emerging problems in biomolecular systems. CHARMM is freely available for academic and nonprofit research at https://academiccharmm.org/program.

Thermodynamic consequences of stapling side-chains on a peptide ligand using a lactam-bridge: A theoretical study on anti-angiogenic peptides targeting VEGF

Peptides are promising therapeutic agents for various biological targets due to their high efficacy and low toxicity, and the design of peptide ligands with high binding affinity to the target of interest is of utmost importance in peptide-based drug design. Introducing a conformational constraint to a flexible peptide ligand using a side-chain lactam-bridge is a convenient and efficient method to improve its binding affinity to the target. However, in general, such a small structural modification to a flexible ligand made with the intent of lowering the configurational entropic penalty for binding may have unintended consequences in different components of the binding enthalpy and entropy, including the configurational entropy component, which are still not clearly understood. Toward probing this, we examine different components of the binding enthalpy and entropy as well as the underlying structure and dynamics, for a side-chain lactam-bridged peptide inhibitor and its flexible analog forming complexes with vascular endothelial growth factor (VEGF), using all-atom molecular dynamics simulations. It is found that introducing a side-chain lactam-bridge constraint into the flexible peptide analog led to a gain in configurational entropy change but losses in solvation entropy, solute internal energy, and solvation energy changes upon binding, pinpointing the opportunities and challenges in drug design. The present study features an interplay between configurational and solvation entropy changes, as well as the one between binding enthalpy and entropy, in ligand-target binding upon imposing a conformational constraint into a flexible ligand.

Physicochemical properties and density functional theory calculation of octahedral UiO-66 with Bis(Trifluoromethanesulfonyl)imide ionic liquids

In this study, the physicochemical properties and molecular interactions between zirconium-based metal-organic framework (UiO-66) and three different ionic liquids based on bis(trifluoromethanesulfonyl)imide anion (EMIM+, BMIM+ and OMIM+) was performed via a combined experimental and computational approach. The ionic liquid loaded UiO-66 or IL@UiO-66 was synthesized and characterized to understand the host-guest interaction. Density functional theory calculation was performed to analyse the electronic structure of IL@UiO-66 to provide molecular insight on the dominant interactions occurred in the hybrid material. Results showed that all ILs were successfully incorporated into the micropores of UiO-66. The 3D framework was retained even after loaded with ILs as analyzed from XRD pattern. FTIR spectrum reveals that interactions of ILs with UiO-66 influenced by the alkyl chain length of the cation. The anion has a profound affinity with the UiO-66 due to the presence of electronegative atoms. Phase transition study from DSC suggested that the incorporation of ILs has stabilized the framework of UiO-66 by shifting the endothermic peak to a higher state. These findings were further elaborated with DFT calculation. Geometrical optimizations confirmed the structural parameter changes of UiO-66 when loaded with ILs. These was mainly contributed by the non-covalent interactions which was confirmed by the reduced density gradient scattered plot. Another important findings are the strength of hydrogen bonding at the host-guest interface was influenced by the alkyl chain length. The molecular orbital analysis also shows that the size of alkyl chain influence the reactivity of the hybrid material. The present study provides fundamental insights on the molecular interaction of UiO-66 and ILs as a hybrid material, which can open new possibilities for advanced material for metal-organic framework applications in energy storage system, catalysis, gas storage and medicinal chemistry.

DOX_BDW: Incorporating Solvation and Desolvation Effects of Cavity Water into Nonfitting Protein-Ligand Binding Affinity Prediction

Accurate prediction of the protein-ligand binding affinity (PLBA) with an affordable cost is one of the ultimate goals in the field of structure-based drug design (SBDD), as well as a great challenge in the computational and theoretical chemistry. Herein, we have systematically addressed the complicated solvation and desolvation effects on the PLBA brought by the difference of the explicit water in the protein cavity before and after ligands bind to the protein-binding site. Based on the new solvation model, a nonfitting method at the first-principles level for the PLBA prediction was developed by taking the bridging and displaced water (BDW) molecules into account simultaneously. The newly developed method, DOX_BDW, was validated against a total of 765 noncovalent and covalent protein-ligand binding pairs, including the CASF2016 core set, Cov_2022 covalent binding testing set, and six testing sets for the hit and lead compound optimization (HLO) simulation. In all of the testing sets, the DOX_BDW method was able to produce PLBA predictions that were strongly correlated with the corresponding experimental data (R = 0.66-0.85). The overall performance of DOX_BDW is better than the current empirical scoring functions that are heavily parameterized. DOX_BDW is particularly outstanding for the covalent binding situation, implying the need for considering an electronic structure in covalent drug design. Furthermore, the method is especially recommended to be used in the HLO scenario of SBDD, where hundreds of similar derivatives need to be screened and refined. The computational cost of DOX_BDW is affordable, and its accuracy is remarkable.