Insights into the DNA binding induced thermal stabilization of transcription factor FOXP3

Mannose 2, 3, 4, 5, 6-O-pentasulfate (MPS): a partial activator of human heparin cofactor II with anticoagulation potential

Thromboembolic diseases are a major cause of mortality in human and the currently available anticoagulants are associated with various drawbacks, therefore the search for anticoagulants that have better safety profile is highly desirable. Compounds that are part of the dietary routine can be modified to possibly increase their anticoagulant potential. We show mannose 2,3,4,5,6-O-pentasulfate (MPS) as a synthetically modified form of mannose that has appreciable anticoagulation properties. An in silico study identified that mannose in sulfated form can bind effectively to the heparin-binding site of antithrombin (ATIII) and heparin cofactor II (HCII). Mannose was sulfated using a simple sulfation strategy-involving triethylamine-sulfur trioxide adduct. HCII and ATIII were purified from human plasma and the binding analysis using fluorometer and isothermal calorimetry showed that MPS binds at a unique site. A thrombin inhibition analysis using the chromogenic substrate showed that MPS partially enhances the activity of HCII. Further an assessment of in vitro blood coagulation assays using human plasma showed that the activated partial thromboplastin time (APTT) and prothrombin time (PT) were prolonged in the presence of MPS. A molecular dynamics simulation analysis of the HCII-MPS complex showed fluctuations in a N-terminal loop and the cofactor binding site of HCII. The results indicate that MPS is a promising lead due to its effect on the in vitro coagulation rate.Communicated by Ramaswamy H. Sarma.

Repurposing of phyto-ligand molecules from the honey bee products for Alzheimer's disease as novel inhibitors of BACE-1: small molecule bioinformatics strategies as amyloid-based therapy

Alzheimer's disease (AD) is one of the neurodegenerative diseases, manifesting dementia, spatial disorientation, language, cognitive, and functional impairment, mainly affects the elderly population with a growing concern about the financial burden on society. Repurposing can improve the traditional progress of drug design applications and could speed up the identification of innovative remedies for AD. The pursuit of potent anti-BACE-1 drugs for AD treatment has become a pot boiler topic in the recent past and to instigate the design of novel improved inhibitors from the bee products. Drug-likeness characteristics (ADMET: absorption, distribution, metabolism, excretion, and toxicity), docking (AutoDock Vina), simulation (GROMACS), and free energy interaction (MM-PBSA, molecular mechanics Poisson-Boltzmann surface area) analyses were performed to identify the lead candidates from the bee products (500 bioactives from the honey, royal jelly, propolis, bee bread, bee wax, and bee venom) for Alzheimer's disease as novel inhibitors of BACE-1 (beta-site amyloid precursor protein cleaving enzyme (1) receptor using appropriate bioinformatics tools. Forty-four bioactive lead compounds were screened from the bee products through high throughput virtual screening on the basis of their pharmacokinetic and pharmacodynamics characteristics, showing favorable intestinal and oral absorption, bioavailability, blood brain barrier penetration, less skin permeability, and no inhibition of cytochrome P450 inhibitors. The docking score of the forty-four ligand molecules was found to be between -4 and -10.3 kcal/mol, respectively, exhibiting strong binding affinity to BACE1 receptor. The highest binding affinity was observed in the rutin (-10.3 kcal/mol), 3,4-dicaffeoylquinic acid (-9.5 kcal/mol), nemorosone (-9.5 kcal/mol), and luteolin (-8.9 kcal/mol). Furthermore, these compounds demonstrated high total binding energy -73.20 to -105.85 kJ/mol), and low root mean square deviation (0.194-0.202 nm), root mean square fluctuation (0.0985-0.1136 nm), radius of gyration (2.12 nm), number of H-bonds (0.778-5.436), and eigenvector values (2.39-3.54 nm²) in the molecular dynamic simulation, signifying restricted motion of Cα atoms, proper folding and flexibility, and highly stable with compact of the BACE1 receptor with the ligands. Docking and simulation studies concluded that rutin, 3,4-dicaffeoylquinic acid, nemorosone, and luteolin are plausibly used as novel inhibitors of BACE1 to combat AD, but further in-depth experimental investigations are warranted to prove these in silico findings.

Targeting SARS-CoV-2 nucleocapsid oligomerization: Insights from molecular docking and molecular dynamics simulations

The outbreak of COVID-19 caused by SARS-CoV-2 virus continually led to infect a large population worldwide. Currently, there is no specific viral protein-targeted therapeutics. The Nucleocapsid (N) protein of the SARS-CoV-2 virus is necessary for viral RNA replication and transcription. The C-terminal domain of N protein (CTD) involves in the self-assembly of N protein into a filament that is packaged into new virions. In this study, the CTD (PDB ID: 6WJI) was targeted for the identification of possible inhibitors of oligomerization of N protein. Herein, multiple computational approaches were employed to explore the potential mechanisms of binding and inhibitor activity of five antiviral drugs toward CTD. The five anti-N drugs studied in this work are 4E1RCat, Silmitasertib, TMCB, Sapanisertib, and Rapamycin. Among the five drugs, 4E1RCat displayed highest binding affinity (-10.95 kcal/mol), followed by rapamycin (-8.91 kcal/mol), silmitasertib (-7.89 kcal/mol), TMCB (-7.05 kcal/mol), and sapanisertib (-6.14 kcal/mol). Subsequently, stability and dynamics of the protein-drug complex were examined with molecular dynamics (MD) simulations. Overall, drug binding increases the stability of the complex with maximum stability observed in the case of 4E1RCat. The CTD-drug complex systems behave differently in terms of the free energy landscape and showed differences in population distribution. Overall, the MD simulation parameters like RMSD, RMSF, Rg, hydrogen bonds analysis, PCA, FEL, and DCCM analysis indicated that 4E1RCat and TMCB complexes were more stable as compared to silmitasertib and sapanisertib and thus could act as effective drug compounds against CTD.Communicated by Ramaswamy H. Sarma.

Delineating the conformational dynamics of intermediate structures on the unfolding pathway of β-lactoglobulin in aqueous urea and dimethyl sulfoxide

The funnel shaped energy landscape model of the protein folding suggests that progression of folding proceeds through multiple pathways, having the multiple intermediates which leads to multidimensional free-energy surface. Herein, we applied all-atom MD simulation to conduct a comparative study on the structure of β-lactoglobulin (β-LgA) in aqueous mixture of 8 M urea and 8 M dimethyl sulfoxide (DMSO), at different temperatures. The cumulative results of multiple simulations suggest a common unfolding pathway of β-LgA, occurred through the stable and meta-stable intermediates (I), in both urea and DMSO. However, the free-energy landscape (FEL) analyses show that the structural transitions of I-states are energetically different. In urea, FEL shows distinct ensemble of intermediates, I1 and I2, separated by the energy barrier of ∼3.0 kcal mol^-1. Similarly, we find the population of two distinct I1 and I2 states in DMSO, however, the I1 appeared transiently around ∼30-35 ns and is short-lived. But, the I2 ensemble is observed structurally compact and long-lived (∼50-150 ns) as compared to unfolding in urea. Furthermore, the I1 and I2 are separated through a high energy barrier of ∼6.0 kcal mol^-1. Thus, our results provide the structural insights of intermediates which essentially bear the signature of a different unfolding pathway of β-LgA in urea and DMSO.Abbreviationsβ-LgAβ-lactoglobulinDMSOdimethyl sulfoxideFELfree-energy landscapeGdmClguanidinium chlorideIintermediate stateMGmolten globule statePMEparticle mesh EwaldQfraction of native contactsRMSDroot mean square deviationRMSFroot mean square fluctuationR_gradius of gyrationSASAsolvent Accessible Surface AreascSASAthe side chain SASATrptryptophanCommunicated by Ramaswamy H. Sarma.

Structure-based virtual screening, molecular dynamics simulation and MM-PBSA toward identifying the inhibitors for two-component regulatory system protein NarL of Mycobacterium Tuberculosis

The nitrate/nitrite response regulatory protein NarL belongs to the two-component regulatory system of Mycobacterium tuberculosis (MTB), plays a crucial role in anaerobic survival of mycobacteria in host. The absence of this protein in humans, makes it an attractive drug target for MTB treatment. However, the specific drug molecules targeting NarL are yet to be identified. In this study, we identified the promising drug candidates using structure based virtual screening of compounds from chemical libraries (ChEMBL and ZINC), followed by the extensive physicochemical properties analyses and molecular dynamics (MD) simulation. As the initial results, we obtained 4,754 bioactive compounds from ChEMBL having anti-tuberculosis activity which is finally narrowed down to the best 10 hits. A similar approach was applied to search for structurally similar compounds from ZINC data, corresponding to the top hits obtained from ChEMBL. Our collective results show that two compounds, ChEMBL509609 (Gscore - 5.054 kcal/mol, Xscore - 6.47 kcal/mol) and ZINC01843143 (Gscore - 5.114 kcal/mol, Xscore - 6.46 kcal/mol) having the best docking score and ADMET profile. The structural stability and dynamics of lead molecules at active site of NarL were examined using MD simulation and the binding free energies were estimated with MM-PBSA. Essential dynamics and MM-PBSA demonstrated that NarL-ChEMBL509609 complex remains the most stable during simulation of 100 ns with the higher binding free energy which may be a suitable candidate for further experimental analysis. AbbreviationsADMEAbsorption, Distribution, Metabolism, And ExcretionBCGBacillus Calmette-GuerinCNSCentral nervous systemDOTSDirectly observed treatment, short courseEDEssential dynamicsHIVHuman immunodeficiency virusHKHistidine kinaseHOAHuman oral absorptionHTVSHigh throughput virtual screeningIRRIIrritationMDMolecular dynamicsMDRMultidrug resistantMTBMycobacterium tuberculosisMUTMutagenicityMWMolecular weightPHOAPercentage of human oral absorptionREPReproductive developmentRgRadius of gyrationRMSDRoot mean square deviationRMSFRoot mean square fluctuationRO5Lipinski's rule of fiveRRResponse regulatorSPStandard precisionSPGStandard precision glideTBTuberculosisTCSTwo-component regulatory systemTDRTotally drug-resistantTUMOTumorigenicityWHOWorld health organizationXDRExtensively drug-resistantXPExtra precisionCommunicated by Ramaswamy H. Sarma.

Characterization of heterogeneous intermediate ensembles on the guanidinium chloride-induced unfolding pathway of β-lactoglobulin

Folding pathway of β-LgA (β-lactoglobulin) evolves through the conformational α→β transition. The α→β transition is a molecular hallmark of various neurodegenerative diseases. Thus, β-LgA may serve as a good model for understanding molecular mechanism of protein aggregation involved in neurodegenerative diseases. Here, we studied the conformational dynamics of β-LgA in 6 M GdmCl at different temperatures using MD simulations. Structural order parameters such as RMSD, R_g, SASA, native contacts (Q), hydrophobic distal-matrix and free-energy landscape (FEL) were used to investigate the conformational transitions. Our results show that GdmCl destabilizes secondary and tertiary structure of β-LgA by weakening the hydrophobic interactions and hydrogen bond network. Multidimensional FEL shows the presence of different unfolding intermediates at 400 K. I1 is long-lived intermediate which has mostly intact native secondary structure, but loose tertiary structure. I2 is structurally compact intermediate formed after the partial loss of secondary structure. The transiently and infrequently buried evolution of W19 shows that intermediate conformational ensembles are structurally heterogeneous. We observed that the intermediate conformations are largely stabilized by non-native H-bonds. The outcome of this work provides the molecular details of intermediates trapped due to non-native interactions that may be regarded as pathogenic conformations involved in neurodegenerative diseases.Communicated by Ramaswamy H. Sarma.