A recently identified variant of SARS-CoV-2 virus, known as the United Kingdom (UK) variant (lineage B.1.1.7), has an N501Y mutation on its spike protein. SARS-CoV-2 spike protein binds with angiotensin-converting enzyme 2 (ACE2), a key protein for the viral entry into the host cells. Here, we report an efficient computational approach, including the simple energy minimizations and binding free energy calculations, starting from an experimental structure of the binding complex along with experimental calibration of the calculated binding free energies, to rapidly and reliably predict the binding affinities of the N501Y mutant with human ACE2 (hACE2) and recently reported miniprotein and hACE2 decoy (CTC-445.2) drug candidates. It has been demonstrated that the N501Y mutation markedly increases the ACE2-spike protein binding affinity (K d) from 22 to 0.44 nM, which could partially explain why the UK variant is more infectious. The miniproteins are predicted to have ∼10,000- to 100,000-fold diminished binding affinities with the N501Y mutant, creating a need for design of novel therapeutic candidates to overcome the N501Y mutation-induced drug resistance. The N501Y mutation is also predicted to decrease the binding affinity of a hACE2 decoy (CTC-445.2) binding with the spike protein by ∼200-fold. This convenient computational approach along with experimental calibration may be similarly used in the future to predict the binding affinities of potential new variants of the spike protein.

Variants of the SARS-CoV-2 virus continue to remain a threat 2 years from the beginning of the pandemic. As more variants arise, and the B.1.1.529 (Omicron) variant threatens to create another wave of infections, a method is needed to predict the binding affinity of the spike protein quickly and accurately with human angiotensin-converting enzyme II (ACE2). We present an accurate and convenient energy minimization/molecular mechanics Poisson–Boltzmann surface area methodology previously used with engineered ACE2 therapeutics to predict the binding affinity of the Omicron variant. Without any additional data from the variants discovered after the publication of our first model, the methodology can accurately predict the binding of the spike/ACE2 variant complexes. From this methodology, we predicted that the Omicron variant spike has a K d of ∼22.69 nM (which is very close to the experimental K d of 20.63 nM published during the review process of the current report) and that spike protein of the new “Stealth” Omicron variant (BA.2) will display a K d of ∼12.9 nM with the wild-type ACE2 protein. This methodology can be used with as-yet discovered variants, allowing for quick determinations regarding the variant’s infectivity versus either the wild-type virus or its variants.

The cannabinoid (CB) receptors (CB 1 R and CB 2 R) represent a promising therapeutic target for several indications such as nociception and obesity. The ligands with nonselectivity can be traced to the high similarity in the binding sites of both cannabinoid receptors. Therefore, the need for selectivity, potency, and G-protein coupling bias has further complicated the design of desired compounds. The bias of currently studied cannabinoid agonists is seldom investigated, and agonists that do exhibit bias are typically nonselective. However, certain long-chain endocannabinoids represent a class of selective and potent CB 1 R agonists. The binding mode for this class of compounds has remained elusive, limiting the implementation of its binding features to currently studied agonists. Hence, in the present study, the binding poses for these long-chain cannabinoids, along with other interesting ligands, with the receptors have been determined, by using a combination of molecular docking and molecular dynamics (MD) simulations along with molecular mechanics-Poisson−Boltzmann surface area (MM-PBSA) binding free energy calculations. The binding poses for the long-chain cannabinoids implicate that a site surrounded by the transmembrane (TM)2, TM7, and extracellular loop (ECL)2 is vital for providing the long-chain ligands with the selectivity for CB 1 R, especially I267 of CB 1 R (corresponding to L182 of CB 2 R). Based on the obtained binding modes, the calculated relative binding free energies and selectivity are all in good agreement with the corresponding experimental data, suggesting that the determined binding poses are reasonable. The computational strategy used in this study may also prove fruitful in applications with other GPCRs or membrane-bound proteins.

No abstract

The central relaxin-3/RXFP3 system plays important roles in stress responses, feeding, and motivation for reward. However, exploration of its therapeutic applications has been hampered by the lack of small molecule ligands and the cross-activation of RXFP1 in the brain and RXFP4 in the periphery. Herein, we report the first structure–activity relationship studies of a series of novel nonpeptide amidinohydrazone-based agonists, which were characterized by RXFP3 functional and radioligand binding assays. Several potent and efficacious RXFP3 agonists (e.g., 10d) were identified with EC50 values <10 nM. These compounds also had high potency at RXFP4 but no agonist activity at RXFP1, demonstrating > 100-fold selectivity for RXFP3/4 over RXFP1. In vitro ADME and pharmacokinetic assessments revealed that the amidinohydrazone derivatives may have limited brain permeability. Collectively, our findings provide the basis for further optimization of lead compounds to develop a suitable agonist to probe RXFP3 functions in the brain.