Effective tailor-made force field parameterization of the several Zn coordination environments in the puzzling FTase enzyme: opening the door to the full understanding of its elusive catalytic mechanism

Metallicious: Automated Force-Field Parameterization of Covalently Bound Metals for Supramolecular Structures

Metal ions play a central, functional, and structural role in many molecular structures, from small catalysts to metal-organic frameworks (MOFs) and proteins. Computational studies of these systems typically employ classical or quantum mechanical approaches or a combination of both. Among classical models, only the covalent metal model reproduces both geometries and charge transfer effects but requires time-consuming parameterization, especially for supramolecular systems containing repetitive units. To streamline this process, we introduce metallicious, a Python tool designed for efficient force-field parameterization of supramolecular structures. Metallicious has been tested on diverse systems including supramolecular cages, knots, and MOFs. Our benchmarks demonstrate that parameters accurately reproduce the reference properties obtained from quantum calculations and crystal structures. Molecular dynamics simulations of the generated structures consistently yield stable simulations in explicit solvent, in contrast to similar simulations performed with nonbonded and cationic dummy models. Overall, metallicious facilitates the atomistic modeling of supramolecular systems, key for understanding their dynamic properties and host-guest interactions. The tool is freely available on GitHub (https://github.com/duartegroup/metallicious).

The milk-derived lactoferrin inhibits V-ATPase activity by targeting its V1 domain

Lactoferrin (Lf), a bioactive milk protein, exhibits strong anticancer and antifungal activities. The search for Lf targets and mechanisms of action is of utmost importance to enhance its effective applications. A common feature among Lf-treated cancer and fungal cells is the inhibition of a proton pump called V-ATPase. Lf-driven V-ATPase inhibition leads to cytosolic acidification, ultimately causing cell death of cancer and fungal cells. Given that a detailed elucidation of how Lf and V-ATPase interact is still missing, herein we aimed to fill this gap by employing a five-stage computational approach. Molecular dynamics simulations of both proteins were performed to obtain a robust sampling of their conformational landscape, followed by clustering, which allowed retrieving representative structures, to then perform protein-protein docking. Subsequently, molecular dynamics simulations of the docked complexes and free binding energy calculations were carried out to evaluate the dynamic binding process and build a final ranking based on the binding affinities. Detailed atomist analysis of the top ranked complexes clearly indicates that Lf binds to the V₁ cytosolic domain of V-ATPase. Particularly, our data suggest that Lf binds to the interfaces between A/B subunits, where the ATP hydrolysis occurs, thus inhibiting this process. The free energy decomposition analysis further identified key binding residues that will certainly aid in the rational design of follow-up experimental studies, hence bridging computational and experimental biochemistry.

Molecular Dynamics Studies of Matrix Metalloproteases

Matrix metalloproteases are multidomain enzymes with a remarkable proteolytic activity located in the extracellular environment. Their catalytic activity and structural properties have been intensively studied during the last few decades using both experimental and theoretical approaches, but many open questions still remain. Extensive molecular dynamics simulations enable the sampling of the configurational space of a molecular system, thus contributing to the characterization of the structure, dynamics, and ligand binding properties of a particular MMP. Based on previous computational experience, we provide in this chapter technical and methodological guidelines that may be useful to and stimulate other researchers to perform molecular dynamics simulations to help address unresolved questions concerning the molecular mode of action of MMPs.

Molecular dynamic simulations and structure-based pharmacophore development for farnesyltransferase inhibitors discovery

Farnesyltransferase is one of the enzyme targets for the development of drugs for diseases, including cancer, malaria, progeria, etc. In the present study, the structure-based pharmacophore models have been developed from five complex structures (1LD7, 1NI1, 2IEJ, 2ZIR and 2ZIS) obtained from the protein data bank. Initially, molecular dynamic (MD) simulations were performed for the complexes for 10 ns using AMBER 12 software. The conformers of the complexes (75) generated from the equilibrated protein were undergone protein-ligand interaction fingerprint (PLIF) analysis. The results showed that some important residues, such as LeuB96, TrpB102, TrpB106, ArgB202, TyrB300, AspB359 and TyrB361, are predominantly present in most of the complexes for interactions. These residues form side chain acceptor and surface (hydrophobic or π-π) kind of interactions with the ligands present in the complexes. The structure-based pharmacophore models were generated from the fingerprint bits obtained from PLIF analysis. The pharmacophore models have 3-4 pharmacophore contours consist of acceptor and metal ligation (Acc & ML), hydrophobic (HydA) and extended acceptor (Acc2) features with the radius ranging between 1-3 Å for Acc & ML and 1-2 Å for HydA. The excluded volumes of the pharmacophore contours radius are between 1-2 Å. Further, the distance between the interacting groups, root mean square deviation (RMSD), root mean square fluctuation (RMSF) and radial distribution function (RDF) analysis were performed for the MD-simulated proteins using PTRAJ module. The generated pharmacophore models were used to screen a set of natural compounds and database compounds to select significant HITs. We conclude that the developed pharmacophore model can be a significant model for the identification of HITs as FTase inhibitors.

Re(I) and Tc(I) complexes for targeting nitric oxide synthase: influence of the chelator in the affinity for the enzyme

Aiming to design (99m) Tc complexes for probing nitric oxide synthase (NOS) by SPECT, we synthesized conjugates (L4-L6) comprising a NOS-recognizing moiety connected to a diamino-propionic acid (dap) chelating unit. The conjugates led to complexes of the type fac-[M(CO)3 (ĸ(3) -L)] (M = Re/(99m) Tc; Re4/Tc4: L = L4; Re5/Tc5: L = L5; Re6/Tc6: L = L6). Enzymatic studies showed that L4 and L5, but not L6, gave complexes (Re4 and Re5) that are less potent than the conjugates. To rationalize these results, we performed docking and molecular dynamics simulations. The high affinity of L4 and L5 is due to the strong interactions between the dap chelator and polar residues of the binding cavity. These interactions are hampered by metallation resulting in complexes with lower affinity. The higher potency of Re5 compared to Re4 was assigned to the increased bulkiness of Re5 and the presence of additional anchoring groups that better fit the active site and provide more extensive contacts. In turn, Re6 is too bulky and its organometallic tail is oriented toward the peripheral pocket of iNOS, leading to loss of contacts and a lower affinity. These results were compared with our previous results obtained with analogue complexes stabilized by a pyrazolyl-diamine chelating unit.

Relevant interactions of antimicrobial iron chelators and membrane models revealed by nuclear magnetic resonance and molecular dynamics simulations

The dynamics and interaction of 3-hydroxy-4-pyridinone fluorescent iron chelators, exhibiting antimicrobial properties, with biological membranes were evaluated through NMR and molecular dynamics simulations. Both NMR and MD simulation results support a strong interaction of the chelators with the lipid bilayers that seems to be strengthened for the rhodamine containing compounds, in particular for compounds that include ethyl groups and a thiourea link. For the latter type of compounds the interaction reaches the hydrophobic core of the lipid bilayer. The molecular docking and MD simulations performed for the potential interaction of the chelators with DC-SIGN receptors provide valuable information regarding the cellular uptake of these compounds since the results show that the fluorophore fragment of the molecular framework is essential for an efficient binding. Putting together our previous and present results, we put forward the hypothesis that all the studied fluorescent chelators have access to the cell, their uptake occurs through different pathways and their permeation properties correlate with a better access to the cell and its compartments and, consequently, with the chelators antimicrobial properties.

Comparative analysis of the performance of commonly available density functionals in the determination of geometrical parameters for zinc complexes

A set of 44 Zinc-ligand bond-lengths and of 60 ligand-metal-ligand bond angles from 10 diverse transition-metal complexes, representative of the coordination spheres of typical biological Zn systems, were used to evaluate the performance of a total of 18 commonly available density functionals in geometry determination. Five different basis sets were considered for each density functional, namely two all-electron basis sets (a double-zeta and triple-zeta formulation) and three basis sets including popular types of effective-core potentials: Los Alamos, Steven-Basch-Krauss, and Stuttgart-Dresden. The results show that there are presently several better alternatives to the popular B3LYP density functional for the determination of Zn-ligand bond-lengths and angles. BB1K, MPWB1K, MPW1K, B97-2 and TPSS are suggested as the strongest alternatives for this effect presently available in most computational chemistry software packages. In addition, the results show that the use of effective-core potentials (in particular Stuttgart-Dresden) has a very limited impact, in terms of accuracy, in the determination of metal-ligand bond-lengths and angles in Zinc-complexes, and is a good and safe alternative to the use of an all-electron basis set such as 6-31G(d) or 6-311G(d,p).

Comparative assessment of theoretical methods for the determination of geometrical properties in biological zinc complexes

In the present study, we have compared the performance of the density functional theory (DFT) functionals B1B95, B3LYP, B97-2, BP86, and BPW91 with MP2 for geometry determination in biological mononuclear Zn complexes. A total of 15 different basis sets, of rather diverse complexity, were tested, several which included also three different types of common effective-core potentials: Los Alamos, Steven-Basch-Krauss, and Stuttgart-Dresden. In addition, the ability to describe mononuclear Zn biological systems using relatively simple models of the metal coordination sphere, comprising only the metal atom and a simplified representation of the ligands at the first coordination sphere, starting from a set of high-resolution X-ray crystallographic structures, is evaluated for 90 combinations of method/basis set. The results show that the use of such models allows for a relatively accurate description of the Zn-ligand bond lengths, although failing to correctly represent the topology of the metal coordination sphere (namely, the angles involving the metal atom) if constraints at the Calpha atoms are not considered. Globally, B3LYP had the best average performance in the test, closely followed by MP2, whereas B1B95 was the least accurate method. The study also points out B3LYP/CEP-121G and B3LYP/SDD, which use, respectively, the Steven-Basch-Krauss and the Stuttgart-Dresden effective-core potentials, as the best compromise between accuracy and CPU time for the geometrical characterization of metal-ligand bond lengths in Zn biological systems.

Theoretical studies on farnesyl transferase: Evidence for thioether product coordination to the active-site zinc sphere

Farnesyltransferase (FTase), an interesting zinc metaloenzyme, has been the subject of great attention in anticancer research over the last decade. However, despite the major accomplishments in the field, some very pungent questions on the farnesylation mechanism still persist. In this study, the authors have analyzed a mechanistic paradox that arises from the existence of several contradicting and inconclusive experimental evidence regarding the existence of direct coordination between the active-site zinc cation and the thioether from the farnesylated peptide product, which include UV-vis spectroscopy data on a Co(2+)-substituted FTase, two X-ray crystallographic structures of the FTase-product complex, and extended X-ray absorption fine structure results. Using high-level theoretical calculations on two models of different sizes, and QM/MM calculations on the full enzyme, the authors have shown that the farnesylated product is Zn coordinated, and that a subsequent step where this Zn bond is broken is coherent with the available kinetic results. Furthermore, an explanation for the contradicting experimental evidence is suggested.