Theoretical characterization of the shikimate 5-dehydrogenase reaction from Mycobacterium tuberculosis by hybrid QC/MM simulations and quantum chemical descriptors

Exploring the electronic structure of knotted proteins: the case of two ornithine transcarbamylase family

CONTEXT

Geometrical knots are rare structural arrangements in proteins in which the polypeptide chain ties itself into a knot, which is very intriguing due to the uncertainty of their impact on the protein properties. Presently, classical molecular dynamics is the most employed technique in the few studies found on this topic, so any information on how the presence of knots affects the reactivity and electronic properties of proteins is even scarcer. Using the electronic structure methods and quantum chemical descriptors analysis, we found that the same amino-acid residues in the knot core have statistically larger values for the unknotted protein, for both hard-hard and soft-soft interaction descriptors. In addition, we present a computationally feasible protocol, where we show it is possible to separate the contribution of the geometrical knot to the reactivity and other electronic structure properties.

METHODS

In order to investigate these systems, we used PRIMoRDiA, a new software developed by our research group, to explore the electronic structure of biological macromolecules. We evaluated several local quantum chemical descriptors to unveil relevant patterns potentially originating from the presence of the geometrical knot in two proteins, belonging to the ornithine transcarbamylase family. We compared several sampled structures from these two enzymes that are highly similar in both tertiary structure and function, but one of them has a knot whereas the other does not. The sampling was carried out through molecular dynamics simulations using ff14SB force field along 50 ns, and the semiempirical convergence was performed with PM7 Hamiltonian.

Thermochemical and Quantum Descriptor Calculations for Gaining Insight into Ricin Toxin A (RTA) Inhibitors

In this work, we performed a study to assess the interactions between the ricin toxin A (RTA) subunit of ricin and some of its inhibitors using modern semiempirical quantum chemistry and ONIOM quantum mechanics/molecular mechanics (QM/MM) methods. Two approaches were followed (calculation of binding enthalpies, ΔH _bind, and reactivity quantum chemical descriptors) and compared with the respective half-maximal inhibitory concentration (IC₅₀) experimental data, to gain insight into RTA inhibitors and verify which quantum chemical method would better describe RTA-ligand interactions. The geometries for all RTA-ligand complexes were obtained after running classical molecular dynamics simulations in aqueous media. We found that single-point energy calculations of ΔH _bind with the PM6-DH+, PM6-D3H4, and PM7 semiempirical methods and ONIOM QM/MM presented a good correlation with the IC₅₀ data. We also observed, however, that the correlation decreased significantly when we calculated ΔH _bind after full-atom geometry optimization with all semiempirical methods. Based on the results from reactivity descriptors calculations for the cases studied, we noted that both types of interactions, molecular overlap and electrostatic interactions, play significant roles in the overall affinity of these ligands for the RTA binding pocket.

PRIMoRDiA: A Software to Explore Reactivity and Electronic Structure in Large Biomolecules

Plenty of enzymes with structural data do not have their mechanism of catalysis elucidated. Reactivity descriptors, theoretical quantities generated from resolved electronic structure, provide a way to predict and rationalize chemical processes of such systems. In this Application Note, we present PRIMoRDiA (PRIMoRDiA Macromolecular Reactivity Descriptors Access), a software built to calculate the reactivity descriptors of large biosystems by employing an efficient and accurate treatment of the large output files produced by quantum chemistry packages. Here, we show the general implementation details and the software main features. Calculated descriptors were applied for a set of enzymatic systems in order to show their relevance for biological studies and the software potential for use in large scale. Also, we test PRIMoRDiA to aid in the interaction depiction between the SARS-CoV-2 main protease and a potential inhibitor.