Computational drug discovery of an inhibitor of APOBEC3B as a treatment for epithelial cancers.

Cancer is one of the leading causes of death in the U.S., and tumorous cancers such as cervical, lung, breast, and ovarian cancers are the most common types. APOBEC3B is a nonessential cytidine deaminase found in humans and theorized to defend against viral infection. However, overexpression of APOBEC3B is linked to cancer in humans, which makes APOBEC3B a potential cancer treatment target through competitive inhibition for several tumorous cancers. Computational studies can help reveal a small molecule inhibitor using high-throughput virtual screening of millions of candidates with relatively little cost. This study aims to narrow the field of potential APOBEC3B inhibition candidates for future in vitro assays and provide an effective scaffold for drug design studies. Another goal of this project is to provide critical amino acid targets in the active site for future drug design studies. This study simulated 7.8 million drug candidates using high-throughput virtual screening and further processed the top scoring 241 molecules from AutoDock Vina, DOCK 6, and de novo design. Using virtual screening, de novo design, and molecular dynamics simulations, a competitive inhibitor candidate was discovered with an average binding free energy score of -46.03 kcal/mol, more than 10 kcal/mol better than the substrate control (dCMP). These results indicate that this molecule (or a structural derivative) may be an effective inhibitor of APOBEC3B and prevent host genome mutagenesis resulting from protein overexpression. Another important finding is the confirmation of essential amino acid targets, such as Tyr250 and Gln213 within the active site of APOBEC3B. Therefore, study used novel computational methods to provide a theoretical scaffold for future drug design studies that may prove useful as a treatment for epithelial cancers.Communicated by Ramaswamy H. Sarma.

Exploring the disruption of SARS-CoV-2 RBD binding to hACE2.

The COVID-19 pandemic was declared due to the spread of the novel coronavirus, SARS-CoV-2. Viral infection is caused by the interaction between the SARS-CoV-2 receptor binding domain (RBD) and the human ACE2 receptor (hACE2). Previous computational studies have identified repurposed small molecules that target the RBD, but very few have screened drugs in the RBD-hACE2 interface. When studies focus solely on the binding affinity between the drug and the RBD, they ignore the effect of hACE2, resulting in an incomplete analysis. We screened ACE inhibitors and previously identified SARS-CoV-2 inhibitors for binding to the RBD-hACE2 interface, and then conducted 500 ns of unrestrained molecular dynamics (MD) simulations of fosinopril, fosinoprilat, lisinopril, emodin, diquafosol, and physcion bound to the interface to assess the binding characteristics of these ligands. Based on MM-GBSA analysis, all six ligands bind favorably in the interface and inhibit the RBD-hACE2 interaction. However, when we repeat our simulation by first binding the drug to the RBD before interacting with hACE2, we find that fosinopril, fosinoprilat, and lisinopril result in a strongly interacting trimeric complex (RBD-drug-hACE2). Hydrogen bonding and pairwise decomposition analyses further suggest that fosinopril is the best RBD inhibitor. However, when lisinopril is bound, it stabilizes the trimeric complex and, therefore, is not an ideal potential drug candidate. Overall, these results reveal important atomistic interactions critical to the binding of the RBD to hACE2 and highlight the significance of including all protein partners in the evaluation of a potential drug candidate.

Variations on the Bergman Cyclization Theme: Electrocyclizations of Ionic Penta-, Hepta-, and Octadiynes.

The Bergman cyclization of (Z)-hexa-3-ene-1,5-diyne to form the aromatic diradical p-benzyne has garnered attention as a potential antitumor agent due to its relatively low cyclization barrier and the stability of the resulting diradical. Here, we present a theoretical investigation of several ionic extensions of the fundamental Bergman cyclization: electrocyclizations of the penta-1,4-diyne anion, hepta-1,6-diyne cation, and octa-1,7-diyne dication, leveraging the spin-flip formulation of the equation-of-motion coupled cluster theory with single and double substitutions (EOM-SF-CCSD). Though the penta-1,4-diyne anion exhibits a large cyclization barrier of +66 kcal mol-1, cyclization of both the hepta-1,6-diyne cation and octa-1,7-diyne dication along a previously unreported triplet pathway requires relatively low energy. We also identified the presence of significant aromaticity in the triplet diradical products of these two cationic cyclizations.

Atomistic insights into the binding of SARS-CoV-2 spike receptor binding domain with the human ACE2 receptor: The importance of residue 493.

SARS-CoV-2 is a coronavirus that has created a global pandemic. The virus contains a spike protein which has been shown to bind to the ACE2 receptor on the surface of human cells. Vaccines have been developed that recognize elements of the SARS-CoV-2 spike protein and they have been successful in preventing infection. Recently, the Omicron variant of the SARS-CoV-2 virus was reported and quickly became a variant of concern due to its transmissibility. This variant contained an unusually large number (32) of point mutations, of which 15 of those mutations are in the receptor binding domain of the spike protein. While several computational and experimental investigations comparing the binding of the Omicron and wild type RBD to the human ACE2 receptor have been conducted, many of these report contradictory findings. In order to assess the differential binding ability, we conducted 2 µs of classical molecular dynamics (cMD) simulation to estimate the binding affinities and behaviors. Based upon MM-GBSA binding affinity, per-residue energy decomposition analysis, center of mass distance measurements, ensemble clustering, pairwise residue decomposition and hydrogen bonding analysis, our results suggest that a single point mutation is responsible for the enhanced binding of the Omicron mutant relative to the WT. While the 15-point mutations in the receptor binding domain contribute positively and negatively to the affinity of the spike protein for the human ACE2 receptor, it is the point mutation Q493R that confers enhanced binding while the Q493K mutation results in similar binding. The MM-GBSA binding estimations over a 2 µs trajectory, suggest that the wild type binds to ACE2 with a value of -29.69 kcal/mol while the Q493K and Q493R Omicron mutants bind with energy values of -26.67 and -34.56 kcal/mol, respectively. These values are significantly different, given the error estimates associated with the MM-GBSA method. In general, while some mutations increase binding, more mutations diminish binding, leading to an overall similar picture of binding for Q493K and enhanced binding for Q493R.

Experimental and Theoretical Study of Oxolan-3-one Thermal Decomposition.

The thermal decomposition of oxolan-3-one, a common component of the bio-oil formed during biomass pyrolysis, has been studied using ab initio calculations and experiments employing pulsed gas-phase pyrolysis with matrix-isolation FTIR product detection. Four pathways for unimolecular decomposition were predicted using computational methods. The dominant reaction channel led to carbon monoxide, formaldehyde, and ethylene, all of which were observed experimentally. The other channels led to an assortment of products including ketene, water, propyne, and acetylene, which were all confirmed in the matrix-isolation FTIR spectra. There is also evidence for the production of substituted ketenes in pyrolysis, most likely hydroxyketene and methylketene.

The Impacts of the Molecular Education and Research Consortium in Undergraduate Computational Chemistry on the Careers of Women in Computational Chemistry.

The Molecular Education and Research Consortium in Undergraduate Computational Chemistry (MERCURY) has supported a diverse group of faculty and students for over 20 years by providing computational resources as well as networking opportunities and professional support. The consortium comprises 38 faculty (42% women) at 34 different institutions, who have trained nearly 900 undergraduate students, more than two-thirds of whom identify as women and one-quarter identify as students of color. MERCURY provides a model for the support necessary for faculty to achieve professional advancement and career satisfaction. The range of experiences and expertise of the consortium members provides excellent networking opportunities that allow MERCURY faculty to support each other's teaching, research, and service needs, including generating meaningful scientific advancements and outcomes with undergraduate researchers as well as being leaders at the departmental, institutional, and national levels. While all MERCURY faculty benefit from these supports, the disproportionate number of women in the consortium, relative to their representation in computational sciences generally, produces a sizable impact on advancing women in the computational sciences. In this report, the women of MERCURY share how the consortium has benefited their careers and the careers of their students.

New 1,2,3-Triazoles from (R)-Carvone: Synthesis, DFT Mechanistic Study and In Vitro Cytotoxic Evaluation.

Aseries of novel 1,4-disubstituted 1,2,3-triazoles were synthesized from an (R)-carvone terminal alkyne derivative via a Cu (I)-catalyzed azide-alkyne cycloaddition reaction using CuSO4,5H2O as the copper (II) source and sodium ascorbate as a reducing agent which reduces Cu (II) into Cu (I). All the newly synthesized 1,2,3-triazoles 9a-h were fully identified on the basis of their HRMS and NMR spectral data and then evaluated for their cell growth inhibition potential by MTS assay against HT-1080 fibrosarcoma, A-549 lung carcinoma, and two breast adenocarcinoma (MCF-7 and MDA-MB-231) cell lines. Compound 9d showed notable cytotoxic effects against the HT-1080 and MCF-7 cells with IC50 values of 25.77 and 27.89 µM, respectively, while compound 9c displayed significant activity against MCF-7 cells with an IC50 value of 25.03 µM. Density functional calculations at the B3LYP/6-31G* level of theory were used to confirm the high reactivity of the terminal alkyne as a dipolarophile. Quantum calculations were also used to investigate the mechanism of both the uncatalyzed and copper (I)-catalyzed azide-alkyne cycloaddition reaction (CuAAC). The catalyzed reaction gives complete regioselectivity via a stepwise mechanism streamlining experimental observations. The calculated free-energy barriers 4.33 kcal/mol and 29.35 kcal/mol for the 1,4- and 1,5-regioisomers, respectively, explain the marked regioselectivity of the CuAAC reaction.

Investigating novel thiazolyl-indazole derivatives as scaffolds for SARS-CoV-2 MPro inhibitors.

COVID-19 is a global pandemic caused by infection with the SARS-CoV-2 virus. Remdesivir, a SARS-CoV-2 RNA polymerase inhibitor, is the only drug to have received widespread approval for treatment of COVID-19. The SARS-CoV-2 main protease enzyme (MPro), essential for viral replication and transcription, remains an active target in the search for new treatments. In this study, the ability of novel thiazolyl-indazole derivatives to inhibit MPro is evaluated. These compounds were synthesized via the heterocyclization of phenacyl bromide with (R)-carvone, (R)-pulegone and (R)-menthone thiosemicarbazones. The binding affinity and binding interactions of each compound were evaluated through Schrödinger Glide docking, AMBER molecular dynamics simulations, and MM-GBSA free energy estimation, and these results were compared with similar calculations of MPro binding various 5-mer substrates (VKLQA, VKLQS, VKLQG) and a previously identified MPro tight-binder X77. From these simulations, we can see that binding is driven by residue specific interactions such as π-stacking with His41, and S/π interactions with Met49 and Met165. The compounds were also experimentally evaluated in a MPro biochemical assay and the most potent compound containing a phenylthiazole moiety inhibited protease activity with an IC50 of 92.9 â€‹µM. This suggests that the phenylthiazole scaffold is a promising candidate for the development of future MPro inhibitors.

Evaluating Halogen-Bond Strength as a Function of Molecular Structure Using Nuclear Magnetic Resonance Spectroscopy and Computational Analysis.

Halogen bonding (XB) is a highly directional, non-covalent intermolecular interaction between a molecule (XB donor) presenting a halogen with an electron-deficient region or sigma hole (σ-hole) and an electron-rich or Lewis-base molecule (XB acceptor). A systematic, experimental, and theoretical study of solution-phase XB strength as a function of the molecular structure for both XB donor and acceptor molecules is presented. The impact of specific structural features is assessed using 19F and 1H nuclear magnetic resonance (NMR) titrations to determine association constants, density functional theory calculations for interaction energies and bond lengths, as well as 19F-1H HOESY NMR measurements of intermolecular cross-relaxation between the interacting XB donor-acceptor adducts. For XB donor molecules (perfluoro-halogenated benzenes), results indicate the critical importance of iodine coupled with electron-withdrawing entities. Prominent structural components of XB acceptor molecules include a central atom working in conjunction with a Lewis-base atom to present high electron density directed at the σ-hole (e.g., tributylphosphine oxide). Additionally, larger surrounding aliphatic R groups (e.g., butyl and octyl) were found to significantly stabilize strong XB, particularly in solvents that promote the interaction. With a more thorough understanding of structure-optimized XB, one can envision harnessing XB interactions more strategically for specific design of optimal materials and chemical applications.

A comparison of the chemical bonding and reactivity of Si8H8O12 and Ge8H8O12: A theoretical study.

We have analyzed the chemical bonding and reactivity in the cubic molecule octahydridosilsesquioxane, Si8H8O12, and its counterpart Ge8H8O12 by means of ab initio quantum chemical methods and group theory. Density functional theory and MP2 methods combined with the basis sets 6-311+G(d) and 6-311++G(2d,p) were used for geometry optimization and vibrational frequency analysis. The geometries of Si8H8O12 and Ge8H8O12 are unstable under Oh symmetry and distort to the rare Th molecular symmetry. The energy gained from this pseudo-Jahn-Teller distortion ranges from 0.78 to 6.14 kcal mol-1 depending on methodological treatment. The Fukui functions and the molecular electrostatic potential were both used as DFT-based reactivity descriptors. Our study shows that Si8H8O12 and Ge8H8O12 are both hard amphoteric molecules. The cavity within each cage is acidic and able to encapsulate hard small bases such as F-. The exterior of the cages is basic and can form stable exohedral complexes with hard acids, as in the case of H+. The insertion of F- in Si8H8O12 and Ge8H8O12 cages gives the most stable endohedral complexes of the series studied, characterized by formation energies of -3.50 and -3.45 eV at CAM-B3LYP/6-311+G(d) and -3.61 and -3.68 eV at the MP2/6-311++G(d,p) level, respectively. The calculated formation energies of the exohedral and endohedral complexes align with the DFT reactivity descriptor analysis.

Identification of Phenazine-Based MEMO1 Small-Molecule Inhibitors: Virtual Screening, Fluorescence Polarization Validation, and Inhibition of Breast Cancer Migration.

Phosphorylation-dependent protein-protein interactions play a significant role in biological signaling pathways; therefore, small molecules that are capable of influencing these interactions can be valuable research tools and have potential as pharmaceutical agents. MEMO1 (mediator of ErbB2-cell driven motility) is a phosphotyrosine-binding protein that interacts with a variety of protein partners and has been found to be upregulated in breast cancer patients. Herein, we report the first small-molecule inhibitors of MEMO1 interactions identified through a virtual screening platform and validated in a competitive fluorescence polarization assay. Initial structure-activity relationships have been investigated for these phenazine-core inhibitors and the binding sites have been postulated using molecular dynamics simulations. The most potent biochemical inhibitor is capable of disrupting the large protein interface with a KI of 2.7 µm. In addition, the most promising phenazine core compounds slow the migration of breast cancer cell lines in a scratch assay.

Understanding the Structure and Apo Dynamics of the Functionally Active JIP1 Fragment.

Recent experiments indicate that the C-Jun amino-terminal kinase-interacting protein 1 (JIP1) binds to and activates the c-Jun N-terminal kinase (JNK) protein. JNK is an integral part of cell apoptosis, and misregulation of this process is a causative factor in diseases such as Alzheimer's disease (AD), obesity, and cancer. It has also been shown that JIP1 may increase the phosphorylation of tau by facilitating the interaction between the tau protein and JNK, which could also be a causative factor in AD. Very little is known about the structure and dynamics of JIP1; however, the amino acid composition of the first 350 residues suggests that it contains an intrinsically disordered region. Molecular dynamics (MD) simulations using AMBER 14 were used to study the structure and dynamics of a functionally active JIP1 10mer fragment to better understand the solution behavior of the fragment. Two microseconds of unbiased MD was performed on the JIP1 10mer fragment in 10 different seeds for a total of 20 µs of simulation time, and from this, seven structurally stable conformations of the 10mer fragment were identified via classical clustering. The 10mer ensemble was also used to build a Markov state model (MSM) that identified four metastable states that encompassed six of the seven conformational families identified by classical dimensional reduction. Based on this MSM, conformational interconversions between the four states occur via two dominant pathways with probability fluxes of 55 and 44% for each individual pathway. Transitions between the initial and final states occur with mean first passage times of 31 (forward) and 16 (reverse) µs.

Isofunctional Clustering and Conformational Analysis of the Arsenate Reductase Superfamily Reveals Nine Distinct Clusters.

Arsenate reductase (ArsC) is a superfamily of enzymes that reduce arsenate. Due to active site similarities, some ArsC can function as low-molecular weight protein tyrosine phosphatases (LMW-PTPs). Broad superfamily classifications align with redox partners (Trx- or Grx-linked). To understand this superfamily's mechanistic diversity, the ArsC superfamily is classified on the basis of active site features utilizing the tools TuLIP (two-level iterative clustering process) and autoMISST (automated multilevel iterative sequence searching technique). This approach identified nine functionally relevant (perhaps isofunctional) protein groups. Five groups exhibit distinct ArsC mechanisms. Three are Grx-linked: group 4AA (classical ArsC), group 3AAA (YffB-like), and group 5BAA. Two are Trx-linked: groups 6AAAAA and 7AAAAAAAA. One is an Spx-like transcriptional regulatory group, group 5AAA. Three are potential LMW-PTP groups: groups 7BAAAA, and 7AAAABAA, which have not been previously identified, and the well-studied LMW-PTP family group 8AAA. Molecular dynamics simulations were utilized to explore functional site details. In several families, we confirm and add detail to literature-based mechanistic information. Mechanistic roles are hypothesized for conserved active site residues in several families. In three families, simulations of the unliganded structure sample specific conformational ensembles, which are proposed to represent either a more ligand-binding-competent conformation or a pathway toward a more binding-competent state; these active sites may be designed to traverse high-energy barriers to the lower-energy conformations necessary to more readily bind ligands. This more detailed biochemical understanding of ArsC and ArsC-like PTP mechanisms opens possibilities for further understanding of arsenate bioremediation and the LMW-PTP mechanism.

The generality of the GUGA MRCI approach in COLUMBUS for treating complex quantum chemistry.

The core part of the program system COLUMBUS allows highly efficient calculations using variational multireference (MR) methods in the framework of configuration interaction with single and double excitations (MR-CISD) and averaged quadratic coupled-cluster calculations (MR-AQCC), based on uncontracted sets of configurations and the graphical unitary group approach (GUGA). The availability of analytic MR-CISD and MR-AQCC energy gradients and analytic nonadiabatic couplings for MR-CISD enables exciting applications including, e.g., investigations of π-conjugated biradicaloid compounds, calculations of multitudes of excited states, development of diabatization procedures, and furnishing the electronic structure information for on-the-fly surface nonadiabatic dynamics. With fully variational uncontracted spin-orbit MRCI, COLUMBUS provides a unique possibility of performing high-level calculations on compounds containing heavy atoms up to lanthanides and actinides. Crucial for carrying out all of these calculations effectively is the availability of an efficient parallel code for the CI step. Configuration spaces of several billion in size now can be treated quite routinely on standard parallel computer clusters. Emerging developments in COLUMBUS, including the all configuration mean energy multiconfiguration self-consistent field method and the graphically contracted function method, promise to allow practically unlimited configuration space dimensions. Spin density based on the GUGA approach, analytic spin-orbit energy gradients, possibilities for local electron correlation MR calculations, development of general interfaces for nonadiabatic dynamics, and MRCI linear vibronic coupling models conclude this overview.

A Molecular Dynamics Investigation of the Thermostability of Cold-Sensitive I707L KlenTaq1 DNA Polymerase and Its Wild-Type Counterpart.

DNA polymerase I from Thermus aquaticus ( Taq DNA polymerase) is useful for polymerase chain reactions because of its exceptional thermostability; however, its activity at low temperatures can cause amplification of unintended products. Mutation of isoleucine 707 to leucine (I707L) slows Taq DNA polymerase at low temperatures, which decreases unwanted amplification due to mispriming. In this work, unrestrained molecular dynamics (MD) simulations were performed on I707L and wild-type (WT) Taq DNA polymerase at 341 and 298 K to determine how the mutation affects the dynamic nature of the protein. The results suggest that I707L Taq DNA polymerase remains relatively immobile at room temperature and becomes more flexible at the higher temperature, while the WT Taq DNA polymerase demonstrates less substantial differences in dynamics at high and low temperatures. These results are in agreement with previous experimental results on the I707L mutant Taq DNA polymerase that show dynamic differences at high and low temperatures. The decreased mobility of the mutant at low temperature suggests that the mutant remains longer in the blocked conformation, and this may lead to reduced activity relative to the WT at 298 K. Principal component analysis revealed that the mutation results in decoupled movements of the Q helix and fingers domain. This decoupled nature of the mutant gives way to an increasingly flexible N-terminal end of the Q helix at 341 K, a characteristic not seen for WT Taq DNA polymerase.

A Multireference Ab Initio Study of the Diradical Isomers of Pyrazine.

Three diradical pyrazine isomers were characterized using highly correlated, multireference methods. The lowest lying singlet and triplet state geometries of 2,3-didehydropyrazine ( ortho), 2,5-didehydropyrazine ( para), and 2,6-didehydropyrazine ( meta) were determined. Two active reference spaces were utilized. The complete active space (CAS) (8,8) includes the σ and σ* orbitals on the dehydrocarbon atoms as well as the valence π and π* orbitals. The CAS (12,10) reference space includes two additional orbitals corresponding to the in-phase and out-of-phase nitrogen lone pair orbitals. Adiabatic and vertical gaps between the lowest lying singlet and triplet states, optimized geometries, canonicalized orbital energies, unpaired electron densities, and spin polarization effects were compared. We find that the singlet states of each diradical isomer contain two significantly weighted configurations, and the larger active space is necessary for the proper physical characterization of both the singlet and triplet states. The singlet-triplet splitting is very small for the 2,3-didehydropyrazine ( ortho) and 2,6-didehydropyrazine ( meta) isomers (+1.8 and -1.4 kcal/mol, respectively) and significant for the 2,5-didehydropyrazine ( para) isomer (+28.2 kcal/mol). Singlet geometries show through-space interactions between the dehydocarbon atoms in the 2,3-didehydropyrazine ( ortho) and 2,6-didehydropyrazine ( meta) isomers. An analysis of the effectively unpaired electrons suggests that the 2,5-didehydropyrazine ( para) isomer also displays through-bond coupling between the diradical electrons.

Halogen Bonding Interactions for Aromatic and Nonaromatic Explosive Detection.

Improved sensing strategies are needed for facile, accurate, and rapid detection of aromatic and nonaromatic explosives. Density functional theory was used to evaluate the relative binding interaction energies between halogen-containing sensor model molecules and nitro-containing explosives. Interaction energies ranged from -18 to -14 kJ/mol and highly directional halogen bonding interactions were observed with bond distances ranging between 3.0 and 3.4 Å. In all geometry optimized structures, the sigma-hole of electropositive potential on the halogen aligned with a lone pair of electrons on the nitro-moiety of the explosive. The computational results predict that the strongest interactions will occur with iodine-based sensors as, of all the halogens studied, iodine is the largest, most polarizable halogen with the smallest electronegativity. Based on these promising proof-of-concept results, synthetically accessible sensors were designed using 1,4-dihalobenzene (X = Cl, Br, and I) with and without tetra-fluoro electron withdrawing groups attached to the benzene ring. These sensing molecules were embedded onto single walled carbon nanotubes that were mechanically abraded onto interdigitated array electrodes, and these were used to measure the responses to explosive model compounds cyclohexanone and dimethyl-dinitro-benzene in nitrogen gas. Amperometric current-time curves for selectors and control molecules, including concentration correlated signal enhancement, as well as response and recovery times, indicate selector responsiveness to these model compounds, with the largest response observed for iodo-substituted sensors.

Analysis of MEMO1 Binding Specificity for ErbB2 Using Fluorescence Polarization and Molecular Dynamics Simulations.

ErbB2 signaling pathways are linked to breast cancer formation, growth, and aggression; therefore, understanding the behavior of proteins associated with these pathways as well as regulatory factors that influence ErbB2 function is essential. MEMO1 is a redox active protein that is shown to associate with phosphorylated ErbB2 and mediate cell motility. We have developed a fluorescence polarization assay to probe the interaction between MEMO1 and an ErbB2-derived peptide containing a phosphorylated tyrosine residue. This interaction is shown to be pH-dependent and stronger with longer peptides as would be expected for protein-protein interactions. We have quantitatively mapped the binding interface of MEMO1 to the peptide using the fluorescence polarization assay and molecular dynamics simulations. We have confirmed that phosphorylation of the peptide is essential for binding and through mutagenesis have identified residues that contribute to favorable interactions. Our results highlight the importance of the protein-protein interactions of MEMO1 that complement the oxidase activity. In the future, these studies will provide a method for screening for selective modulators of MEMO1, which will allow for additional biological investigations.