Quantum Biochemistry and MM-PBSA Description of the ZIKV NS2B-NS3 Protease: Insights into the Binding Interactions beyond the Catalytic Triad Pocket.

The Zika virus protease NS2B-NS3 has a binding site formed with the participation of a H51-D75-S135 triad presenting two forms, active and inactive. Studies suggest that the inactive conformation is a good target for the design of inhibitors. In this paper, we evaluated the co-crystallized structures of the protease with the inhibitors benzoic acid (5YOD) and benzimidazole-1-ylmethanol (5H4I). We applied a protocol consisting of two steps: first, classical molecular mechanics energy minimization followed by classical molecular dynamics were performed, obtaining stabilized molecular geometries; second, the optimized/relaxed geometries were used in quantum biochemistry and molecular mechanics/Poisson-Boltzmann surface area (MM-PBSA) calculations to estimate the ligand interactions with each amino acid residue of the binding pocket. We show that the quantum-level results identified essential residues for the stabilization of the 5YOD and 5H4I complexes after classical energy minimization, matching previously published experimental data. The same success, however, was not observed for the MM-PBSA simulations. The application of quantum biochemistry methods seems to be more promising for the design of novel inhibitors acting on NS2B-NS3.

The urokinase plasminogen activator binding to its receptor: a quantum biochemistry description within an in/homogeneous dielectric function framework with application to uPA-uPAR peptide inhibitors.

Despite being recognized as a therapeutic target in the processes of cancer cell proliferation and metastasis for over 50 years, the interaction of the urokinase plasminogen activator uPA with its receptor uPAR still needs an improved understanding. High resolution crystallographic data (PDB ) of the uPA-uPAR binding geometry was used to perform quantum biochemistry computations within the density functional theory (DFT) framework. A divide to conquer methodology considering a mixed homogeneous/inhomogeneous dielectric model and explicitly taking water molecules into account was employed to obtain a large set of uPA-uPAR residue-residue interaction energies. In order of importance, not only were Phe25 > Tyr24 > Trp30 > Ile28 shown to be the most relevant uPA residues binding it to uPAR, but the residues Lys98 > His87 > Gln40 > Asn22 > Lys23 > Val20 also had significant interaction energies, which helps to explain published experimental mutational data. Furthermore, the results obtained with the uPA-uPAR in/homogeneous dielectric function show that a high dielectric constant value ε = 40 is adequate to take into account the electrostatic environment at the interface between the proteins, while using a smaller value of ε (<10) leads to an overestimation of the uPA-uPAR binding energy. Hot spots of the uPA-uPAR binding domain were identified and a quantum biochemistry description of the uPAR blockers uPA21-30 and cyclo21,29uPA21-29[(S21C;H29C)] was performed, demonstrating that cyclization improves the stability of mimetic peptides without compromising their binding energies to uPAR.