A New Approach for Investigating the Molecular Recognition of Protein: Toward Structure-Based Drug Design Based on the 3D-RISM Theory

RISMiCal: A software package to perform fast RISM/3D-RISM calculations

Solvent plays an essential role in a variety of chemical, physical, and biological processes that occur in the solution phase. The reference interaction site model (RISM) and its three-dimensional extension (3D-RISM) serve as powerful computational tools for modeling solvation effects in chemical reactions, biological functions, and structure formations. We present the RISM integrated calculator (RISMiCal) program package, which is based on RISM and 3D-RISM theories with fast GPU code. RISMiCal has been developed as an integrated RISM/3D-RISM program that has interfaces with external programs such as Gaussian16, GAMESS, and Tinker. Fast 3D-RISM programs for single- and multi-GPU codes written in CUDA would enhance the availability of these hybrid methods because they require the performance of many computationally expensive 3D-RISM calculations. We expect that our package can be widely applied for chemical and biological processes in solvent. The RISMiCal package is available at https://rismical-dev.github.io.

Free Energy Density of a Fluid and Its Role in Solvation and Binding

The concept that a fluid has a position-dependent free energy density appears in the literature but has not been fully developed or accepted. We set this concept on an unambiguous theoretical footing via the following strategy. First, we set forth four desiderata that should be satisfied by any definition of the position-dependent free energy density, f(R), in a system comprising only a fluid and a rigid solute: its volume integral, plus the fixed internal energy of the solute, should be the system free energy; it deviates from its bulk value, fbulk, near a solute but should asymptotically approach fbulk with increasing distance from the solute; it should go to zero where the solvent density goes to zero; and it should be well-defined in the most general case of a fluid made up of flexible molecules with an arbitrary interaction potential. Second, we use statistical thermodynamics to formulate a definition of the free energy density that satisfies these desiderata. Third, we show how any free energy density satisfying the desiderata may be used to analyze molecular processes in solution. In particular, because the spatial integral of f(R) equals the free energy of the system, it can be used to compute free energy changes that result from the rearrangement of solutes as well as the forces exerted on the solutes by the solvent. This enables the use of a thermodynamic analysis of water in protein binding sites to inform ligand design. Finally, we discuss related literature and address published concerns regarding the thermodynamic plausibility of a position-dependent free energy density. The theory presented here has applications in theoretical and computational chemistry and may be further generalizable beyond fluids, such as to solids and macromolecules.

RNA-ligand molecular docking: advances and challenges

With rapid advances in computer algorithms and hardware, fast and accurate virtual screening has led to a drastic acceleration in selecting potent small molecules as drug candidates. Computational modeling of RNA-small molecule interactions has become an indispensable tool for RNA-targeted drug discovery. The current models for RNA-ligand binding have mainly focused on the docking-and-scoring method. Accurate docking and scoring should tackle four crucial problems: (1) conformational flexibility of ligand, (2) conformational flexibility of RNA, (3) efficient sampling of binding sites and binding poses, and (4) accurate scoring of different binding modes. Moreover, compared with the problem of protein-ligand docking, predicting ligand binding to RNA, a negatively charged polymer, is further complicated by additional effects such as metal ion effects. Thermodynamic models based on physics-based and knowledge-based scoring functions have shown highly encouraging success in predicting ligand binding poses and binding affinities. Recently, kinetic models for ligand binding have further suggested that including dissociation kinetics (residence time) in ligand docking would result in improved performance in estimating in vivo drug efficacy. More recently, the rise of deep-learning approaches has led to new tools for predicting RNA-small molecule binding. In this review, we present an overview of the recently developed computational methods for RNA-ligand docking and their advantages and disadvantages.

Biomolecular Simulations with the Three-Dimensional Reference Interaction Site Model with the Kovalenko-Hirata Closure Molecular Solvation Theory

The statistical mechanics-based 3-dimensional reference interaction site model with the Kovalenko-Hirata closure (3D-RISM-KH) molecular solvation theory has proven to be an essential part of a multiscale modeling framework, covering a vast region of molecular simulation techniques. The successful application ranges from the small molecule solvation energy to the bulk phase behavior of polymers, macromolecules, etc. The 3D-RISM-KH successfully predicts and explains the molecular mechanisms of self-assembly and aggregation of proteins and peptides related to neurodegeneration, protein-ligand binding, and structure-function related solvation properties. Upon coupling the 3D-RISM-KH theory with a novel multiple time-step molecular dynamic (MD) of the solute biomolecule stabilized by the optimized isokinetic Nosé-Hoover chain thermostat driven by effective solvation forces obtained from 3D-RISM-KH and extrapolated forward by generalized solvation force extrapolation (GSFE), gigantic outer time-steps up to picoseconds to accurately calculate equilibrium properties were obtained in this new quasidynamics protocol. The multiscale OIN/GSFE/3D-RISM-KH algorithm was implemented in the Amber package and well documented for fully flexible model of alanine dipeptide, miniprotein 1L2Y, and protein G in aqueous solution, with a solvent sampling rate ~150 times faster than a standard MD simulation in explicit water. Further acceleration in computation can be achieved by modifying the extent of solvation layers considered in the calculation, as well as by modifying existing closure relations. This enhanced simulation technique has proven applications in protein-ligand binding energy calculations, ligand/solvent binding site prediction, molecular solvation energy calculations, etc. Applications of the RISM-KH theory in molecular simulation are discussed in this work.

A computational method to simulate global conformational changes of proteins induced by cosolvent

A computational method to investigate the global conformational change of a protein is proposed by combining the linear response path following (LRPF) method and three-dimensional reference interaction site model (3D-RISM) theory, which is referred to as the LRPF/3D-RISM method. The proposed method makes it possible to efficiently simulate protein conformational changes caused by either solutions of varying concentrations or the presence of cosolvent species by taking advantage of the LRPF and 3D-RISM. The proposed method is applied to the urea-induced denaturation of ubiquitin. The LRPF/3D-RISM trajectories successfully simulate the early stage of the denaturation process within the simulation time of 300 ns, whereas no significant structural change is observed even in the 1 μs standard MD simulation. The obtained LRPF/3D-RISM trajectories reproduce the mechanism of the urea denaturation of ubiquitin reported in previous studies, and demonstrate the high efficiency of the method.

Molecular Recognition and Self-Organization in Life Phenomena Studied by a Statistical Mechanics of Molecular Liquids, the RISM/3D-RISM Theory

There are two molecular processes that are essential for living bodies to maintain their life: the molecular recognition, and the self-organization or self-assembly. Binding of a substrate by an enzyme is an example of the molecular recognition, while the protein folding is a good example of the self-organization process. The two processes are further governed by the other two physicochemical processes: solvation and the structural fluctuation. In the present article, the studies concerning the two molecular processes carried out by Hirata and his coworkers, based on the statistical mechanics of molecular liquids or the RISM/3D-RISM theory, are reviewed.

Modeling the interaction of SARS-CoV-2 binding to the ACE2 receptor via molecular theory of solvation

The angiotensin-converting enzyme 2 (ACE2) protein is a cell gate receptor for the SARS-CoV-2 virus, responsible for the development of symptoms associated with the Covid-19 disease.

First-principles modeling of chemistry in mixed solvents: Where to go from here?

Mixed solvents (i.e., binary or higher order mixtures of ionic or nonionic liquids) play crucial roles in chemical syntheses, separations, and electrochemical devices because they can be tuned for specific reactions and applications. Apart from fully explicit solvation treatments that can be difficult to parameterize or computationally expensive, there is currently no well-established first-principles regimen for reliably modeling atomic-scale chemistry in mixed solvent environments. We offer our perspective on how this process could be achieved in the near future as mixed solvent systems become more explored using theoretical and computational chemistry. We first outline what makes mixed solvent systems far more complex compared to single-component solvents. An overview of current and promising techniques for modeling mixed solvent environments is provided. We focus on so-called hybrid solvation treatments such as the conductor-like screening model for real solvents and the reference interaction site model, which are far less computationally demanding than explicit simulations. We also propose that cluster-continuum approaches rooted in physically rigorous quasi-chemical theory provide a robust, yet practical, route for studying chemical processes in mixed solvents.

New Protocol for Predicting the Ligand-Binding Site and Mode Based on the 3D-RISM/KH Theory

An efficient algorithm to find the binding position and mode of small ligands bound at an active site of protein is proposed based on the spatial distribution function (SDF) obtained from the three-dimensional reference interaction site model (3D-RISM) theory with the Kovalenko-Hirata (KH) closure relation. The ligand examined includes hydrophobic, acidic, and basic molecules and zwitterions. Eighteen different types of proteins, which serve as targets for those ligands, are selected to examine the robustness of the algorithm. An imaginary atom, referred to as an "anchor site", is defined at the center of geometry of a ligand molecule that serves as a center for searching the binding position and mode of the ligand molecule in the translational and rotational spaces. The probable binding sites (PBSs) are identified based on the SDFs of the ligand molecules around the protein, and the PBS is ranked according to the peak height of SDF. The deviations from the mean height of the peak values of SDFs for 50 PBSs are analyzed based on the z-score, which is a measure of prominence of the site. The PBS found at the closest distance from the anchor site of the crystal structure is referred to as the "nearest site". The orientation of the ligand molecule at each PBS is explored by changing the Euler angles, and the most probable binding mode is determined based on the superposition approximation. The binding position of ligand molecules is successfully predicted as one of the distinct peaks in SDF of the anchor site, with a few exceptions. The binding mode of the ligand molecule predicted based on the superposition approximation is consistent with the X-ray crystal structure in nine systems, a half of the systems investigated. The significance of the results is discussed in detail. An application of the new protocol to fragment-based drug discovery is suggested.

Role of Mg2+ Ions in DNA Hydrolysis by EcoRV, Studied by the 3D-Reference Interaction Site Model and Molecular Dynamics

The role of Mg²⁺ ions during precursor formation in DNA hydrolysis by the homodimeric restriction enzyme EcoRV was elucidated based on the 3D-reference interaction site model (RISM) theory and the molecular dynamics (MD) simulation. From an analysis of the spatial distribution of Mg²⁺ in an active site using 3D-RISM, we identified a new position for Mg²⁺ in the X-ray EcoRV-DNA complex structure ( 1rvb ). We refer to the position as site IV^†. Site IV^† is almost at the same position as that of a Ca²⁺ ion in the superimposed X-ray crystallographic active-site structure of the PvuII-DNA complex ( 1f0o ). 3D-RISM was also used to locate the position of water molecules, including the water nucleophile at the active site. MD simulations were carried out with the initial structure having two Mg²⁺ ions at site IV^† and at site I*, experimentally identified by Horton et al., to find a stable complex structure in which the DNA fragment was rearranged to orient the scissile bond direction toward the water nucleophile. The equilibrium active-site structure of the EcoRV-DNA complex obtained from the MD simulation was similar to the superimposed X-ray crystallographic structure of the BamHI-DNA complex ( 2bam ). In the active-site structure, two metal ions have almost the same position (≤1.0 Å) as that of 2bam , and the scissile phosphate is twisted to orient the scissile bond toward the water nucleophile, as is the case in 2bam . We propose the equilibrium active-site structure obtained in this study as a precursor for the hydrolysis reaction of EcoRV.

How Water's Properties Are Encoded in Its Molecular Structure and Energies

How are water's material properties encoded within the structure of the water molecule? This is pertinent to understanding Earth's living systems, its materials, its geochemistry and geophysics, and a broad spectrum of its industrial chemistry. Water has distinctive liquid and solid properties: It is highly cohesive. It has volumetric anomalies-water's solid (ice) floats on its liquid; pressure can melt the solid rather than freezing the liquid; heating can shrink the liquid. It has more solid phases than other materials. Its supercooled liquid has divergent thermodynamic response functions. Its glassy state is neither fragile nor strong. Its component ions-hydroxide and protons-diffuse much faster than other ions. Aqueous solvation of ions or oils entails large entropies and heat capacities. We review how these properties are encoded within water's molecular structure and energies, as understood from theories, simulations, and experiments. Like simpler liquids, water molecules are nearly spherical and interact with each other through van der Waals forces. Unlike simpler liquids, water's orientation-dependent hydrogen bonding leads to open tetrahedral cage-like structuring that contributes to its remarkable volumetric and thermal properties.

Dynamics theory for molecular liquids based on an interaction site model

Dynamics theories for molecular liquids based on an interaction site model have been developed over the past few decades and proved to be powerful tools to investigate various dynamical phenomena.

Development of molecular closures for the reference interaction site model theory with application to square-well and Lennard-Jones homonuclear diatomics

Inspired by significant improvements obtained for the performances of the polymer reference interaction site model (PRISM) theory of the fluid phase when coupled with 'molecular closures' (Schweizer and Yethiraj 1993 J. Chem. Phys. 98 9053), we exploit a matrix generalization of this concept, suitable for the more general RISM framework. We report a preliminary test of the formalism, as applied to prototype square-well homonuclear diatomics. As for the structure, comparison with Monte Carlo shows that molecular closures are slightly more predictive than their 'atomic' counterparts, and thermodynamic properties are equally accurate. We also devise an application of molecular closures to models interacting via continuous, soft-core potentials, by using well established prescriptions in liquid state perturbation theories. In the case of Lennard-Jones dimers, our scheme definitely improves over the atomic one, providing semi-quantitative structural results, and quite good estimates of internal energy, pressure and phase coexistence. Our finding paves the way to a systematic employment of molecular closures within the RISM framework to be applied to more complex systems, such as molecules constituted by several non-equivalent interaction sites.

Designing molecular complexes using free-energy derivatives from liquid-state integral equation theory

Complex formation between molecules in solution is the key process by which molecular interactions are translated into functional systems. These processes are governed by the binding or free energy of association which depends on both direct molecular interactions and the solvation contribution. A design goal frequently addressed in pharmaceutical sciences is the optimization of chemical properties of the complex partners in the sense of minimizing their binding free energy with respect to a change in chemical structure. Here, we demonstrate that liquid-state theory in the form of the solute-solute equation of the reference interaction site model provides all necessary information for such a task with high efficiency. In particular, computing derivatives of the potential of mean force (PMF), which defines the free-energy surface of complex formation, with respect to potential parameters can be viewed as a means to define a direction in chemical space toward better binders. We illustrate the methodology in the benchmark case of alkali ion binding to the crown ether 18-crown-6 in aqueous solution. In order to examine the validity of the underlying solute-solute theory, we first compare PMFs computed by different approaches, including explicit free-energy molecular dynamics simulations as a reference. Predictions of an optimally binding ion radius based on free-energy derivatives are then shown to yield consistent results for different ion parameter sets and to compare well with earlier, orders-of-magnitude more costly explicit simulation results. This proof-of-principle study, therefore, demonstrates the potential of liquid-state theory for molecular design problems.

Massively parallel implementation of 3D-RISM calculation with volumetric 3D-FFT

A new three-dimensional reference interaction site model (3D-RISM) program for massively parallel machines combined with the volumetric 3D fast Fourier transform (3D-FFT) was developed, and tested on the RIKEN K supercomputer. The ordinary parallel 3D-RISM program has a limitation on the number of parallelizations because of the limitations of the slab-type 3D-FFT. The volumetric 3D-FFT relieves this limitation drastically. We tested the 3D-RISM calculation on the large and fine calculation cell (2048(3) grid points) on 16,384 nodes, each having eight CPU cores. The new 3D-RISM program achieved excellent scalability to the parallelization, running on the RIKEN K supercomputer. As a benchmark application, we employed the program, combined with molecular dynamics simulation, to analyze the oligomerization process of chymotrypsin Inhibitor 2 mutant. The results demonstrate that the massive parallel 3D-RISM program is effective to analyze the hydration properties of the large biomolecular systems.

Analysis of biomolecular solvation sites by 3D-RISM theory

We derive, implement, and apply equilibrium solvation site analysis for biomolecules. Our method utilizes 3D-RISM calculations to quickly obtain equilibrium solvent distributions without either necessity of simulation or limits of solvent sampling. Our analysis of these distributions extracts highest likelihood poses of solvent as well as localized entropies, enthalpies, and solvation free energies. We demonstrate our method on a structure of HIV-1 protease where excellent structural and thermodynamic data are available for comparison. Our results, obtained within minutes, show systematic agreement with available experimental data. Further, our results are in good agreement with established simulation-based solvent analysis methods. This method can be used not only for visual analysis of active site solvation but also for virtual screening methods and experimental refinement.

Biomolecular electrostatics and solvation: a computational perspective

An understanding of molecular interactions is essential for insight into biological systems at the molecular scale. Among the various components of molecular interactions, electrostatics are of special importance because of their long-range nature and their influence on polar or charged molecules, including water, aqueous ions, proteins, nucleic acids, carbohydrates, and membrane lipids. In particular, robust models of electrostatic interactions are essential for understanding the solvation properties of biomolecules and the effects of solvation upon biomolecular folding, binding, enzyme catalysis, and dynamics. Electrostatics, therefore, are of central importance to understanding biomolecular structure and modeling interactions within and among biological molecules. This review discusses the solvation of biomolecules with a computational biophysics view toward describing the phenomenon. While our main focus lies on the computational aspect of the models, we provide an overview of the basic elements of biomolecular solvation (e.g. solvent structure, polarization, ion binding, and non-polar behavior) in order to provide a background to understand the different types of solvation models.

3D-RISM-Dock: A New Fragment-Based Drug Design Protocol

We explore a new approach in the rational design of specificity in molecular recognition of small molecules based on statistical-mechanical integral equation theory of molecular liquids in the form of the three-dimensional reference interaction site model with the Kovalenko-Hirata closure (3D-RISM-KH). The numerically stable iterative solution of conventional 3D-RISM equations includes the fragmental decomposition of flexible ligands, which are treated as distinct species in solvent mixtures of arbitrary complexity. The computed density functions for solution (including ligand) molecules are obtained as a set of discrete spatial grids that uniquely describe the continuous solvent-site distribution around the protein solute. Potentials of mean force derived from these distributions define the scoring function interfaced with the AutoDock program for an automated ranking of docked conformations. As a case study in terms of solvent composition, we analyze cooperative interactions encountered in the binding of a flexible thiamine molecule to the prion protein at near-physiological conditions. The predicted location and residency times of computed binding modes are in excellent agreement with the available experimental data.

Placevent: an algorithm for prediction of explicit solvent atom distribution-application to HIV-1 protease and F-ATP synthase

We have created a simple algorithm for automatically predicting the explicit solvent atom distribution of biomolecules. The explicit distribution is coerced from the three-dimensional (3D) continuous distribution resulting from a 3D reference interaction site model (3D-RISM) calculation. This procedure predicts optimal location of solvent molecules and ions given a rigid biomolecular structure and the solvent composition. We show examples of predicting water molecules near the KNI-272 bound form of HIV-1 protease and predicting both sodium ions and water molecules near the rotor ring of F-adenosine triphosphate (ATP) synthase. Our results give excellent agreement with experimental structure with an average prediction error of 0.39-0.65 Å. Further, unlike experimental methods, this method does not suffer from the partial occupancy limit. Our method can be performed directly on 3D-RISM output within minutes. It is extremely useful for examining multiple specific solvent-solute interactions, as a convenient method for generating initial solvent structures for molecular dynamics calculations, and may assist in refinement of experimental structures. © 2012 Wiley Periodicals, Inc.