RESUMO
A comprehensive understanding of molecular interactions and functions is imperative for unraveling the intricacies of viral protein behavior and conformational dynamics during cellular entry. Focusing on the SARS-CoV-2 spike protein (SARS-CoV-2 sp), a Principal Component Analysis (PCA) on a subset comprising 131 A-chain structures in presence of various inhibitors was conducted. Our analyses unveiled a compelling correlation between PCA modes and Anisotropic Network Model (ANM) modes, underscoring the reliability and functional significance of low-frequency modes in adapting to diverse inhibitor binding scenarios. The role of HR1 in viral processing, both linear Normal Mode Analysis (NMA) and Nonlinear NMA were implemented. Linear NMA exhibited substantial inter-structure variability, as evident from a higher Root Mean Square Deviation (RMSD) range (7.30 Å), nonlinear NMA show stability throughout the simulations (RMSD 4.85 Å). Frequency analysis further emphasized that the energy requirements for conformational changes in nonlinear modes are notably lower compared to their linear counterparts. Using simulations of molecular dynamics at constant pH (cpH-MD), we successfully predicted the pKa order of the interconnected residues within the HR1 mutations at lower pH values, suggesting a transition to a post-fusion structure. The pKa determination study illustrates the profound effects of pH variations on protein structure. Key results include pKa values of 9.5179 for lys-921 in the D936H mutant, 9.50 for the D950N mutant, and a slightly higher value of 10.49 for the D936Y variant. To further understand the behavior and physicochemical characteristics of the protein in a biologically relevant setting, we also examine hydrophobic regions in the prefused states of the HR1 protein mutants D950N, D936Y, and D936H in our study. This analysis was conducted to ascertain the hydrophobic moment of the protein within a lipid environment, shedding light on its behavior and physicochemical properties in a biologically relevant context.
Assuntos
COVID-19 , SARS-CoV-2 , Humanos , SARS-CoV-2/genética , Reprodutibilidade dos Testes , Proteínas/química , Simulação de Dinâmica Molecular , Concentração de Íons de HidrogênioRESUMO
Zika virus (ZIKV) is a global health concern and has been linked to severe neurological pathologies. Although no medication is available yet, many efforts to develop antivirals and host cell binding inhibitors led to attractive drug-like scaffolds, mainly targeting the nonstructural NS2B/NS3 protease (NS2B/NS3pro). NS2B/NS3pro active site has several titratable residues susceptible to pH changes and ligand binding; hence, understanding these residues' protonation is essential to drug design efforts targeting the active site. Here we use in silico methods to probe non-covalent binding and its effect on pKa shifts of the active site residues on a ligand-free protease and with a non-peptidic competitive inhibitor (Ki=13.5 µM). By applying constant pH molecular dynamics, we found that the catalytic residues of the unbound NS2B/NS3pro achieved the protonation needed for the serine protease mechanism over the pH value of 8.5. Nevertheless, the protease in the holo state achieved this same scenario at lower pH values. Also, non-covalent binding affected the catalytic triad (H51, D75, and S135) by stabilizing their distances and interaction network. Thus, NS2B/NS3pro residues configuration for activity might be both pH-dependent and influenced by ligand binding. However, compound presence within the binding site destabilized the NS2B, interfering with the closed and active conformation necessary for substrate binding and catalysis. Our outcomes provide valuable insights into non-covalent inhibitor behavior and its effect on protease active site residues, impacting optimization and design of novel compounds. Communicated by Ramaswamy H. Sarma.
Assuntos
Antivirais , Inibidores de Proteases , Zika virus , Sítios de Ligação , Concentração de Íons de Hidrogênio , Peptídeo Hidrolases/química , Inibidores de Proteases/farmacologia , Inibidores de Proteases/química , Conformação Proteica , Serina Endopeptidases/química , Proteínas não Estruturais Virais/química , Zika virus/efeitos dos fármacos , Antivirais/química , Antivirais/farmacologiaRESUMO
The vitamin D receptor (VDR) constitutes a promising therapeutic target for the treatment of cancer. Unfortunately, its natural agonist calcitriol does not have clinical utility due to its potential to induce hypercalcemic effects at the concentrations required to display antitumoral activity. For this reason, the search for new calcitriol analogues with adequate therapeutic profiles has been actively pursued by the scientific community. We have previously reported the obtaining and the biological activity evaluation of new calcitriol analogues by modification of its sidechain, which exhibited relevant antiproliferative and selectivity profiles against tumoral and normal cells. In this work we conducted molecular modeling studies (i.e. molecular docking, molecular dynamics, constant pH molecular dynamics (CpHMD) and free energy of binding analysis) to elucidate at an atomistic level the molecular basis related to the potential of the new calcitriol analogues to achieve selectivity between tumoral and normal cells. Two histidine residues (His305 and His397) were found to exhibit a particular tautomeric configuration that produces the observed bioactivity. Also, different acid-based properties were observed for His305 and His307 with His305 showing an increased acidity (pKa 5.2) compared to His397 (pKa 6.8) and to the typical histidine residue. This behavior favored the pharmacodynamic interaction of the calcitriol analogues exhibiting selectivity for tumoral cells when VDR was modeled at the more acidic tumoral environment (pH â 6) compared to the case when VDR was modeled at pH 7.4 (normal cell environment). On the other hand, non-selective compounds, including calcitriol, exhibited a similar interaction pattern with VDR when the receptor was modeled at both pH conditions. The results presented constitute the first evidence on the properties of the VDR receptor in different physicochemical environments and thus represent a significant contribution to the in silico screening and design of new calcitriol analogues.