Intrinsically smooth discretisation of Connolly's solvent-excluded molecular surface

Investigation and drug design for novel molecules from natural products as inhibitors for controlling multiple myeloma disease using in-silico tools

Multiple myeloma (MM) is a disease that causes plasma cell growth in the bone marrow and immune globulin buildup in blood and urine. Despite recent advances in MM therapy, many still die due to its high mortality rate. A study using computational simulations analyzed 100 natural ingredients from the SANC database to determine if they inhibited the IgH domain, a known cause of multiple myeloma. Natural component Diospyrin inhibited the IgH enzyme with the best binding energy of -10.3 kcal/mol and three carbon-hydrogen bonds, followed by Parviflorone F complex with a binding energy of -10.1 kcal/mol and two conventional-hydrogen bonds. As a result, the Molecular Dynamic simulation was used to test the stability of the two complexes. During the simulation, the Diospyrin molecule dissociated from the protein at roughly 67.5 ns, whereas the Parviflorone F molecule stayed attached to the protein throughout. The latter was the subject of the investigation. The analysis of the production run data revealed that the Parviflorone F molecule exhibits a variety of conformations within the binding pocket while keeping a relatively constant distance from the protein's center of mass. The analysis of the production run data revealed that the Parviflorone F molecule exhibited a variety of conformations within the binding pocket while keeping a relatively constant distance from the protein's center of mass. The root mean square deviation (RMSD) plots for both the protein and complex showed a stable and steady average value of 4.4 Å for the first 82 nanoseconds of manufacture. As a result, the average value increased to 8.3 Å. Furthermore, the components of the binding free energy, as computed by MM-GBSA, revealed that the mean binding energy of the Parviflorone F molecule was -23.88 kcal/mol. Finally, after analyzing all of the examination data, Parviflorone F was identified as a powerful inhibitor of the IgH domain and hence of the MM disease, which requires further in-vivo conformation.Communicated by Ramaswamy H. Sarma.

Characteristic Binding Landscape of Estrogen Receptor-α36 Protein Enhances Promising Cancer Drug Design

Breast cancer (BC) remains the most common cancer among women worldwide, and estrogen receptor-α expression is a critical diagnostic factor for BC. Estrogen receptor (ER-α36) is a dominant-negative effector of ER-α66-mediated estrogen-responsive gene pathways. ER-α36 is a novel target that mediates the non-genomic estrogen signaling pathway. However, the crystallized structure of ER-α36 remains unavailable for molecular studies. ER-positive and triple-negative BC tumors aggressively resist the FDA-approved drugs; therefore, highly potent structure-based inhibitors with preeminent benefits over toxicity will preferably replace the current BC treatment. Broussoflanol B (BFB), a B. papyrifera bark compound, exhibits potent growth inhibitory activity in ER-negative BC cells by inducing cell cycle arrest. For the first time, we unravel the comparative dynamic events of the enzymes' structures and the binding mechanisms of BFB when bound to the ER-α36 and ER-α66 ligand-binding domain using an all-atom molecular dynamics simulations approach and MM/PBSA-binding-free energy calculations. The dynamic findings have revealed that ER-α36 and ER-α66 LBD undergo timescale "coiling", opening and closing conformations favoring the high-affinity BFB-bound ER-α36 (ΔG = -52.57 kcal/mol) compared to the BFB-bound ER-α66 (ΔG = -42.41 kcal/mol). Moreover, the unbound (1.260 Å) and bound ER-α36 (1.182 Å) exhibit the highest flexibilities and atomistic motions relative to the ER-α66 systems. The RMSF (Å) of the unbound ER-α36 and ER-α66 exhibit lesser stabilities than the BFB-bound systems, resulting in higher structural flexibilities and atomistic motions than the bound variants. These findings present a model that describes the mechanisms by which the BFB compound induces downregulation-accompanied cell cycle arrest at the Gap0 and Gap1 phases.

Nuclear-Electronic Orbital QM/MM Approach: Geometry Optimizations and Molecular Dynamics

Hybrid quantum mechanical/molecular mechanical (QM/MM) methods allow simulations of chemical reactions in atomistic solvent and heterogeneous environments such as proteins. Herein, the nuclear-electronic orbital (NEO) QM/MM approach is introduced to enable the quantization of specified nuclei, typically protons, in the QM region using a method such as NEO-density functional theory (NEO-DFT). This approach includes proton delocalization, polarization, anharmonicity, and zero-point energy in geometry optimizations and dynamics. Expressions for the energies and analytical gradients associated with the NEO-QM/MM method, as well as the previously developed polarizable continuum model (NEO-PCM), are provided. Geometry optimizations of small organic molecules hydrogen bonded to water in either dielectric continuum solvent or explicit atomistic solvent illustrate that aqueous solvation can strengthen hydrogen-bonding interactions for the systems studied, as indicated by shorter intermolecular distances at the hydrogen-bond interface. We then performed a real-time direct dynamics simulation of a phenol molecule in explicit water using the NEO-QM/MM method. These developments and initial examples provide the foundation for future studies of nuclear-electronic quantum dynamics in complex chemical and biological environments.

A Statistical Perspective on Microsolvation

The lack of a procedure to determine equilibrium thermodynamic properties of a small system interacting with a bath is frequently seen as a weakness of conventional statistical mechanics. A typical example for such a small system is a solute surrounded by an explicit solvation shell. One way to approach this problem is to enclose the small system of interest in a large bath of explicit solvent molecules, considerably larger than the system itself. The explicit inclusion of the solvent degrees of freedom is obviously limited by the available computational resources. A potential remedy to this problem is a microsolvation approach where only a few explicit solvent molecules are considered and surrounded by an implicit solvent bath. Still, the sampling of the solvent degrees of freedom is challenging with conventional grand canonical Monte Carlo methods, since no single chemical potential for the solvent molecules can be defined in the realm of small-system thermodynamics. In this work, a statistical thermodynamic model based on the grand canonical ensemble is proposed that avoids the conventional system size limitations and accurately characterizes the properties of the system of interest subject to the thermodynamic constraints of the bath. We extend an existing microsolvation approach to a generalized multibath "microstatistical" model and show that the previously derived approaches result as a limit of our model. The framework described here is universal and we validate our method numerically for a Lennard-Jones model fluid.

State-specific solvation for restricted active space spin-flip (RAS-SF) wave functions based on the polarizable continuum formalism

The restricted active space spin-flip (RAS-SF) formalism is a particular form of single-reference configuration interaction that can describe some forms of strong correlation at a relatively low cost and which has recently been formulated for the description of charge-transfer excited states. Here, we introduce both equilibrium and nonequilibrium versions of a state-specific solvation correction for vertical transition energies computed using RAS-SF wave functions, based on the framework of a polarizable continuum model (PCM). Ground-state polarization is described using the solvent's static dielectric constant and in the nonequilibrium solvation approach that polarization is modified upon vertical excitation using the solvent's optical dielectric constant. Benchmark calculations are reported for well-studied models of photo-induced charge transfer, including naphthalene dimer, C₂H₄⋯C₂F₄, pentacene dimer, and perylene diimide (PDI) dimer, several of which are important in organic photovoltaic applications. For the PDI dimer, we demonstrate that the charge-transfer character of the excited states is enhanced in the presence of a low-dielectric medium (static dielectric constant ɛ₀ = 3) as compared to a gas-phase calculation (ɛ₀ = 1). This stabilizes mechanistic traps for singlet fission and helps to explain experimental singlet fission rates. We also examine the effects of nonequilibrium solvation on charge-separated states in an intramolecular singlet fission chromophore, where we demonstrate that the energetic ordering of the states changes as a function of solvent polarity. The RAS-SF + PCM methodology that is reported here provides a framework to study charge-separated states in solution and in photovoltaic materials.

Software for the frontiers of quantum chemistry: An overview of developments in the Q-Chem 5 package

This article summarizes technical advances contained in the fifth major release of the Q-Chem quantum chemistry program package, covering developments since 2015. A comprehensive library of exchange-correlation functionals, along with a suite of correlated many-body methods, continues to be a hallmark of the Q-Chem software. The many-body methods include novel variants of both coupled-cluster and configuration-interaction approaches along with methods based on the algebraic diagrammatic construction and variational reduced density-matrix methods. Methods highlighted in Q-Chem 5 include a suite of tools for modeling core-level spectroscopy, methods for describing metastable resonances, methods for computing vibronic spectra, the nuclear-electronic orbital method, and several different energy decomposition analysis techniques. High-performance capabilities including multithreaded parallelism and support for calculations on graphics processing units are described. Q-Chem boasts a community of well over 100 active academic developers, and the continuing evolution of the software is supported by an "open teamware" model and an increasingly modular design.

Probing Interfacial Effects on Ionization Energies: The Surprising Banality of Anion-Water Hydrogen Bonding at the Air/Water Interface

Liquid microjet photoelectron spectroscopy is an increasingly common technique to measure vertical ionization energies (VIEs) of aqueous solutes, but the interpretation of these experiments is subject to questions regarding sensitivity to bulk versus interfacial solvation environments. We have computed aqueous-phase VIEs for a set of inorganic anions, using a combination of molecular dynamics simulations and electronic structure calculations, with results that are in excellent agreement with experiment regardless of whether the simulation data are restricted to ions at the air/water interface or to those in bulk aqueous solution. Although the computed VIEs are sensitive to ion-water hydrogen bonding, we find that the short-range solvation structure is sufficiently similar in both environments that it proves impossible to discriminate between the two on the basis of the VIE, a conclusion that has important implications for the interpretation of liquid-phase photoelectron spectroscopy. More generally, analysis of the simulation data suggests that the surface activity of soft anions is largely a second or third solvation shell effect, arising from disruption of water-water hydrogen bonds and not from significant changes in first-shell anion-water hydrogen bonding.

Benzene Radical Anion Microsolvated in Ammonia Clusters: Modeling the Transition from an Unbound Resonance to a Bound Species

The benzene radical anion, well-known in organic chemistry as the first intermediate in the Birch reduction of benzene in liquid ammonia, exhibits intriguing properties from the point of view of quantum chemistry. Notably, it has the character of a metastable shape resonance in the gas phase, while measurements in solution find it to be experimentally detectable and stable. In this light, our previous calculations performed in bulk liquid ammonia explicitly reveal that solvation leads to stabilization. Here, we focus on the transition of the benzene radical anion from an unstable gas-phase ion to a fully solvated bound species by explicit ionization calculations of the radical anion solvated in molecular clusters of increasing size. The computational cost of the largest systems is mitigated by combining density functional theory with auxiliary methods including effective fragment potentials or approximating the bulk by polarizable continuum models. Using this methodology, we obtain the cluster size dependence of the vertical binding energy of the benzene radical anion converging to the value of -2.3 eV at a modest computational cost.

How Well Does a Solvated Octa-acid Capsule Shield the Embedded Chromophore? A Computational Analysis Based on an Anisotropic Dielectric Continuum Model

The optical properties of chromophores embedded in a water-solvated dimer of octa-acid that forms a molecular-shaped capsule are investigated. In particular, we address the anisotropic dielectric environment that appears to blue-shift excitation energies compared to the free aqueous chromophores. Recently we reported that using an effective scalar dielectric constant ε ≈ 3 appears to reproduce the measured spectra of the embedded coumarins, suggesting that the capsule provides a significant, albeit not perfect, screening of the aqueous dielectric environment. Here, we report absorption energies using a theoretical treatment that includes continuum solvation affected by an anisotropic dielectric function reflecting the high-dielectric environment outside of the capsule and the low-dielectric region within. We report time-dependent density functional theory calculations using a range-separated functional with the Poisson boundary conditions that model the anisotropic dielectric environment. Our calculations find that the anisotropic environment due to the water-solvated hydrophobic capsule is equivalent to a homogeneous effective dielectric constant of ≈3. The calculated values also appear to reproduce measured absorption of the embedded coumarin, where we study the effect of the hydrophobic capsule on the excited state.