Computational design of small molecular modulators of protein-protein interactions with a novel thermodynamic cycle: Allosteric inhibitors of HIV-1 integrase

Adenine Methylation Enhances the Conformational Flexibility of an RNA Hairpin Tetraloop

The N6-methyladenosine modification is one of the most abundant post-transcriptional modifications in ribonucleic acid (RNA) molecules. Using molecular dynamics simulations and alchemical free-energy calculations, we studied the structural and energetic implications of incorporating this modification in an adenine mononucleotide and an RNA hairpin structure. At the mononucleotide level, we found that the syn configuration is more favorable than the anti configuration by 2.05 ± 0.15 kcal/mol. The unfavorable effect of methylation was due to the steric overlap between the methyl group and a nitrogen atom in the purine ring. We then probed the effect of methylation in an RNA hairpin structure containing an AUCG tetraloop, which is recognized by a "reader" protein (YTHDC1) to promote transcriptional silencing of long noncoding RNAs. While methylation had no significant conformational effect on the hairpin stem, the methylated tetraloop showed enhanced conformational flexibility compared to the unmethylated tetraloop. The increased flexibility was associated with the outward flipping of two bases (A6 and U7) which formed stacking interactions with each other and with the C8 and G9 bases in the tetraloop, leading to a conformation similar to that in the RNA/reader protein complex. Therefore, methylation-induced conformational flexibility likely facilitates RNA recognition by the reader protein.

The Drug-Induced Interface That Drives HIV-1 Integrase Hypermultimerization and Loss of Function

Allosteric HIV-1 integrase (IN) inhibitors (ALLINIs) are an emerging class of small molecules that disrupt viral maturation by inducing the aberrant multimerization of IN. Here, we present cocrystal structures of HIV-1 IN with two potent ALLINIs, namely, BI-D and the drug candidate Pirmitegravir. The structures reveal atomistic details of the ALLINI-induced interface between the HIV-1 IN catalytic core and carboxyl-terminal domains (CCD and CTD). Projecting from their principal binding pocket on the IN CCD dimer, the compounds act as molecular glue by engaging a triad of invariant HIV-1 IN CTD residues, namely, Tyr226, Trp235, and Lys266, to nucleate the CTD-CCD interaction. The drug-induced interface involves the CTD SH3-like fold and extends to the beginning of the IN carboxyl-terminal tail region. We show that mutations of HIV-1 IN CTD residues that participate in the interface with the CCD greatly reduce the IN-aggregation properties of Pirmitegravir. Our results explain the mechanism of the ALLINI-induced condensation of HIV-1 IN and provide a reliable template for the rational development of this series of antiretrovirals through the optimization of their key contacts with the viral target. IMPORTANCE Despite the remarkable success of combination antiretroviral therapy, HIV-1 remains among the major causes of human suffering and loss of life in poor and developing nations. To prevail in this drawn-out battle with the pandemic, it is essential to continue developing advanced antiviral agents to fight drug resistant HIV-1 variants. Allosteric integrase inhibitors (ALLINIs) are an emerging class of HIV-1 antagonists that are orthogonal to the current antiretroviral drugs. These small molecules act as highly specific molecular glue, which triggers the aggregation of HIV-1 integrase. In this work, we present high-resolution crystal structures that reveal the crucial interactions made by two potent ALLINIs, namely, BI-D and Pirmitegravir, with HIV-1 integrase. Our results explain the mechanism of drug action and will inform the development of this promising class of small molecules for future use in antiretroviral regimens.

Exploring PROTAC Cooperativity with Coarse-Grained Alchemical Methods

Proteolysis targeting chimera (PROTAC) is a novel drug modality that facilitates the degradation of a target protein by inducing proximity with an E3 ligase. In this work, we present a new computational framework to model the cooperativity between PROTAC-E3 binding and PROTAC-target binding principally through protein-protein interactions (PPIs) induced by the PROTAC. Due to the scarcity and low resolution of experimental measurements, the physical and chemical drivers of these non-native PPIs remain to be elucidated. We develop a coarse-grained (CG) approach to model interactions in the target-PROTAC-E3 complexes, which enables converged thermodynamic estimations using alchemical free energy calculation methods despite an unconventional scale of perturbations. With minimal parametrization, we successfully capture fundamental principles of cooperativity, including the optimality of intermediate PROTAC linker lengths that originates from configurational entropy. We qualitatively characterize the dependency of cooperativity on PROTAC linker lengths and protein charges and shapes. Minimal inclusion of sequence- and conformation-specific features in our current force field, however, limits quantitative modeling to reproduce experimental measurements, but further development of the CG model may allow for efficient computational screening to optimize PROTAC cooperativity.

Identification and Inhibition of the Druggable Allosteric Site of SARS-CoV-2 NSP10/NSP16 Methyltransferase through Computational Approaches

Since its emergence in early 2019, the respiratory infectious virus, SARS-CoV-2, has ravaged the health of millions of people globally and has affected almost every sphere of life. Many efforts are being made to combat the COVID-19 pandemic’s emerging and recurrent waves caused by its evolving and more infectious variants. As a result, novel and unexpected targets for SARS-CoV-2 have been considered for drug discovery. 2′-O-Methyltransferase (nsp10/nsp16) is a significant and appealing target in the SARS-CoV-2 life cycle because it protects viral RNA from the host degradative enzymes via a cap formation process. In this work, we propose prospective allosteric inhibitors that target the allosteric site, SARS-CoV-2 MTase. Four drug libraries containing ~119,483 compounds were screened against the allosteric site of SARS-CoV-2 MTase identified in our research. The identified best compounds exhibited robust molecular interactions and alloscore-score rankings with the allosteric site of SARS-CoV-2 MTase. Moreover, to further assess the dynamic stability of these compounds (CHEMBL2229121, ZINC000009464451, SPECS AK-91811684151, NCI-ID = 715319), a 100 ns molecular dynamics simulation, along with its holo-form, was performed to provide insights on the dynamic nature of these allosteric inhibitors at the allosteric site of the SARS-CoV-2 MTase. Additionally, investigations of MM-GBSA binding free energies revealed a good perspective for these allosteric inhibitor–enzyme complexes, indicating their robust antagonistic action on SARS-CoV-2 (nsp10/nsp16) methyltransferase. We conclude that these allosteric repressive agents should be further evaluated through investigational assessments in order to combat the proliferation of SARS-CoV-2.

Structure-based virtual screening workflow to identify antivirals targeting HIV-1 capsid

We have identified novel HIV-1 capsid inhibitors targeting the PF74 binding site. Acting as the building block of the HIV-1 capsid core, the HIV-1 capsid protein plays an important role in the viral life cycle and is an attractive target for antiviral development. A structure-based virtual screening workflow for hit identification was employed, which includes docking 1.6 million commercially-available drug-like compounds from the ZINC database to the capsid dimer, followed by applying two absolute binding free energy (ABFE) filters on the 500 top-ranked molecules from docking. The first employs the Binding Energy Distribution Analysis Method (BEDAM) in implicit solvent. The top-ranked compounds are then refined using the Double Decoupling method in explicit solvent. Both docking and BEDAM refinement were carried out on the IBM World Community Grid as part of the FightAIDS@Home project. Using this virtual screening workflow, we identified 24 molecules with calculated binding free energies between − 6 and − 12 kcal/mol. We performed thermal shift assays on these molecules to examine their potential effects on the stability of HIV-1 capsid hexamer and found that two compounds, ZINC520357473 and ZINC4119064 increased the melting point of the latter by 14.8 °C and 33 °C, respectively. These results support the conclusion that the two ZINC compounds are primary hits targeting the capsid dimer interface. Our simulations also suggest that the two hit molecules may bind at the capsid dimer interface by occupying a new sub-pocket that has not been exploited by existing CA inhibitors. The possible causes for why other top-scored compounds suggested by ABFE filters failed to show measurable activity are discussed.

Molecular Dynamics Free Energy Simulations Reveal the Mechanism for the Antiviral Resistance of the M66I HIV-1 Capsid Mutation

While drug resistance mutations can often be attributed to the loss of direct or solvent-mediated protein-ligand interactions in the drug-mutant complex, in this study we show that a resistance mutation for the picomolar HIV-1 capsid (CA)-targeting antiviral (GS-6207) is mainly due to the free energy cost of the drug-induced protein side chain reorganization in the mutant protein. Among several mutations, M66I causes the most suppression of the GS-6207 antiviral activity (up to ~84,000-fold), and only 83- and 68-fold reductions for PF74 and ZW-1261, respectively. To understand the molecular basis of this drug resistance, we conducted molecular dynamics free energy simulations to study the structures, energetics, and conformational free energy landscapes involved in the inhibitors binding at the interface of two CA monomers. To minimize the protein-ligand steric clash, the I66 side chain in the M66I-GS-6207 complex switches to a higher free energy conformation from the one adopted in the apo M66I. In contrast, the binding of GS-6207 to the wild-type CA does not lead to any significant M66 conformational change. Based on an analysis that decomposes the absolute binding free energy into contributions from two receptor conformational states, it appears that it is the free energy cost of side chain reorganization rather than the reduced protein-ligand interaction that is largely responsible for the drug resistance against GS-6207.