The quality of the 1.14*CM1A and 1.20*CM5 charge models was evaluated for calculations of free energies of hydration. For a set of 426 neutral molecules, 1.14*CM1A and 1.20*CM5 yield MADs of 1.26 and 1.20 kcal/mol, respectively. The 1.14*CM1A charges, which can be readily obtained for large systems, exhibit large deviations only for a subset of functional groups. The results for these cases were systematically improved using Localized Bond Charge Corrections (LBCC) by which off-setting adjustments are made to the partial charges for atoms in specified bond types. Only 19 LBCCs were needed to yield 1.14*CM1A-LBCC charges that reduce the errors for the 426 ΔGhyd values to only 0.61 kcal/mol. The modified charge method was also tested in computation of heats of vaporization and densities for pure organic liquids, yielding average errors of 1.40 kcal/mol and 0.024 g/cm3, similar to those for 1.14*CM1A.

We present the first receptor-based pharmacophore model for HIV-1 integrase. The development of "dynamic" pharmacophore models is a new method that accounts for the inherent flexibility of the active site and aims to reduce the entropic penalties associated with binding a ligand. Furthermore, this new drug discovery method overcomes the limitation of an incomplete crystal structure of the target protein. A molecular dynamics (MD) simulation describes the flexibility of the uncomplexed protein. Many conformational models of the protein are saved from the MD simulations and used in a series of multi-unit search for interacting conformers (MUSIC) simulations. MUSIC is a multiple-copy minimization method, available in the BOSS program; it is used to determine binding regions for probe molecules containing functional groups that complement the active site. All protein conformations from the MD are overlaid, and conserved binding regions for the probe molecules are identified. Those conserved binding regions define the dynamic pharmacophore model. Here, the dynamic model is compared to known inhibitors of the integrase as well as a three-point, ligand-based pharmacophore model from the literature. Also, a "static" pharmacophore model was determined in the standard fashion, using a single crystal structure. Inhibitors thought to bind in the active site of HIV-1 integrase fit the dynamic model but not the static model. Finally, we have identified a set of compounds from the Available Chemicals Directory that fit the dynamic pharmacophore model, and experimental testing of the compounds has confirmed several new inhibitors.

A consensus virtual screening protocol has been applied to ca. 2000 approved drugs to seek inhibitors of the main protease (M pro ) of SARS-CoV-2, the virus responsible for COVID-19. 42 drugs emerged as top candidates, and after visual analyses of the predicted structures of their complexes with M pro , 17 were chosen for evaluation in a kinetic assay for M pro inhibition. Remarkably 14 of the compounds at 100-μM concentration were found to reduce the enzymatic activity and 5 provided IC 50 values below 40 μM: manidipine (4.8 μM), boceprevir (5.4 μM), lercanidipine (16.2 μM), bedaquiline (18.7 μM), and efonidipine (38.5 μM). Structural analyses reveal a common cloverleaf pattern for the binding of the active compounds to the P1, P1′, and P2 pockets of M pro . Further study of the most active compounds in the context of COVID-19 therapy is warranted, while all of the active compounds may provide a foundation for lead optimization to deliver valuable chemotherapeutics to combat the pandemic.

A consensus virtual screening protocol has been applied to ca. 2000 approved drugs to seek inhibitors of the main protease (Mpro) of SARS-CoV-2, the virus responsible for COVID-19. 42 drugs emerged as top candidates, and after visual analyses of the predicted structures of their complexes with Mpro, 17 were chosen for evaluation in a kinetic assay for Mpro inhibition. Remarkably 14 of the compounds at 100-μM concentration were found to reduce the enzymatic activity and 5 provided IC50 values below 40 μM: manidipine (4.8 μM), boceprevir (5.4 μM), lercanidipine (16.2 μM), bedaquiline (18.7 μM), and efonidipine (38.5 μM). Structural analyses reveal a common cloverleaf pattern for the binding of the active compounds to the P1, P1’, and P2 pockets of Mpro. Further study of the most active compounds in the context of COVID-19 therapy is warranted, while all of the active compounds may provide a foundation for lead optimization to deliver valuable chemotherapeutics to combat the pandemic.

Development of resistance remains a major challenge for drugs to treat HIV-1 infections, including those targeting the essential viral polymerase, HIV-1 reverse transcriptase (RT). Resistance associated with the Tyr181Cys mutation in HIV-1 RT has been a key roadblock in the discovery of nonnucleoside RT inhibitors (NNRTIs). It is the principal point mutation that arises from treatment of HIV-infected patients with nevirapine, the first-in-class drug still widely used, especially in developing countries. We report covalent inhibitors of Tyr181Cys RT (CRTIs) that can completely knock out activity of the resistant mutant and of the particularly challenging Lys103Asn/Tyr181Cys variant. Conclusive evidence for the covalent modification of Cys181 is provided from enzyme inhibition kinetics, mass spectrometry, protein crystallography, and antiviral activity in infected human T-cell assays. The CRTIs are also shown to be selective for Cys181 and have lower cytotoxicity than the approved NNRTI drugs efavirenz and rilpivirine.