On the Quest of Dioxygen by Monomeric Sarcosine Oxidase. A Molecular Dynamics Investigation

It is reported here on random acceleration molecular dynamics (RAMD) simulations with the 2GF3 bacterial monomeric sarcosine oxidase (MSOX), O₂, and furoic acid in place of sarcosine, solvated by TIP3 H₂O in a periodic box. An external tiny force, acting randomly on O₂, accelerated its relocation, from the center of activation between residue K265 and the si face of the flavin ring of the flavin adenine dinucleotide cofactor, to the surrounding solvent. Only three of the four O₂ gates previously described for this system along a composite method technique were identified, while two more major O₂ gates were found. The RAMD simulations also revealed that the same gate can be reached by O₂ along different pathways, often involving traps for O₂. Both the residence time of O₂ in the traps, and the total trajectory time for O₂ getting to the solvent, could be evaluated. The new quick pathways discovered here suggest that O₂ exploits all nearby interstices created by the thermal fluctuations of the protein, not having necessarily to look for the permanent large channel used for uptake of the FADH cofactor. To this regard, MSOX resembles closely KijD3 N‐oxygenase. These observations solicit experimental substantiation, in a long term aim at discovering whether gates and pathways for the small gaseous ligands inside the proteins are under Darwinian functional evolution or merely stochastic control operates.     

On the Permeation by Dioxygen of the Cofactor‐Independent Unusual Oxygenase RhCC,in Complex with Substrate 4‐Hydroxyphenylenolpyruvate. A Molecular Dynamics Investigation

This work deals with a trimeric bacterial protein, RhCC, which, although belonging to the tautomerase superfamily, shows oxygenase activity. A model of the complex from RhCC and substrate 4‐hydroxyphenylenolpyruvate (4HPP), fitting the observation of extra electron densities from X‐ray diffraction of the crystal, could be built by autodocking. When subjected to molecular dynamics (MD) aided by an external random force applied to a O₂ molecule placed above 4HPP, this model evolved with O₂ egressing toward the bulk solvent from two nearly opposite gates. These were located between the nearly parallel helices 75 – 91 and 15 – 33 of either chain C (gate SE) or chain B (gate FL). Alternatively, with four O₂ molecules in the bulk solvent, unbiased MD led to O₂ entering the protein from gate SE and getting to 4HPP, while forming a stabilizing salt bridge between the 4HPP carboxylate and P1.C ⁺NH₂, thus providing scientific ground for a refined model of the complex.     

On the Pathways of Biologically Relevant Diatomic Gases through Proteins. Dioxygen and Heme Oxygenase from the Perspective of Molecular Dynamics

This work deals with dioxygen (O₂) binding sites and pathways through inducible human heme oxygenase (HO‐1). The experimentally known distal binding site 1, and sites 2–3 above it, could be reproduced by means of non‐deterministic random‐acceleration molecular‐dynamics (RAMD) simulations. In addition, RAMD revealed the proximal binding site 5, a deeply‐seated binding site 4, which lies behind heme, as well as a few gates communicating with the external medium. In getting from site 1 to the main gate, which lies on the protein front opposed to site 4, O₂ follows chiefly the shortest direct pathway. Less frequently, O₂ visits intermediate sites 2, 4, or 5 along longer pathways. A similarity between HO‐1, myoglobin, and cytoglobin in using, for diatomic gas delivery, the direct shortest pathway from the heme center to the surrounding medium, is emphasized. Otherwise, comparing other proteins and diatomic gases, each system reveals its peculiarities as to sites, gates, and pathways. Thus, relating these properties to the physiological functions of the proteins remains in general a challenge for future studies.     

On Dioxygen and Substrate Access to Soluble Methane Monooxygenases: An all‐Atom Molecular Dynamics Investigation in Water Solution

In a preliminary exploration of the dummy model for diiron proteins, random‐acceleration molecular dynamics (RAMD) revealed that a pure four‐helix bundle structure, like hemerythrin, constitutes an efficient cage for dioxygen (O₂), which can only leave from defined, albeit very broad, gates. However, this well ordered structure does not constitute an archetype on which to compare O₂ permeation of other diiron proteins, like the complex of soluble methane monooxygenase hydroxylase with the regulatory protein (sMMOH‐MMOB). The reason is that with this complex, unlike hemerythrin, the four helices of the four‐helix bundle are heavily bent, and RAMD showed that most traps for O₂ lie outside them. It was also observed that, in spite of a nearly identical van der Waals radius for O₂ and the natural substrate CH₄, the latter behaves under RAMD as a bulkier molecule than O₂, requiring a higher external force to be brought out of sMMOH‐MMOB along trajectories of viable length. All that determined with sMMOH‐MMOB multiple gates and multiple pathways to each of them through several binding pockets, for both O₂ and CH₄. Of the two equally preferred pathways for O₂, at right angle with one another, one proved to be in accordance with the Xe‐atom mapping for sMMOH. In contrast, none of the pathways identified for CH₄ proved to be in accordance with such mapping, CH₄ looking for more open avenues instead.     

Binding Pockets and Pathways for Dioxygen through the KijD3 N‐Oxygenase in Complex with Flavin Mononucleotide Cofactor and a 3‐Aminoglucose Substrate: Predictions from Molecular Dynamics Simulations

In this work, two protein systems, Kij3D? FMN? AKM? O₂ and Kij3D? FMN? O₂, made of KijD3 N‐oxygenase, flavin mononucleotide (FMN) cofactor, dTDP‐3‐amino‐2,3,6‐trideoxy‐4‐keto‐3‐methyl‐D ‐glucose (AKM) substrate, and dioxygen (O₂), have been assembled by adding a molecule of O₂, and removing (or not) AKM, to crystal data for the Kij3D? FMN? AKM complex. Egress of AKM and O₂ from these systems was then investigated by applying a tiny external random force, in turn, to their center of mass in the course of molecular dynamics in explicit H₂O. It turned out that the wide AKM channel, even when emptied, does not constitute the main route for O₂ egress. Other routes appear to be also viable, while various binding pockets (BPs) outside the active center are prone to trap O₂. By reversing the reasoning, these can also be considered as routes for uptake of O₂ by the protein, before or after AKM uptake, while BPs may serve as reservoirs of O₂. This shows that the small molecule O₂ is capable of permeating the protein by exploiting all nearby interstices that are created on thermal fluctuations of the protein, rather than having necessarily to look for farther, permanent channels.     

On Dioxygen Permeation through a Dehydrogenase‐Pyrroloquinoline Quinone Complex. A Molecular‐Dynamics Investigation

In this work, an all atom model of the quinoprotein dehydrogenase PqqC in complex with the PQQ (=4,5‐dihydro‐4,5‐dioxo‐1H‐pyrrolo[2,3‐f]quinoline‐2,7,9‐tricarboxylic acid) cofactor and dioxygen (O₂), solvated with TIP3 water in periodic boxes, was subjected to random‐acceleration molecular dynamics (RAMD). It was found that O₂ leaves the active binding pocket, in front of PQQ, to get to the solvent, as easily as with a variety of other O₂‐activating enzymes, O₂ carriers, and gas‐sensing proteins. The shortest pathway, orthogonal to the center of the mean plane of PQQ, was largely preferred by O₂ over pathways slightly deviating from this line. These observations challenge the interpretation of an impermeable active binding pocket of PqqC‐PQQ, as drawn from both X‐ray diffraction data of the crystal at low temperature and physiological experimentation.     

Unveiling the Pathways of Dioxygen Through the C2 Component of the Environmentally Relevant Monooxygenase p‐Hydroxyphenylacetate Hydroxylase from Acinetobacter baumannii: A Molecular Dynamics Investigation

In this work, models of the homotetrameric C2 component of the monooxygenase p‐hydroxyphenylacetate hydroxylase from Acinetobacter baumannii, in complex with dioxygen (O₂) and, or not, the substrate p‐hydroxyphenylacetate (HPA) were built. Both models proved to be amenable to random‐acceleration molecular dynamics (RAMD) simulations, whereby a tiny randomly oriented external force, acting on O₂ at the active site in front of flavin mononucleotide (FMNH^?), accelerated displacement of O₂ toward the bulk solvent. This allowed us to carry out a sufficiently large number of RAMD simulations to be of statistical significance. The two systems behaved very similarly under RAMD, except for O₂ leaving the active site more easily in the absence of HPA, but then finding similar obstacles in getting to the gate as when the active site was sheltered by HPA. This challenges previous conclusions that HPA can only reach the active center after that the C4aOOH derivative of FMNH^? is formed, requiring uptake of O₂ at the active site before HPA. According to these RAMD simulations, O₂ could well get to FMNH^? also in the presence of the substrate at the active site.     

Molecular Dynamics Simulations of a Calmodulin-peptide Complex in Solution

Abstract

We have performed an 4-ns MD simulation of calmodulin complexed with a target peptide in explicit water, under realistic conditions of constant temperature and pressure, in the presence of a physiological concentration of counterions and using Ewald summation to avoid truncation of long-range electrostatic forces. During the simulation the system tended to perform small fluctuations around a structure similar to, but somewhat looser than the starting crystal structure. The calmodulin-peptide complex was quite rigid and did not exhibit any large amplitude domain motions such as previously seen in apo- and calcium-bound calmodulin. We analyzed the calmodulin-peptide interactions by calculating buried surface areas, CHARMM interaction energies and continuum model interaction free energies. In the trajectory, the protein surface area buried by contact with the peptide is 1373 Ų, approximately evenly divided between the calmodulin N-terminal, C-terminal and central linker regions. A majority of this buried surface, 803 ·A², comes from nonpolar residues, in contrast to the protein as a whole, for which the surface is made up of mostly polar and charged groups. Our continuum calculations indicate that the largest favorable contribution to pep- tide binding comes from burial of molecular surface upon complex formation. Electrostatic contributions are favorable but smaller in the trajectory structures, and actually unfavorable for binding in the crystal structure. Since nonpolar groups make up most of buried surface of the protein, our calculations suggest that the hydrophobic effect is the main driving force for binding the helical peptide to calmodulin, consistent with thermodynamic analysis of experimental data. Besides the burial of nonpolar surface area, secondary contributions to peptide binding come from burial of polar surface and electrostatic interactions. In the nonpolar interactions a crucial role is played by the nine methionines of calmodulin. In the electrostatic interactions the negatively charged protein residues and positively charged peptide residues play a dominant role.     

On the Atomic Velocities in Molecular and Langevin Dynamics Simulations of Soft-Sphere Systems

Abstract

Time dependent probability distributions of the changes of direction of atomic velocities are considered in order to examine in detail the shape of the trajectories obtained through molecular simulations. We have analysed the atomic motions obtained from molecular dynamics simulations of soft-sphere systems at three very different states, i.e. a dilute fluid, a liquid at high density, and a solid. The methodology has also been used to check the reliability of the velocity evolution obtained when it is assumed that a single particle obeys the generalized Langevin equation and the effect of the other particles is represented by friction and random forces.     

Molecular Dynamics Simulations of Proteins in Water Without the Truncation of Long-range Coulomb Interactions

A new program package (COSMOS90) for molecular dynamics simulations was developed to simulate large molecular systems consisting of more than tens of thousands of atoms without the truncation of long-range coulomb interactions. This program package was based on a new approximation scheme (PPPC) for calculating efficiently the coulomb interactions without sacrificing accuracy. In this approximation scheme, the group of charges at a long distance from each atom was represented by a total charge and total dipole moment of the group. In order to assess the accuracy of PPPC and the ability of COSMOS90, molecular dynamics simulations were carried out for a large system consisting of 16108 atoms (human lysozyme in water) for 50 ps using this program package. The coulomb energy per solute atom was calculated with only five percent of the error found in the 10 Å cut-off approximation (about 0.9 kcal/mol versus 18 kcal/mol, respectively). The molecular dynamics simulations using COSMOS90 require no more CPU time than the simulations based on the 10 Å cut-off approximation of the conventional programs for macromolecular simulations.     

On Dioxygen Permeation of MaL Laccase from the Thermophilic Fungus Melanocarpus albomyces: An all‐Atom Molecular Dynamics Investigation

In this article, biased molecular dynamics (MD) simulations of O₂ egress from the active center of MaL laccase toward the bulk solvent were described. Parameterization of the set of four Cu(II) ions, in the framework of CHARMM‐36 FF, was carried out on a recent dummy‐atom model that takes into account the Jahn–Teller effect. By carrying out a number of statistically relevant MD simulations, under a tiny randomly oriented external force applied to the molecule O₂, three preferred gates for O₂ egress from the enzyme and a few intermediate binding pockets (BPs) for O₂ were visible; all the gates and pockets were located on two of the three domains. This wide distribution of preferred gates notwithstanding the molecule O₂ was seen to follow specific pathways, exploiting consistently the interstices created by the enzyme thermal fluctuations. These are features that can be imagined to have evolved to make MaL laccase extremely efficient as catalysts in various reactions that require O₂.     

Reactivity of an Excess Electron with Monovalent Cations in Bulk Water by Mixed Quantum Classical Molecular Dynamics Simulations

A mixed quantum classical molecular dynamics (QCMD) simulation of the silver and sodium cations in presence of an excess electron is reported. The silver cation is shown to be reduced by the hydrated electron and to form a stable, highly polarized, neutral atom. On the contrary, the sodium cation is not reduced and a metastable contact ion pair is observed. The resulting absorption spectra of both species are compared with experiments and shown to be in good agreement. Furthermore, the free energy curve for the charge separation was calculated and rationalized in terms of a thermodynamic cycle. Finally, a direct, reactive, molecular dynamics trajectory provides some useful informations on the reduction mechanism.     

Molecular Dynamics Simulation of Dioxygen Pathways through Mini Singlet Oxygen Generator (miniSOG), a Genetically Encoded Marker and Killer Protein

In this work, molecular dynamics (MD) simulations of the permeation of proteins by small gases of biological significance have been extended from gas carrier, sensor, and enzymatic proteins to genetically encoded tags and killer proteins. To this end, miniSOG was taken as an example of current high interest, using a biased form of MD, called random‐acceleration MD. Various egress gates and binding pockets for dioxygen, as an indistinguishable mimic of singlet dioxygen, were found on both above and below the isoalloxazine plane of the flavin mononucleotide cofactor in miniSOG. Of such gates and binding pockets, those lying within two opposite cones, coaxial with a line normal to the isoalloxazine plane, and with the vertex at the center of such a plane are those most visited by the escaping gas molecule. Out of residues most capable of quenching ¹O₂, Y30, lying near the base of one such a cone, and H85, near the base of the opposite cone, are held to be most responsible for the reduced quantum yield of ¹O₂ with folded miniSOG with respect to free flavin mononucleotide in solution.     

Knowledge-based Homology Modeling and Experimental Determination of Amino Acid Side Chain Accessibility by the Laser Photo CIDNP (Chemically Induced Dynamic Nuclear Polarization) Approach in Solution: Lessons from the Small Sialidase of Clostridium perfringens

The success of knowledge-based homology modelling is critically dependent on the predictive potency of the program structure-based calculations, which attempt to translate homologous sequences into three-dimensional structures, and on the actual relevance of the crystal structure for the protein topology. As quality control, experimental data for selected parameters of the proteins conformation are required. Using the crystal structure of the sialidase of Salmonella typhimurium as framework for model building of the homologous enzyme from Clostridium perfringens, a set of energy-minimised conformers is derived. These proteins present e.g. Tyr, Trp and His residues with an assessable area on the surface, since the side chains of these amino acid residues are responsive to chemically induced dynamic nuclear polarization (CIDNP), monitored by NMR. Hence, as first lesson, a comparative analysis for model-derived and experimentally determined values can be performed. The second lesson of this study concerns the notable impact of single amino acid substitutions (Tyr/Phe, Cys/Ser) on the surface accessibility of the CIDNP-reactive amino acid side chains in mutant forms of the sialidase. Corroborating the predictions from the theoretical calculations, the spectra of the engineered mutants reveal marked and non-uniform alterations. Thus, the effect of apparently rather conservative amino acid substitutions on a distinct conformational aspect of this protein, even at distant sites, should not be underestimated.     

Consequences of Sequential Ca2+ Occupancy for the Structure and Dynamics of CalbindinD9K: Computational Simulations and Comparison to Experimental Determinations in Solution

Abstract

Molecular dynamics (MD) simulations of the structures of calbindin_D9K (CAB) with different occupancies of the two Ca²⁺ binding sites were carried out to gain insight into structural and energetic consequences of sequential Ca²⁺ binding. The aim of the study is to identify effects of Ca-binding site occupancy that relate to the properties and functions of Ca-binding proteins containing EF-hand motifs. Two different models of solvation were employed, one defined by a linear, distance dependent dielectric permittivity (ε = r) and inclusion only of the 36 crystallographically observed water molecules, and the other with the protein immersed in a 9Å shell of explicit waters and ε = 1. Experimental results from x-ray crystallography, and insights from NMR and from measurements of hydrogen exchange rates in these systems served as tests and guides for assessing the quality, validity and mechanistic interpretation of the results from the computational study. The results of the MD simulations agree very well with earlier experimental observations that the structure of calbindin_D9k is rather insensitive to removal of Ca²⁺, and indicate that this insensitivity is not dependent on the order in which the ions are removed. The calculated values of the electrostatic potentials at the Ca²⁺ binding sites are very similar, in agreement with the small differences in the measured microscopic binding constants. Details of the dynamic mechanisms of molecular flexibility revealed by the MD simulations are also in good agreement with experimental findings, including the local properties identified from comparisons of hydrogen exchange rates in various parts of the structures of sequentially occupied forms of CAB. Estimation of the changes in configurational entropy from the rms fluctuations in the structures of CAB at various levels of Ca²⁺ occupancy in the EF-hands, supports earlier suggestions relating the dynamic properties of the protein to the observed cooperativity in the binding of Ca²⁺. The computational approaches and the results of the MD simulations are evaluated in relation to the study of effects of Ca²⁺ occupancy in calmodulin and troponin C where ion binding determines function and is known to trigger significant changes in structural and dynamic properties.     

Disruption Mechanism in the Helix of SPF Peptide by Interchanging E5 and K10 Residues: Inference from Molecular Dynamics Study

Abstract

Seminalplasmin (SPLN) is a 47-residue peptide (SDEKASPDKHHRFSLSRYAKLANR LANPKLLETFLSKWIGDRGNRSV) from bovine seminal plasma. It has broad spectrum antimicrobial activity, without any hemolytic activity. The 28–40 segment of SPLN with the sequence PKLLETFLSKWIG, designated as SPF, is the most hydrophobic stretch of SPLN and primarily responsible for the membrane-perturbing activity of SPLN. It was reported that SPF has a helical structure and the interchange of E5 and K10 residues disrupted the helical structure. The present paper reports a possible mechanism of disruption of the helical structure of SPF peptide during the interchange of E5 and K10 residues. The result is based on simulated annealing and molecular dynamics simulation studies on SPF and its four analogues with K10E, K10D, E5K, and E5K & K10E substitutions. It showed that K10 residue has a critical role in maintaining the highest helical content and the positions of charged residues are also very important for maintaining the helical structure of the SPF peptide. Formation of some new long-range hydrogen bonds and the rupture of some shortrange hydrogen bonds involving the tenth residue led to the disruption of helical structure of SPF peptide when E5 and K10 residues are interchanged.     

On CO Egress from,and Re‐Uptake by,the Enzyme MauG,as a Mimic of the Acquisition of Oxidizing Agents by the pre‐MADH_MauG System. A Molecular Mechanics Approach

In this work, by applying a non‐deterministic, randomly‐oriented minimal force to the dissociated CO ligand of the MauG‐CO system, the molecular‐dynamics (MD) behavior of this system could be quickly unraveled. It turned out that CO has no marked directional egress from the high‐spin c‐heme iron distal pocket. Rather, CO is able to exploit all interstices created during the protein fluctuations. Nonetheless, no steady route toward the surrounding solvent was ever observed: CO jumped first into other binding pockets before being able to escape the protein. In a few cases, on hitting the surrounding H₂O molecules, CO was observed to reverse direction, re‐entering the protein. A contention that conformational inversion of the P107 ring provides a gate to the iron ion is not supported by the present simulations.     

Binding Pockets and Permeation Channels for Dioxygen through Cofactorless 3‐Hydroxy‐2‐methylquinolin‐4‐one 2,4‐Dioxygenase in Association with Its Natural Substrate, 3‐Hydroxy‐2‐methylquinolin‐4(1H)‐one. A Perspective from Molecular Dynamics Simulations

This work describes an investigation of pathways and binging pockets (BPs) for dioxygen (O₂) through the cofactorless oxygenase 3‐hydroxy‐2‐methylquinolin‐4‐one 2,4‐dioxygenase in complex with its natural substrate, 3‐hydroxy‐2‐methylquinolin‐4(1H)‐one, in aqueous solution. The investigation tool was random‐acceleration molecular dynamics (RAMD), whereby a tiny, randomly oriented external force is applied to O₂ in order to accelerate its movements. In doing that, care was taken that the external force only continues, if O₂ moves along a direction for a given period of time, otherwise the force changed direction randomly. Gates for expulsion of O₂ from the protein, which can also be taken as gates for O₂ uptake, were found throughout almost the whole external surface of the protein, alongside a variety of BPs for O₂. The most exploited gates and BPs were not found to correspond to the single gate and BP proposed previously from the examination of the static model from X‐ray diffraction analysis of this system. Therefore, experimental investigations of this system that go beyond the static model are urgently needed.     

Molecular Dynamics Simulation in Solvent of the Estrogen Receptor Protein DNA Binding Domain in Complex with a Non-Consensus Estrogen Response Element DNA Sequence

Abstract

We investigated protein/DNA interactions, using molecular dynamics simulations computed between a 10 Angstom water layer model of the estrogen receptor (ER) protein DNA binding domain (DBD) amino acids and DNA of a non-consensus estrogen response element (ERE) consisting of 29 nucleotide base pairs. This ERE nucleotide sequence occurs naturally upstream of the Xenopus laevis Vitelligenin AI gene. The ER DBD is encoded by three exons. Namely, exons 2 and 3 which encode the two zinc binding motifs and a sequence of exon 4 which encodes a predicted alpha helix. We generated a computer model of the ER DBD using atomic coordinates derived from the average of 30 nuclear magnetic resonance (NMR) spectroscopy coordinate sets. Amino acids on the carboxyl end of the ER DBD were disordered in both X-ray crystallography and NMR determinations and no coordinates were reported. This disordered region includes 10 amino acids of a predicted alpha helix encoded in exon 4 at the exon 3/4 splice junction. These amino acids are known to be important in DNA binding and are also believed to function as a nuclear translocation signal sequence for the ER protein. We generated a computer model of the predicted alpha helix consisting of the 10 amino acids encoded in exon 4 and attached this helix to the carboxyl end of the ER DBD at the exon 3/4 splice junction site. We docked the ER DBD model within the DNA major groove halfsites of the 29 base pair non-consensus ERE and flanking nucleotides. We constructed a solvated model with the ER DBD/ERE complex surrounded by a ten Angstrom water layer and conducted molecular dynamics simulations. Hydrogen bonding interactions were monitored. In addition, van der Waals and electrostatic interaction energies were calculated. Amino acids of the ER DBD DNA recognition helix formed both direct and water mediated hydrogen bonds at cognate codon-anticodon nucleotide base and backbone sites within the ERE DNA right major groove halfsite. Amino acids of the ER DBD exon 4 encoded predicted alpha helix formed direct and water mediated H-bonds with base and backbone sites of their cognate codon-anticodon nucleotides within the minor grooves flanking the ERE DNA major groove halfsites. These interactions together induced bending of the DNA into the protein.