A computational analysis of the binding model of MDM2 with inhibitors

Binding mechanism of inhibitors to SARS-CoV-2 main protease deciphered by multiple replica molecular dynamics simulations

The outbreak caused by SARS-CoV-2 has received extensive worldwide attention. As the main protease (M^pro) in SARS-CoV-2 has no human homologues, it is feasible to reduce the possibility of targeting the host protein by accidental drugs. Thus, M^pro has been an attractive target of efficient drug design for anti-SARS-CoV-2 treatment. In this work, multiple replica molecular dynamics (MRMD) simulations, principal component analysis (PCA), free energy landscapes (FELs), and the molecular mechanics-generalized Born surface area (MM-GBSA) method were integrated together to decipher the binding mechanism of four inhibitors masitinib, O6K, FJC and GQU to M^pro. The results indicate that the binding of four inhibitors clearly affects the structural flexibility and internal dynamics of M^pro along with dihedral angle changes of key residues. The analysis of FELs unveils that the stability in the relative orientation and geometric position of inhibitors to M^pro is favorable for inhibitor binding. Residue-based free energy decomposition reveals that the inhibitor-M^pro interaction networks involving hydrogen bonding interactions and hydrophobic interactions provide significant information for the design of potent inhibitors against M^pro. The hot spot residues including H41, M49, F140, N142, G143, C145, H163, H164, M165, E166 and Q189 identified by computational alanine scanning are considered as reliable targets of clinically available inhibitors inhibiting the activities of M^pro.

RBCK1 promotes p53 degradation via ubiquitination in renal cell carcinoma

Renal cell carcinoma (RCC) accounts for approximately 3% of adult malignancies, and the incidence of RCC continues to rise worldwide. Although RCC can be treated with surgery at an early stages, the five-year survival rates have been observed to decline dramatically in patients with advanced disease. Most patients with RCC treated with cytotoxic or targeted drugs will develop resistance at some point during therapy. Thus, it is necessary to identify novel therapeutic targets for RCC. Here, we found that RANBP2-type and C3HC4-type zinc finger-containing 1 (RBCK1) expression was upregulated in human RCC samples. Analysis of multiple public databases revealed the correlation between RBCK1 expression and poor prognosis in RCC patients. Subsequently, we performed RBCK1 depletion experiments in RCC cells that severely affected the in vivo and in vitro proliferation of renal cancer cells. The effects of RBCK1 on cell proliferation could be rescued with p53 expression knockdown in two cell lines expressing wild-type p53. Further experiments demonstrated that RBCK1 could facilitate p53 poly-ubiquitination and degradation by direct interaction with p53. Together, our results show that RBCK1 may serve as a promising target for RCC therapy by restoring p53 functions.

Performance of a docking/molecular dynamics protocol for virtual screening of nutlin-class inhibitors of Mdmx

A virtual screening protocol involving docking and molecular dynamics has been tested against the results of fluorescence polarization assays testing the potency of a series of compounds of the nutlin class for inhibition of the interaction between p53 and Mdmx, an interaction identified as a driver of certain cancers. The protocol uses a standard docking method (AutoDock) with a cutoff based on the AutoDock score (ADscore), followed by molecular dynamics simulation with a cutoff based on root-mean-square-deviation (RMSD) from the docked pose. An analysis of the experimental and computational results shows modest performance of ADscore alone, but dramatically improved performance when RMSD is also used.

The conformational feasibility for the formation of reaching dimer in ASV and HIV integrase: a molecular dynamics study

Retroviral integrases are reported to form alternate dimer assemblies like the core-core dimer and reaching dimer. The core-core dimer is stabilized predominantly by an extensive interface between two catalytic core domains. The reaching dimer is stabilized by N-terminal domains that reach to form intermolecular interfaces with the other subunit's core and C-terminal domains (CTD), as well as CTD-CTD interactions. In this study, molecular dynamics (MD), Brownian dynamics (BD) simulations, and free energy analyses, were performed to elucidate determinants for the stability of the reaching dimer forms of full-length Avian Sarcoma Virus (ASV) and Human Immunodeficiency Virus (HIV) IN, and to examine the role of the C-tails (the last ~16-18 residues at the C-termini) in their structural dynamics. The dynamics of an HIV reaching dimer derived from small angle X-ray scattering and protein crosslinking data, was compared with the dynamics of a core-core dimer model derived from combining the crystal structures of two-domain fragments. The results showed that the core domains in the ASV reaching dimer express free dynamics, whereas those in the HIV reaching dimer are highly stable. BD simulations suggest a higher rate of association for the HIV core-core dimer than the reaching dimer. The predicted stability of these dimers was therefore ranked in the following order: ASV reaching dimer < HIV reaching dimer < composite core-core dimer. Analyses of MD trajectories have suggested residues that are critical for intermolecular contacts in each reaching dimer. Tests of these predictions and insights gained from these analyses could reveal a potential pathway for the association and dissociation of full-length IN multimers.

Medicinal Chemistry Strategies to Disrupt the p53-MDM2/MDMX Interaction

The growth inhibitory activity of p53 tumor suppressor is tightly regulated by interaction with two negative regulatory proteins, murine double minute 2 (MDM2) and X (MDMX), which are overexpressed in about half of all human tumors. The elucidation of crystallographic structures of MDM2/MDMX complexes with p53 has been pivotal for the identification of several classes of inhibitors of the p53-MDM2/MDMX interaction. The present review provides in silico strategies and screening approaches used in drug discovery as well as an overview of the most relevant classes of small-molecule inhibitors of the p53-MDM2/MDMX interaction, their progress in pipeline, and highlights particularities of each class of inhibitors. Most of the progress made with high-throughput screening has led to the development of inhibitors belonging to the cis-imidazoline, piperidinone, and spiro-oxindole series. However, novel potent and selective classes of inhibitors of the p53-MDM2 interaction with promising antitumor activity are emerging. Even with the discovery of the 3D structure of complex p53-MDMX, only two small molecules were reported as selective p53-MDMX antagonists, WK298 and SJ-172550. Dual inhibition of the p53-MDM2/MDMX interaction has shown to be an alternative approach since it results in full activation of the p53-dependent pathway. The knowledge of structural requirements crucial to the development of small-molecule inhibitors of the p53-MDMs interactions has enabled the identification of novel antitumor agents with improved in vivo efficacy.

Energetic Landscape of MDM2-p53 Interactions by Computational Mutagenesis of the MDM2-p53 Interaction

The ubiquitin ligase MDM2, a principle regulator of the tumor suppressor p53, plays an integral role in regulating cellular levels of p53 and thus a prominent role in current cancer research. Computational analysis used MUMBO to rotamerize the MDM2-p53 crystal structure 1YCR to obtain an exhaustive search of point mutations, resulting in the calculation of the ΔΔG comprehensive energy landscape for the p53-bound regulator. The results herein have revealed a set of residues R65-E69 on MDM2 proximal to the p53 hydrophobic binding pocket that exhibited an energetic profile deviating significantly from similar residues elsewhere in the protein. In light of the continued search for novel competitive inhibitors for MDM2, we discuss possible implications of our findings on the drug discovery field.

Probing Origin of Binding Difference of inhibitors to MDM2 and MDMX by Polarizable Molecular Dynamics Simulation and QM/MM-GBSA Calculation

Binding abilities of current inhibitors to MDMX are weaker than to MDM2. Polarizable molecular dynamics simulations (MD) followed by Quantum mechanics/molecular mechanics generalized Born surface area (QM//MM-GBSA) calculations were performed to investigate the binding difference of inhibitors to MDM2 and MDMX. The predicted binding free energies not only agree well with the experimental results, but also show that the decrease in van der Walls interactions of inhibitors with MDMX relative to MDM2 is a main factor of weaker bindings of inhibitors to MDMX. The analyses of dihedral angles based on MD trajectories suggest that the closed conformation formed by the residues M53 and Y99 in MDMX leads to a potential steric clash with inhibitors and prevents inhibitors from arriving in the deep of MDMX binding cleft, which reduces the van der Waals contacts of inhibitors with M53, V92, P95 and L98. The calculated results using the residue-based free energy decomposition method further prove that the interaction strength of inhibitors with M53, V92, P95 and L98 from MDMX are obviously reduced compared to MDM2. We expect that this study can provide significant theoretical guidance for designs of potent dual inhibitors to block the p53-MDM2/MDMX interactions.

Probing Difference in Binding Modes of Inhibitors to MDMX by Molecular Dynamics Simulations and Different Free Energy Methods

The p53-MDMX interaction has attracted extensive attention of anti-cancer drug development in recent years. This current work adopted molecular dynamics (MD) simulations and cross-correlation analysis to investigate conformation changes of MDMX caused by inhibitor bindings. The obtained information indicates that the binding cleft of MDMX undergoes a large conformational change and the dynamic behavior of residues obviously change by the presence of different structural inhibitors. Two different methods of binding free energy predictions were employed to carry out a comparable insight into binding mechanisms of four inhibitors PMI, pDI, WK23 and WW8 to MDMX. The data show that the main factor controlling the inhibitor bindings to MDMX arises from van der Waals interactions. The binding free energies were further divided into contribution of each residue and the derived information gives a conclusion that the hydrophobic interactions, such as CH-CH, CH-π and π-π interactions, are responsible for the inhibitor associations with MDMX.

Characterizing the Free-Energy Landscape of MDM2 Protein-Ligand Interactions by Steered Molecular Dynamics Simulations

Inhibition of p53-MDM2 interaction by small molecules is considered to be a promising approach to re-activate wild-type p53 for tumor suppression. Several inhibitors of the MDM2-p53 interaction were designed and studied by the experimental methods and the molecular dynamics simulation. However, the unbinding mechanism was still unclear. The steered molecular dynamics simulations combined with Brownian dynamics fluctuation-dissipation theorem were employed to obtain the free-energy landscape of unbinding between MDM2 and their four ligands. It was shown that compounds 4 and 8 dissociate faster than compounds 5 and 7. The absolute binding free energies for these four ligands are in close agreement with experimental results. The open movement of helix II and helix IV in the MDM2 protein-binding pocket upon unbinding is also consistent with experimental MDM2-unbound conformation. We further found that different binding mechanisms among different ligands are associated with H-bond with Lys51 and Glu25. These mechanistic results may be useful for improving ligand design.

Binding modes of three inhibitors 8CA, F8A and I4A to A-FABP studied based on molecular dynamics simulation

Adipocyte fatty-acid binding protein (A-FABP) is an important target of drug designs treating some diseases related to lipid-mediated biology. Molecular dynamics (MD) simulations coupled with solvated interaction energy method (SIE) were carried out to study the binding modes of three inhibitors 8CA, F8A and I4A to A-FABP. The rank of our predicted binding affinities is in accordance with experimental data. The results show that the substitution in the position 5 of N-benzyl and the seven-membered ring of N-benzyl-indole carboxylic acids strengthen the I4A binding, while the substitution in the position 2 of N-benzyl weakens the F8A binding. Computational alanine scanning and dynamics analyses were performed and the results suggest that the polar interactions of the positively charged residue R126 with the three inhibitors provide a significant contribution to inhibitor bindings. This polar interaction induces the disappearance of the correlated motion of the C terminus of A-FABP relative to the N terminus and favors the stability of the binding complex. This study is helpful for the rational design of potent inhibitors within the fields of metabolic disease, inflammation and atherosclerosis.

Practical Aspects of Free-Energy Calculations: A Review

Free-energy calculations in the framework of classical molecular dynamics simulations are nowadays used in a wide range of research areas including solvation thermodynamics, molecular recognition, and protein folding. The basic components of a free-energy calculation, that is, a suitable model Hamiltonian, a sampling protocol, and an estimator for the free energy, are independent of the specific application. However, the attention that one has to pay to these components depends considerably on the specific application. Here, we review six different areas of application and discuss the relative importance of the three main components to provide the reader with an organigram and to make nonexperts aware of the many pitfalls present in free energy calculations.

A computational analysis of binding modes and conformation changes of MDM2 induced by p53 and inhibitor bindings

Molecular dynamics (MD) simulations followed by principal component analysis were performed to study the conformational change of MDM2 induced by p53 and two inhibitor (P4 and MI63a) bindings. The results show that the hydrophobic cleft of MDM2 is very flexible and adaptive to different structural binding partners. The cleft tends to become wider and more stable as MDM2 binds to the three binding partners, while unbound MDM2 shows a narrower and pretty flexible cleft, which agrees with recent experimental data and theoretical studies. It was also found that the binding of P4 and p53 stabilizes the motion of the loop L2 linking the helix α2 and β strand (β3), but the presence of MI63a makes the motion of L2 disordered. In addition, the binding free energies of the three partners to MDM2 were calculated using molecular mechanics generalized Born surface area to explain the binding modes of these three partners to MDM2. This study will be helpful not only for better understanding the functional, concerted motion of MDM2, but also for the rational design of potent anticancer drugs targeting the p53-MDM2 interaction.

Pharmacophore based 3D-QSAR modeling and free energy analysis of VEGFR-2 inhibitors

VEGFR-2, a transmembrane tyrosine kinase receptor is responsible for angiogenesis and has been an attractive target in treating cancers. The inhibition mechanism of structurally diverse urea derivatives, reported as VEGFR-2 inhibitors, was explored by pharmacophore modeling, QSAR, and molecular dynamics based free energy analysis.The pharmacophore hypothesis AADRR, resulted in a highly significant atom based 3D-QSAR model (r(2) = 0.94 and q(2) = 0.84). Binding free energy analysis of the docked complexes of highly active and inactive compounds, after 7 ns MD simulation, revealed the importance of van der Waals interaction in VEGFR-2 inhibition. The decomposition of binding free energy on a per residue basis disclosed that the residues in hinge region and hydrophobic pocket play a role in discriminating the active and inactive inhibitors. Thus, the present study proposes a pharmacophore hypothesis representing the identified interactions pattern and its further application as a template in screening databases to identify novel VEGFR-2 inhibitor scaffolds.

Molecular dynamic simulation insights into the normal state and restoration of p53 function

As a tumor suppressor protein, p53 plays a crucial role in the cell cycle and in cancer prevention. Almost 50 percent of all human malignant tumors are closely related to a deletion or mutation in p53. The activity of p53 is inhibited by over-active celluar antagonists, especially by the over-expression of the negative regulators MDM2 and MDMX. Protein-protein interactions, or post-translational modifications of the C-terminal negative regulatory domain of p53, also regulate its tumor suppressor activity. Restoration of p53 function through peptide and small molecular inhibitors has become a promising strategy for novel anti-cancer drug design and development. Molecular dynamics simulations have been extensively applied to investigate the conformation changes of p53 induced by protein-protein interactions and protein-ligand interactions, including peptide and small molecular inhibitors. This review focuses on the latest MD simulation research, to provide an overview of the current understanding of interactions between p53 and its partners at an atomic level.

Simulating molecular mechanisms of the MDM2-mediated regulatory interactions: a conformational selection model of the MDM2 lid dynamics

Diversity and complexity of MDM2 mechanisms govern its principal function as the cellular antagonist of the p53 tumor suppressor. Structural and biophysical studies have demonstrated that MDM2 binding could be regulated by the dynamics of a pseudo-substrate lid motif. However, these experiments and subsequent computational studies have produced conflicting mechanistic models of MDM2 function and dynamics. We propose a unifying conformational selection model that can reconcile experimental findings and reveal a fundamental role of the lid as a dynamic regulator of MDM2-mediated binding. In this work, structure, dynamics and energetics of apo-MDM2 are studied as a function of posttranslational modifications and length of the lid. We found that the dynamic equilibrium between "closed" and "semi-closed" lid forms may be a fundamental characteristic of MDM2 regulatory interactions, which can be modulated by phosphorylation, phosphomimetic mutation as well as by the lid size. Our results revealed that these factors may regulate p53-MDM2 binding by fine-tuning the thermodynamic equilibrium between preexisting conformational states of apo-MDM2. In agreement with NMR studies, the effect of phosphorylation on MDM2 interactions was more pronounced with the truncated lid variant that favored the thermodynamically dominant closed form. The phosphomimetic mutation S17D may alter the lid dynamics by shifting the thermodynamic equilibrium towards the ensemble of "semi-closed" conformations. The dominant "semi-closed" lid form and weakened dependence on the phosphorylation seen in simulations with the complete lid can provide a rationale for binding of small p53-based mimetics and inhibitors without a direct competition with the lid dynamics. The results suggested that a conformational selection model of preexisting MDM2 states may provide a robust theoretical framework for understanding MDM2 dynamics. Probing biological functions and mechanisms of MDM2 regulation would require further integration of computational and experimental studies and may help to guide drug design of novel anti-cancer therapeutics.

Computational studies of difference in binding modes of peptide and non-peptide inhibitors to MDM2/MDMX based on molecular dynamics simulations

Inhibition of p53-MDM2/MDMX interaction is considered to be a promising strategy for anticancer drug design to activate wild-type p53 in tumors. We carry out molecular dynamics (MD) simulations to study the binding mechanisms of peptide and non-peptide inhibitors to MDM2/MDMX. The rank of binding free energies calculated by molecular mechanics generalized Born surface area (MM-GBSA) method agrees with one of the experimental values. The results suggest that van der Waals energy drives two kinds of inhibitors to MDM2/MDMX. We also find that the peptide inhibitors can produce more interaction contacts with MDM2/MDMX than the non-peptide inhibitors. Binding mode predictions based on the inhibitor-residue interactions show that the π–π, CH–π and CH–CH interactions dominated by shape complimentarity, govern the binding of the inhibitors in the hydrophobic cleft of MDM2/MDMX. Our studies confirm the residue Tyr99 in MDMX can generate a steric clash with the inhibitors due to energy and structure. This finding may theoretically provide help to develop potent dual-specific or MDMX inhibitors.

Hot spots and transient pockets: predicting the determinants of small-molecule binding to a protein-protein interface

Protein-protein interfaces are considered difficult targets for small-molecule protein-protein interaction modulators (PPIMs ). Here, we present for the first time a computational strategy that simultaneously considers aspects of energetics and plasticity in the context of PPIM binding to a protein interface. The strategy aims at identifying the determinants of small-molecule binding, hot spots, and transient pockets, in a protein-protein interface in order to make use of this knowledge for predicting binding modes of and ranking PPIMs with respect to their affinity. When applied to interleukin-2 (IL-2), the computationally inexpensive constrained geometric simulation method FRODA outperforms molecular dynamics simulations in sampling hydrophobic transient pockets. We introduce the PPIAnalyzer approach for identifying transient pockets on the basis of geometrical criteria only. A sequence of docking to identified transient pockets, starting structure selection based on hot spot information, RMSD clustering and intermolecular docking energies, and MM-PBSA calculations allows one to enrich IL-2 PPIMs from a set of decoys and to discriminate between subgroups of IL-2 PPIMs with low and high affinity. Our strategy will be applicable in a prospective manner where nothing else than a protein-protein complex structure is known; hence, it can well be the first step in a structure-based endeavor to identify PPIMs.

Insights into scFv:drug binding using the molecular dynamics simulation and free energy calculation

Molecular dynamics simulations and free energy calculation have been performed to study how the single-chain variable fragment (scFv) binds methamphetamine (METH) and amphetamine (AMP). The structures of the scFv:METH and the scFv:AMP complexes are analyzed by examining the time-dependence of their RMSDs, by analyzing the distance between some key atoms of the selected residues, and by comparing the averaged structures with their corresponding crystallographic structures. It is observed that binding an AMP to the scFv does not cause significant changes to the binding pocket of the scFv:ligand complex. The binding free energy of scFv:AMP without introducing an extra water into the binding pocket is much stronger than scFv:METH. This is against the first of the two scenarios postulated in the experimental work of Celikel et al. (Protein Science 18, 2336 (2009)). However, adding a water to the AMP (at the position of the methyl group of METH), the binding free energy of the scFv:AMP-H2O complex, is found to be significantly weaker than scFv:METH. This is consistent with the second of the two scenarios given by Celikel et al. Decomposition of the binding energy into ligand-residue pair interactions shows that two residues (Tyr175 and Tyr177) have nearly-zero interactions with AMP in the scFv:AMP-H2O complex, whereas their interactions with METH in the scFv:METH complex are as large as -0.8 and -0.74 kcal mol(-1). The insights gained from this study may be helpful in designing more potent antibodies in treating METH abuse.

Molecular modeling and molecular dynamics simulation studies on pyrrolopyrimidine-based α-helix mimetic as dual inhibitors of MDM2 and MDMX

Inhibition of the interactions between the tumor suppressor protein p53 and its negative regulators, the MDM2 and MDMX oncogenic proteins, is increasingly gaining interest in cancer therapy and drug design. In this study, we carry out molecular docking, molecular dynamics (MD) simulations, and molecular mechanics Poisson-Boltzmann and generalized Born/surface area (MM-PB/GBSA) binding free energy calculations on an active compound 3a and an inactive compound NC-1, which share a common pyrrolopyrimidine-based scaffold. MD simulations and MM-PB/GBSA calculations show that the compound NC-1 may not bind to MDM2 and MDMX, in agreement with the experimental results. Detailed MM-PB/GBSA calculations on the MDM2-3a and MDMX-3a complexes unravel that the binding free energies are similar for the two complexes. Furthermore, the van der Waals energy is the largest component of the binding free energy for both complexes, which indicates that the interactions between the compound 3a and MDM2 and MDMX are dominated by shape complementarity. In addition, the analysis of individual residue contribution and protein-ligand binding mode show that the three functional groups on R₁, R₂, and R₃ of the compound 3a can mimic the spatial orientation of the side chains of Phe19, Trp23, and Leu26 of p53, respectively. The obtained computational results suggest that the compound 3a can act as a dual inhibitor of MDM2-p53 and MDMX-p53 interactions, consistent with the experimental results.

Free energy calculations of protein-ligand interactions

In the calculation of free energies of binding for protein-ligand complexes, we distinguish endpoint methods, methods involving alchemical modifications and methods that physically displace the ligand from the protein. Most methodological advances seem to come from a clever combination of multiple existing methods to enhance the sampling or to utilize specific advantages of various approaches. The coupling parameters common in thermodynamic integration and in Hamiltonian replica exchange are for instance combined to yield replica exchange thermodynamic integration. As new methods mostly aim to improve efficiency or to attain more complete sampling, there are good prospects to understand and tackle the sampling problem better and to shift the focus towards the scoring problem in the context of more robust and accurate force fields.