Structural Fluctuations of Aromatic Residues in an Apo-Form Reveal Cryptic Binding Sites: Implications for Fragment-Based Drug Design

Lipid Trolling to Optimize A3 Adenosine Receptor-Positive Allosteric Modulators (PAMs)

A3 adenosine receptor (A3AR) positive allosteric modulators (PAMs) (2,4-disubstituted-1H-imidazo[4,5-c]quinolin-4-amines) allosterically increase the Emax of A3AR agonists, but not potency, due to concurrent orthosteric antagonism. Following mutagenesis/homology modeling of the proposed lipid-exposed allosteric binding site on the cytosolic side, we functionalized the scaffold, including heteroatom substitutions and exocyclic phenylamine extensions, to increase allosteric binding. Strategically appended linear alkyl-alkynyl chains with terminal amino/guanidino groups improved allosteric effects at both human and mouse A3ARs. The chain length, functionality, and attachment position were varied to modulate A3AR PAM activity. For example, 26 (MRS8247, p-alkyne-linked 8 methylenes) and homologues increased agonist Cl-IB-MECA's Emax and potency ([35S]GTPγS binding). The putative mechanism involves a flexible, terminally cationic chain penetrating the lipid environment for stable electrostatic anchoring to cytosolic phospholipid head groups, suggesting "lipid trolling", supported by molecular dynamic simulation of the active-state model. Thus, we have improved A3AR PAM activity through rational design based on an extrahelical, lipidic binding site.

Binding Mechanism of Riboswitch to Natural Ligand Elucidated by McMD-Based Dynamic Docking Simulations

Flavin mononucleotide riboswitches are common among many pathogenic bacteria and are therefore considered to be an attractive target for antibiotics development. The riboswitch binds riboflavin (RBF, also known as vitamin B2), and although an experimental structure of their complex has been solved with the ligand bound deep inside the RNA molecule in a seemingly unreachable state, the binding mechanism between these molecules is not yet known. We have therefore used our Multicanonical Molecular Dynamics (McMD)-based dynamic docking protocol to analyze their binding mechanism by simulating the binding process between the riboswitch aptamer domain and the RBF, starting from the apo state of the riboswitch. Here, the refinement stage was crucial to identify the native binding configuration, as several other binding configurations were also found by McMD-based docking simulations. RBF initially binds the interface between P4 and P6 including U61 and G62, which forms a gateway where the ligand lingers until this gateway opens sufficiently to allow the ligand to pass through and slip into the hidden binding site including A48, A49, and A85.

Grasping cryptic binding sites to neutralize drug resistance in the field of anticancer

Drug resistance is a significant obstacle to successful cancer treatment. The utilization and development of cryptic binding sites (CBSs) in proteins involved in cancer-related drug-resistance (CRDR) could help to overcome that drug resistance. However, there is no comprehensive review of the successful use of CBSs in addressing CRDR. Here, we have systematically summarized and analyzed the opportunities and challenges of using CBSs in addressing CRDR and revealed the key role that CBSs have in targeting CRDR. First, we have identified the CRDR targets and the corresponding CBSs. Second, we discuss the mechanisms by which CBSs can overcome CRDR. Finally, we have provided examples of successful CBS applications in addressing CRDR. We hope that this approach will provide guidance to biologists and chemists in effectively utilizing CBSs for the development of new drugs to alleviate CRDR.

Mutual induced-fit mechanism drives binding between intrinsically disordered Bim and cryptic binding site of Bcl-xL

The intrinsically disordered region (IDR) of Bim binds to the flexible cryptic site of Bcl-xL, a pro-survival protein involved in cancer progression that plays an important role in initiating apoptosis. However, their binding mechanism has not yet been elucidated. We have applied our dynamic docking protocol, which correctly reproduced both the IDR properties of Bim and the native bound configuration, as well as suggesting other stable/meta-stable binding configurations and revealed the binding pathway. Although the cryptic site of Bcl-xL is predominantly in a closed conformation, initial binding of Bim in an encounter configuration leads to mutual induced-fit binding, where both molecules adapt to each other; Bcl-xL transitions to an open state as Bim folds from a disordered to an α-helical conformation while the two molecules bind each other. Finally, our data provides new avenues to develop novel drugs by targeting newly discovered stable conformations of Bcl-xL.

Free energy and kinetic rate calculation via non-equilibrium molecular simulation: application to biomolecules

Non-equilibrium molecular dynamics (NEMD) simulation has been recognized as a powerful tool for examining biomolecules and provides fruitful insights into not only non-equilibrium but also equilibrium processes. We review recent advances in NEMD simulation and relevant, fundamental results of non-equilibrium statistical mechanics. We first introduce Crooks fluctuation theorem and Jarzynski equality that relate free energy difference to work done on a physical system during a non-equilibrium process. The theorems are beneficial for the analysis of NEMD trajectories. We then describe rate theory, a framework to calculate molecular kinetics from a non-equilibrium process; this theoretical framework enables us to calculate a reaction time-mean-first passage time-from NEMD trajectories. We, in turn, present recent NEMD techniques that apply an external force to a system to enhance molecular dissociation and introduce their application to biomolecules. Lastly, we show the current status of an appropriate selection of reaction coordinates for NEMD simulation.

Computer simulation of molecular recognition in biomolecular system: from in silico screening to generalized ensembles

Prediction of ligand-receptor complex structure is important in both the basic science and the industry such as drug discovery. We report various computation molecular docking methods: fundamental in silico (virtual) screening, ensemble docking, enhanced sampling (generalized ensemble) methods, and other methods to improve the accuracy of the complex structure. We explain not only the merits of these methods but also their limits of application and discuss some interaction terms which are not considered in the in silico methods. In silico screening and ensemble docking are useful when one focuses on obtaining the native complex structure (the most thermodynamically stable complex). Generalized ensemble method provides a free-energy landscape, which shows the distribution of the most stable complex structure and semi-stable ones in a conformational space. Also, barriers separating those stable structures are identified. A researcher should select one of the methods according to the research aim and depending on complexity of the molecular system to be studied.

Molecular Structure of the Surface-Immobilized Super Uranyl Binding Protein

Recently, a super uranyl binding protein (SUP) was developed, which exhibits excellent sensitivity/selectivity to bind uranyl ions. It can be immobilized onto a surface in sensing devices to detect uranyl ions. Here, sum frequency generation (SFG) vibrational spectroscopy was applied to probe the interfacial structures of surface-immobilized SUP. The collected SFG spectra were compared to the calculated orientation-dependent SUP SFG spectra using a one-excitonic Hamiltonian approach based on the SUP crystal structures to deduce the most likely surface-immobilized SUP orientation(s). Furthermore, discrete molecular dynamics (DMD) simulation was applied to refine the surface-immobilized SUP conformations and orientations. The immobilized SUP structures calculated from DMD simulations confirmed the SUP orientations obtained from SFG data analyzed based on the crystal structures and were then used for a new round of SFG orientation analysis to more accurately determine the interfacial orientations and conformations of immobilized SUP before and after uranyl ion binding, providing an in-depth understanding of molecular interactions between SUP and the surface and the effect of uranyl ion binding on the SUP interfacial structures. We believe that the developed method of combining SFG measurements, DMD simulation, and Hamiltonian data analysis approach is widely applicable to study biomolecules at solid/liquid interfaces.

Cryptic-site binding mechanism of medium-sized Bcl-xL inhibiting compounds elucidated by McMD-based dynamic docking simulations

We have performed multicanonical molecular dynamics (McMD) based dynamic docking simulations to study and compare the binding mechanism between two medium-sized inhibitors (ABT-737 and WEHI-539) that bind to the cryptic site of Bcl-xL, by exhaustively sampling the conformational and configurational space. Cryptic sites are binding pockets that are transiently formed in the apo state or are induced upon ligand binding. Bcl-xL, a pro-survival protein involved in cancer progression, is known to have a cryptic site, whereby the shape of the pocket depends on which ligand is bound to it. Starting from the apo-structure, we have performed two independent McMD-based dynamic docking simulations for each ligand, and were able to obtain near-native complex structures in both cases. In addition, we have also studied their interactions along their respective binding pathways by using path sampling simulations, which showed that the ligands form stable binding configurations via predominantly hydrophobic interactions. Although the protein started from the apo state, both ligands modulated the pocket in different ways, shifting the conformational preference of the sub-pockets of Bcl-xL. We demonstrate that McMD-based dynamic docking is a powerful tool that can be effectively used to study binding mechanisms involving a cryptic site, where ligand binding requires a large conformational change in the protein to occur.

Identification of Cryptic Binding Sites Using MixMD with Standard and Accelerated Molecular Dynamics

Protein dynamics play an important role in small molecule binding and can pose a significant challenge in the identification of potential binding sites. Cryptic binding sites have been defined as sites which require significant rearrangement of the protein structure to become physically accessible to a ligand. Mixed-solvent MD (MixMD) is a computational protocol which maps the surface of the protein using molecular dynamics (MD) of the unbound protein solvated in a 5% box of probe molecules with explicit water. This method has successfully identified known active and allosteric sites which did not require reorganization. In this study, we apply the MixMD protocol to identify known cryptic sites of 12 proteins characterized by a wide range of conformational changes. Of these 12 proteins, three require reorganization of side chains, five require loop movements, and four require movement of more significant structures such as whole helices. In five cases, we find that standard MixMD simulations are able to map the cryptic binding sites with at least one probe type. In two cases (guanylate kinase and TIE-2), accelerated MD, which increases sampling of torsional angles, was necessary to achieve mapping of portions of the cryptic binding site missed by standard MixMD. For more complex systems where movement of a helix or domain is necessary, MixMD was unable to map the binding site even with accelerated dynamics, possibly due to the limited timescale (100 ns for individual simulations). In general, similar conformational dynamics are observed in water-only simulations and those with probe molecules. This could imply that the probes are not driving opening events but rather take advantage of mapping sites that spontaneously open as part of their inherent conformational behavior. Finally, we show that docking to an ensemble of conformations from the standard MixMD simulations performs better than docking the apo crystal structure in nine cases and even better than half of the bound crystal structures. Poorer performance was seen in docking to ensembles of conformations from the accelerated MixMD simulations.

Druggable Transient Pockets in Protein Kinases

Drug discovery using small molecule inhibitors is reaching a stalemate due to low selectivity, adverse off-target effects and inevitable failures in clinical trials. Conventional chemical screening methods may miss potent small molecules because of their use of simple but outdated kits composed of recombinant enzyme proteins. Non-canonical inhibitors targeting a hidden pocket in a protein have received considerable research attention. Kii and colleagues identified an inhibitor targeting a transient pocket in the kinase DYRK1A during its folding process and termed it FINDY. FINDY exhibits a unique inhibitory profile; that is, FINDY does not inhibit the fully folded form of DYRK1A, indicating that the FINDY-binding pocket is hidden in the folded form. This intriguing pocket opens during the folding process and then closes upon completion of folding. In this review, we discuss previously established kinase inhibitors and their inhibitory mechanisms in comparison with FINDY. We also compare the inhibitory mechanisms with the growing concept of “cryptic inhibitor-binding sites.” These sites are buried on the inhibitor-unbound surface but become apparent when the inhibitor is bound. In addition, an alternative method based on cell-free protein synthesis of protein kinases may allow the discovery of small molecules that occupy these mysterious binding sites. Transitional folding intermediates would become alternative targets in drug discovery, enabling the efficient development of potent kinase inhibitors.

Asymmetric dynamics of dimeric SARS-CoV-2 and SARS-CoV main proteases in an apo form: Molecular dynamics study on fluctuations of active site, catalytic dyad, and hydration water

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has been widely spread around the world. It is necessary to examine the viral proteins that play a notorious role in the invasion of our body. The main protease (3CL^pro) facilitates the maturation of the coronavirus. It is thought that the dimerization of 3CL^pro leads to its catalytic activity; the detailed mechanism has, however, not been suggested. Furthermore, the structural differences between the predecessor SARS-CoV 3CL^pro and SARS-CoV-2 3CL^pro have not been fully understood. Here, we show the structural and dynamical differences between the two main proteases, and demonstrate the relationship between the dimerization and the activity via atomistic molecular dynamics simulations. Simulating monomeric and dimeric 3CL^pro systems for each protease, we show that (i) global dynamics between the two different proteases are not conserved, (ii) the dimerization stabilizes the catalytic dyad and hydration water molecules behind the dyad, and (iii) the substrate-binding site (active site) and hydration water molecules in each protomer fluctuate asymmetrically. We then speculate the roles of hydration water molecules in their catalytic activity.