Zinc binding in proteins and solution: a simple but accurate nonbonded representation

Coarse-grained force field for ZIF-8: A study on adsorption, diffusion, and structural properties

Metal-organic frameworks (MOFs) are revolutionizing a spectrum of industries, from groundbreaking gas storage solutions to transformative biological system applications. The intricate architecture of these materials necessitates the use of advanced computational techniques for a comprehensive understanding of their molecular structure and prediction of their physical properties. Coarse-grained (CG) simulations shine a spotlight on the often-neglected influences of defects, pressure effects, and spatial disorders on the performance of MOFs. These simulations are not just beneficial but indispensable for high-demand applications, such as mixed matrix membranes and intricate biological system interfaces. In this work, we propose an optimized CG force field tailored for ZIF-8. Our work provides a deep dive into sorption isotherms and diffusion coefficients of small molecules. We demonstrate the structural dynamics of ZIF-8, particularly how it responds to pressurization, which affects its crystal structure and leads to local changes in aperture size and area. Emphasizing the game-changing potential of CG simulations, we explore the characteristics of amorphization in ZIF-8. Through computational exploration, we aim to bridge the knowledge gap, enhancing the potential applications of nanoporous materials for various applications.

Co-Assembly of Cancer Drugs with Cyclo-HH Peptides: Insights from Simulations and Experiments

Peptide-based nanomaterials can serve as promising drug delivery agents, facilitating the release of active pharmaceutical ingredients while reducing the risk of adverse reactions. We previously demonstrated that Cyclo-Histidine-Histidine (Cyclo-HH), co-assembled with cancer drug Epirubicin, zinc, and nitrate ions, can constitute an attractive drug delivery system, combining drug self-encapsulation, enhanced fluorescence, and the ability to transport the drug into cells. Here, we investigated both computationally and experimentally whether Cyclo-HH could co-assemble, in the presence of zinc and nitrate ions, with other cancer drugs with different physicochemical properties. Our studies indicated that Methotrexate, in addition to Epirubicin and its epimer Doxorubicin, and to a lesser extent Mitomycin-C and 5-Fluorouracil, have the capacity to co-assemble with Cyclo-HH, zinc, and nitrate ions, while a significantly lower propensity was observed for Cisplatin. Epirubicin, Doxorubicin, and Methorexate showed improved drug encapsulation and drug release properties, compared to Mitomycin-C and 5-Fluorouracil. We demonstrated the biocompatibility of the co-assembled systems, as well as their ability to intracellularly release the drugs, particularly for Epirubicin, Doxorubicin, and Methorexate. Zinc and nitrate were shown to be important in the co-assembly, coordinating with drugs and/or Cyclo-HH, thereby enabling drug-peptide as well as drug-drug interactions in successfully formed nanocarriers. The insights could be used in the future design of advanced cancer therapeutic systems with improved properties.

Comparison of the Backbone Dynamics of Dehaloperoxidase-Hemoglobin Isoenzymes

Dehaloperoxidase (DHP) is a multifunctional hemeprotein with a functional switch generally regulated by the chemical class of the substrate. Its two isoforms, DHP-A and DHP-B, differ by only five amino acids and have an almost identical protein fold. However, the catalytic efficiency of DHP-B for oxidation by a peroxidase mechanism ranges from 2- to 6-fold greater than that of DHP-A depending on the conditions. X-ray crystallography has shown that many substrates and ligands have nearly identical binding in the two isoenzymes, suggesting that the difference in catalytic efficiency could be due to differences in the conformational dynamics. We compared the backbone dynamics of the DHP isoenzymes at pH 7 through heteronuclear relaxation dynamics at 11.75, 16.45, and 19.97 T in combination with four 300 ns MD simulations. While the overall dynamics of the isoenzymes are similar, there are specific local differences in functional regions of each protein. In DHP-A, Phe35 undergoes a slow chemical exchange between two conformational states likely coupled to a swinging motion of Tyr34. Moreover, Asn37 undergoes fast chemical exchange in DHP-A. Given that Phe35 and Asn37 are adjacent to Tyr34 and Tyr38, it is possible that their dynamics modulate the formation and migration of the active tyrosyl radicals in DHP-A at pH 7. Another significant difference is that both distal and proximal histidines have a 15-18% smaller S2 value in DHP-B, thus their greater flexibility could account for the higher catalytic activity. The distal histidine grants substrate access to the distal pocket. The greater flexibility of the proximal histidine could also accelerate H2O2 activation at the heme Fe by increased coupling of an amino acid charge relay to stabilize the ferryl Fe(IV) oxidation state in a Poulos-Kraut "push-pull"-type peroxidase mechanism.

Impact of Li, Na and Zn metal cation concentration in EMIM-TFSI ionic liquids on ion clustering, structure and dynamics

We use molecular dynamics calculations to investigate the behavior of metal cations (Li, Na and Zn) within ionic liquids (ILs), specifically EMIM-TFSI, and their impact on key properties, particularly focusing on ion-ion correlations and their influence on diffusion and conductivity. The study explores the competition between metal cations and EMIM ions for binding to TFSI and analyzes ion pair dynamics, revealing that metal cation-TFSI pairs exhibit significantly longer lifetimes compared to TFSI-EMIM pairs. This competitive interaction and the increased stability of metal cation-TFSI pairs at higher concentrations leads to reduced ion exchange, resulting in decreased diffusion and conductivity. The observations underscore the importance of ion size and charge in determining their behavior regarding IL dynamics. Overall, this work provides valuable insights for designing ILs with customized properties, particularly in the context of optimizing conductivity and addressing energy storage challenges.

DP/MM: A Hybrid Model for Zinc-Protein Interactions in Molecular Dynamics

Zinc-containing proteins are vital for many biological processes, yet accurately modeling them using classical force fields is hindered by complicated polarization and charge transfer effects. This study introduces DP/MM, a hybrid force field scheme that utilizes a deep potential model to correct the atomic forces of zinc ions and their coordinated atoms, elevating them from MM to QM levels of accuracy. Trained on the difference between MM and QM atomic forces across diverse zinc coordination groups, the DP/MM model faithfully reproduces structural characteristics of zinc coordination during simulations, such as the tetrahedral coordination of Cys4 and Cys3His1 groups. Furthermore, DP/MM allows water exchange in the zinc coordination environment. With its unique blend of accuracy, efficiency, flexibility, and transferability, DP/MM serves as a valuable tool for studying structures and dynamics of zinc-containing proteins and also represents a pioneering approach in the evolving landscape of machine learning potentials for molecular modeling.

Atomistic molecular simulations of Aβ-Zn conformational ensembles

The amyloid-forming Aβ peptide is able to interact with metal cations to form very stable complexes that influence fibril formation and contribute to the onset of Alzheimer's disease. Multiple structures of peptides derived from Aβ in complex with different metals have been resolved experimentally to provide an atomic-level description of the metal-protein interactions. However, Aβ is intrinsically disordered, and hence more amenable to an ensemble description. Molecular dynamics simulations can now reach the timescales needed to generate ensembles for these type of complexes. However, this requires accurate force fields both for the protein and the protein-metal interactions. Here we use state-of-the-art methods to generate force field parameters for the Zn(II) cations in a set of complexes formed by different Aβ variants and combine them with the Amber99SB*-ILDN optimized force field. Upon comparison of NMR experiments with the simulation results, further optimized with a Bayesian/Maximum entropy approach, we provide an accurate description of the molecular ensembles for most Aβ-metal complexes. We find that the resulting conformational ensembles are more heterogeneous than the NMR models deposited in the Protein Data Bank.

Molecular dynamics of JUNO-IZUMO1 complexation suggests biologically relevant mechanisms in fertilization

JUNO-IZUMO1 binding is the first known physical link created between the sperm and egg membranes in fertilization, however, how this initiates sperm-egg fusion remains elusive. As advanced structural insights will help to combat the infertility crisis, or advance fertility control, we employed all-atom Molecular Dynamics (MD) to derive dynamic structural insights that are difficult to obtain experimentally. We found that the hydrated JUNO-IZUMO1 interface is composed of a large set of short-lived non-covalent interactions. The contact interface is destabilized by strategically located point mutations, as well as by Zn2+ ions, which shift IZUMO1 into the non-binding "boomerang" conformation. We hypothesize that the latter might explain how the transient zinc spark, as released after sperm entry into the oocyte, might contribute to block polyspermy. To address a second mystery, we performed another set of simulations, as it was previously suggested that JUNO in solution is unable to bind to folate despite it belonging to the folate receptor family. MD now suggests that JUNO complexation with IZUMO1 opens up the binding pocket thereby enabling folate insertion. Our MD simulations thus provide crucial new hypotheses how the dynamics of the JUNO-IZUMO1 complex upon solvation might regulate fertility.

Computed Protein-Protein Enthalpy Signatures as a Tool for Identifying Conformation Sampling Problems

Understanding the thermodynamic signature of protein-peptide binding events is a major challenge in computational chemistry. The complexity generated by both components possessing many degrees of freedom poses a significant issue for methods that attempt to directly compute the enthalpic contribution to binding. Indeed, the prevailing assumption has been that the errors associated with such approaches would be too large for them to be meaningful. Nevertheless, we currently have no indication of how well the present methods would perform in terms of predicting the enthalpy of binding for protein-peptide complexes. To that end, we carefully assembled and curated a set of 11 protein-peptide complexes where there is structural and isothermal titration calorimetry data available and then computed the absolute enthalpy of binding. The initial "out of the box" calculations were, as expected, very modest in terms of agreement with the experiment. However, careful inspection of the outliers allows for the identification of key sampling problems such as distinct conformations of peptide termini not being sampled or suboptimal cofactor parameters. Additional simulations guided by these aspects can lead to a respectable correlation with isothermal titration calorimetry (ITC) experiments (R2 of 0.88 and an RMSE of 1.48 kcal/mol overall). Although one cannot know prospectively whether computed ITC values will be correct or not, this work shows that if experimental ITC data are available, then this in conjunction with computed ITC, can be used as a tool to know if the ensemble being simulated is representative of the true ensemble or not. That is important for allowing the correct interpretation of the detailed dynamics of the system with respect to the measured enthalpy. The results also suggest that computational calorimetry is becoming increasingly feasible. We provide the data set as a resource for the community, which could be used as a benchmark to help further progress in this area.

Comparison of force fields to study the zinc-finger containing protein NPL4, a target for disulfiram in cancer therapy

Molecular dynamics (MD) simulations are a powerful approach to studying the structure and dynamics of proteins related to health and disease. Advances in the MD field allow modeling proteins with high accuracy. However, modeling metal ions and their interactions with proteins is still challenging. NPL4 is a zinc-binding protein and works as a cofactor for p97 to regulate protein homeostasis. NPL4 is of biomedical importance and has been proposed as the target of disulfiram, a drug recently repurposed for cancer treatment. Experimental studies proposed that the disulfiram metabolites, bis-(diethyldithiocarbamate)‑copper and cupric ions, induce NPL4 misfolding and aggregation. However, the molecular details of their interactions with NPL4 and consequent structural effects are still elusive. Here, biomolecular simulations can help to shed light on the related structural details. To apply MD simulations to NPL4 and its interaction with copper the first important step is identifying a suitable force field to describe the protein in its zinc-bound states. We examined different sets of non-bonded parameters because we want to study the misfolding mechanism and cannot rule out that the zinc may detach from the protein during the process and copper replaces it. We investigated the force-field ability to model the coordination geometry of the metal ions by comparing the results from MD simulations with optimized geometries from quantum mechanics (QM) calculations using model systems of NPL4. Furthermore, we investigated the performance of a force field including bonded parameters to treat copper ions in NPL4 that we obtained based on QM calculations.

Surfactant Additives Containing Hydrophobic Fluorocarbon Chains and Hydrophilic Sulfonate Anion for Highly Reversible Zn Anode

Aqueous zinc-ion batteries (AZIBs) show enormous potential as a large-scale energy storage technique. However, the growth of Zn dendrites and serious side reactions occurring at the Zn anode hinder the practical application of AZIBs. For the first time, we reported a fluorine-containing surfactant, i.e., potassium perfluoro-1-butanesulfonate (PPFBS), as an additive to the 2 M ZnSO₄ electrolyte. Benefitting from its hydrophilic sulfonate anion and hydrophobic long fluorocarbon chain, PPFBS can promote the uniform distribution of Zn²⁺ flux at the anode/electrolyte interface, allowing the Zn/Zn cell to cycle for 2200 h. Furthermore, PPFBS could inhibit side reactions due to the existence of the perfluorobutyl sulfonate (C₄F₉SO₃^-) adsorption layer and the presence of C₄F₉SO₃^- in the solvation structure of Zn²⁺. The former can reduce the amount of H₂O molecules and SO₄^2- in contact with the Zn anode and C₄F₉SO₃^- entering the Zn²⁺-solvation structure by replacing SO₄^2-. The Zn/Cu cell exhibits a superior average CE of 99.47% over 500 cycles. When coupled with the V₂O₅ cathode, the full cell shows impressive cycle stability. This work provides a simple, effective, and economical solution to the common issues of the Zn anode.

Varying molecular interactions explain aspects of crowder-dependent enzyme function of a viral protease

Biochemical processes in cells, including enzyme-catalyzed reactions, occur in crowded conditions with various background macromolecules occupying up to 40% of cytoplasm's volume. Viral enzymes in the host cell also encounter such crowded conditions as they often function at the endoplasmic reticulum membranes. We focus on an enzyme encoded by the hepatitis C virus, the NS3/4A protease, which is crucial for viral replication. We have previously found experimentally that synthetic crowders, polyethylene glycol (PEG) and branched polysucrose (Ficoll), differently affect the kinetic parameters of peptide hydrolysis catalyzed by NS3/4A. To gain understanding of the reasons for such behavior, we perform atomistic molecular dynamics simulations of NS3/4A in the presence of either PEG or Ficoll crowders and with and without the peptide substrates. We find that both crowder types make nanosecond long contacts with the protease and slow down its diffusion. However, they also affect the enzyme structural dynamics; crowders induce functionally relevant helical structures in the disordered parts of the protease cofactor, NS4A, with the PEG effect being more pronounced. Overall, PEG interactions with NS3/4A are slightly stronger but Ficoll forms more hydrogen bonds with NS3. The crowders also interact with substrates; we find that the substrate diffusion is reduced much more in the presence of PEG than Ficoll. However, contrary to NS3, the substrate interacts more strongly with Ficoll than with PEG crowders, with the substrate diffusion being similar to crowder diffusion. Importantly, crowders also affect the substrate-enzyme interactions. We observe that both PEG and Ficoll enhance the presence of substrates near the active site, especially near catalytic H57 but Ficoll crowders increase substrate binding more than PEG molecules.

MM/GBSA prediction of relative binding affinities of carbonic anhydrase inhibitors: effect of atomic charges and comparison with Autodock4Zn

Carbonic anhydrase is an attractive drug target for the treatment of many diseases. This paper examines the ability of end-state MM/GBSA methods to rank inhibitors of carbonic anhydrase in terms of their binding affinities. The MM/GBSA binding energies were evaluated using different atomic charge schemes (Mulliken, ESP and NPA) at different levels of theories, including Hartree-Fock, B3LYP-D3(BJ), and M06-2X with the 6-31G(d,p) basis set. For a large test set of 32 diverse inhibitors, the use of B3LYP-D3(BJ) ESP atomic charges yielded the strongest correlation with experiment (R² = 0.77). The use of the recently enhanced Autodock Vina and zinc optimised AD4_Zn force field also predicted ligand binding affinities with moderately strong correlation (R² = 0.64) at significantly lower computational cost. However, the docked poses deviate significantly from crystal structures. Overall, this study demonstrates the applicability of docking to estimate ligand binding affinities for a diverse range of CA inhibitors, and indicates that more theoretically robust MM/GBSA simulations show promise for improving the accuracy of predicted binding affinities, as long as a validated set of parameters is used.

Modeling of minimal systems based on ATP-Zn coordination for chemically fueled self-assembly

Following nature's example, there is currently strong interest in using adenosine 5'-triphosphate (ATP) as a fuel for the self-assembly of functional materials with transient/non-equilibrium behaviours. These hold great promise for applications, e.g. in catalysis and drug delivery. In a recent seminal work [Maiti et al., Nat. Chem., 2016, 8, 725], binding of ATP to the metallosurfactant zinc hexadecyl-1,4,7-triazacyclononane ([ZnC₁₆ TACN]²⁺) was exploited to produce ATP-fueled transient vesicles. Crucial to the complex formation is the ability of ATP to bind to the metal ion. As a first step to unveil the key elements underlying this process, we investigate the interaction of ATP with Zn²⁺ and with methyl-1,4,7-triazacyclononane ([ZnCH₃ TACN]²⁺), using all-atom molecular dynamics simulations. The free energy landscape of the complex formation is sampled using well-tempered metadynamics with three collective variables, corresponding to the coordination numbers of Zn²⁺ with the oxygen atoms of the three phosphate groups. We find that the structure of the ternary complex is controlled by direct triphosphate coordination to zinc, with a minor role played by the interactions between ATP and CH₃ TACN which, however, may be important for the build-up of supramolecular assemblies.

Ligands-induced open-close conformational change during DapE catalysis: Insights from molecular dynamics simulations

The microbial enzyme DapE plays a critical role in the lysine biosynthetic pathway and is considered as a potentially safe antibiotic target. In this study, atomistic simulations are employed to identify the modes of essential dynamics that define the conformational response of the enzyme to ligand binding and unbinding. The binding modes and the binding affinities of the products to the DapE enzyme are estimated from the MM-PBSA method, and the residues contributing to the ligand binding are identified. Various structural analyses and the principal component analysis of the molecular dynamics trajectories reveal that the removal of products from the active site causes a significant change in the overall enzyme structure. Both Cartesian and dihedral principal component analyses are used to characterize the structural changes in terms of domain unfolding and domain twisting motions. In the most dominant mode, that is, the domain unfolding motion, the two catalytic domains move away from the two dimerization domains of the dimeric enzyme, representing a closed-to-open conformational change. The conformational changes are initiated by the coordinated movement of three loops (Asp75-Pro82, Gly240-Asn244, and Thr347-Glu353) that trigger a domain-level movement. From multiple short trajectories, the time constant associated with the domain opening motion is estimated as 43.6 ns. Physiologically, this close-to-open conformational change is essential for the regeneration of the initial state of the enzyme for the subsequent cycle of catalytic action and provides the apo enzyme enough flexibility for efficient substrate binding.

Non-bonding energy directed designing of HDAC2 inhibitors through molecular dynamics simulation

Designing an inhibitor having strong affinity in the active site pocket is the cherished goal of structure based drug designing. To achieve this, it is considerably important to predict which structural scaffold is better suited for change to increase affinity. We have explored five HDAC2 co-crystals having PDB ligand code-SHH (vorinostat), LLX, 20Y, IWX (BRD4884) and 6EZ (BRD7232). For analyzing protein-ligand interaction at an atomistic level, we have employed the NAMD molecular dynamics (MD) package. The obtained 100 ns long MD trajectories were subjected to quantitative estimations of non-bonding energies (NBEs) for inferring their interactions with the whole protein or its composite active site (CAS). In addition, relative ΔG_bind was calculated to rank the inhibitors. These inhibitors' NBEs reveal that the phenyl moieties are the major structural scaffold where modifications should be attempted. We designed new compounds (NCs) via introducing hydroxyl groups at 4,5 position of the phenyl moiety of 6EZ, called NC1. Improvement in NC1 further encouraged us for CAP modification by isochromane and isoindoline moieties in place of oxabicyclooctane in NC1, resulting in NC2 and NC3. We also explored trifluoromethyl oxadiazole in 6EZ (NC4 and NC5) and SHH (NC6 and NC7). This moiety acts as a ZBG in NC4 while acting as a part of the foot-pocket in the rest. NC2 and NC6 have highest favorable NBEs among all studied ligands due increased favorable electrostatic contribution. We expect these NBEs data will provide atomistic level insights and benefit in designing new and improved HDAC2 inhibitors. Communicated by Ramaswamy H. Sarma.

Interplay between the Enamine and Imine Forms of the Hydrolyzed Imipenem in the Active Sites of Metallo-β-lactamases and in Water Solution

Deactivation of the β-lactam antibiotics in the active sites of the β-lactamases is among the main mechanisms of bacterial antibiotic resistance. As drugs of last resort, carbapenems are efficiently hydrolyzed by metallo-β-lactamases, presenting a serious threat to human health. Our study reveals mechanistic aspects of the imipenem hydrolysis by bizinc metallo-β-lactamases, NDM-1 and L1, belonging to the B1 and the B3 subclasses, respectively. The results of QM(PBE0-D3/6-31G**)/MM simulations show that the enamine product with the protonated nitrogen atom is formed as the major product in NDM-1 and as the only product in the L1 active site. In NDM-1, there is also another reaction pathway that leads to the formation of the (S)-enantiomer of the imine form of the hydrolyzed imipenem; this process occurs with the higher energy barriers. The absence of the second pathway in L1 is due to the different amino acid composition of the active site loop. In L1, the hydrophobic Pro226 residue is located above the pyrroline ring of imipenem that blocks protonation of the carbon atom. Electron density analysis is performed at the stationary points to compare reaction pathways in L1 and NDM-1. Tautomerization from the enamine to the imine form likely happens in solution after the dissociation of the hydrolyzed imipenem from the active site of the enzyme. Classical molecular dynamics simulations of the hydrolyzed imipenem in solution, both with the neutral enamine and the negatively charged N-C₂-C₃ fragment, demonstrate a huge diversity of conformations. The vast majority of conformations blocks the C₃-atom from the side required for the (S)-imine formation upon tautomerization. Thus, according to our calculations, formation of the (R)-imine is more likely. QM(PBE0-D3/6-31G**)/MM molecular dynamics simulations of the hydrolyzed imipenem with the negatively charged N-C₂-C₃ fragment followed by the Laplacian bond order analysis demonstrate that the N═C₂-C₃^- resonance structure is the most pronounced that facilitates formation of the imine form. The proposed mechanism of the enzymatic enamine formation and its subsequent tautomerization to the imine form in solution is in agreement with the recent spectroscopic and NMR studies.

Adhesion of Bis-Salphen-Based Coordination Polymers to Graphene: Insights from Free Energy Perturbation Study

Manipulation of nanoscale objects using molecular self-assembly is a potent tool to achieve large scale nanopatterning with small effort. Coordination polymers of bis-salphen compounds based on zinc have demonstrated their ability to align carbon nanotubes into micro-scale networks with an unusual “rings-and-rods” pattern. This paper investigates how the compounds interact with pristine and functionalized graphene using density functional theory calculations and molecular dynamic simulations. Using the free energy perturbation method we will show how the addition of phenyl side groups to the core compound and functionalization of graphene affect the stability, mobility and conformation adopted by a dimer of bis-(Zn)salphen compound adsorbed on graphene surface and what it can reveal about the arrangement of chains of bis-(Zn)salphen polymer around carbon nanotubes during the self-assembly of microscale networks.

Absolute binding free energies of mucroporin and its analog mucroporin-M1 with the heptad repeat 1 domain and RNA-dependent RNA polymerase of SARS-CoV-2

The peptide Mucroporin and its analog Mucroporin-M1 were studied using the molecular docking and molecular dynamics simulation of their complexation with two protein targets, the Heptad Repeat 1 (HR1) domain and RNA-dependent RNA polymerase (RdRp) of SARS-CoV-2. The molecular docking of the peptide-protein complexes was performed using the glowworm swarm optimization algorithm. The lowest energy poses were submitted to molecular dynamics simulation. Then, the binding free energies of Mucroporin and its analog Mucroporin-M1 with these two protein targets were calculated using the Multistate Bennett Acceptance Ratio (MBAR) method. It was verified that the peptides/HR1 domain complex showed stability in the interaction site determined by molecular docking. It was also found that Mucroporin-M1 has a much higher affinity than Mucroporin to the HR1 protein target. The peptides showed similar stability and affinity at the NTP binding site in the RdRp protein. Additional experimental studies are needed to confirm the antiviral activity of Mucroporin-M1 and a possible mechanism of action against SARS-CoV-2. However, here we indicate that Mucroporin-M1 may have potential antiviral activity against the HR1 domain with the possibility for further peptide optimization.Communicated by Ramaswamy H. Sarma.

Implications of conformational flexibility, lipid binding, and regulatory domains in cell traversal-protein CelTOS for apicomplexan migration

Malaria and other apicomplexan-caused diseases affect millions of humans, agricultural animals, and pets. Cell traversal is a common feature used by multiple apicomplexan parasites to migrate through host cells and can be exploited to develop therapeutics against these deadly parasites. Here, we provide insights into the mechanism of the Cell-traversal protein for ookinetes and sporozoites (CelTOS), a conserved cell-traversal protein in apicomplexan parasites and malaria vaccine candidate. CelTOS has previously been shown to form pores in cell membranes to enable traversal of parasites through cells. We establish roles for the distinct protein regions of Plasmodium vivax CelTOS and examine the mechanism of pore formation. We further demonstrate that CelTOS dimer dissociation is required for pore formation, as disulfide bridging between monomers inhibits pore formation, and this inhibition is rescued by disulfide-bridge reduction. We also show that a helix-destabilizing amino acid, Pro127, allows CelTOS to undergo significant conformational changes to assemble into pores. The flexible C terminus of CelTOS is a negative regulator that limits pore formation. Finally, we highlight that lipid binding is a prerequisite for pore assembly as mutation of a phospholipids-binding site in CelTOS resulted in loss of lipid binding and abrogated pore formation. These findings identify critical regions in CelTOS and will aid in understanding the egress mechanism of malaria and other apicomplexan parasites as well as have implications for studying the function of other essential pore-forming proteins.

Aromatic Schiff bases confer inhibitory efficacy against New Delhi metallo-β-lactamase-1 (NDM-1)

The irregular use of antibiotics has created a natural selection pressure for bacteria to adapt resistance. Bacterial resistance caused by metallo-β-lactamases (MβLs) has been the most prevalent in terms of posing a threat to human health. The New Delhi metallo-β-lactamase-1 (NDM-1) has been shown to be capable of hydrolyzing almost all β-lactams. In this work, eight aromatic Schiff bases 1-8 were prepared and identified by enzyme kinetic assays to be the potent inhibitors of NDM-1 (except 4). These molecules exhibited a more than 95 % inhibition, and an IC₅₀ value in the range of 0.13-19 μM on the target enzyme, and 3 was found to be the most effective inhibitor (IC₅₀ = 130 nM). Analysis of structure-activity relationship revealed that the o-hydroxy phenyl improved the inhibitory activity of Schiff bases on NDM-1. The inhibition mode assays including isothermal titration calorimetry (ITC) disclosed that both compounds 3 and 5 exhibited a reversibly mixed inhibition on NDM-1, with a K_i value of 1.9 and 10.8 μM, respectively. Antibacterial activity tests indicated that a dose of 64 μg·mL^-1 Schiff bases resulted in 2-128-fold reduction in MICs of cefazolin on E. coli producing NDM-1 (except 4). Cytotoxicity assays showed that both Schiff bases 3 and 5 have low cytotoxicity on the mouse fibroblast (L929) cells at a concentration of up to 400 μM. Docking studies suggested that the hydroxyl group interacts with Gln123 and Glu152 of NDM-1, and the amino groups interact with the backbone amide groups of Glu152 and Asp223. This study provided a novel scaffold for the development of NDM-1 inhibitors.

Force Field Parameterization of Actinyl Molecular Cations Using the 12-6-4 Model

In this work, a set of 12-6-4 force fields (FFs) parameters were developed for the actinyl molecular cations, AnO2n+ (n = 1, 2), from uranium to plutonium for classical molecular dynamics (MD) for four water models: TIP3P, SPC/E, OPC3, and TIP4Pew. Such a non-bonded potential model taking into account the induced dipole between the metallic center and the surrounding molecules has shown better performances for various cations than the classic 12-6 non-bonded potentials. The parametrization method proposed elsewhere for metallic cations has been extended to these molecular cations. In contrast to the actinyl 12-6 FFs from the literature, the new models reproduce correctly both solvation and thermodynamic properties, thanks to the inclusion of the induced dipole term (C4). The transferability of such force fields was assessed by performing MD simulations of carbonato actinyl species, which are highly implicated in actinide migration or actinide extraction from seawater. A highly satisfying agreement was found when comparing the EXAFS signals computed from our MD simulation to the experimental ones. The set of FFs developed here opens new possibilities for the study of actinide chemistry.

Structure-activity relationship and molecular modelling studies of quinazolinedione derivatives MMV665916 as potential antimalarial agent

A series of new quinazolinedione derivatives have been readily synthesized and evaluated for their in vitro antiplasmodial growth inhibition activity. Most of the compounds inhibited P. falciparum FcB1 strain in the low to medium micromolar concentration. The 2-ethoxy 8ag', 2-trifluoromethoxy 8ai' and 4-fluoro-2-methoxy 8ak' showed the best inhibitory activity with EC₅₀ values around 5 µM and were non-toxic to the primary human fibroblast cell line AB943. However, these compounds were less potent than the original hit MMV665916, which showed remarkable growth inhibition with EC₅₀ value of 0.4 µM and presented the highest selectivity index (SI > 250). In addition, a novel approach for determining the docking poses of these quinazolinedione derivatives with their potential protein target, the P. falciparum farnesyltransferase PfFT, was investigated.

Machine Learning-Enabled Pipeline for Large-Scale Virtual Drug Screening

Virtual screening is receiving renewed attention in drug discovery, but progress is hampered by challenges on two fronts: handling the ever-increasing sizes of libraries of drug-like compounds and separating true positives from false positives. Here, we developed a machine learning-enabled pipeline for large-scale virtual screening that promises breakthroughs on both fronts. By clustering compounds according to molecular properties and limited docking against a drug target, the full library was trimmed by 10-fold; the remaining compounds were then screened individually by docking; and finally, a dense neural network was trained to classify the hits into true and false positives. As illustration, we screened for inhibitors against RPN11, the deubiquitinase subunit of the proteasome, and a drug target for breast cancer.

An Iron Shield to Protect Epigallocatehin-3-Gallate from Degradation: Multifunctional Self-Assembled Iron Oxide Nanocarrier Enhances Protein Kinase CK2 Intracellular Targeting and Inhibition

Protein kinase CK2 is largely involved in cell proliferation and apoptosis and is generally recognized as an Achilles’ heel of cancer, being overexpressed in several malignancies. The beneficial effects of (−)-epigallocatechin-3-gallate (EGCG) in the prevention and treatment of several diseases, including cancer, have been widely reported. However, poor stability and limited bioavailability hinder the development of EGCG as an effective therapeutic agent. The combination of innovative nanomaterials and bioactive compounds into nanoparticle-based systems demonstrates the synergistic advantages of nanocomplexes as compared to the individual components. In the present study, we developed a self-assembled core-shell nanohybrid (SAMN@EGCG) combining EGCG and intrinsic dual-signal iron oxide nanoparticles (Surface Active Maghemite Nanoparticles). Interestingly, nano-immobilization on SAMNs protects EGCG from degradation, preventing its auto-oxidation. Most importantly, the nanohybrid was able to successfully deliver EGCG into cancer cells, displaying impressive protein kinase CK2 inhibition comparable to that obtained with the most specific CK2 inhibitor, CX-4945 (5.5 vs. 3 µM), thus promoting the phytochemical exploitation as a valuable alternative for cancer therapy. Finally, to assess the advantages offered by nano-immobilization, we tested SAMN@EGCG against Pseudomonas aeruginosa, a Gram-negative bacterium involved in severe lung infections. An improved antimicrobial effect with a drastic drop of MIC from 500 to 32.7 μM was shown.

Metal-ion promiscuity of microbial enzyme DapE at its second metal-binding site

Metalloenzymes are ubiquitous in nature catalyzing a number of crucial biochemical processes in animal and plant kingdoms. For better adaptation to the relative abundance of different metal ions in different cellular fluids, many of these enzymes exhibit metal promiscuity. The microbial enzyme DapE, an essential enzyme for bacterial growth and survival and a potentially safe target for antibiotics, continues to show enzyme activity when the two zinc ions in its active site are replaced by other transition metal ions. The effect of metal-ion substitution at the second metal-binding site of DapE on its substrate affinity and catalytic efficiency is investigated by QM/MM treatment of the enzyme-substrate complex, by modelling the enzyme with Mn(II), Co(II), Ni(II), or Cu(II) ion in place of Zn(II) at its second metal-binding site, while retaining Zn(II) ion at the first metal-binding site. On the basis of substrate binding energy and activation energy barrier for the chemical catalysis, it is found that Zn-Mn DapE shows poor binding affinity as well as inefficient chemical catalysis. Although Zn-Cu and Zn-Ni DapEs show activation energy barriers comparable to that of wild-type Zn-Zn DapE, their weaker substrate affinity renders these mixed-metal enzymes less efficient. On the other hand, Zn-Co DapE is found to outperform the naturally occurring Zn-Zn DapE, both in terms of substrate affinity and chemical catalysis. The observed metal promiscuity may have played an important role in the survival of bacteria even in those cellular media where Zn ions are in limited supply. Metal nonspecificity in the catalysis of DapE enzyme allows bacteria to thrive in different cellular media.

Force Field Parameters for Fe2+4S2-4 Clusters of Dihydropyrimidine Dehydrogenase, the 5-Fluorouracil Cancer Drug Deactivation Protein: A Step towards In Silico Pharmacogenomics Studies

The dimeric dihydropyrimidine dehydrogenase (DPD), metalloenzyme, an adjunct anti-cancer drug target, contains highly specialized 4 × Fe²⁺₄S²⁻₄ clusters per chain. These clusters facilitate the catalysis of the rate-limiting step in the pyrimidine degradation pathway through a harmonized electron transfer cascade that triggers a redox catabolic reaction. In the process, the bulk of the administered 5-fluorouracil (5-FU) cancer drug is inactivated, while a small proportion is activated to nucleic acid antimetabolites. The occurrence of missense mutations in DPD protein within the general population, including those of African descent, has adverse toxicity effects due to altered 5-FU metabolism. Thus, deciphering mutation effects on protein structure and function is vital, especially for precision medicine purposes. We previously proposed combining molecular dynamics (MD) and dynamic residue network (DRN) analysis to decipher the molecular mechanisms of missense mutations in other proteins. However, the presence of Fe²⁺₄S²⁻₄ clusters in DPD poses a challenge for such in silico studies. The existing AMBER force field parameters cannot accurately describe the Fe²⁺ center coordination exhibited by this enzyme. Therefore, this study aimed to derive AMBER force field parameters for DPD enzyme Fe²⁺ centers, using the original Seminario method and the collation features Visual Force Field Derivation Toolkit as a supportive approach. All-atom MD simulations were performed to validate the results. Both approaches generated similar force field parameters, which accurately described the human DPD protein Fe²⁺₄S²⁻₄ cluster architecture. This information is crucial and opens new avenues for in silico cancer pharmacogenomics and drug discovery related research on 5-FU drug efficacy and toxicity issues.

Parametrization of Trivalent and Tetravalent Metal Ions for the OPC3, OPC, TIP3P-FB, and TIP4P-FB Water Models

Commonly seen in rare-earth chemistry and materials science, highly charged metal ions play key roles in many chemical processes. Computer simulations have become an important tool for scientific research nowadays. Meaningful simulations require reliable parameters. In the present work, we parametrized 18 M(III) and 6 M(IV) metal ions for four new water models (OPC3, OPC, TIP3P-FB, TIP4P-FB) in conjunction with each of the 12-6 and 12-6-4 nonbonded models. Similar to what was observed previously, issues with the 12-6 model can be fixed by using the 12-6-4 model. Moreover, the four new water models showed comparable performance or considerable improvement over the previous water models (TIP3P, SPC/E, and TIP4P_EW) in the same category (3-point or 4-point water models, respectively). Finally, we reported a study of a metalloprotein system demonstrating the capability of the 12-6-4 model to model metalloproteins. The reported parameters will facilitate accurate simulations of highly charged metal ions in aqueous solution.

Benzimidazole and Benzoxazole Zinc Chelators as Inhibitors of Metallo-β-Lactamase NDM-1

Bacterial expression of β-lactamases, which hydrolyze β-lactam antibiotics, contributes to the growing threat of antibacterial drug resistance. Metallo-β-lactamases, such as NDM-1, use catalytic zinc ions in their active sites and hydrolyze nearly all clinically available β-lactam antibiotics. Inhibitors of metallo-β-lactamases are urgently needed to overcome this resistance mechanism. Zinc-binding compounds are promising leads for inhibitor development, as many NDM-1 inhibitors contain zinc-binding pharmacophores. Here, we evaluated 13 chelating agents containing benzimidazole and benzoxazole scaffolds as NDM-1 inhibitors. Six of the compounds showed potent inhibitory activity with IC50 values as low as 0.38 μM, and several compounds restored the meropenem susceptibility of NDM-1-expressing E. coli. Spectroscopic and docking studies suggest ternary complex formation as the mechanism of inhibition, making these compounds promising for development as NDM-1 inhibitors.

Determination of van der Waals Parameters Using a Double Exponential Potential for Nonbonded Divalent Metal Cations in TIP3P Solvent

A double exponential (DE) functional form for Lennard-Jones (LJ) interactions, proposed in our previous study, has many advantages over LJ potentials including a natural softcore characteristic for the convenience of the pathway-based free-energy calculations, fast convergence, and flexibility in use. In this work, we put the first step on the application of the DE functional form by identifying a DE potential, coined DE-TIP3P, for molecular simulations using the TIP3P water model. The developed DE-TIP3 potential was better than LJ potential in reproducing the experimental water properties. Afterward, we developed the nonbonded models of 15 divalent metal ions, which frequently appear and play vital roles in biological systems, to be consistent with the DE-TIP3P potential and TIP3P water model. Our nonbonded models were as good as the complicated nonbonded dummy cationic models by Jiang et al. and the nonbonded 12-6-4 LJ models by Li and Merz in reproducing the experimental properties of those ions. Moreover, our nonbonded models achieved a better performance than the compromise (CM) LJ models and 12-6-4 LJ models, developed by Li and Merz, in reproducing the properties of MgCl₂ in aqueous solution.

Parameterization of Monovalent Ions for the OPC3, OPC, TIP3P-FB, and TIP4P-FB Water Models

Monovalent ions play significant roles in various biological and material systems. Recently, four new water models (OPC3, OPC, TIP3P-FB, and TIP4P-FB), with significantly improved descriptions of condensed phase water, have been developed. The pairwise interaction between the metal ion and water necessitates the development of ion parameters specifically for these water models. Herein, we parameterized the 12-6 and the 12-6-4 nonbonded models for 12 monovalent ions with the respective four new water models. These monovalent ions contain eight cations including alkali metal ions (Li⁺, Na⁺, K⁺, Rb⁺, Cs⁺), transition-metal ions (Cu⁺ and Ag⁺), and Tl⁺ from the boron family, along with four halide anions (F^-, Cl^-, Br^-, I^-). Our parameters were designed to reproduce the target hydration free energies (the 12-6 hydration free energy (HFE) set), the ion-oxygen distances (the 12-6 ion-oxygen distance (IOD) set), or both of them (the 12-6-4 set). The 12-6-4 parameter set provides highly accurate structural features overcoming the limitations of the routinely used 12-6 nonbonded model for ions. Specifically, we note that the 12-6-4 parameter set is able to reproduce experimental hydration free energies within 1 kcal/mol and experimental ion-oxygen distances within 0.01 Å simultaneously. We further reproduced the experimentally determined activity derivatives for salt solutions, validating the ion parameters for simulations of ion pairs. The improved performance of the present water models over our previous parameter sets for the TIP3P, TIP4P, and SPC/E water models (Li, P. et al J. Chem. Theory Comput. 2015 11 1645 1657) highlights the importance of the choice of water model in conjunction with the metal ion parameter set.

Exploring of paritaprevir and glecaprevir resistance due to A156T mutation of HCV NS3/4A protease: molecular dynamics simulation study

Hepatitis C virus (HCV) NS3/4A serine protease is a promising drug target for the discovery of anti-HCV drugs. However, its amino acid mutations, particularly A156T, commonly lead to rapid emergence of drug resistance. Paritaprevir and glecaprevir, the newly FDA-approved HCV drugs, exhibit distinct resistance profiles against the A156T mutation of HCV NS3/4A serine protease. To illustrate their different molecular resistance mechanisms, molecular dynamics simulations and binding free energy calculations were carried out on the two compounds complexed with both wild-type (WT) and A156T variants of HCV NS3/4A protease. QM/MM-GBSA-based binding free energy calculations revealed that the binding affinities of paritaprevir and glecaprevir towards A156T NS3/4A were significantly reduced by ∼4 kcal/mol with respect to their WT complexes, which were in line with the experimental resistance folds. Moreover, the relatively weak intermolecular interactions with amino acids such as H57, R155, and T156 of NS3 protein, the steric effect and the destabilized protein binding surface, which is caused by the loss of salt bridge between R123 and D168, are the main contributions for the higher fold-loss in potency of glecaprevir due to A156T mutation. An insight into the difference of molecular mechanism of drug resistance against the A156T substitution among the two classes of serine protease inhibitors could be useful for further optimization of new generation HCV NS3/4A inhibitors with enhanced inhibitory potency.Communicated by Ramaswamy H. Sarma.

Origin of the Phosphoprotein Phosphatase (PPP) sequence family in Bacteria: Critical ancestral sequence changes, radiation patterns and substrate binding features

Background

Phosphoprotein phosphatases (PPP) belong to the PPP Sequence family, which in turn belongs to the broader metallophosphoesterase (MPE) superfamily. The relationship between the PPP Sequence family and other members of the MPE superfamily remains unresolved, in particular what transitions took place in an ancestral MPE to ultimately produce the phosphoprotein specific phosphatases (PPPs).

Methods

We use structural and sequence alignment data, phylogenetic tree analysis, sequence signature (Weblogo) analysis, in silico protein-peptide modeling data, and in silico mutagenesis to trace a likely route of evolution from MPEs to the PPP Sequence family. Hidden Markov Model (HMM) based iterative database search strategies were utilized to identify PPP Sequence Family members from numerous bacterial groups.

Results

Using Mre11 as proxy for an ancestral nuclease-like MPE we trace a possible evolutionary route that alters a single active site substrate binding His-residue to yield a new substrate binding accessory, the "2-Arg-Clamp". The 2-Arg-Clamp is not found in MPEs, but is present in all PPP Sequence family members, where the phosphomonesterase reaction predominates. Variation in position of the clamp arginines and a supplemental sequence loop likely provide substrate specificity for each PPP Sequence family group.

Conclusions

Loss of a key substrate binding His-in MPEs opened the path to bind novel substrates and evolution of the 2-Arg-Clamp, a sequence change seen in both bacterial and eukaryotic phosphoprotein phosphatases.General significance: We establish a likely evolutionary route from nuclease-like MPE to PPP Sequence family enzymes, that includes the phosphoprotein phosphatases.

Combining Ramachandran plot and molecular dynamics simulation for structural-based variant classification: Using TP53 variants as model

The wide application of new DNA sequencing technologies is generating vast quantities of genetic variation data at unprecedented speed. Developing methodologies to decode the pathogenicity of the variants is imperatively demanding. We hypothesized that as deleterious variants may function through disturbing structural stability of their affected proteins, information from structural change caused by genetic variants can be used to identify the variants with deleterious effects. In order to measure the structural change for proteins with large size, we designed a method named RP-MDS composed of Ramachandran plot (RP) and Molecular Dynamics Simulation (MDS). Ramachandran plot captures the variant-caused secondary structural change, whereas MDS provides a quantitative measure for the variant-caused globular structural change. We tested the method using variants in TP53 DNA binding domain of 219 residues as the model. In total, RP-MDS identified 23 of 38 (60.5%) TP53 known Pathogenic variants and 17 of 42 (41%) TP53 VUS that caused significant changes of P53 structure. Our study demonstrates that RP-MDS method provides a powerful protein structure-based tool to screen deleterious genetic variants affecting large-size proteins.

Role of carboxylic group pattern on protein surface in the recognition of iron oxide nanoparticles: A key for protein corona formation

The knowledge of protein-nanoparticle interplay is of crucial importance to predict the fate of nanomaterials in biological environments. Indeed, protein corona on nanomaterials is responsible for the physiological response of the organism, influencing cell processes, from transport to accumulation and toxicity. Herein, a comparison using four different proteins reveals the existence of patterned regions of carboxylic groups acting as recognition sites for naked iron oxide nanoparticles. Readily interacting proteins display a distinctive surface distribution of carboxylic groups, recalling the geometric shape of an ellipse. This is morphologically complementary to nanoparticles curvature and compatible with the topography of exposed Fe^III sites laying on the nanomaterial surface. The recognition site, absent in non-interacting proteins, promotes the nanoparticle harboring and allows the formation of functional protein coronas. The present work envisages the possibility of predicting the composition and the biological properties of protein corona on metal oxide nanoparticles.

AMBER force field parameters for the Zn (II) ions of the tunneling-fold enzymes GTP cyclohydrolase I and 6-pyruvoyl tetrahydropterin synthase

The folate biosynthesis pathway is an essential pathway for cell growth and survival. Folate derivatives serve as a source of the one-carbon units in several intracellular metabolic reactions. Rapidly dividing cells rely heavily on the availability of folate derivatives for their proliferation. As a result, drugs targeting this pathway have shown to be effective against tumor cells and pathogens, but drug resistance against the available antifolate drugs emerged quickly. Therefore, there is a need to develop new treatment strategies and identify alternative metabolic targets. The two de novo folate biosynthesis pathway enzymes, GTP cyclohydrolase I (GCH1) and 6-pyruvoyl tetrahydropterin synthase (PTPS), can provide an alternative strategy to overcome the drug resistance that emerged in the two primary targeted enzymes dihydrofolate reductase and dihydropteroate synthase. Both GCH1 and PTPS enzymes contain Zn²⁺ ions in their active sites, and to accurately study their dynamic behaviors using all-atom molecular dynamics (MD) simulations, appropriate parameters that can describe their metal sites should be developed and validated. In this study, force field parameters of the GCH1 and PTPS metal centers were generated using quantum mechanics (QM) calculations and then validated through MD simulations to ensure their accuracy in describing and maintaining the Zn²⁺ ion coordination environment. The derived force field parameters will provide accurate and reliable MD simulations involving these two enzymes for future in-silico identification of drug candidates against the GCH1 and PTPS enzymes. Communicated by Ramaswamy H. Sarma.

Investigating the role of N-terminal domain in phosphodiesterase 4B-inhibition by molecular dynamics simulation

Phosphodiesterase 4B (PDE4B) is a potential therapeutic target for the inflammatory respiratory diseases such as congestive obstructive pulmonary disease (COPD) and asthma. The sequence identity of ∼88% with its isoform PDE4D is the key barrier in developing selective PDE4B inhibitors which may help to overcome associated side effects. Despite high sequence identity, both isoforms differ in few residues present in N-terminal (UCR2) and C-terminal (CR3) involved in catalytic site formation. Previously, we designed and tested specific PDE4B inhibitors considering N-terminal residues as a part of the catalytic cavity. In continuation, current work thoroughly presents an MD simulation-based analysis of N-terminal residues and their role in ligand binding. The various parameters viz. root mean square deviation (RMSD), radius of gyration (Rg), root mean square fluctuation (RMSF), principal component analysis (PCA), dynamical cross-correlation matrix (DCCM) analysis, secondary structure analysis and residue interaction mapping were investigated to establish rational. Results showed that UCR2 reduced RMSF values for the metal binding pocket (31.5 ± 11 to 13.12 ± 6 Ų) and the substrate-binding pocket (38.8 ± 32 to 17.3 ± 11 Ų). UCR2 enhanced anti-correlated motion at the active site region that led to the improved ligand-binding affinity of PDE4B from -24.57 ± 3 to -35.54 ± 2 kcal/mol. Further, the atomic-level analysis indicated that T-π and π-π interactions between inhibitors and residues are vital forces that regulate inhibitor association to PDE4B with high affinity. In conclusion, UCR2, the N-terminal domain, embraces the dynamics of PDE4B active site and stabilizes PDE4B inhibitor interactions. Therefore the N-terminal domain needs to be considered while designing next-generation, selective PDE4B-inhibitors as potential anti-inflammatory drugs. Communicated by Ramaswamy H. Sarma.

Discovery and Optimization of Chromone Derivatives as Novel Selective Phosphodiesterase 10 Inhibitors

Phosphodiesterase 10 (PDE10) inhibitors have received much attention as promising therapeutic agents for central nervous system (CNS) disorders such as schizophrenia and Huntington's disease. Recently, a hit compound 1 with a novel chromone scaffold has shown moderate inhibitory activity against PDE10A (IC₅₀ = 500 nM). Hit-to-lead optimization has resulted in compound 3e with an improved inhibitory activity (IC₅₀ = 6.5 nM), remarkable selectivity (>95-fold over other PDEs), and good metabolic stability (RLM t_1/2 = 105 min) by using an integrated strategy (molecular modeling, chemical synthesis, bioassay, and cocrystal structure). The cocrystal structural information provides insights into the binding pattern of 3e in the PDE10A catalytic domain to highlight the key role of the halogen and hydrogen bonds toward Tyr524 and Tyr693, respectively, thereby resulting in high selectivity against other PDEs. These new observations are of benefit for the rational design of the next generation PDE10 inhibitors for CNS disorders.

Improving the binding affinity estimations of protein-ligand complexes using machine-learning facilitated force field method

Scoring functions are routinely deployed in structure-based drug design to quantify the potential for protein-ligand (PL) complex formation. Here, we present a new scoring function Bappl+ that is designed to predict the binding affinities of non-metallo and metallo PL complexes. Bappl+ outperforms other state-of-the-art scoring functions, achieving a high Pearson correlation coefficient of up to ~ 0.76 with low standard deviations. The biggest contributors to the increased performance are the use of a machine-learning model and the enlarged training dataset. We have also evaluated the performance of Bappl+ on target-specific proteins, which highlighted the limitations of our function and provides a way for further improvements. We believe that Bappl+ methodology could prove valuable in ranking candidate molecules against a target metallo or non-metallo protein by reliably predicting their binding affinities, thus helping in the drug discovery process.

AMBER Force Field Parameters for Cobalt-Containing Biological Systems: A Systematic Derivation Study

In the present work, the parameterization of a set of cobalt-containing systems has been performed to create a comprehensive library for bonding parameters of biological Co-containing systems. A standard process for the extraction and validation of parameters was employed, which could be used to create force field parameters for the other metal-containing systems. All protein data banks were searched to extract common chemical groups in bonding with cobalt, and finally, 16 structures were designed to represent the binding model of the chemical moieties with cobalt. The Hessian matrix of each structure was computed at the B3LYP/6-311++G(2d,2p) level of theory and the Seminario method was employed to compute cobalt bond stretching and angle bending parameters. Validation of the derived parameters was performed using structural minimization and molecular dynamics (MD) simulations of four models. Further validation was performed using an extensive MD simulation on carbonic anhydrase II as a common cobalt-containing metalloprotein. The results demonstrated that among models, the bonded model in combination with the RESP charges can produce the most reliable and accurate structural conformations for the metal site of cobalt-containing systems.

Assessing Structural Determinants of Zn2+ Binding to Human HV1 via Multiple MD Simulations

Voltage-gated proton channels (H_V1) are essential for various physiological tasks but are strongly inhibited by Zn²⁺ cations. Some determinants of Zn²⁺ binding have been elucidated experimentally and in computational studies. However, the results have always been interpreted under the assumption that Zn²⁺ binds to monomeric H_V1 despite evidence that H_V1 expresses as a dimer and that the dimer has a higher affinity for zinc than the monomer and experimental data that suggest coordination in the dimer interface. The results of former studies are also controversial, e.g., supporting either one single or two binding sites. Some structural determinants of the binding are still elusive. We performed a series of molecular dynamics simulations to address different structures of the human proton channel, the monomer and two plausible dimer conformations, to compare their respective potential to interact with and bind Zn²⁺ via the essential histidines. The series consisted of several copies of the system to generate independent trajectories and increase the significance compared to a single simulation. The amount of time simulated totals 29.9 μs for 126 simulations of systems comprising ∼59,000 to ∼187,000 atoms. Our approach confirms the existence of two binding sites in monomeric and dimeric human H_V1. The dimer interface is more efficient for attracting and binding Zn²⁺ via the essential histidines than the monomer or a dimer with the histidines in the periphery. The higher affinity is due to the residues in the dimer interface that create an attractive electrostatic potential funneling the zinc cations toward the binding sites.

Binding of divalent cations to acetate: molecular simulations guided by Raman spectroscopy

We combine Raman-MCR vibrational spectroscopy experiments with ab initio and classical MD simulations to gain molecular insights into carboxylate–cation binding.

Understanding of the drug resistance mechanism of hepatitis C virus NS3/4A to paritaprevir due to D168N/Y mutations: A molecular dynamics simulation perspective

Hepatitis C virus (HCV) NS3/4A protease is an attractive target for the development of antiviral therapy. However, the evolution of antiviral drug resistance is a major problem for treatment of HCV infected patients. Understanding of drug-resistance mechanisms at molecular level is therefore very important for the guidance of further design of antiviral drugs with high efficiency and specificity. Paritaprevir is a potent inhibitor against HCV NS3/4A protease genotype 1a. Unfortunately, this compound is highly susceptible to the substitution at D168 in the protease. In this work, molecular dynamics simulations of paritaprevir complexed with wild-type (WT) and two mutated strains (D168 N and D168Y) were carried out. Due to such mutations, the inhibitor-protein hydrogen bonding between them was weakened and the salt-bridge network among residues R123, R155 and D168 responsible for inhibitor binding was disrupted. Moreover, the per-residue free energy decomposition suggested that the main contributions from key residues such as Q80, V132, K136, G137 and R155 were lost in the D168 N/Y mutations. These lead to a lower binding affinity of paritaprevir for D168 N/Y variants of the HCV NS3/4A protease, consistent with the experimental data. This detailed information could be useful for further design of high potency anti-HCV NS3/4A inhibitors.

Dithiocarbamate as a Valuable Scaffold for the Inhibition of Metallo-β-Lactmases

The 'superbug' infection caused by metallo-β-lactamases (MβLs) has grown into an emergent health threat. Given the clinical importance of MβLs, a novel scaffold, dithiocarbamate, was constructed. The obtained molecules, DC1, DC8 and DC10, inhibited MβLs NDM-1, VIM-2, IMP-1, ImiS and L1 from all three subclasses, exhibiting an IC₅₀ < 26 μM. DC1 was found to be the best inhibitor of ImiS (IC₅₀ < 0.22 μM). DC1-2, DC4, DC8 and DC10 restored antimicrobial effects of cefazolin and imipenem against E. coli-BL21, producing NDM-1, ImiS or L1, and DC1 showed the best inhibition of E. coli cells, expressing the three MβLs, resulting in a 2-16-fold reduction in the minimum inhibitory concentrations (MICs) of both antibiotics. Kinetics and isothermal titration calorimetry (ITC) assays showed that DC1 exhibited a reversible, and partially mixed inhibition, of NDM-1, ImiS and L1, with Ki values of 0.29, 0.14 and 5.06 µM, respectively. Docking studies suggest that the hydroxyl and carbonyl groups of DC1 form coordinate bonds with the Zn (II) ions, in the active center of NDM-1, ImiS and L1, thereby inhibiting the activity of the enzymes. Cytotoxicity assays showed that DC1, DC3, DC7 and DC9 have low toxicity in L929 mouse fibroblastic cells, at a dose of up to 250 μM. These studies revealed that the dithiocarbamate is a valuable scaffold for the development of MβLs inhibitors.

Conformation and Domain Movement Analysis of Human Matrix Metalloproteinase-2: Role of Associated Zn2+ and Ca2+ Ions

Matrix metaloproteinase-2 (MMP-2) is an extracellular Zn²⁺ protease specific to type I and IV collagens. Its expression is associated with several inflammatory, degenerative, and malignant diseases. Conformational properties, domain movements, and interactions between MMP-2 and its associated metal ions were characterized using a 1.0 µs molecular dynamics simulation. Dihedral principle component analysis revealed ten families of conformations with the greatest degree of variability occurring in the link region connecting the catalytic and hemopexin domains. Dynamic cross-correlation analysis indicated domain movements corresponding to the opening and closing of the hemopexin domain in relation to the fibronectin and catalytic domains facilitated by the link region. Interaction energies were calculated using the molecular mechanics Poisson Boltzman surface area-interaction entropy (MMPBSA-IE) analysis method and revealed strong binding energies for the catalytic Zn²⁺ ion 1, Ca²⁺ ion 1, and Ca²⁺ ion 3 with significant conformational stability at the binding sites of Zn²⁺ ion 1 and Ca²⁺ ion 1. Ca²⁺ ion 2 diffuses freely away from its crystallographically defined binding site. Zn²⁺ ion 2 plays a minor role in conformational stability of the catalytic domain while Ca²⁺ ion 3 is strongly attracted to the highly electronegative sidechains of the Asp residues around the central β-sheet core of the hemopexin domain; however, the interacting residue sidechain carboxyl groups are outside of Ca²⁺ ion 3′s coordination sphere.