Elucidating structural variability in p53 conformers using combinatorial refinement strategies and molecular dynamics

Low molecular weight proteins and protein assemblies can now be investigated using cryo-electron microscopy (EM) as a complement to traditional structural biology techniques. It is important, however, to not lose sight of the dynamic information inherent in macromolecules that give rise to their exquisite functionality. As computational methods continue to advance the field of biomedical imaging, so must strategies to resolve the minute details of disease-related entities. Here, we employed combinatorial modeling approaches to assess flexible properties among low molecular weight proteins (~100 kDa or less). Through a blend of rigid body refinement and simulated annealing, we determined new hidden conformations for wild type p53 monomer and dimer forms. Structures for both states converged to yield new conformers, each revealing good stereochemistry and dynamic information about the protein. Based on these insights, we identified fluid parts of p53 that complement the stable central core of the protein responsible for engaging DNA. Molecular dynamics simulations corroborated the modeling results and helped pinpoint the more flexible residues in wild type p53. Overall, the new computational methods may be used to shed light on other small protein features in a vast ensemble of structural data that cannot be easily delineated by other algorithms.

Structural and functional effects of the L84S mutant in the SARS-COV-2 ORF8 dimer based on microsecond molecular dynamics study

The L84S mutation has been observed frequently in the ORF8 protein of SARS-CoV-2, which is an accessory protein involved in various important functions such as virus propagation, pathogenesis, and evading the immune response. However, the specific effects of this mutation on the dimeric structure of ORF8 and its impacts on interactions with host components and immune responses are not well understood. In this study, we performed one microsecond molecular dynamics (MD) simulation and analyzed the dimeric behavior of the L84S and L84A mutants in comparison to the native protein. The MD simulations revealed that both mutations caused changes in the conformation of the ORF8 dimer, influenced protein folding mechanisms, and affected the overall structural stability. In particular, the ⁷³YIDI⁷⁶ motif has found to be significantly affected by the L84S mutation, leading to structural flexibility in the region connecting the C-terminal β4 and β5 strands. This flexibility might be responsible for virus immune modulation.  The free energy landscape (FEL) and principle component analysis (PCA) have also supported our investigation. Overall, the L84S and L84A mutations affect the ORF8 dimeric interfaces by reducing the frequency of protein-protein interacting residues (Arg52, Lys53, Arg98, Ile104, Arg115, Val117, Asp119, Phe120, and Ile121) in the ORF8 dimer.  Our findings provide detail insights for further research in designing structure-based therapeutics against the SARS-CoV-2.Communicated by Ramaswamy H. Sarma.

Repurposing fusion inhibitor peptide against SARS-CoV-2

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is continuously evolving. Although several vaccines were approved, this pandemic is still a major threat to public life. Till date, no established therapies are available against SARS-CoV-2. Peptide inhibitors hold great promise for this viral pathogen due to their efficacy, safety, and specificity. In this study, seventeen antiviral peptides which were known to inhibit SARS-CoV-1 are collected and computationally screened against heptad repeat 1 (HR1) of the SARS-CoV-2 spike protein (S2). Out of 17 peptides, Fp13 and Fp14 showed better binding affinity toward HR1 compared to a control peptide EK1 (a modified pan-coronavirus fusion inhibitor) in molecular docking. To explore the time-dependent interactions of the fusion peptide with HR1, molecular dynamics simulation was performed incorporating lipid membrane. During 100 ns MD simulation, structural and energy parameters of Fp13-HR1 and Fp14-HR1 complexes demonstrated lower fluctuations compared to the control EK1-HR1 complex. Furthermore, principal component analysis and free energy landscape study revealed that these two peptides (Fp13 and Fp14) strongly bind to the HR1 with higher affinity than that of control EK1. Tyr917, Asn919, Gln926, lys933, and Gln949 residues in HR1 protein were found to be crucial residues for peptide interaction. Notably, Fp13, Fp14 showed reasonably better binding free energy and hydrogen bond contribution than that of EK1. Taken together, Fp13 and Fp14 peptides may be highly specific for HR1 which can potentially prevent the formation of the fusion core and could be further developed as therapeutics for treatment or prophylaxis of SARS-CoV-2 infection.

Designing potent inhibitors against the multidrug resistance P-glycoprotein

The multidrug transporter P-glycoprotein is an ATP binding cassette (ABC) exporter responsible for resistance to tumor cells during chemotherapy. This study was designed with computational approaches aimed at identifying the best potent inhibitors of P-glycoprotein. Although many compounds have been suggested to inhibit P-glycoprotein, however, their information on bioavailability, selectivity, ADMET properties, and molecular interactions has not been revealed. Molecular docking, ADMET analysis, molecular dynamics, Principal component analysis (PCA), and binding free energy calculations were performed. Two compounds D1 and D2 showed the best docking score against P-glycoprotein and both compounds have 4-thiazolidinone derivatives containing indolin-3 one moiety are novel anti-tumor compounds. ADMET calculation analysis predicted D1 and D2 to have acceptable pharmacokinetic properties. The MD simulation discloses that D1-P-glycoprotein and D2-P-glycoprotein complexes are in stable conformation as apo-form. Hydrophobic amino acid such as phenylalanine plays significant on the interactions of inhibitors. Principal component analysis shows that both complexes are relatively similar variables as apo-form except planarity and Columbo energy profile. In addition, Quantitative Structural Activity Relationship (QSAR) of the ligand candidates were subjected to the principal component analysis (PCA) for pattern recognition. Partial-least-square (PLS) regression analysis was further utilized to model drug candidates' QSAR for subsequent prediction of the binding energy of validated drug candidates. PCA revealed groupings of the drug candidates based on the similarity or differences in drug candidates QSAR. Moreover, the developed PLS regression accurately predicted the values of the binding energy of drug candidates, with low residual error of prediction.Communicated by Ramaswamy H. Sarma.

Cysteine focused covalent inhibitors against the main protease of SARS-CoV-2

In viral replication and transcription, the main protease (Mpro) of SARS-CoV-2 plays an important role and appears to be a vital target for drug design. In Mpro, there is a Cys-His catalytic dyad, and ligands that interact with the Cys145 assumed to be an effective approach to inhibit the Mpro. In this study, approximately 1400 cysteine-focused ligands were screened to identify the best candidates that can act as potent inhibitors against Mpro. Our results show that the selected ligands strongly interact with the key Cys145 and His41 residues. Covalent docking was performed for the selected candidates containing the acrylonitrile group, which can form a covalent bond with Cys145. All atoms molecular dynamics (MD) simulation was performed on the selected four inhibitors including L1, L2, L3 and L4 to validate the docking interactions. Our results were also compared with a control ligand, α-ketoamide (11r). Principal component analysis on structural and energy data obtained from the MD trajectories shows that L1, L3, L4 and α-ketoamide (11r) have structural similarity with the apo-form of the Mpro. Quantitative structure-activity relationship method was employed for pattern recognition of the best ligands, which discloses that ligands containing acrylonitrile and amide warheads can show better performance. ADMET analysis displays that our selected candidates appear to be safer inhibitors. Our combined studies suggest that the best cysteine focused ligands can help to design an effective lead drug for COVID-19 treatment. Communicated by Ramaswamy H. Sarma.

Prediction of Deleterious Non-synonymous SNPs of Human STK11 Gene by Combining Algorithms, Molecular Docking, and Molecular Dynamics Simulation

Serine-threonine kinase11 (STK11) is a tumor suppressor gene which plays a key role in regulating cell growth and apoptosis. It is widely known as a multitasking kinase and engaged in cell polarity, cell cycle arrest, chromatin remodeling, energy metabolism, and Wnt signaling. The substitutions of single amino acids in highly conserved regions of the STK11 protein are associated with Peutz-Jeghers syndrome (PJS), which is an autosomal dominant inherited disorder. The abnormal function of the STK11 protein is still not well understood. In this study, we classified disease susceptible single nucleotide polymorphisms (SNPs) in STK11 by using different computational algorithms. We identified the deleterious nsSNPs, constructed mutant protein structures, and evaluated the impact of mutation by employing molecular docking and molecular dynamics analysis. Our results show that W239R and W308C variants are likely to be highly deleterious mutations found in the catalytic kinase domain, which may destabilize structure and disrupt the activation of the STK11 protein as well as reduce its catalytic efficiency. The W239R mutant is likely to have a greater impact on destabilizing the protein structure compared to the W308C mutant. In conclusion, these mutants can help to further realize the large pool of disease susceptibilities linked with catalytic kinase domain activation of STK11 and assist to develop an effective drug for associated diseases.

Exploring the potent inhibitors and binding modes of phospholipase A2 through in silico investigation

Snake venom of Naja naja comprises of several types of enzymes, and among them, water-soluble proteolytic enzyme, phospholipase A2 (PLA₂), is noteworthy for its numerous adverse effects, such as cytotoxicity, cardiotoxicity, hemolytic, anti-coagulant, and hypotensive effects, including being highly potent as a neurotoxin. Limited anti-venom therapy (with their lower efficacy) has attracted considerable pharmacological interest to develop potent inhibitors of PLA₂. Thus, 34 experimentally proven and diverse synthetic inhibitors of PLA₂ were screened primarily on the basis of Glide extra precision docking and MM-GBSA rescoring function. Then, ten potential hits were subjected to induced fit docking, in which top three potential inhibitors were considered, and those were found to interact with Ca²⁺, disulfide binding site, and phosphatidylcholine activation sites, thereby, possibly disrupting the catalytic activity of Ca²⁺ as well as the inflammatory functions of PLA₂. These compounds showed positive remarks on various physiochemical properties and pharmacologically relevant descriptors. Gap energy and thermodynamic properties were investigated by employing density functional theory for all compounds to understand their chemical reactivity and thermodynamic stability. Molecular dynamics simulation was performed for 100 ns in order to evaluate the stability and binding modes of docked complexes, and the energy of binding was calculated through MM-PBSA analysis. On the whole, the proposed compounds could be used for targeted inhibition. Communicated by Ramaswamy H. Sarma.

Quantum chemical calculation and binding modes of H1R; a combined study of molecular docking and DFT for suggesting therapeutically potent H1R antagonist

Histamine-1 receptor (H1R) belongs to the family of rhodopsin-like G-protein-coupled receptors expressed in cells that mediates allergies and other pathophysiological diseases. For alleviation of allergic symptoms, H1R antagonists are therapeutic drugs; of which the most frequently prescribed are second generation drugs, such as; Cetirizine, Loratadine, Hydroxyzine, Desloratadine, Bepotastine, Acrivastine and Rupatadine. To understand their potency, binding affinity and interaction; we have employed molecular docking and quantum chemical study such as; Induced-fit docking and calculation of quantum chemical descriptors. This study also introduces the binding site characterization of H1R, with its known antagonists and Curcumin (our proposed alternative H1R antagonist); useful for future drug target site. The interactive binding site residues of H1R are found to be; Lys-191, Tyr-108, Asp-107, Tyr-100, Lys-179, Lys-191, Thr-194, Trp-428, Phe-432, Tyr-458, Hie-450, with most of these shown to be inhibited by naturally-occurring compound curcumin. Amongst the FDA approved drugs, Hydroxyzine showed best ligand binding affinity, calculated as - 141.491 kcal/mol and naturally occurring compound, Curcumin showed binding affinity of - 87.046 kcal/mol. The known antagonists of H1R has been used for hypothesizing curcumin as naturally occurring lead compound for the target using accurate molecular docking simulation study. Though the pharmacological action of known inhibitor is already established, they could differ from their reactivity, which we have also focused in our study for predicting drug reactivity.