Parameters for molecular dynamics simulations of iron-sulfur proteins

Multiscale modeling of unfolding and bond dissociation of rubredoxin metalloprotein

Mechanical properties of proteins that have a crucial effect on their operation. This study used a molecular dynamics simulation package to investigate rubredoxin unfolding on the atomic scale. Different simulation techniques were applied, and due to the dissociation of covalent/hydrogen bonds, this protein demonstrates several intermediate states in force-extension behavior. A conceptual model based on the cohesive finite element method was developed to consider the intermediate damages that occur during unfolding. This model is based on force-displacement curves derived from molecular dynamics results. The proposed conceptual model is designed to accurately identify bond rupture points and determine the associated forces. This is achieved by conducting a thorough comparison between molecular dynamics and cohesive finite element results. The utilization of a viscoelastic cohesive zone model allows for the consideration of loading rate effects. This rate-dependent model can be further developed and integrated into the multiscale modeling of large assemblies of metalloproteins, providing a comprehensive understanding of mechanical behavior while maintaining a reduced computational cost.

Computational Characterization of the Inhibition Mechanism of Xanthine Oxidoreductase by Topiroxostat

Xanthine oxidase (XO) is a member of the molybdopterin-containing enzyme family. It interconverts xanthine to uric acid as the last step of purine catabolism in the human body. The high uric acid concentration in the blood directly leads to human diseases like gout and hyperuricemia. Therefore, drugs that inhibit the biosynthesis of uric acid by human XO have been clinically used for many years to decrease the concentration of uric acid in the blood. In this study, the inhibition mechanism of XO and a new promising drug, topiroxostat (code: FYX-051), is investigated by employing molecular dynamics (MD) and quantum mechanics/molecular mechanics (QM/MM) calculations. This drug has been reported to act as both a noncovalent and covalent inhibitor and undergoes a stepwise inhibition by all its hydroxylated metabolites, which include 2-hydroxy-FYX-051, dihydroxy-FYX-051, and trihydroxy-FYX-051. However, the detailed mechanism of inhibition of each metabolite remains elusive and can be useful for designing more effective drugs with similar inhibition functions. Hence, herein we present the computational investigation of the structural and dynamical effects of FYX-051 and the calculated reaction mechanism for all of the oxidation steps catalyzed by the molybdopterin center in the active site. Calculated results for the proposed reaction mechanisms for each metabolite's inhibition reaction in the enzyme's active site, binding affinities, and the noncovalent interactions with the surrounding amino acid residues are consistent with previously reported experimental findings. Analysis of the noncovalent interactions via energy decomposition analysis (EDA) and noncovalent interaction (NCI) techniques suggests that residues L648, K771, E802, R839, L873, R880, R912, F914, F1009, L1014, and A1079 can be used as key interacting residues for further hybrid-type inhibitor development.

Effects of Active-Center Reduction of Plant-Type Ferredoxin on Its Structure and Dynamics: Computational Analysis Using Molecular Dynamics Simulations

"Plant-type" ferredoxins (Fds) in the thylakoid membranes of plants, algae, and cyanobacteria possess a single [2Fe-2S] cluster in active sites and mediate light-induced electron transfer from Photosystem I reaction centers to various Fd-dependent enzymes. Structural knowledge of plant-type Fds is relatively limited to static structures, and the detailed behavior of oxidized and reduced Fds has not been fully elucidated. It is important that the investigations of the effects of active-center reduction on the structures and dynamics for elucidating electron-transfer mechanisms. In this study, model systems of oxidized and reduced Fds were constructed from the high-resolution crystal structure of Chlamydomonas reinhardtii Fd1, and three 200 ns molecular dynamics simulations were performed for each system. The force field parameters of the oxidized and reduced active centers were independently obtained using quantum chemical calculations. There were no substantial differences in the global conformations of the oxidized and reduced forms. In contrast, active-center reduction affected the hydrogen-bond network and compactness of the surrounding residues, leading to the increased flexibility of the side chain of Phe61, which is essential for the interaction between Fd and the target protein. These computational results will provide insight into the electron-transfer mechanisms in the Fds.

Cosolvent Sites-Based Discovery of Mycobacterium Tuberculosis Protein Kinase G Inhibitors

Computer-aided drug discovery methods play a major role in the development of therapeutically important small molecules, but their performance needs to be improved. Molecular dynamics simulations in mixed solvents are useful in understanding protein–ligand recognition and improving molecular docking predictions. In this work, we used ethanol as a cosolvent to find relevant interactions for ligands toward protein kinase G, an essential protein of Mycobacterium tuberculosis (Mtb). We validated the hot spots by screening a database of fragment-like compounds and another one of known kinase inhibitors. Next, we performed a pharmacophore-guided docking simulation and found three low micromolar inhibitors, including one with a novel chemical scaffold that we expanded to four derivative compounds. Binding affinities were characterized by intrinsic fluorescence quenching assays, isothermal titration calorimetry, and the analysis of melting curves. The predicted binding mode was confirmed by X-ray crystallography. Finally, the compounds significantly inhibited the viability of Mtb in infected THP-1 macrophages.

Specificity and Reactivity of Mycobacterium tuberculosis Serine/Threonine Kinases PknG and PknB

Mycobacterium tuberculosis (Mtb), the causative agent of Tuberculosis, has 11 eukaryotic-like serine/threonine protein kinases, which play essential roles in cell growth, signal transduction, and pathogenesis. Protein kinase G (PknG) regulates the carbon and nitrogen metabolism by phosphorylation of the glycogen accumulation regulator (GarA) protein at Thr21. Protein kinase B (PknB) is involved in cell wall synthesis and cell shape, as well as phosphorylates GarA but at Thr22. While PknG seems to be constitutively activated and recognition of GarA requires phosphorylation in its unstructured tail, PknB activation is triggered by phosphorylation of its activation loop, which allows binding of the forkhead-associated domain of GarA. In the present work, we used molecular dynamics and quantum-mechanics/molecular mechanics simulations of the catalytically competent complex and kinase activity assays to understand PknG/PknB specificity and reactivity toward GarA. Two hydrophobic residues in GarA, Val24 and Phe25, seem essential for PknG binding and allow specificity for Thr21 phosphorylation. On the other hand, phosphorylated residues in PknB bind Arg26 in GarA and regulate its specificity for Thr22. We also provide a detailed analysis of the free energy profile for the phospho-transfer reaction and show why PknG has a constitutively active conformation not requiring priming phosphorylation in contrast to PknB. Our results provide new insights into these two key enzymes relevant for Mtb and the mechanisms of serine/threonine phosphorylation in bacteria.

The two redox states of the human NEET proteins' [2Fe-2S] clusters

The NEET proteins constitute a unique class of [2Fe-2S] proteins. The metal ions bind to three cysteines and one histidine. The proteins' clusters exist in two redox states; the oxidized protein (containing two FeIII ions) can transfer the cluster to apo-acceptor protein(s), while the reduced form (containing one ferrous ion) remains bound to the protein frame. Here, we perform in silico and in vitro studies on human NEET proteins in both reduced and oxidized forms. Quantum chemical calculations on all available human NEET proteins structures suggest that reducing the cluster weakens the Fe-NHis and Fe-SCys bonds, similar to what is seen in other Fe-S proteins (e.g., ferredoxin and Rieske protein). We further show that the extra electron in the [2Fe-2S]+ clusters of one of the NEET proteins (mNT) is localized on the His-bound iron ion, consistently with our previous spectroscopic studies. Kinetic measurements demonstrate that the mNT [2Fe-2S]+ is released only by an increase in temperature. Thus, the reduced state of human NEET proteins [2Fe-2S] cluster is kinetically inert. This previously unrecognized kinetic inertness of the reduced state, along with the reactivity of the oxidized state, is unique across all [2Fe-2S] proteins. Finally, using a coevolutionary analysis, along with molecular dynamics simulations, we provide insight on the observed allostery between the loop L2 and the cluster region. Specifically, we show that W75, R76, K78, K79, F82 and G85 in the latter region share similar allosteric characteristics in both redox states.

Electrostatics and water occlusion regulate covalently-bound flavin mononucleotide cofactors of Vibrio cholerae respiratory complex NQR

Proteins like NADH:ubiquinone oxidoreductase (NQR), an essential enzyme and ion pump in the physiology of several pathogenic bacteria, tightly regulate the redox properties of their cofactors. Although flavin mononucleotide (FMN) is fully reduced in aqueous solution, FMN in subunits B and C of NQR exclusively undergo one-electron transitions during its catalytic cycle. Here, we perform ab initio calculations and molecular dynamics simulations to elucidate the mechanisms that regulate the redox state of FMN in NQR. QM/MM calculations show that binding site electrostatics disfavor anionic forms of FMNH₂ , but permit a neutral form of the fully reduced flavin. The potential energy surface is unaffected by covalent bonding between FMN and threonine. Molecular dynamics simulations show that the FMN binding sites are inaccessible by water, suggesting that further reductions of the cofactors are limited or prohibited by the availability of water and other proton donors. These findings provide a deeper understanding of the mechanisms used by NQR to regulate electron transfer through the cofactors and perform its physiologic role. They also provide the first, to our knowledge, evidence of the simple concept that proteins regulate flavin redox states via water occlusion.

Proteins maintain hydration at high [KCl] concentration regardless of content in acidic amino acids

Proteins of halophilic organisms, which accumulate molar concentrations of KCl in their cytoplasm, have a much higher content in acidic amino acids than proteins of mesophilic organisms. It has been proposed that this excess is necessary to maintain proteins hydrated in an environment with low water activity, either via direct interactions between water and the carboxylate groups of acidic amino acids or via cooperative interactions between acidic amino acids and hydrated cations. Our simulation study of five halophilic proteins and five mesophilic counterparts does not support either possibility. The simulations use the AMBER ff14SB force field with newly optimized Lennard-Jones parameters for the interactions between carboxylate groups and potassium ions. We find that proteins with a larger fraction of acidic amino acids indeed have higher hydration levels, as measured by the concentration of water in their hydration shell and the number of water/protein hydrogen bonds. However, the hydration level of each protein is identical at low (b_KCl = 0.15 mol/kg) and high (b_KCl = 2 mol/kg) KCl concentrations; excess acidic amino acids are clearly not necessary to maintain proteins hydrated at high salt concentration. It has also been proposed that cooperative interactions between acidic amino acids in halophilic proteins and hydrated cations stabilize the folded protein structure and would lead to slower dynamics of the solvation shell. We find that the translational dynamics of the solvation shell is barely distinguishable between halophilic and mesophilic proteins; if such a cooperative effect exists, it does not have that entropic signature.

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.

Development of parameters compatible with the CHARMM36 force field for [Fe4S4]2+ clusters and molecular dynamics simulations of adenosine-5'-phosphosulfate reductase in GROMACS 2019

DFT calculations were used to obtain parameters compatible with the CHARMM36 force field for iron-sulfur clusters (Fe-S) of the type [Fe₄S₄]²⁺ that are coordinated to dissimilatory adenosine-5'-phosphosulfate reductase (APSrAB). Classical molecular dynamics (MD) simulations were performed on two APSrAB systems to validate the parameters and verify the stability of the studied systems. The time analysis of the parameters inserted into the force field was in reasonable agreement with the experimental X-ray diffraction data. The analysis of the time evolution of the studied systems indicated that these systems and, in particular, the clusters in their respective cavities had a good stability and were in agreement with what was observed in previous works. The parameters obtained provide the basis for the study of APSrAB as well as other systems that contain [Fe₄S₄]²⁺ through the CHARMM36 force field.

Novel complexes with ONNO tetradentate coumarin schiff-base donor ligands: x-ray structures, DFT calculations, molecular dynamics and potential anticarcinogenic activity

The synthesis of eight novel Zn(II), Co(II), Cu(II), Ni(II) and Pt(II) complexes (2-9) derived from the ONNO tetradentate coumarin Schiff-Base donor ligands, L₁ and the novel L₂, was performed. All compounds were characterized by analytical, spectrometry and spectroscopy techniques. Complexes 2-4 were also characterized by DFT calculations and the structures of 5 and 6 were determined by single-crystal X-ray diffraction analysis. A cytotoxicity study was carried out through an MTT assay in the carcinogenic cell line HeLa and the noncarcinogenic cell lines HFF-1 and HaCaT. The results indicated that among all the evaluated compounds, 2 and 6 presented the best anticarcinogenic potential against HeLa cells with an IC₅₀ of 3.5 and 4.1 µM, respectively. In addition, classical molecular dynamics simulations were performed on the synthesized coordination compounds bound to G4 DNA architectures in the scope of shedding light on their inhibition mode and the most conserved interactions that may lead to the biological activity of the compounds.

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.

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.

Design and Mechanism of (S)-3-Amino-4-(difluoromethylenyl)cyclopent-1-ene-1-carboxylic Acid, a Highly Potent γ-Aminobutyric Acid Aminotransferase Inactivator for the Treatment of Addiction

γ-Aminobutyric acid (GABA) is the major inhibitory neurotransmitter in the central nervous system. Inhibition of GABA aminotransferase (GABA-AT), a pyridoxal 5'-phosphate (PLP)-dependent enzyme that degrades GABA, has been established as a possible strategy for the treatment of substance abuse. The raised GABA levels that occur as a consequence of this inhibition have been found to antagonize the rapid release of dopamine in the ventral striatum (nucleus accumbens) that follows an acute challenge by an addictive substance. In addition, increased GABA levels are also known to elicit an anticonvulsant effect in patients with epilepsy. We previously designed the mechanism-based inactivator (1S,3S)-3-amino-4-difluoromethylenyl-1-cyclopentanoic acid (2), now called CPP-115, that is 186 times more efficient in inactivating GABA-AT than vigabatrin, the only FDA-approved drug that is an inactivator of GABA-AT. CPP-115 was found to have high therapeutic potential for the treatment of cocaine addiction and for a variety of epilepsies, has successfully completed a Phase I safety clinical trial, and was found to be effective in the treatment of infantile spasms (West syndrome). Herein we report the design, using molecular dynamics simulations, synthesis, and biological evaluation of a new mechanism-based inactivator, (S)-3-amino-4-(difluoromethylenyl)cyclopent-1-ene-1-carboxylic acid (5), which was found to be almost 10 times more efficient as an inactivator of GABA-AT than CPP-115. We also present the unexpected crystal structure of 5 bound to GABA-AT, as well as computational analyses used to assist the structure elucidation process. Furthermore, 5 was found to have favorable pharmacokinetic properties and low off-target activities. In vivo studies in freely moving rats showed that 5 was dramatically superior to CPP-115 in suppressing the release of dopamine in the corpus striatum, which occurs subsequent to either an acute cocaine or nicotine challenge. Compound 5 also attenuated increased metabolic demands (neuronal glucose metabolism) in the hippocampus, a brain region that encodes spatial information concerning the environment in which an animal receives a reinforcing or aversive drug. This multidisciplinary computational design to preclinical efficacy approach should be applicable to the design and improvement of mechanism-based inhibitors of other enzymes whose crystal structures and inactivation mechanisms are known.

Molecular Dynamics Simulations of the [2Fe-2S] Cluster-Binding Domain of NEET Proteins Reveal Key Molecular Determinants That Induce Their Cluster Transfer/Release

The NEET proteins are a novel family of iron-sulfur proteins characterized by an unusual three cysteine and one histidine coordinated [2Fe-2S] cluster. Aberrant cluster release, facilitated by the breakage of the Fe-N bond, is implicated in a variety of human diseases, including cancer. Here, the molecular dynamics in the multi-microsecond timescale, along with quantum chemical calculations, on two representative members of the family (the human NAF-1 and mitoNEET proteins), show that the loss of the cluster is associated with a dramatic decrease in secondary and tertiary structure. In addition, the calculations provide a mechanism for cluster release and clarify, for the first time, crucial differences existing between the two proteins, which are reflected in the experimentally observed difference in the pH-dependent cluster reactivity. The reliability of our conclusions is established by an extensive comparison with the NMR data of the solution proteins, in part measured in this work.

Topology and parameter data of thirteen non-natural amino acids for molecular simulations with CHARMM22

In this article we provide a data package containing the topology files and parameters compatible with the CHARMM22 force field for thirteen non-natural amino acids. The force field parameters were derived based on quantum mechanical (QM) calculations involving geometry optimization and potential energy surface scanning at the HF 6-31G(d) and HF 6-311G(d,p) levels of theory. The resulting energy data points were fitted to mathematical functions representing each component of the CHARMM22 force field. Further fine-tuning of the parameters utilized molecular mechanics energies, which were iteratively calculated and compared to the corresponding QM values until the latter were satisfactorily reproduced. The final force field data were validated with molecular dynamics simulations in explicit solvent conditions.

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.

Effect of including torsional parameters for histidine-metal interactions in classical force fields for metalloproteins

Classical force-field parameters of the metal site of metalloproteins usually comprise only the partial charges of the involved atoms, as well as the bond-stretching and bending parameters of the metal-ligand interactions. Although for certain metal ligands such as histidine residues, the torsional motions at the metal site play an important role for the dynamics of the protein, no such terms have been considered to be crucial in the parametrization of the force fields, and they have therefore been omitted in the parametrization. In this work, we have optimized AMBER-compatible force-field parameters for the reduced state of the metal site of copper, zinc superoxide dismutase (SOD1) and assessed the effect of including torsional parameters for the histidine-metal interactions in molecular dynamics simulations. On the basis of the obtained results, we recommend that torsion parameters of the metal site are included when processes at the metal site are investigated or when free-energy calculations are performed. As the torsion parameters mainly affect the structure of the metal site, other kinds of structural studies can be performed without considering the torsional parameters of the metal site.

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.

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.