Parameters for Molecular Dynamics Simulations of Manganese-Containing Metalloproteins

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).

Origin of Enantioselectivity in Engineered Cytochrome c-Catalyzed Carbon-Radical FePP Hydrolysis Revealed Using QM/MM (ABEEM Polarizable Force Field) and MD Simulations

The origin of highly efficient asymmetric aminohydroxylation of styrene catalyzed by engineered cytochrome c is investigated by the developed Atom-Bond Electronegativity Equalization Method polarizable force field (ABEEM PFF), which is a combined outcome of electronic and steric effects. Model molecules were used to establish the charge parameters of the ABEEM PFF, for which the bond-stretching and angle-bending parameters were obtained by using a combination of modified Seminario and scan methods. The interactions between carbon-radical Fe-porphyrin (FePP) and waters are simulated by molecular dynamics, which shows a clear preference for the pre-R over the pre-S. This preference is attributed to the hydrogen-bond between the mutated 100S and 101P residues as well as van der Waals interactions, enforcing a specific conformation of the carbon-radical FePP complex within the binding pocket. Meanwhile, the hydrogen-bond between water and the nitrogen atom in the active intermediate dictates the stereochemical outcome. Quantum mechanics/molecular mechanics (QM/MM (ABEEM PFF)) and free-energy perturbation calculations elucidate that the 3RTS is characterized by sandwich-like structure among adjacent amino acid residues, which exhibits greater stability than crowed arrangement in 3STS and enables the R enantiomer to form more favorably. Thus, this study provides mechanistic insight into the catalytic reaction of hemoproteins.

Insights into the action of the pharmaceutical metformin: Targeted inhibition of the gut microbial enzyme agmatinase

Metformin is the first-line treatment for type 2 diabetes, yet its mechanism of action is not fully understood. Recent studies suggest metformin's interactions with gut microbiota are responsible for exerting therapeutic effects. In this study, we report that metformin targets the gut microbial enzyme agmatinase, as a competitive inhibitor, which may impair gut agmatine catabolism. The metformin inhibition constant (Ki) of E. coli agmatinase is 1 mM and relevant in the gut where the drug concentration is 1-10 mM. Metformin analogs phenformin, buformin, and galegine are even more potent inhibitors of E. coli agmatinase (Ki = 0.6, 0.1, and 0.007 mM, respectively) suggesting a shared mechanism. Agmatine is a known effector of human host metabolism and has been reported to augment metformin's therapeutic effects for type 2 diabetes. This gut-derived inhibition mechanism gives new insights on metformin's action in the gut and may lead to significant discoveries in improving metformin therapy.

The investigation into the effect of the length of RGD peptides and temperature on the interaction with the αIIbβ3 integrin: a molecular dynamic study

The tripeptide Arg-Gly-Asp acid (RGD) is a protein sequence in the binding of proteins to cell surfaces, and is involved in various biological processes such as cell adhesion to the extracellular matrix, platelet activation, hemostasis, etc. The C2 domain of the Von Willebrand Factor (VWF), containing the RGD motif, plays an important role in the initial homeostasis process. It binds to the αIIbβ3 integrin and stimulates platelet aggregation. We have investigated, using the molecular Dynamic (MD) simulation method, the effect of the RGD-peptide length, and temperature variation, on the binding to the αIIbβ3 integrin receptor. We examined 10 different structural modes of the αIIbβ3 at three different temperatures; 237 K, 310 K and 318 K. Our findings show that the amino acids that form a binding pocket include Asp₂₂₄, Tyr₂₃₄, Ser₂₂₆, Tyr₁₉₀, Tyr₁₈₉, Trp₂₆₀, Trp₂₆₂, Asp_259, Lys₂₅₃, Arg₂₁₄, Asp₂₁₇, Ser₁₆₁ and Ala₂₁₈ and that the ligand-receptor interaction was increased at higher temperatures. It was also found that the increase in the number of ligands' amino acids and their types (% glycine) plays an important role in the stability, conformation, and ligand-receptor interaction.Communicated by Ramaswamy H. Sarma.

Divalent-Metal-Ion Selectivity of the CRISPR-Cas System-Associated Cas1 Protein: Insights from Classical Molecular Dynamics Simulations and Electronic Structure Calculations

CRISPR-associated protein 1 (Cas1) is a universally conserved essential metalloenzyme of the clustered regularly interspaced short palindromic repeat (CRISPR) immune system of prokaryotes (bacteria, archaea) that can cut and integrate a part of viral DNA to its host genome with the help of other proteins. The integrated DNA acts as a memory of viral infection, which can be transcribed to RNA and stop future infection by recognition (based on the RNA/DNA complementarity principle) followed by protein-mediated degradation of the viral DNA. It has been proposed that the presence of a single manganese (Mn²⁺) ion in a conserved divalent-metal-ion binding pocket (key residues: E190, H254, D265, D268) of Cas1 is crucial for its function. Cas1-mediated DNA degradation was proposed to be hindered by metal substitution, metal chelation, or mutation of the binding pocket residues. Cas1 is active toward dsDNA degradation with both Mn²⁺ and Mg²⁺. X-ray structures of Cas1 revealed an intricate atomic interaction network of the divalent-metal-ion binding pocket and opened up the possibility of modeling related metal ions (viz., Mg²⁺, Ca²⁺) in the binding pocket of wild-type (WT) and mutated Cas1 proteins for computational analysis, which includes (1) quantitative estimation of the energetics of the divalent-metal-ion preference and (2) exploring the structural and dynamical aspects of the protein in response to divalent-metal-ion substitution or amino acid mutation. Using the X-ray structure of the Cas1 protein from Pseudomonas aeruginosa as a template (PDB 3GOD), we performed (∼2.23 μs) classical molecular dynamics (MD) simulations to compare structural and dynamical differences between Mg²⁺- and Ca²⁺-bound binding pockets of wild-type (WT) and mutant (E190A, H254A, D265A, D268A) Cas1. Furthermore, reduced binding pocket models were generated from X-ray and molecular dynamics (MD) trajectories, and the resulting structures were subjected to quantum chemical calculations. Results suggest that Cas1 prefers Mg²⁺ binding relative to Ca²⁺ and the preference is the strongest for WT and the weakest for the D268A mutant. Quantum chemical calculations indicate that Mn²⁺ is the most preferred relative to both Mg²⁺ and Ca²⁺ in the wild-type and mutant Cas1. Substitution of Mg²⁺ by Ca²⁺ does not alter the interaction network between Cas1 and the divalent metal ion but increases the wetness of the binding pocket by introducing a single water molecule in the first coordination shell of the latter. The strength of metal-ion preference (Mg²⁺ versus Ca²⁺) seems to be dependent on the solvent accessibility of the divalent-metal-ion binding pocket, strongest for wild-type Cas1 (in which the metal-ion binding pocket is dry, which includes two water molecules) and the weakest for the D268A mutant (in which the metal-ion binding pocket is wet, which includes four water molecules).

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.

Developing and Assessing Nonbonded Dummy Models of Magnesium Ion with Different Hydration Free Energy References

A large diversity in the targeted hydration free energies (HFEs) during model parameterization of metal ions was reported in the literature with a difference by dozens of kcal/mol. Here, we developed a series of nonbonded dummy models of the Mg²⁺ ion targeting different HFE references in TIP3P water, followed by assessments of the designed models in the simulations of MgCl₂ solution and biological systems. Together with the comparison of existing models, we conclude that the difference in the targeted HFEs has a limited influence on the model performance, while the usability of these models differs from case to case. The feasibility of reproducing more properties of Mg²⁺ such as diffusion constants and water exchange rates using a nonbonded dummy model is demonstrated. Underestimated activity derivative and osmotic coefficient of MgCl₂ solutions in high concentration reveal a necessity for further optimization of ion-pair interactions. The developed dummy models are applicable to metal coordination with Asp, Glu, and His residues in metalloenzymes, and the performance in predicting monodentate or bidentate binding modes of Asp/Glu residues depends on the complexity of metal centers and the choice of protein force fields. When both the binding modes coexist, the nonbonded dummy models outperform point charge models, probably in need of considering polarization of metal-binding residues by, for instance, charge calibration in classical force fields. This work is valuable for the use and further development of magnesium ion models for simulations of metal-containing systems with good accuracy.

Passive Diffusion of Ciprofloxacin and its Metalloantibiotic: A Computational and Experimental study

Fluoroquinolones (FQ) are antibiotics widely used in clinical practise, but the development of bacterial resistance to these drugs is currently a critical public health problem. In this context, ternary copper complexes of FQ (CuFQPhen) have been studied as a potential alternative. In this study, we compared the passive diffusion across the lipid bilayer of one of the most used FQ, ciprofloxacin (Cpx), and its ternary copper complex, CuCpxPhen, that has shown previous promising results regarding antibacterial activity and membrane partition. A combination of spectroscopic studies and molecular dynamics simulations were used and two different model membranes tested: one composed of anionic phospholipids, and the other composed of zwitterionic phospholipids. The obtained results showed a significantly higher membrane permeabilization activity, larger partition, and a more favourable free energy landscape for the permeation of CuCpxPhen across the membrane, when compared to Cpx. Furthermore, the computational results indicated a more favourable translocation of CuCpxPhen across the anionic membrane, when compared to the zwitterionic one, suggesting a higher specificity towards the former. These findings are important to decipher the influx mechanism of CuFQPhen in bacterial cells, which is crucial for the ultimate use of CuFQPhen complexes as an alternative to FQ to tackle multidrug-resistant bacteria.

Systematic Parameterization and Simulation of Boronic Acid-β-Lactamase Aqueous Solution in Developing the ABEEMσπ Polarizable Force Field

Boronic acid, an inhibitor of β-lactamase, has begun to be applied to the treatment of biological infections and tumors. Scientists are working to develop new and more effective boronic acid. Molecular dynamics (MD) simulation provides a powerful auxiliary tool for drug design. However, the current force fields have no boron-related parameters. In this work, an atom-bond electronegativity equalization method at the σπ level (ABEEMσπ) polarizable force field (ABEEMσπ PFF) of boronic acid and β-lactamase has been developed to determine the potential functions and parameters. The interaction between boron and serine in β-lactamase is regarded as a bonded mode. The interaction between them is simulated by the Morse potential energy function, which is close to the experimental change of the stretching potential energy in a large range. The potential energy surfaces of the bond length, bond angle, and dihedral angle of boronic acid-β-lactamase have the same stability point and change trend as M06-2X/6-311G**. For 47 boronic acid-β-lactamase training molecules, the linear correlation coefficient (R) of the charge distribution between the ABEEMσπ PFF and HF/STO-3G is greater than 0.96. Attributed to the fact that the charge distribution of the ABEEMσπ PFF can fluctuate with the change of geometry and environment, the polarization effect and charge-transfer effect are well reflected. The binding ability of different boronic acids with the same β-lactamase is different. A total of 10 boronic acid-β-lactamase model molecules and 10 boronic acid-β-lactamase and water complexes are simulated. The order of binding energy of five large model molecules calculated by the ABEEMσπ PFF is consistent with that of the MP2 method. The binding energies of boronic acid-β-lactamase and water complexes are close to those of the MP2 method. The results of MD simulation of five aqueous boronic acid-β-lactamase complexes in the NVT ensemble verify the rationality of boron-related parameters of the ABEEMσπ PFF, which have a good application prospect. This study lays a solid theoretical foundation for further study of the inhibition of boronic acid on β-lactamase.

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.

Hydration of Closely Related Manganese and Magnesium Porphyrins in Aqueous Solutions: Ab Initio Quantum Mechanical Charge Field Molecular Dynamics Simulation Study

To the best of our knowledge, the current study based on ab initio quantum mechanical charge field molecular dynamics (QMCF-MD) is the first to explore the difference in the hydration behavior between Mn(II)- and Mg(II)-associated porphyrins (Mn(II)-POR and Mg(II)-POR) in aqueous solution. The simulation study highlights similar and dissimilar characteristics of the structural, dynamical, and thermodynamical properties of these closely related metals bound to porphyrins in aqueous solution. The structural analysis is based on radial and angular distribution functions, coordination number distributions, and angular-radial distributions. Both hydrated systems demonstrate similar pentacoordinated structures formed via the axial coordination of one water molecule to the metal ion in addition to the four nitrogen atoms of the porphyrin ring. However, in the case of Mn(II)-POR, the formation of a distorted square pyramidal geometry was observed. It was envisaged as a weak coordination of the water molecule to the Mn(II) atom and thus higher atomic fluctuation for all atoms in contrast to that for the hydrated Mg(II)-POR. The dynamical data in terms of the mean residence times, velocity autocorrelation function, free energy, and other parameters revealed the difference in the metal binding effect because the Mn(II) atom was observed to inhibit H-bond formation more than the presence of Mg(II) atoms in the core of the porphyrin. The current study thus highlights the significant differences in the structural and dynamical properties of Mn(II)- and Mg(II)-associated porphyrin systems.

The N-glycan structures of the antigenic variants of chlorovirus PBCV-1 major capsid protein help to identify the virus-encoded glycosyltransferases

The chlorovirus Paramecium bursaria chlorella virus 1 (PBCV-1) is a large dsDNA virus that infects the microalga Chlorella variabilis NC64A. Unlike most other viruses, PBCV-1 encodes most, if not all, of the machinery required to glycosylate its major capsid protein (MCP). The structures of the four N-linked glycans from the PBCV-1 MCP consist of nonasaccharides, and similar glycans are not found elsewhere in the three domains of life. Here, we identified the roles of three virus-encoded glycosyltransferases (GTs) that have four distinct GT activities in glycan synthesis. Two of the three GTs were previously annotated as GTs, but the third GT was identified in this study. We determined the GT functions by comparing the WT glycan structures from PBCV-1 with those from a set of PBCV-1 spontaneous GT gene mutants resulting in antigenic variants having truncated glycan structures. According to our working model, the virus gene a064r encodes a GT with three domains: domain 1 has a β-l-rhamnosyltransferase activity, domain 2 has an α-l-rhamnosyltransferase activity, and domain 3 is a methyltransferase that decorates two positions in the terminal α-l-rhamnose (Rha) unit. The a075l gene encodes a β-xylosyltransferase that attaches the distal d-xylose (Xyl) unit to the l-fucose (Fuc) that is part of the conserved N-glycan core region. Last, gene a071r encodes a GT that is involved in the attachment of a semiconserved element, α-d-Rha, to the same l-Fuc in the core region. Our results uncover GT activities that assemble four of the nine residues of the PBCV-1 MCP N-glycans.

New insights about the monomer and homodimer structures of the human AOX1

We conducted MD simulations to provide a comprehensive study on the human aldehyde oxidase and on the impact that the allosteric inhibitor thioridazine and malonate ions have on its structure, particularly on the catalytic tunnel.

Membrane partition of bis-(3-hydroxy-4-pyridinonato) zinc(ii) complexes revealed by molecular dynamics simulations

The class of 3-hydroxy-4-pyridinone ligands is widely known and valuable for biomedical and pharmaceutical purposes. Their chelating properties towards biologically-relevant transition metal ions highlight their potential biomedical utility. A set of 3-hydroxy-4-pyridinone Zn(ii) complexes at different concentrations was studied for their ability to interact with lipid phases. We employed umbrella sampling simulations to attain the potential-of-mean force for a set of ligands and one Zn(ii) complex, as these permeated a 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) hydrated bilayer system. In addition, we used conventional molecular dynamics simulations to study the behavior of various Zn(ii) complexes in hydrated bilayer systems. This work discusses: (i) the partition of 3-hydroxy-4-pyridinone ligands to bilayer phases; (ii) self-aggregation in crowded environments of Zn(ii) complexes; and (iii) possible mechanisms for the membrane translocation of Zn(ii) complexes. We observed distinct interactions for the studied complexes, and distinct membrane partition coefficients (K_mem) depending on the considered ligand. The more hydrophobic ligand, 1-hexyl-3-hydroxy-2-methyl-4(1H)-pyridinone, partitioned more favorably to lipid phases (at least two orders of magnitude higher K_mem when compared to the other ligands), and the corresponding Zn(ii) complex was also prone to self-aggregation when an increased concentration of the complex was employed. We also observed that the inclusion of a coordinated water molecule in the parameterization of the Zn(ii) coordination sphere, as proposed in the available crystallographic structure of the complex, decreased the partition coefficient and membrane permeability for the tested complex.

The membrane partition of hydroxypyridinones and of zinc complexes explored by molecular dynamics.

Substrate-analog binding and electrostatic surfaces of human manganese superoxide dismutase

Superoxide dismutases (SODs) are enzymes that play a key role in protecting cells from toxic oxygen metabolites by disproportionation of two molecules of superoxide into molecular oxygen and hydrogen peroxide via cyclic reduction and oxidation at the active site metal. The azide anion is a potent competitive inhibitor that binds directly to the metal and is used as a substrate analog to superoxide in studies of SOD. The crystal structure of human MnSOD-azide complex was solved and shows the putative binding position of superoxide, providing a model for binding to the active site. Azide is bound end-on at the sixth coordinate position of the manganese ion. Tetrameric electrostatic surfaces were calculated incorporating accurate partial charges for the active site in three states, including a state with superoxide coordinated to the metal using the position of azide as a model. These show facilitation of the anionic ligand to the active site pit via a 'valley' of positively-charged surface patches. Surrounding ridges of negative charge help guide the superoxide anion. Within the active site pit, Arg173 and Glu162 further guide and align superoxide for efficient catalysis. Superoxide coordination at the sixth position causes the electrostatic surface of the active site pit to become nearly neutral. A model for electrostatic-mediated diffusion, and efficient binding of superoxide for catalysis is presented.

Metal Ion Modeling Using Classical Mechanics

Metal ions play significant roles in numerous fields including chemistry, geochemistry, biochemistry, and materials science. With computational tools increasingly becoming important in chemical research, methods have emerged to effectively face the challenge of modeling metal ions in the gas, aqueous, and solid phases. Herein, we review both quantum and classical modeling strategies for metal ion-containing systems that have been developed over the past few decades. This Review focuses on classical metal ion modeling based on unpolarized models (including the nonbonded, bonded, cationic dummy atom, and combined models), polarizable models (e.g., the fluctuating charge, Drude oscillator, and the induced dipole models), the angular overlap model, and valence bond-based models. Quantum mechanical studies of metal ion-containing systems at the semiempirical, ab initio, and density functional levels of theory are reviewed as well with a particular focus on how these methods inform classical modeling efforts. Finally, conclusions and future prospects and directions are offered that will further enhance the classical modeling of metal ion-containing systems.

VFFDT: A New Software for Preparing AMBER Force Field Parameters for Metal-Containing Molecular Systems

Force fields are fundamental to molecular dynamics simulations. However, the incompleteness of force field parameters has been a long-standing problem, especially for metal-related systems. In our previous work, we adopted the Seminario method based on the Hessian matrix to systematically derive the zinc-related force field parameters for AMBER. In this work, in order to further simplify the whole protocol, we have implemented a user-friendly Visual Force Field Derivation Toolkit (VFFDT) to derive the force field parameters via simply clicking on the bond or angle in the 3D viewer, and we have further extended our previous program to support the Hessian matrix output from a variety of quantum mechanics (QM) packages, including Gaussian 03/09, ORCA 3.0, QChem, GAMESS-US, and MOPAC 2009/2012. In this toolkit, a universal VFFDT XYZ file format containing the raw Hessian matrix is available for all of the QM packages, and an instant force field parametrization protocol based on a semiempirical quantum mechanics (SQM) method is introduced. The new function that can automatically obtain the relevant parameters for zinc, copper, iron, etc., which can be exported in AMBER Frcmod format, has been added. Furthermore, our VFFDT program can read and write files in AMBER Prepc, AMBER Frcmod, and AMBER Mol2 format and can also be used to customize, view, copy, and paste the force field parameters in the context of the 3D viewer, which provides utilities complementary to ANTECHAMBER, MCPB, and MCPB.py in the AmberTools.

CHARMM Force-Fields with Modified Polyphosphate Parameters Allow Stable Simulation of the ATP-Bound Structure of Ca(2+)-ATPase

Adenosine triphosphate (ATP) is an indispensable energy source in cells. In a wide variety of biological phenomena like glycolysis, muscle contraction/relaxation, and active ion transport, chemical energy released from ATP hydrolysis is converted to mechanical forces to bring about large-scale conformational changes in proteins. Investigation of structure-function relationships in these proteins by molecular dynamics (MD) simulations requires modeling of ATP in solution and ATP bound to proteins with accurate force-field parameters. In this study, we derived new force-field parameters for the triphosphate moiety of ATP based on the high-precision quantum calculations of methyl triphosphate. We tested our new parameters on membrane-embedded sarcoplasmic reticulum Ca(2+)-ATPase and four soluble proteins. The ATP-bound structure of Ca(2+)-ATPase remains stable during MD simulations, contrary to the outcome in shorter simulations using original parameters. Similar results were obtained with the four ATP-bound soluble proteins. The new force-field parameters were also tested by investigating the range of conformations sampled during replica-exchange MD simulations of ATP in explicit water. Modified parameters allowed a much wider range of conformational sampling compared with the bias toward extended forms with original parameters. A diverse range of structures agrees with the broad distribution of ATP conformations in proteins deposited in the Protein Data Bank. These simulations suggest that the modified parameters will be useful in studies of ATP in solution and of the many ATP-utilizing proteins.

Insights into the key interactions between human protein phosphatase 5 and cantharidin using molecular dynamics and site-directed mutagenesis bioassays

Serine/threonine protein phosphatase 5 (PP5) is a promising novel target for anticancer therapies. This work aims to uncover the key interactions at the atomic level between PP5 and three inhibitors (cantharidin, norcantharidin and endothall). We found that, unlike previous report, Arg 100 contributes less to PP5-inhibitor binding, and the residues His 69, Asn 128, His 129, Arg 225, His 252 and Arg 250 are of importance to PP5-inhibitor binding. The hydrophobic interactions established between the residues Val 254, Phe 271 and Tyr 276, especially Glu 253, are very important to enhance the inhibitive interaction. We suggested that, to increase the inhibitory activity, the interactions of inhibitor with three negatively charged unfavorable interaction residues, Asp 99, Glu 130 and Asp 213, should be avoided. However, the interactions of inhibitor with favorable interaction residue Arg 250 could enhance the inhibitory activity. The Manganese ion 2 (MN2) unfavorably contribute to the total interaction free energies. The coordination between MN2 and chemical group of inhibitor should be eliminated. This work provides insight into how cantharidin and its analogs bind to PP5c at the atomic level and will facilitate modification of cantharidin-like chemicals to rationally develop more specific and less cytotoxic anti-cancer drugs.

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.

Insights into the structural determinants of substrate specificity and activity in mouse aldehyde oxidases

In this work, a combination of homology modeling and molecular dynamics (MD) simulations was used to investigate the factors that modulate substrate specificity and activity of the mouse AOX isoforms: mAOX1, mAOX2 (previously mAOX3l1), mAOX3 and mAOX4. The results indicate that the AOX isoform structures are highly preserved and even more conserved than the corresponding amino acid sequences. The only differences are at the protein surface and substrate-binding site region. The substrate-binding site of all isoforms consists of two regions: the active site, which is highly conserved among all isoforms, and a isoform-specific region located above. We predict that mAOX1 accepts a broader range of substrates of different shape, size and nature relative to the other isoforms. In contrast, mAOX4 appears to accept a more restricted range of substrates. Its narrow and hydrophobic binding site indicates that it only accepts small hydrophobic substrates. Although mAOX2 and mAOX3 are very similar to each other, we propose the following pairs of overlapping substrate specificities: mAOX2/mAOX4 and mAOX3/mAXO1. Based on these considerations, we propose that the catalytic activity between all isoforms should be similar but the differences observed in the binding site might influence the substrate specificity of each enzyme. These results also suggest that the presence of several AOX isoforms in mouse allows them to oxidize more efficiently a wider range of substrates. This contrasts with the same or other organisms that only express one isoform and are less efficient or incapable of oxidizing the same type of substrates.

Force field independent metal parameters using a nonbonded dummy model

The cationic dummy atom approach provides a powerful nonbonded description for a range of alkaline-earth and transition-metal centers, capturing both structural and electrostatic effects. In this work we refine existing literature parameters for octahedrally coordinated Mn²⁺, Zn²⁺, Mg²⁺, and Ca²⁺, as well as providing new parameters for Ni²⁺, Co²⁺, and Fe²⁺. In all the cases, we are able to reproduce both M²⁺–O distances and experimental solvation free energies, which has not been achieved to date for transition metals using any other model. The parameters have also been tested using two different water models and show consistent performance. Therefore, our parameters are easily transferable to any force field that describes nonbonded interactions using Coulomb and Lennard-Jones potentials. Finally, we demonstrate the stability of our parameters in both the human and Escherichia coli variants of the enzyme glyoxalase I as showcase systems, as both enzymes are active with a range of transition metals. The parameters presented in this work provide a valuable resource for the molecular simulation community, as they extend the range of metal ions that can be studied using classical approaches, while also providing a starting point for subsequent parametrization of new metal centers.

Identification of crucial amino acids in mouse aldehyde oxidase 3 that determine substrate specificity

In order to elucidate factors that determine substrate specificity and activity of mammalian molybdo-flavoproteins we performed site directed mutagenesis of mouse aldehyde oxidase 3 (mAOX3). The sequence alignment of different aldehyde oxidase (AOX) isoforms identified variations in the active site of mAOX3 in comparison to other AOX proteins and xanthine oxidoreductases (XOR). Based on the structural alignment of mAOX3 and bovine XOR, differences in amino acid residues involved in substrate binding in XORs in comparison to AOXs were identified. We exchanged several residues in the active site to the ones found in other AOX homologues in mouse or to residues present in bovine XOR in order to examine their influence on substrate selectivity and catalytic activity. Additionally we analyzed the influence of the [2Fe-2S] domains of mAOX3 on its kinetic properties and cofactor saturation. We applied UV-VIS and EPR monitored redox-titrations to determine the redox potentials of wild type mAOX3 and mAOX3 variants containing the iron-sulfur centers of mAOX1. In addition, a combination of molecular docking and molecular dynamic simulations (MD) was used to investigate factors that modulate the substrate specificity and activity of wild type and AOX variants. The successful conversion of an AOX enzyme to an XOR enzyme was achieved exchanging eight residues in the active site of mAOX3. It was observed that the absence of the K889H exchange substantially decreased the activity of the enzyme towards all substrates analyzed, revealing that this residue has an important role in catalysis.