Dipankar Roy

3D-RISM-KH molecular solvation theory based on statistical mechanics has been an engine of the multiscale methods framework, which also includes molecular simulation techniques. Its applications range from the solvation energy of small molecules to the phase behavior of polymers and biomolecules. Molecular solvation theory predicts and explains the molecular mechanisms and functioning of a variety of chemical and biomolecular systems. This includes the self-assembly and conformational stability of synthetic organic rosette nanotubes (RNTs), the aggregation of peptides and proteins related to neurodegeneration, the binding of ligands to proteins, and the solvation properties of biomolecules related to their functions. The replica RISM-KH-VM molecular solvation theory predicts and explains the structure, thermodynamics, and electrochemistry of electrolyte solutions sorbed in nanoporous carbon supercapacitor electrodes, and is part of recent research and development efforts. A new quasidynamics protocol couples multiple time step molecular dynamics (MTS-MD) stabilized with an optimized isokinetic Nosé–Hoover (OIN) thermostat driven by 3D-RISM-KH mean solvation forces at gigantic outer time steps of picoseconds, which are extrapolated forward at short inner time steps of femtoseconds with generalized solvation force extrapolation (GSFE). The OIN/3D-RISM-KH/GSFE quasidynamics is implemented in the Amber Molecular Dynamics package. It is validated on miniprotein 1L2Y and protein G in ambient aqueous solution, and shows the rate of sampling 150 times faster than in standard MD simulations on these biomolecules in explicit water. The self-consistent field version of Kohn–Sham DFT in 3D-RISM-KH mean solvation forces is implemented in the Amsterdam Density Functional (ADF) package. Its applications range from solvation thermochemistry, conformational equilibria, and photochemistry to activation barriers of different nanosystems in solutions and ionic liquids.

The partition coefficients of drug and drug-like molecules between an aqueous and organic phase are an important property for developing new therapeutics. The predictive power of computational methods is used extensively to predict partition coefficients of molecules. The application of quantum chemical calculations is used to develop methods to develop structure–activity relationship models for such prediction, either based on molecular fragment methods, or via direct calculation of solvation free energy in solvent continuum. The applicability, merits, and shortcomings of these developments are revisited here.

The 3D-reference interaction site model (3D-RISM) molecular solvation theory in combination with the Kovalenko–Hirata (KH) closure is extended to seven heterocyclic liquids to understand their liquid states and to test the performance of the theory in solvation free energy (SFE) calculations of solutes in select solvents. The computed solvent site distribution profiles were compared with the all-atom molecular dynamics (MD) simulations, showing comparable performances. The computational results were compared against the structural parameters for liquids, whenever available, as well as against the experimental SFEs. The liquids are found to have local ordered structures held together via weak interactions in both the RISM and MD simulations. The 3D-RISM-KH computed SFEs are in good agreement with the benchmark values for the tetrahydrothiophene-S,S-dioxide, and showed comparatively larger deviations in the case of the SFEs in the tetrahydrofuran continuum.

In this work, we have examined the efficiency of continuum solvation models, used with density functional theory methods, in calculating octanol-water partition coefficients (logP) of 56 fluorine containing drug molecules. The solvation model based on density model computed logP values that are in good agreement with the benchmark values. The conductor-like polarizable continuum models computed results have issues in predicting correct trend, often with reversal of sign from benchmark. The choice of basis set does not show significant effect, and the selection of atomic radii affects geometry convergence during calculations.

The molecular mechanism of chloroquine resistance by the chloroquine resistance transporter protein of Plasmodium sp. is explored using molecular modeling and computational methods. The key mutation, lysine(K)-76 to threonine(T) (LYS76THR) in the transporter protein pertains to increased recognition of the protonated forms of the antimalarial drug. Such enhanced affinity can promote drug efflux from host digestive vacuole, rendering aminoquinoline-based treatment ineffective.

The three-dimensional reference interaction site model of the molecular solvation theory with the Kovalenko–Hirata closure is used to calculate the free energy of solvation of organic solutes in liquid aliphatic ketones. The ketone solvent sites were modeled using a modified united-atom force field. The successful application of these solvation models in calculating ketone–water partition coefficients of a large number of solutes supports the validation and benchmarking reported here.

The Three-Dimensional Reference Interaction Site Model (3D-RISM) with Kovalenko−Hirata (KH) closure is applied to calculate the 1,9-Decadiene/Water partition coefficients for a diverse class of compounds. The liquid state of 1,9-Decadiene is represented with the united atom TraPPE force field parameters. The 3D-RISM-KH computed partition functions are in good agreement with the experimental results. Our computational scheme can be used for a quantitative structure partitioning prediction for decadiene-water system, which has been used in membrane-mimicking of the egg-lecithin/water permeability experiments.

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.

The parallel artificial membrane permeability assay (PAMPA), a non-cellular lab-based assay, is extensively used to measure the permeability of pharmaceutical compounds. PAMPA experiments provide a working mimic of a molecule passing through cells and PAMPA values are widely used to estimate drug absorption parameters. There is an increased interest in developing computational methods to predict PAMPA permeability values. We developed an in silico model to predict the permeability of compounds based on the PAMPA assay. We used the three-dimensional reference interaction site model (3D-RISM) theory with the Kovalenko–Hirata (KH) closure to calculate the excess chemical potentials of a large set of compounds and predicted their apparent permeability with good accuracy (mean absolute error or MAE = 0.69 units) when compared to a published experimental data set. Furthermore, our in silico PAMPA protocol performed very well in the binary prediction of 288 compounds as being permeable or impermeable (precision = 94%, accuracy = 93%). This suggests that our in silico protocol can mimic the PAMPA assay and could aid in the rapid discovery or screening of potentially therapeutic drug leads that can be delivered to a desired tissue.