Molecular mechanism of regulation of RhoA GTPase by phosphorylation of RhoGDI.

Rho-specific guanine nucleotide dissociation inhibitors (RhoGDIs) play a crucial role in the regulation of Rho family GTPases. They act as negative regulators that prevent the activation of Rho GTPases by forming complexes with the inactive GDP-bound state of GTPase. Release of Rho GTPase from the RhoGDI-bound complex is necessary for Rho GTPase activation. Biochemical studies provide evidence of a "phosphorylation code," where phosphorylation of some specific residues of RhoGDI selectively releases its GTPase partner (RhoA, Rac1, Cdc42, etc.). This work attempts to understand the molecular mechanism behind this specific phosphorylation-induced reduction in binding affinity. Using several microseconds long atomistic molecular dynamics simulations of the wild-type and phosphorylated states of the RhoA-RhoGDI complex, we propose a molecular-interaction-based mechanistic model for the dissociation of the complex. Phosphorylation induces major structural changes, particularly in the positively charged polybasic region (PBR) of RhoA and the negatively charged N-terminal region of RhoGDI that contribute most to the binding affinity. Molecular mechanics Poisson-Boltzmann surface area binding energy calculations show a significant weakening of interaction on phosphorylation at the RhoA-specific site of RhoGDI. In contrast, phosphorylation at a Rac1-specific site does not affect the overall binding affinity significantly, which confirms the presence of a phosphorylation code. RhoA-specific phosphorylation leads to a reduction in the number of contacts between the PBR of RhoA and the N-terminal region of RhoGDI, which manifests a reduction of the binding affinity. Using hydrogen bond occupancy analysis and energetic perturbation network, we propose a mechanistic model for the allosteric response, i.e., long-range signal propagation from the site of phosphorylation to the PBR and buried geranylgeranyl group in the form of rearrangement and rewiring of hydrogen bonds and salt bridges. Our results highlight the crucial role of specific electrostatic interactions in manifestation of the phosphorylation code.

hERG stereoselective modulation by mexiletine-derived ureas: Molecular docking study, synthesis, and biological evaluation.

Long QT syndrome (LQTS) is a disorder of cardiac electrophysiology resulting in life-threatening arrhythmias; nowadays, only a few drugs are available for the management of LQTS. Focusing our attention on LQT2, one of the most common subtypes of LQTS caused by mutations in the human ether-à-go-go-related gene (hERG), in the present work, the stereoselectivity of the recently discovered mexiletine-derived urea 8 was investigated on the hERG potassium channel. According to preliminary in silico predictions, in vitro studies revealed a stereoselective behavior, with the meso form showing the greatest hERG opening activity. In addition, functional studies on guinea pig isolated left atria, aorta, and ileum demonstrated that 8 does not present any cardiac or intestinal liability in our ex vivo studies. Due to its overall profile, (R,S)-8 paves the way for the design and development of a new series of compounds potentially useful in the treatment of both congenital and drug-induced forms of LQTS.

Integrated docking and enhanced sampling-based selection of repurposing drugs for SARS-CoV-2 by targeting host dependent factors.

Since the onset of global pandemic, the most focused research currently in progress is the development of potential drug candidates and clinical trials of existing FDA approved drugs for other relevant diseases, in order to repurpose them for the COVID-19. At the same time, several high throughput screenings of drugs have been reported to inhibit the viral components during the early course of infection but with little proven efficacies. Here, we investigate the drug repurposing strategies to counteract the coronavirus infection which involves several potential targetable host proteins involved in viral replication and disease progression. We report the high throughput analysis of literature-derived repurposing drug candidates that can be used to target the genetic regulators known to interact with viral proteins based on experimental and interactome studies. In this work we have performed integrated molecular docking followed by molecular dynamics (MD) simulations and free energy calculations through an expedite in silico process where the number of screened candidates reduces sequentially at every step based on physicochemical interactions. We elucidate that in addition to the pre-clinical and FDA approved drugs that targets specific regulatory proteins, a range of chemical compounds (Nafamostat, Chloramphenicol, Ponatinib) binds to the other gene transcription and translation regulatory proteins with higher affinity and may harbour potential for therapeutic uses. There is a rapid growing interest in the development of combination therapy for COVID-19 to target multiple enzymes/pathways. Our in silico approach would be useful in generating leads for experimental screening for rapid drug repurposing against SARS-CoV-2 interacting host proteins.Communicated by Ramaswamy H. Sarma.

Organoselenium Compounds as Acetylcholinesterase Inhibitors: Evidence and Mechanism of Mixed Inhibition.

Acetylcholinesterase (AChE) inhibitors are actively used for the effective treatment of Alzheimer's disease. In recent years, the neuroprotective effects of organoselenium compounds such as ebselen and diselenides on the AChE activity have been investigated as potential therapeutic agents. In this work, we have carried out systematic kinetic and intrinsic fluorescence assays in combination with docking and molecular dynamics (MD) simulations to elucidate the molecular mechanism of the mixed inhibition of AChE by ebselen and diphenyl diselenide (DPDSe) molecules. Our MD simulations demonstrate significant heterogeneity in the binding modes and allosteric hotspots for DPDSe on AChE due to non-specific interactions. We have further identified that both ebselen and DPDSe can strongly bind around the peripheral anionic site (PAS), leading to non-competitive inhibition similar to other PAS-binding inhibitors. We also illustrate the entry of the DPDSe molecule into the gorge through a "side door", which offers an alternate entry point for AChE inhibitors as compared to the usual substrate entry point of the gorge. Together with results from experiments, these simulations provide mechanistic insights into the mixed type of inhibition for AChE using DPDSe as a promising inhibitor for AChE.

Protonation-Induced Dynamic Allostery in PDZ Domain: Evidence of Perturbation-Independent Universal Response Network.

Dynamic allostery is a relatively new paradigm where certain external perturbations may lead to modulation of conformational dynamics at a distant part of a protein without significant changes in the overall structure. While most well-characterized examples of dynamic allostery involve binding with other entities like small molecules, peptides, or nucleic acids, in this work we demonstrate that chemical modifications like protonation may lead to significant dynamical allosteric response in a PDZ domain protein. Tuning the protonation states of two histidine residues (H317 and H372), we identify the allosteric pathways responsible for the dynamic response. Interestingly, the same set of residues that constitute the allosteric response network upon ligand binding seem to be responsible for protonation-induced dynamic allostery. Thus, we propose the existence of an inherent universal response network in signaling proteins, where the same set of residues can respond to varying types of external perturbations in terms of rearrangement of hydrogen-bonded network and redistribution of electrostatic interaction energies.

Hidden electrostatic basis of dynamic allostery in a PDZ domain.

Allosteric effect implies ligand binding at one site leading to structural and/or dynamical changes at a distant site. PDZ domains are classic examples of dynamic allostery without conformational changes, where distal side-chain dynamics is modulated on ligand binding and the origin has been attributed to entropic effects. In this work, we unearth the energetic basis of the observed dynamic allostery in a PDZ3 domain protein using molecular dynamics simulations. We demonstrate that electrostatic interaction provides a highly sensitive yardstick to probe the allosteric modulation in contrast to the traditionally used structure-based parameters. There is a significant population shift in the hydrogen-bonded network and salt bridges involving side chains on ligand binding. The ligand creates a local energetic perturbation that propagates in the form of dominolike changes in interresidue interaction pattern. There are significant changes in the nature of specific interactions (nonpolar/polar) between interresidue contacts and accompanied side-chain reorientations that drive the major redistribution of energy. Interestingly, this internal redistribution and rewiring of side-chain interactions led to large cancellations resulting in small change in the overall enthalpy of the protein, thus making it difficult to detect experimentally. In contrast to the prevailing focus on the entropic or dynamic effects, we show that the internal redistribution and population shift in specific electrostatic interactions drive the allosteric modulation in the PDZ3 domain protein.

Nucleotide Dependent Switching in Rho GTPase: Conformational Heterogeneity and Competing Molecular Interactions.

Ras superfamily of GTPases regulate myriad cellular processes through a conserved nucleotide (GTP/GDP) dependent switching mechanism. Unlike Ras family of GTPases, for the Rho GTPases, there is no clear evidence for the existence of "sub-states" such as state 1 &state 2 in the GTP bound form. To explore the nucleotide dependent conformational space of the Switch I loop and also to look for existence of state 1 like conformations in Rho GTPases, atomistic molecular dynamics and metadynamics simulations on RhoA were performed. These studies demonstrate that both the nucleotide-free state and the GDP bound "OFF" state have very similar conformations, whereas the GTP bound "ON" state has unique conformations with signatures of two intermediate states. The conformational free energy landscape for these systems suggests the presence of multiple intermediate states. Interestingly, the energetic penalty of exposing the non-polar residues in the GTP bound form is counter balanced by the favourable hydrogen bonded interactions between the γ-phosphate group of GTP with the highly conserved Tyr34 and Thr37 residues. These competing molecular interactions lead to a tuneable energy landscape of the Switch I conformation, which can undergo significant changes based on the local environment including changes upon binding to effectors.