The N501Y and K417N mutations in spike protein of SARS-CoV-2 and their combination arise questions but the data about their mechanism of action at molecular level is limited. Here, we present Free energy perturbation (FEP) calculations for the interactions of the spike S1 receptor binding domain (RBD) with both the ACE2 receptor and an antibody derived from COVID-19 patients. Our results shown that the S1 RBD-ACE2 interactions were significantly increased whereas those with the STE90-C11 antibody dramatically decreased; about over 100 times. The K417N mutation had much more pronounced effect and in a combination with N501Y fully abolished the antibody effect. This may explain the observed in UK and South Africa more spread of the virus but also raise an important question about the possible human immune response and the success of already available vaccines.

The N501Y and K417N mutations in the spike protein of SARS-CoV-2 and their combination gave rise to questions, but the data on their mechanism of action at the molecular level were limited. In this study, we present free energy perturbation (FEP) calculations, performed at the end of December 2020, for the interactions of the spike S1 receptor-binding domain (RBD) with both the ACE2 receptor and an antibody derived from COVID-19 patients. Our results showed that the S1 RBD-ACE2 interactions were significantly increased whereas those with the STE90-C11 antibody dramatically decreased. The K417N mutation in a combination with N501Y fully abolished the antibody effect. However, Lys417Asn seems to have a compensatory mechanism of action increasing the S1 RBD-ACE2 free energy of binding. This may explain the increased spread of the virus observed in the U.K. and South Africa and also gives rise to an important question regarding the possible human immune response and the success of the already available vaccines. Notably, when the experimental data became available confirming our calculations, it was demonstrated that protein−protein FEP can be a useful tool for providing urgent data to the scientific community.

Recent improvements to the free energy perturbation (FEP) calculations, especially FEP+ , established their utility for pharmaceutical lead optimization. Herein, we propose a modified version of the FEP/REST (i.e., replica exchange with solute tempering) sampling protocol, based on detail studies on several targets by probing a large number of perturbations with different sampling schemes. Improved FEP+ binding affinity predictions for regular flexible-loop motions and considerable structural changes can be obtained by extending the prior to REST (pre-REST) sampling time from 0.24 ns/λ to 5 ns/λ and 2 × 10 ns/λ, respectively. With this new protocol, much more precise ∆∆G values of the individual perturbations, including the sign of the transformations and decreased error were obtained. We extended the REST simulations from 5 ns to 8 ns to achieve reasonable free energy convergence. Implementing REST to the entire ligand as opposed to solely the perturbed region, and also some important flexible protein residues (pREST region) in the ligand binding domain (LBD) has considerably improved the FEP+ results in most of the studied cases. Preliminary molecular dynamics (MD) runs were useful for establishing the correct binding mode of the compounds and thus precise alignment for FEP+ . Our improved protocol may further increase the FEP+ accuracy.

No abstract

In recent years, the nuclear receptors (NR) dynamics have been studied extensively by various approaches. However, the transition path of helix 12 (H12) to an agonist or an antagonist conformation and the exchange pathway between these states is not clear yet. A number of accelerated molecular dynamics (aMD) runs were performed on both an ERα monomer and a homodimer with a total length of 2.2 μs. We have been able to sample reasonably well the H12 conformational landscape to reproduce precisely both the agonist and the antagonist conformations, starting from an unfolded position, and to describe the transition path between them, even in the presence of an agonist ligand. These conformations were the most prevalent, suggesting that the extended H12 state is not likely to exist and that the natural ERα H12 position might exist in both the agonist and antagonist states. Remarkably, the H12 transition occurs and is regulated only in a dimer form and the proper agonist or antagonist H12 conformation can be achieved solely in one of the dimer subunits. These results clearly demonstrate that clusters of the two well-known H12 states exist by themselves in the protein free energy landscape, i.e. they are not constituted directly by the ligands, and dimerization favors the switch between them. Conversely, in a monomer, no transitions have been observed. Thus, the dimer formation helps the constitution of populations of discrete H12 conformational states and reshapes the conformational landscape. Further analyses have shown that these observations can be explained by specific interface and long range protein-protein interactions, resulting in conformational fluctuations in helices 5 and 11. Based on these results, a new ERα activation/deactivation mechanism and a sequence of binding events during receptor activity modulation have been suggested according to which ligands control the H12 conformation via alterations of the inter-dimer interactions. These findings agree with the HDX and fluorescence experiments and provide an explanation on a structural basis of these data, demonstrating that the dynamics of H12 are not altered greatly upon ligand binding and large fluctuations at the end of H11 are present.

The B-RAF kinase plays an important role both in tumor induction and maintenance in several cancers. The molecular basis of the inactive B-RAF(WT) and B-RAF(V600E) inhibition and selectivity of a series of inhibitors was examined with a combination of molecular dynamics (MD), free energy MM-PBSA and local-binding energy (LBE) approaches. The conformational stability of the unbounded kinases and in particular the processes of the B-RAF (V600E) mutant activation were analyzed. A unique salt bridge network formed mainly by the catalytic residues was identified in the unbounded B-RAFs. The reorganization of this network and the restriction of the active segment flexibility upon ligand binding inhibit both B-RAF(WT) and B-RAF (V600E), thus appearing as an important factor for ligand selectivity. A significant correlation between the binding energies of the compounds in B-RAF(WT) and their inhibition effects on B-RAF (V600E) was revealed, which can explain the low mutant selectivity observed for numerous inhibitors. Our results suggest that the interactions between the activation segment and the alpha C-helix, as well as between the residues in the salt bridge network, are the major mechanism of the B-RAF (V600E) activation. Overall data revealed the important role of Lys601 for ligand activity, selectivity and protein stabilization, proposing an explanation of the observed strong kinase activation in the K601E mutated form.

The structural and dynamical properties of the peroxisome proliferator-activated receptor γ (PPARγ) nuclear receptor have been broadly studied in its agonist state but little is known about the key features required for the receptor antagonistic activity. Here we report a series of molecular dynamics (MD) simulations in combination with free energy estimation of the recently discovered class of non-covalent PPARγ antagonists. Their binding modes and dynamical behavior are described in details. Two key interactions have been detected within the cavity between helices H3, H11 and the activation helix H12, as well as with H12. The strength of the ligand-amino acid residues interactions has been analyzed in relation to the specificity of the ligand dynamical and antagonistic features. According to our results, the PPARγ activation helix does not undergo dramatic conformational changes, as seen in other nuclear receptors, but rather perturbations that occur through a significant ligand-induced reshaping of the ligand-receptor and the receptor-coactivator binding pockets. The H12 residue Tyr473 and the charge clamp residue Glu471 play a central role for the receptor transformations. Our results also demonstrate that MD can be a helpful tool for the compound phenotype characterization (full agonists, partial agonists or antagonists) when insufficient experimental data are available.

The structural and dynamical properties of PPARγ receptor in a complex with either partial or full agonists have been intensively studied but little is known about the receptor antagonistic conformation. A composition of microsecond accelerated molecular dynamics (aMD) simulation show that like partial agonists a non-covalent PPARγ full antagonist can bind in different modes of similar population size and free energies of binding. Four different and periodically exchanging ligand conformations are detected and described. The studied antagonist interacts with different receptor substructures and affects both the co-activator and the Cdk5 phosphorylation sites and, presumably, the natural complex with the DNA. However, no significant changes in the conformational states of the activation helix 12, and in particular an antagonist orientation, have been recorded. Finally, our results show also that the aMD approach can be successfully used in recovering the possible binding modes, considering fully the receptor flexibility, and is not dependent on the starting conformation.