QM modeling
Dotaz
Zobrazit nápovědu
Since its early days in the 19th century, medicinal chemistry has concentrated its efforts on the treatment of diseases, using tools from areas such as chemistry, pharmacology, and molecular biology. The understanding of biological mechanisms and signaling pathways is crucial information for the development of potential agents for the treatment of diseases mainly because they are such complex processes. Given the limitations that the experimental approach presents, computational chemistry is a valuable alternative for the study of these systems and their behavior. Thus, classical molecular dynamics, based on Newton's laws, is considered a technique of great accuracy, when appropriated force fields are used, and provides satisfactory contributions to the scientific community. However, as many configurations are generated in a large MD simulation, methods such as Statistical Inefficiency and Optimal Wavelet Signal Compression Algorithm are great tools that can reduce the number of subsequent QM calculations. Accordingly, this review aims to briefly discuss the importance and relevance of medicinal chemistry allied to computational chemistry as well as to present a case study where, through a molecular dynamics simulation of AMPK protein (50 ns) and explicit solvent (TIP3P model), a minimum number of snapshots necessary to describe the oscillation profile of the protein behavior was proposed. For this purpose, the RMSD calculation, together with the sophisticated OWSCA method was used to propose the minimum number of snapshots.
- MeSH
- algoritmy MeSH
- farmaceutická chemie MeSH
- kvantová teorie MeSH
- lidé MeSH
- proteinkinasy aktivované AMP metabolismus chemie MeSH
- simulace molekulární dynamiky * MeSH
- Check Tag
- lidé MeSH
- Publikační typ
- časopisecké články MeSH
- přehledy MeSH
The ISWI family protein SMARCA5 contains the ATP-binding pocket that coordinates the catalytic Mg2+ ion and water molecules for ATP hydrolysis. In this study, we demonstrate that SMARCA5 can also possess an alternative metal-binding ability. First, we isolated SMARCA5 on the cobalt column (IMAC) to near homogeneity. Examination of the interactions of SMARCA5 with metal-chelating supports showed that, apart from Co2+, it binds to Cu2+, Zn2+ and Ni2+. The efficiency of the binding to the last-listed metal was influenced by the chelating ligand, resulting in a strong preference for Ni-NTA over the Ni-CM-Asp equivalent. To gain insight in the preferential affinity for the Ni-NTA ligand, QM calculations were performed on model systems and metal-ligand complexes with a limited protein fragment of SMARCA5 containing the double-histidine (dHis) motif. The calculations correlated the observed affinity with the relative stability of the d-block metals to tetradentate ligand coordination over tridentate, as well as their overall octahedral coordination capacity. Likewise, binding free energies derived from model imidazole complexes mirrored the observed Ni-NTA/Ni-CM-Asp preferential affinity. Finally, similar calculations on complexes with a SMARCA5 peptide fragment derived from the AlphaFold structural prediction, captured almost accurately the expected relative stability of the TM complexes, and produced a large energetic separation (~10 kcal∙mol-1) between Ni-NTA and Ni-CM-Asp in favour of the former.
- MeSH
- adenosintrifosfatasy MeSH
- chromozomální proteiny, nehistonové metabolismus chemie MeSH
- kovy chemie metabolismus MeSH
- kvantová teorie MeSH
- lidé MeSH
- ligandy MeSH
- molekulární modely MeSH
- restrukturace chromatinu MeSH
- vazba proteinů * MeSH
- Check Tag
- lidé MeSH
- Publikační typ
- časopisecké články MeSH
Pseudouridine is the most frequently naturally occurring RNA modification, found in all classes of biologically functional RNAs. Compared to uridine, pseudouridine contains an additional hydrogen bond donor group and is therefore widely regarded as a structure stabilizing modification. However, the effects of pseudouridine modifications on the structure and dynamics of RNAs have so far only been investigated in a limited number of different structural contexts. Here, we introduced pseudouridine modifications into the U-turn motif and the adjacent U:U closing base pair of the neomycin-sensing riboswitch (NSR)-an extensively characterized model system for RNA structure, ligand binding, and dynamics. We show that the effects of replacing specific uridines with pseudouridines on RNA dynamics crucially depend on the exact location of the replacement site and can range from destabilizing to locally or even globally stabilizing. By using a combination of NMR spectroscopy, MD simulations and QM calculations, we rationalize the observed effects on a structural and dynamical level. Our results will help to better understand and predict the consequences of pseudouridine modifications on the structure and function of biologically important RNAs.
Phosphate···π, also called anion···π, contacts occur between nucleobases and anionic phosphate oxygens (OP2) in r(GNRA) and r(UNNN) U-turn motifs (N = A,G,C,U; R = A,G). These contacts were investigated using state-of-the-art quantum-chemical methods (QM) to characterize their physicochemical properties and to serve as a reference to evaluate AMBER force field (AFF) performance. We found that phosphate···π interaction energies calculated with the AFF for dimethyl phosphate···nucleobase model systems are less stabilizing in comparison with double-hybrid DFT and that minimum contact distances are larger for all nucleobases. These distance stretches are also observed in large-scale AFF vs QM/MM computations and classical molecular dynamics (MD) simulations on several r(gcGNRAgc) tetraloop hairpins when compared to experimental data extracted from X-ray/cryo-EM structures (res. ≤ 2.5 Å) using the WebFR3D bioinformatic tool. MD simulations further revealed shifted OP2/nucleobase positions. We propose that discrepancies between the QM and AFF result from a combination of missing polarization in the AFF combined with too large AFF Lennard-Jones (LJ) radii of nucleobase carbon atoms in addition to an exaggerated short-range repulsion of the r-12 LJ repulsive term. We compared these results with earlier data gathered on lone pair···π contacts in CpG Z-steps occurring in r(UNCG) tetraloops. In both instances, charge transfer calculations do not support any significant n → π* donation effects. We also investigated thiophosphate···π contacts that showed reduced stabilizing interaction energies when compared to phosphate···π contacts. Thus, we challenge suggestions that the experimentally observed enhanced thermodynamic stability of phosphorothioated r(GNRA) tetraloops can be explained by larger London dispersion.
- MeSH
- fosfáty MeSH
- RNA * chemie MeSH
- simulace molekulární dynamiky * MeSH
- termodynamika MeSH
- výpočetní biologie MeSH
- Publikační typ
- časopisecké články MeSH
Mangiferin is a glycosylated xanthone widely distributed in nature, which exhibits wide pharmacological activities, highlighting its anti-cancer properties. Mangiferin interferes with inflammation, lipid, and calcium signaling, which selectively inhibits multiple NFkB target genes as interleukin-6, tumor necrosis factor, plasminogen, and matrix metalloproteinase, among others. In this work, the interactions of this polyphenol with MMP-9 and NF-κβ are characterized by using computational chemistry methods. The results show MMP-9 inhibition by mangiferina is characterized for the interact with the catalytic Zn atom through a penta-coordinate structure. It is also demonstrated through a strong charge transfer established between mangiferin and Zn in the QM/MM study. Concerning the mangiferin/NF-κβ system, the 92.3% of interactions between p50 sub-unity and DNA are maintained with a binding energy of - 8.04 kcal/mol. These findings indicate that mangiferin blocks the p50-p65/DNA interaction resulting in the loss of the functions of this hetero-dimeric member and suggesting inhibition of the cancer progression. Experimental results concerning the anti-cancer properties of mangiferin show that this natural compound can inhibit selectively MMP-9 and NF-ƙβ. Although the anti-tumor properties of mangiferin are well defined, its molecular mechanisms of actions are not described. In this work, a computational study is carried out to characterize the interactions of mangiferin with these molecular targets. The results obtained corroborate the anti-proliferative and anti-apoptotic activity of mangiferin and provide a depiction of its mechanisms of action.
BACKGROUND AND AIMS: We previously identified subsets of patients with NAFLD with different metabolic phenotypes. Here we align metabolomic signatures with cardiovascular disease (CVD) and genetic risk factors. APPROACH AND RESULTS: We analyzed serum metabolome from 1154 individuals with biopsy-proven NAFLD, and from four mouse models of NAFLD with impaired VLDL-triglyceride (TG) secretion, and one with normal VLDL-TG secretion. We identified three metabolic subtypes: A (47%), B (27%), and C (26%). Subtype A phenocopied the metabolome of mice with impaired VLDL-TG secretion; subtype C phenocopied the metabolome of mice with normal VLDL-TG; and subtype B showed an intermediate signature. The percent of patients with NASH and fibrosis was comparable among subtypes, although subtypes B and C exhibited higher liver enzymes. Serum VLDL-TG levels and secretion rate were lower among subtype A compared with subtypes B and C. Subtype A VLDL-TG and VLDL-apolipoprotein B concentrations were independent of steatosis, whereas subtypes B and C showed an association with these parameters. Serum TG, cholesterol, VLDL, small dense LDL5,6 , and remnant lipoprotein cholesterol were lower among subtype A compared with subtypes B and C. The 10-year high risk of CVD, measured with the Framingham risk score, and the frequency of patatin-like phospholipase domain-containing protein 3 NAFLD risk allele were lower in subtype A. CONCLUSIONS: Metabolomic signatures identify three NAFLD subgroups, independent of histological disease severity. These signatures align with known CVD and genetic risk factors, with subtype A exhibiting a lower CVD risk profile. This may account for the variation in hepatic versus cardiovascular outcomes, offering clinically relevant risk stratification.
- MeSH
- apolipoproteiny B MeSH
- fosfolipasy metabolismus MeSH
- játra patologie MeSH
- kardiovaskulární nemoci * epidemiologie etiologie MeSH
- lipoproteiny VLDL MeSH
- myši MeSH
- nealkoholová steatóza jater * patologie MeSH
- rizikové faktory kardiovaskulárních chorob MeSH
- rizikové faktory MeSH
- triglyceridy metabolismus MeSH
- VLDL-cholesterol metabolismus MeSH
- zvířata MeSH
- Check Tag
- myši MeSH
- zvířata MeSH
- Publikační typ
- časopisecké články MeSH
- práce podpořená grantem MeSH
- Research Support, N.I.H., Extramural MeSH
Inhibition of the human O-linked β-N-acetylglucosaminidase (hOGA, GH84) enzyme is pharmacologically relevant in several diseases such as neurodegenerative and cardiovascular disorders, type 2 diabetes, and cancer. Human lysosomal hexosaminidases (hHexA and hHexB, GH20) are mechanistically related enzymes; therefore, selective inhibition of these enzymes is crucial in terms of potential applications. In order to extend the structure-activity relationships of OGA inhibitors, a series of 2-acetamido-2-deoxy-d-glucono-1,5-lactone sulfonylhydrazones was prepared from d-glucosamine. The synthetic sequence involved condensation of N-acetyl-3,4,6-tri-O-acetyl-d-glucosamine with arenesulfonylhydrazines, followed by MnO2 oxidation to the corresponding glucono-1,5-lactone sulfonylhydrazones. Removal of the O-acetyl protecting groups by NH3/MeOH furnished the test compounds. Evaluation of these compounds by enzyme kinetic methods against hOGA and hHexB revealed potent nanomolar competitive inhibition of both enzymes, with no significant selectivity towards either. The most efficient inhibitor of hOGA was 2-acetamido-2-deoxy-d-glucono-1,5-lactone 1-naphthalenesulfonylhydrazone (5f, Ki = 27 nM). This compound had a Ki of 6.8 nM towards hHexB. To assess the binding mode of these inhibitors to hOGA, computational studies (Prime protein-ligand refinement and QM/MM optimizations) were performed, which suggested the binding preference of the glucono-1,5-lactone sulfonylhydrazones in an s-cis conformation for all test compounds.
- MeSH
- antigeny nádorové chemie metabolismus MeSH
- beta-hexosaminidasa, beta řetězec chemie metabolismus MeSH
- histonacetyltransferasy chemie metabolismus MeSH
- hyaluronoglukosaminidasa chemie metabolismus MeSH
- hydrazony chemická syntéza chemie farmakologie MeSH
- inhibitory enzymů chemická syntéza chemie farmakologie MeSH
- laktony chemie MeSH
- lidé MeSH
- molekulární konformace MeSH
- molekulární modely MeSH
- oxidy chemie MeSH
- sloučeniny manganu chemie MeSH
- vztahy mezi strukturou a aktivitou MeSH
- Check Tag
- lidé MeSH
- Publikační typ
- časopisecké články MeSH
This work aims to synthesize akaganeite nanoparticles (AKNPs) by using microwave and use them to adsorb Congo red dye (CR) from the aqueous solution. The AKNPs with an average particle size of about 50 nm in width and 100 nm in length could be fabricated in 20 min. The effects of pH, CR initial concentration, adsorption time, and adsorbent dosage on the adsorption process were investigated and the artificial neural network (ANN) was used to analyze the adsorption data. The various ANN structures were examined in training the data to find the optimal model. The structure with training function, TRAINLM; adaptation learning function, LARNGDM; transfer function, LOGSIG (in hidden layer) and PURELIN (in output layer); and 10 neutrons in hidden layer having the highest correlation (R2 = 0.996) and the lowest MSE (4.405) is the optimal ANN structure. The consistency between the experimental data and the data predicted by the ANN model showed that the behavior of the adsorption process of CR onto AKNPs under different conditions can be estimated by the ANN model. The adsorption kinetics was studied by fitting the data into pseudo-first-order, pseudo-second-order, Elovich, and intraparticle diffusion models. The results showed that the adsorption kinetics obeyed the pseudo-second-order model and governed by several steps. The adsorption isotherms at the different temperatures were studied by fitting the data to Langmuir, Freundlich, and Temkin isotherm models. The R2 obtained from the Langmuir model was above 0.9 and the highest value in three of four temperatures, suggesting that the adsorption isotherms were the best fit to the Langmuir model and the maximum adsorption capacity was estimated to be more than 150 mg/g. Thermodynamic studies suggested that the adsorption of CR onto AKNPs was a spontaneous and endothermic process and physicochemical adsorption. The obtained results indicated the potential application of microwave-synthesize AKNPs for removing organic dyes from aqueous solutions.
Photosystem II (PSII) is a multi-subunit pigment-protein complex and is one of several protein assemblies that function cooperatively in photosynthesis in plants and cyanobacteria. As more structural data on PSII become available, new questions arise concerning the nature of the charge separation in PSII reaction center (RC). The crystal structure of PSII RC from cyanobacteria Thermosynechococcus vulcanus was selected for the computational study of conformational changes in photosystem II associated to the charge separation process. The parameterization of cofactors and lipids for classical MD simulation with Amber force field was performed. The parametrized complex of PSII was embedded in the lipid membrane for MD simulation with Amber in Gromacs. The conformational behavior of protein and the cofactors directly involved in the charge separation were studied by MD simulations and QM/MM calculations. This study identified the most likely mechanism of the proton-coupled reduction of plastoquinone QB. After the charge separation and the first electron transfer to QB, the system undergoes conformational change allowing the first proton transfer to QB- mediated via Ser264. After the second electron transfer to QBH, the system again adopts conformation allowing the second proton transfer to QBH-. The reduced QBH2 would then leave the binding pocket.
Three-dimensional structure models refined using low-resolution data from crystallographic or electron cryo-microscopy experiments can benefit from high-quality restraints derived from quantum-chemical methods. However, nonperiodic atom-centered quantum-chemistry codes do not inherently account for nearest-neighbor interactions of crystallographic symmetry-related copies in a satisfactory way. Here, these nearest-neighbor effects have been included in the model by expanding to a super-cell and then truncating the super-cell to only include residues from neighboring cells that are interacting with the asymmetric unit. In this way, the fragmentation approach can adequately and efficiently include nearest-neighbor effects. It has previously been shown that a moderately sized X-ray structure can be treated using quantum methods if a fragmentation approach is applied. In this study, a target protein (PDB entry 4gif) was partitioned into a number of large fragments. The use of large fragments (typically hundreds of atoms) is tractable when a GPU-based package such as TeraChem is employed or cheaper (semi-empirical) methods are used. The QM calculations were run at the HF-D3/6-31G level. The models refined using a recently developed semi-empirical method (GFN2-xTB) were compared and contrasted. To validate the refinement procedure for a non-P1 structure, a standard set of crystallographic metrics were used. The robustness of the implementation is shown by refining 13 additional protein models across multiple space groups and a summary of the refinement metrics is presented.
