- MeSH
- entropie MeSH
- inteligence MeSH
- kvantová teorie MeSH
- lidé MeSH
- mozek fyziologie MeSH
- systémová teorie MeSH
- teoretické modely MeSH
- vědomí * MeSH
- znalosti MeSH
- Check Tag
- lidé MeSH
- Publikační typ
- přehledy MeSH
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
We designed 0D, 1D, and 2D supramolecular assemblies made of diaryliodonium salts (functioning as double σ-hole donors) and carboxylates (as σ-hole acceptors). The association was based on two charge-supported halogen bonds (XB), which occurred between IIII sites of the iodonium cations and the carboxylate anions. The sequential introduction of the carboxylic groups in the aryl ring of the benzoic acid added a dimension to the 0D supramolecular organization of the benzoate, which furnished 1D-chained and 2D-layered structures when terephthalate and trimesate anions, correspondingly, were applied as XB acceptors. The structure-directing XB were studied using DFT calculations under periodic boundary conditions and were followed by the one-electron-potential analysis and the Bader atoms-in-molecules topological analysis of electron density. These theoretical methods confirmed the existence of the XB and verified the philicities of the interaction partners in the designed solid-state structures.
Guanine radical cation (G•+) is a key intermediate in many oxidative processes occurring in nucleic acids. Here, by combining mixed Quantum Mechanical/Molecular Mechanics calculations and Molecular Dynamics (MD) simulations, we study how the structural behaviour of a tract GGG(TTAGGG)3 (hereafter Tel21) of the human telomeric sequence, folded in an antiparallel quadruple helix, changes when one of the G bases is ionized to G•+ (Tel21+). Once assessed that the electron-hole is localized on a single G, we perform MD simulations of twelve Tel21+ systems, differing in the position of G•+ in the sequence. When G•+ is located in the tetrad adjacent to the diagonal loop, we observe substantial structural rearrangements, which can decrease the electrostatic repulsion with the inner Na+ ions and increase the solvent exposed surface of G•+. Analysis of solvation patterns of G•+ provides new insights on the main reactions of G•+, i.e. the deprotonation at two different sites and hydration at the C8 atom, the first steps of the processes producing 8oxo-Guanine. We suggest the main structural determinants of the relative reactivity of each position and our conclusions, consistent with the available experimental trends, can help rationalizing the reactivity of other G-quadruplex topologies.
- MeSH
- DNA chemie MeSH
- G-kvadruplexy * MeSH
- guanin chemie MeSH
- ionty chemie MeSH
- konformace nukleové kyseliny MeSH
- kvantová teorie * MeSH
- lidé MeSH
- molekulární modely MeSH
- oxidační stres * MeSH
- rozpustnost MeSH
- simulace molekulární dynamiky * MeSH
- telomery chemie MeSH
- Check Tag
- lidé MeSH
- Publikační typ
- časopisecké články MeSH
This work extends the multi-scale computational scheme for the quantum mechanics (QM) calculations of Nuclear Magnetic Resonance (NMR) chemical shifts (CSs) in proteins that lack a well-defined 3D structure. The scheme couples the sampling of an intrinsically disordered protein (IDP) by classical molecular dynamics (MD) with protein fragmentation using the adjustable density matrix assembler (ADMA) and density functional theory (DFT) calculations. In contrast to our early investigation on IDPs (Pavlíková Přecechtělová et al., J. Chem. Theory Comput., 2019, 15, 5642-5658) and the state-of-the art NMR calculations for structured proteins, a partial re-optimization was implemented on the raw MD geometries in vibrational normal mode coordinates to enhance the accuracy of the MD/ADMA/DFT computational scheme. In addition, machine-learning based cluster analysis was performed on the scheme to explore its potential in producing protein structure ensembles (CLUSTER ensembles) that yield accurate CSs at a reduced computational cost. The performance of the cluster-based calculations is validated against results obtained with conventional structural ensembles consisting of MD snapshots extracted from the MD trajectory at regular time intervals (REGULAR ensembles). CS calculations performed with the refined MD/ADMA/DFT framework employed the 6-311++G(d,p) basis set that outperformed IGLO-III calculations with the same density functional approximation (B3LYP) and both explicit and implicit solvation. The partial geometry optimization did not universally improve the agreement of computed CSs with the experiment but substantially decreased errors associated with the ensemble averaging. A CLUSTER ensemble with 50 structures yielded ensemble averages close to those obtained with a REGULAR ensemble consisting of 500 MD frames. The cluster based calculations thus required only a fraction of the computational time.
