Factors Governing Selectivity of Dopamine Receptor Binding Compounds for D2R and D3R Subtypes.

Targeting the D3 dopamine receptor (D3R) is a promising pharmacotherapeutic strategy for the treatment of many disorders. The structure of the D3R is similar to the D2 dopamine receptor (D2R), especially in the transmembrane spanning regions that form the orthosteric binding site, making it difficult to identify D3R selective pharmacotherapeutic agents. Here, we examine the molecular basis for the high affinity D3R binding and D3R vs D2R binding selectivity of substituted phenylpiperazine thiopheneamides. We show that removing the thiophenearylamide portion of the ligand consistently decreases the affinity of these ligands at D3R, while not affecting their affinity at the D2R. Our long (>10 µs) molecular dynamics simulations demonstrated that both dopamine receptor subtypes adopt two major conformations that we refer to as closed or open conformations, with D3R sampling the open conformation more frequently than D2R. The binding of ligands with conjoined orthosteric-allosteric binding moieties causes the closed conformation to populate more often in the trajectories. Also, significant differences were observed in the extracellular loops (ECL) of these two receptor subtypes leading to the identification of several residues that contribute differently to the ligand binding for the two receptors that could potentially contribute to ligand binding selectivity. Our observations also suggest that the displacement of ordered water in the binding pocket of D3R contributes to the affinity of the compounds containing an allosteric binding motif. These studies provide a better understanding of how a bitopic mode of engagement can determine ligands that bind selectively to D2 and D3 dopamine receptor subtypes.

Probing Protein Allostery as a Residue-Specific Concept via Residue Response Maps.

Allosteric regulation is a well-established phenomenon defined as a distal conformational or dynamical change of the protein upon allosteric effector binding. Here, we developed a novel approach to delineate allosteric effects in proteins. In this approach, we applied robust machine learning methods, including deep neural network and random forest, on extensive molecular dynamics (MD) simulations to distinguish otherwise similar allosteric states of proteins. Using the PDZ3 domain of PDS-95 as a model protein, we demonstrated that the allosteric effects could be represented as residue-specific properties through two-dimensional property-residue maps, which we refer to as "residue response maps". These maps were constructed through two machine learning methods and could accurately describe how different properties of various residues are affected upon allosteric perturbation on protein. Based on the "residue response maps", we propose allostery as a residue-specific concept, suggesting that all residues could be considered as allosteric residues because each residue "senses" the allosteric events through changing its single or multiple attributes in a quantitatively unique way. The "residue response maps" could be used to fingerprint a protein based on the unique patterns of residue responses upon binding events, providing a novel way to systematically describe the protein allosteric effects of each residue upon perturbation.

A Genetically Encoded, Phage-Displayed Cyclic-Peptide Library.

Superior to linear peptides in biological activities, cyclic peptides are considered to have great potential as therapeutic agents. To identify cyclic-peptide ligands for therapeutic targets, phage-displayed peptide libraries in which cyclization is achieved by the covalent conjugation of cysteines have been widely used. To resolve drawbacks related to cysteine conjugation, we have invented a phage-display technique in which its displayed peptides are cyclized through a proximity-driven Michael addition reaction between a cysteine and an amber-codon-encoded Nϵ -acryloyl-lysine (AcrK). Using a randomized 6-mer library in which peptides were cyclized at two ends through a cysteine-AcrK linker, we demonstrated the successful selection of potent ligands for TEV protease and HDAC8. All selected cyclic peptide ligands showed 4- to 6-fold stronger affinity to their protein targets than their linear counterparts. We believe this approach will find broad applications in drug discovery.

Stress and interferon signalling-mediated apoptosis contributes to pleiotropic anticancer responses induced by targeting NGLY1.

BACKGROUND: Although NGLY1 is known as a pivotal enzyme that catalyses the deglycosylation of denatured glycoproteins, information regarding the responses of human cancer and normal cells to NGLY1 suppression is limited. METHODS: We examined how NGLY1 expression affects viability, tumour growth, and responses to therapeutic agents in melanoma cells and an animal model. Molecular mechanisms contributing to NGLY1 suppression-induced anticancer responses were revealed by systems biology and chemical biology studies. Using computational and medicinal chemistry-assisted approaches, we established novel NGLY1-inhibitory small molecules. RESULTS: Compared with normal cells, NGLY1 was upregulated in melanoma cell lines and patient tumours. NGLY1 knockdown caused melanoma cell death and tumour growth retardation. Targeting NGLY1 induced pleiotropic responses, predominantly stress signalling-associated apoptosis and cytokine surges, which synergise with the anti-melanoma activity of chemotherapy and targeted therapy agents. Pharmacological and molecular biology tools that inactivate NGLY1 elicited highly similar responses in melanoma cells. Unlike normal cells, melanoma cells presented distinct responses and high vulnerability to NGLY1 suppression. CONCLUSION: Our work demonstrated the significance of NGLY1 in melanoma cells, provided mechanistic insights into how NGLY1 inactivation leads to eradication of melanoma with limited impact on normal cells, and suggested that targeting NGLY1 represents a novel anti-melanoma strategy.

Structural and energetic analysis of 2-aminobenzimidazole inhibitors in complex with the hepatitis C virus IRES RNA using molecular dynamics simulations.

Despite the many biological functions of RNA, very few drugs have been designed or found to target RNA. Here we report the results of molecular dynamics (MD) simulations and binding energy analyses on hepatitis C virus internal ribosome entry site (IRES) RNA in complex with highly charged 2-aminobenzimidazole inhibitors. Initial coordinates were taken from NMR and crystallography studies that had yielded different binding modes. During MD simulations, the RNA-inhibitor complex is stable in the crystal conformation but not in the NMR conformation. Additionally, we found that existing and standard MD trajectory postprocessing free energy methods, such as the MM-GBSA and MM-PBSA approaches available in AMBER, seem unsuitable to properly rank the binding energies of complexes between highly charged molecules. A better correlation with the experimental data was found using a rather simple binding enthalpy calculation based on the explicitly solvated potential energies. In anticipation of further growth in the use of small molecules to target RNA, we include results addressing the impact of charge assignment on docking, the structural role of magnesium in the IRES-inhibitor complex, the entropic contribution to binding energy, and simulations of a plausible scaffold design for new inhibitors.

A Quick Route to Multiple Highly Potent SARS-CoV-2 Main Protease Inhibitors*.

The COVID-19 pathogen, SARS-CoV-2, requires its main protease (SC2MPro ) to digest two of its translated long polypeptides to form a number of mature proteins that are essential for viral replication and pathogenesis. Inhibition of this vital proteolytic process is effective in preventing the virus from replicating in infected cells and therefore provides a potential COVID-19 treatment option. Guided by previous medicinal chemistry studies about SARS-CoV-1 main protease (SC1MPro ), we have designed and synthesized a series of SC2MPro inhibitors that contain ß-(S-2-oxopyrrolidin-3-yl)-alaninal (Opal) for the formation of a reversible covalent bond with the SC2MPro active-site cysteine C145. All inhibitors display high potency with Ki values at or below 100 nM. The most potent compound, MPI3, has as a Ki value of 8.3 nM. Crystallographic analyses of SC2MPro bound to seven inhibitors indicated both formation of a covalent bond with C145 and structural rearrangement from the apoenzyme to accommodate the inhibitors. Virus inhibition assays revealed that several inhibitors have high potency in inhibiting the SARS-CoV-2-induced cytopathogenic effect in both Vero E6 and A549/ACE2 cells. Two inhibitors, MPI5 and MPI8, completely prevented the SARS-CoV-2-induced cytopathogenic effect in Vero E6 cells at 2.5-5 µM and A549/ACE2 cells at 0.16-0.31 µM. Their virus inhibition potency is much higher than that of some existing molecules that are under preclinical and clinical investigations for the treatment of COVID-19. Our study indicates that there is a large chemical space that needs to be explored for the development of SC2MPro inhibitors with ultra-high antiviral potency.

Filtering out Low-Affinity Bitropic Ligands for Dopamine Receptors Based on Ligand Conformation.

We investigated the correlation between the conformations of a set of published 90 bitopic compounds on their affinity for two subtypes of the human dopamine receptor, D2R and D3R. Using molecular dynamics simulations, we showed that the compounds with large populations of compact conformation in the free solution are weak binders to both subtypes of the receptor. Our study provides a computational approach to quickly filter out low-affinity dopamine receptor ligands before their costly chemical synthesis.

A Speedy Route to Multiple Highly Potent SARS-CoV-2 Main Protease Inhibitors.

The COVID-19 pathogen, SARS-CoV-2, requires its main protease (SC2M Pro ) to digest two of its translated polypeptides to form a number of mature proteins that are essential for viral replication and pathogenesis. Inhibition of this vital proteolytic process is effective in preventing the virus from replication in infected cells and therefore provides a potential COVID-19 treatment option. Guided by previous medicinal chemistry studies about SARS-CoV-1 main protease (SC1M Pro ), we have designed and synthesized a series of SC2M Pro inhibitors that contain ß-( S -2-oxopyrrolidin-3-yl)-alaninal (Opal) for the formation of a reversible covalent bond with the SC2M Pro active site cysteine C145. All inhibitors display high potency with IC 50 values at or below 100 nM. The most potent compound MPI3 has as an IC 50 value as 8.5 nM. Crystallographic analyses of SC2M Pro bound to 7 inhibitors indicated both formation of a covalent bond with C145 and structural rearrangement from the apoenzyme to accommodate the inhibitors. Virus inhibition assays revealed that several inhibitors have high potency in inhibiting the SARS-CoV-2-induced cytopathogenic effect in both Vero E6 and A549 cells. Two inhibitors MP5 and MPI8 completely prevented the SARS-CoV-2-induced cytopathogenic effect in Vero E6 cells at 2.5-5 µM and A549 cells at 0.16-0.31 µM. Their virus inhibition potency is much higher than some existing molecules that are under preclinical and clinical investigations for the treatment of COVID-19. Our study indicates that there is a large chemical space that needs to be explored for the development of SC2M Pro inhibitors with extreme potency. Due to the urgent matter of the COVID-19 pandemic, MPI5 and MPI8 may be quickly advanced to preclinical and clinical tests for COVID-19.

Structural and functional insights into the bona fide catalytic state of Streptococcus pyogenes Cas9 HNH nuclease domain.

The CRISPR-associated endonuclease Cas9 from Streptococcus pyogenes (SpyCas9), along with a programmable single-guide RNA (sgRNA), has been exploited as a significant genome-editing tool. Despite the recent advances in determining the SpyCas9 structures and DNA cleavage mechanism, the cleavage-competent conformation of the catalytic HNH nuclease domain of SpyCas9 remains largely elusive and debatable. By integrating computational and experimental approaches, we unveiled and validated the activated Cas9-sgRNA-DNA ternary complex in which the HNH domain is neatly poised for cleaving the target DNA strand. In this catalysis model, the HNH employs the catalytic triad of D839-H840-N863 for cleavage catalysis, rather than previously implicated D839-H840-D861, D837-D839-H840, or D839-H840-D861-N863. Our study contributes critical information to defining the catalytic conformation of the HNH domain and advances the knowledge about the conformational activation underlying Cas9-mediated DNA cleavage.

Consensus Conformations of Dinucleoside Monophosphates Described with Well-Converged Molecular Dynamics Simulations.

Dinucleoside monophosphates (DNMPs) have been described using various experimental approaches as flexible molecules which generate ensembles populating at least a small set of different conformations in solution. However, due to limitations of each approach in its ability to delineate the ensemble of conformations, an accurate and quantitative description of certain conformational features has not been performed for all DNMPs. Here, we apply a temperature replica-exchange molecular dynamics approach to fully and quickly converge conformational distributions of all RNA DNMPs immersed in the TIP3P water model using the AMBER ff14 force field. For a selection of DNMPs, the conformational ensembles were also generated when immersed in the OPC water model using alternative AMBER and CHARMM force fields. The OPC water model and other force field choices did not introduce new conformational classes but shifted the populations among existing conformations. Except for pyrimidine-pyrimidine dinucleosides, all other DNMPs populated four major conformations (which are defined in the main text and labeled A-form, Ladder, Inverted, and Sheared), in addition to an Extended form. Pyrimidine-pyrimidines did not generate the Sheared conformation. Distinguishing features and stabilizing factors of each conformation were identified and assessed based on the known experimental interpretations. The configuration of the glycosidic bond and the nonbonding interactions of hydrogen bond acceptors with the 2'-hydroxyl group were found to play determining roles in stabilizing particular conformations which could serve as a guide for potential force field modifications to improve the accuracy. Additionally, we computed stacking free energies based on the DNMP conformational distributions and found significant discrepancies with a previous study. Our investigation determined that the AMBER force field was incorrectly implemented in the previous study. In the future, this simulation approach can be used to quickly analyze the effects of new force field modifications in shifting the conformational populations of DNMPs, and can can be further applied to foresee such effects in larger RNA motifs including tetranucleotides and tetraloops.

Investigating the ion dependence of the first unfolding step of GTPase-Associating Center ribosomal RNA.

The interactions in the tertiary structure of a ribosomal RNA fragment in the GTPase Associating Center (GAC) have been experimentally studied, but the roles of the bound and diffuse cations in its folding pathway have not yet been fully elucidated. Melting experiments have shown that the temperature of the first of the two distinguishable transitions in the unfolding pathway of the GAC RNA can be regulated by altering the magnesium concentration, yet the physical interpretation of such ion-dependent effects on folding have not been clearly understood in spite of the availability of crystal structures that depict many GAC RNA-ion interactions. Here, we use umbrella sampling and molecular dynamics (MD) simulations to provide a physical description for the first transition in this unfolding pathway, with a focus on the role of a chelated magnesium ion. Our results indicate that the presence of cations mediating the local interaction of two loops stabilizes the folded state relative to the unfolded or partially folded states. Also, our findings suggest that a bridging magnesium ion between the two loops improves the stabilizing effect. This is consistent with the multistep unfolding pathway proposed for the GAC RNA and highlights the importance of ions in the first unfolding step. The results suggest how MD simulations can provide insight into RNA unfolding pathways as a complementary approach to experiments.

Analogues of Arylamide Phenylpiperazine Ligands To Investigate the Factors Influencing D3 Dopamine Receptor Bitropic Binding and Receptor Subtype Selectivity.

We have previously reported on the ability of arylamide phenylpiperazines to bind selectively to the D3 versus the D2 dopamine receptor subtype. For these studies, we used LS-3-134 as the prototypic arylamide phenylpiperazine ligand because it binds with high affinity at D3 dopamine receptor (0.17 nM) and exhibits >150-fold D3 vs D2 receptor binding selectivity. Our goal was to investigate how the composition and size of the nonaromatic ring structure at the piperazine position of substituted phenylpiperazine analogues might influence binding affinity at the human D2 and D3 dopamine receptors. Two factors were identified as being important for determining the binding affinity of bitropic arylamide phenylpiperazines at the dopamine D3 receptor subtype. One factor was the strength of the salt bridge between the highly conserved residue Asp3.32 with the protonated nitrogen of the nonaromatic ring at the piperazine position. The second factor was the configuration of the unbound ligand in an aqueous solution. These two factors were found to be related to the logarithm of the affinities using a simple correlation model, which could be useful when designing high affinity subtype selective bitropic ligands. While this model is based upon the interaction of arylamide phenylpiperazines with the D2 and D3 D2-like dopamine receptor subtypes, it provides insights into the complexity of the factors that define a bitropic mode of the binding at GPCRs.

Computational Assessment of Potassium and Magnesium Ion Binding to a Buried Pocket in GTPase-Associating Center RNA.

An experimentally well-studied model of RNA tertiary structures is a 58mer rRNA fragment, known as GTPase-associating center (GAC) RNA, in which a highly negative pocket walled by phosphate oxygen atoms is stabilized by a chelated cation. Although such deep pockets with more than one direct phosphate to ion chelation site normally include magnesium, as shown in one GAC crystal structure, another GAC crystal structure and solution experiments suggest potassium at this site. Both crystal structures also depict two magnesium ions directly bound to the phosphate groups comprising this controversial pocket. Here, we used classical molecular dynamics simulations as well as umbrella sampling to investigate the possibility of binding of potassium versus magnesium inside the pocket and to better characterize the chelation of one of the binding magnesium ions outside the pocket. The results support the preference of the pocket to accommodate potassium rather than magnesium and suggest that one of the closely binding magnesium ions can only bind at high magnesium concentrations, such as might be present during crystallization. This work illustrates the complementary utility of molecular modeling approaches with atomic-level detail in resolving discrepancies between conflicting experimental results.