Type-II transmembrane serine proteases are effective pharmacological targets for host defence against viral entry and in certain cancer cell progressions. These serine proteases cleave viral spike proteins to expose the fusion peptide for cell entry, which is essential to the life cycle of the virus. TMPRSS2 inhibitors can also fight against respiratory viruses that employ them for cell entry. Our study combining virtual screening, all-atom molecular dynamics, and well-tempered metadynamics simulation identifies vicenin-2, neohesperidin, naringin, and rhoifolin as promising TMPRSS2 antagonists. The binding energies obtained are - 16.3, - 15.4, - 13.6, and - 13.8 kcal/mol for vicenin-2, neohesperidin, naringin, and rhoifolin respectively. The RMSD, RMSF, PCA, DCCM, and binding free energy profiles also correlate with the stable binding of these ligands at the active site of TMPRSS2. The study reveals that these molecules could be promising lead molecules for combating future outbreaks of coronavirus and other respiratory viruses.

Activation of heterotrimeric G-proteins (Gαβγ) downstream to receptor tyrosine kinases (RTKs) is a well-established crosstalk between the signaling pathways mediated by G-protein coupled receptors (GPCRs) and RTKs. While GPCR serves as a guanine exchange factor (GEF) in the canonical activation of Gα that facilitates the exchange of GDP for GTP, the mechanism through which RTK phosphorylations induce Gα activation remains unclear. Recent experimental studies revealed that the epidermal growth factor receptor (EGFR), a well-known RTK, phosphorylates the helical domain tyrosine residues Y154 and Y155 and accelerates the GDP release from the Gαi3, a subtype of Gα-protein. Using well-tempered metadynamics and extensive unbiased molecular dynamics simulations, we captured the GDP release event and identified the intermediates between bound and unbound states through Markov state models. In addition to weakened salt bridges at the domain interface, phosphorylations induced the unfolding of helix αF, which contributed to increased flexibility near the hinge region, facilitating a greater distance between domains in the phosphorylated Gαi3. Although the larger domain separation in the phosphorylated system provided an unobstructed path for the nucleotide, the accelerated release of GDP was attributed to increased fluctuations in several conserved regions like P-loop, switch 1, and switch 2. Overall, this study provides atomistic insights into the activation of G-proteins induced by RTK phosphorylations and identifies the specific structural motifs involved in the process. The knowledge gained from the study could establish a foundation for targeting non-canonical signaling pathways and developing therapeutic strategies against the ailments associated with dysregulated G-protein signaling.

New psychoactive substances (NPS) targeting cannabinoid receptor 1 pose a significant threat to society as recreational abusive drugs that have pronounced physiological side effects. These greater adverse effects compared to classical cannabinoids have been linked to the higher downstream β-arrestin signaling. Thus, understanding the mechanism of differential signaling will reveal important structure-activity relationship essential for identifying and potentially regulating NPS molecules. In this study, we simulate the slow (un)binding process of NPS MDMB-Fubinaca and classical cannabinoid HU-210 from CB1 using multi-ensemble simulation to decipher the effects of ligand binding dynamics on downstream signaling. The transition-based reweighing method is used for the estimation of transition rates and underlying thermodynamics of (un)binding processes of ligands with nanomolar affinities. Our analyses reveal major interaction differences with transmembrane TM7 between NPS and classical cannabinoids. A variational autoencoder-based approach, neural relational inference (NRI), is applied to assess the allosteric effects on intracellular regions attributable to variations in binding pocket interactions. NRI analysis indicate a heightened level of allosteric control of NPxxY motif for NPS-bound receptors, which contributes to the higher probability of formation of a crucial triad interaction (Y7.53-Y5.58-T3.46) necessary for stronger β-arrestin signaling. Hence, in this work, MD simulation, data-driven statistical methods, and deep learning point out the structural basis for the heightened physiological side effects associated with NPS, contributing to efforts aimed at mitigating their public health impact.

Phosphodiesterase-5 (PDE5) is responsible for regulating the concentration of the second messenger molecule cGMP by hydrolyzing it into 5\'-GMP. PDE5 is implicated in erectile dysfunction and cardiovascular diseases. The substrate binding site in the catalytic domain of PDE5 is surrounded by several dynamic structural motifs (including the α14 helix, M-loop, and H-loop) that are known to switch between inactive and active conformational states via currently unresolved structural intermediates. We evaluated the conformational dynamics of these structural motifs in the apo state and upon binding of an allosteric inhibitor (evodiamine) or avanafil, a competitive inhibitor. We employed enhanced sampling-based replica exchange solute scaling (REST2) method, principal component analysis (PCA), time-lagged independent component analysis (tICA), molecular dynamics (MD) simulations, and well-tempered metadynamics simulations to probe the conformational changes in these structural motifs. Our results support a regulatory mechanism for PDE5, where the α14 helix alternates between an inward (lower activity) conformation and an outward (higher activity) conformation that is accompanied by the folding/unfolding of the α8\' and α8″ helices of the H-loop. When the allosteric inhibitor evodiamine is bound to PDE5, the inward (inactive) state of the α14 helix is preferred, thus preventing substrate access to the catalytic site. In contrast, competitive inhibitors of PDE5 block catalysis by occupying the active site accompanied by stabilization of the outward conformation of the α14 helix. Defining the conformational dynamics underlying regulation of PDE5 activation will be helpful in rational design of next-generation small molecules modulators of PDE5 activity.

Pancreatic lipase (PL) is a keen target for anti-obesity therapy that reduces dietary fat absorption. Here, we investigated the binding patterns of 220 PL inhibitors having experimental IC50 values, using molecular docking and binding energy calculations. Screening of these compounds illustrated most of them bound at the catalytic site (S1-S2 channel) and a few compounds are at the non-catalytic site (S2-S3 channel/S1-S3 channel) of PL. This binding pattern could be due to structural uniqueness or bias in conformational search. A strong correlation of pIC50 values with SP/XP docking scores, binding energies (ΔGMMGBSA) assured the binding poses are more true positives. Further, understanding of each class and subclasses of polyphenols indicated tannins preferred non-catalytic site wherein binding energies are underestimated due to huge desolvation energy. In contrast, most of the flavonoids and furan-flavonoids have good binding energies due to strong interactions with catalytic residues. While scoring functions limited the understanding of sub-classes of flavonoids. Hence, focused on 55 potent PL inhibitors of IC50 < 5 µM for better in vivo efficacy. The prediction of bioactivity, drug-likeness properties, led to 14 bioactive compounds. The low root mean square deviation (0.1-0.2 nm) of these potent flavonoids and non-flavonoid/non-polyphenols PL-inhibitor complexes during 100 ns molecular dynamics runs (MD) as well as binding energies obtained from both MD and well-tempered metadynamics, support strong binding to catalytic site. Based on the bioactivity, ADMET properties, and binding affinity data of MD and wt-metaD of potent PL-inhibitors suggests Epiafzelechin 3-O-gallate, Sanggenon C, and Sanggenofuran A shall be promising inhibitors at in vivo conditions.Communicated by Ramaswamy H. Sarma.

Choline geranate (CAGE) ionic liquids (ILs) stabilize insulin, thereby aiding its oral delivery, whereas ethanol (EtOH) affects its stability by disrupting the hydrophobic interactions. In this study, cognizance of the stabilization mechanism of insulin dimer in the presence of both CAGE ILs and EtOH mixtures is achieved through biased and unbiased molecular dynamics (MD) simulations. Here, two order parameters are employed to study the insulin dimer dissociation using well-tempered metadynamics (WT-MetaD). The stability of insulin is found to be strongly maintained until a 0.20 mole fraction of EtOH. Besides, higher concentrations of EtOH marginally affect the insulin stability. Moreover, geranate anions form a higher number of H-bonding interactions with water molecules, which aids insulin stabilization. Conversely, the addition of EtOH minimizes the water-mediated H-bonding interactions of geranate. Additionally, geranate traps the EtOH molecules, thereby preventing the interactions between insulin and EtOH. Furthermore, the free energy landscape (FEL) reveals the absence of dimer dissociation along with noticeable deviations in the distances R and the number of contacts Q. The dimerization free energy of insulin was calculated to be -16.1 kcal/mol at a 0.20 mole fraction of EtOH. Moreover, increments in mole fractions of EtOH effectuate a decrease in the insulin stability. Thus, the present study represents CAGE ILs as efficient insulin dimer stabilizes at low concentrations of EtOH.

Using supramolecules for protein function regulation is an effective strategy in chemical biology and drug discovery. However, due to the presence of multiple binding sites on protein surfaces, protein function regulation via selective binding of supramolecules is challenging. Recently, the functions of 14-3-3 proteins, which play an important role in regulating intracellular signaling pathways via protein-protein interactions, have been modulated using a supramolecular tweezer, CLR01. However, the binding mechanisms of the tweezer molecule to 14-3-3 proteins are still unclear, which has hindered the development of novel supramolecules targeting the 14-3-3 proteins. Herein, the binding mechanisms of the tweezer to the lysine residues on 14-3-3σ (an isoform in 14-3-3 protein family) were explored by well-tempered metadynamics. The results indicated that the inclusion complex formed between the protein and supramolecule is affected by both kinetic and thermodynamic factors. In particular, simulations confirmed that K214 could form a strong binding complex with the tweezer; the binding free energy was calculated to be -10.5 kcal·mol-1 with an association barrier height of 3.7 kcal·mol-1. In addition, several other lysine residues on 14-3-3σ were identified as being well-recognized by the tweezer, which agrees with experimental results, although only K214/tweezer was co-crystallized. Additionally, the binding mechanisms of the tweezer to all lysine residues were analyzed by exploring the representative conformations during the formation of the inclusion complex. This could be helpful for the development of new inhibitors based on tweezers with more functions against 14-3-3 proteins via modifications of CLR01. We also believe that the proposed computational strategies can be extended to understand the binding mechanism of multi-binding sites proteins with supramolecules and will, thus, be useful toward drug design.

The addition of aprotic solvents results in higher reactivities and selectivities in many key aqueous phase biomass reactions, including the acid-catalyzed conversion of fructose to 5-hydroxyl methyl furfural (HMF). The addition of certain co-solvents inhibits the formation of humins via preferential solvation of key functional groups and can alter reaction kinetics. An important factor in this context is the relative stability of the hydronium ion (the catalyst) in the vicinity of the biomass moiety as compared to that in bulk, as it could determine its efficacy in the protonation step. Hence, in the present work, molecular dynamics (MD) simulations of HMF (the model product) and fructose (the model reactant) in acidic water and water-DMSO mixtures are performed to analyze their interaction with the hydronium ions. We show that the presence of DMSO favors the interaction of the hydronium ion with fructose, whereas it has a detrimental effect on the interaction of hydronium ion with HMF. Well-tempered metadynamics (WT-MTD) simulations are performed to determine the relative stability of the hydronium ion in the immediate vicinity of fructose and HMF, as compared to that in the bulk solvent phase, as a function of solvent composition. We find that DMSO improves the stabilization of the hydronium ions in the first solvation shell of fructose compared to that in the bulk solvent. On the other hand, hydronium ions become less stable in the immediate vicinity of HMF, as the concentration of DMSO increases.

Computational drug design is increasingly becoming important with new and unforeseen diseases like COVID-19. In this study, we present a new computational de novo drug design and repurposing method and applied it to find plausible drug candidates for the receptor binding domain (RBD) of SARS-CoV-2 (COVID-19). Our study comprises three steps: atom-by-atom generation of new molecules around a receptor, structural similarity mapping to existing approved and investigational drugs, and validation of their binding strengths to the viral spike proteins based on rigorous all-atom, explicit-water well-tempered metadynamics free energy calculations. By choosing the receptor binding domain of the viral spike protein, we showed that some of our new molecules and some of the repurposable drugs have stronger binding to RBD than hACE2. To validate our approach, we also calculated the free energy of hACE2 and RBD, and found it to be in an excellent agreement with experiments. These pool of drugs will allow strategic repurposing against COVID-19 for a particular prevailing conditions.

Riboswitches are metabolite sensing aptamer domains present in non-coding regions in RNA and act as gene-regulating elements. Thiamine pyrophosphate (TPP) riboswitch is evolved as a new target for developing antibiotics against many pathogenic bacteria. The earlier reports suggest that the modification of the pyrophosphate group in the ligand molecule can enhance gene expression. In this work, we have examined the binding affinity and efficacy of TPP and two recently reported ligands, CH2-TPP, and CF2-TPP, using Well-tempered metadynamics (WT-MtD) simulations. The experimental in vitro assays show that both TPP and CH2-TPP repress the gene expression to the same extent. The calculated binding energies correlate well with the experimental study and show the same trend of binding affinity of ligands for the TPP riboswitch. The root mean square fluctuation profiles suggest that the CH2-TPP and TPP trigger higher fluctuations in P1 and L3 region, and such fluctuations in the P1 region is involved in the gene regulation process. The metal ion mediated contact of TPP ligand with pyrophosphate binding helix is found to be critical in the gene regulation process. The simulation results corroborate the experimental observations that the role of conformational changes occurring in different riboswitch regions upon ligand binding is essential to repress the gene expression process. This work sheds light on the subtle change in the ligand structure that can induce a more considerable impact on binding affinity and efficacy of ligands with riboswitch.