Molecular dynamics simulation of superoxide interacting with superoxide dismutase

MD simulation reveals differential binding of Cm(III) and Th(IV) with serum transferrin at acidic pH

The iron carrier human serum transferrin (sTf) is known to transport other metals, including some actinides (An). Radiotoxic An are routinely involved in the nuclear fuel cycle and the possibility of their accidental exposure cannot be ruled out. Understanding An interaction with sTf assumes a greater significance for the development of safe and efficacious chelators for their removal from the blood stream. Here we report several 100 ns equilibrium MD simulations of Cm(III)- and Th(IV)-loaded sTf at various protonation states of the protein to explore the possibility of the two An ions release and speciation. The results demonstrate variation in protonation state of dilysine pair (K206 and K296) and the tyrosine (Y188) residue is necessary for the opening of Cm(III)-bound protein and the release of the ion. For the tetravalent thorium, protonation of dilysine pair suffices to cause conformational changes of protein. However, in none of the protonation states, Th(IV) releases from sTf because of its strong electrostatic interaction with D63 in the first shell of the sTf binding cleft. Analysis of hydrogen bond, water bridge, and the evaluation of potential of mean forces of the An ions' release from sTf, substantiate the differential behavior of Cm(III) and Th(IV) at endosomal pH. The results provide insight in the regulation of Cm(III) and Th(IV) bioavailability that may prove useful for effective design of their decorporating agents and as well may help the future design of radiotherapy based on tetravalent ions.

Accelerating the weighted histogram analysis method by direct inversion in the iterative subspace

The weighted histogram analysis method (WHAM) for free energy calculations is a valuable tool to produce free energy differences with the minimal errors. Given multiple simulations, WHAM obtains from the distribution overlaps the optimal statistical estimator of the density of states, from which the free energy differences can be computed. The WHAM equations are often solved by an iterative procedure. In this work, we use a well-known linear algebra algorithm which allows for more rapid convergence to the solution. We find that the computational complexity of the iterative solution to WHAM and the closely-related multiple Bennett acceptance ratio (MBAR) method can be improved by using the method of direct inversion in the iterative subspace. We give examples from a lattice model, a simple liquid and an aqueous protein solution.

Molecular dynamics simulations of the enzyme Cu, Zn superoxide dismutase

The enzyme Cu, Zn superoxide dismutase (Cu,Zn-SOD) is a ubiquitous oxireductase, which is responsible for the cellular defense against oxidative stress caused by the high toxicity of the superoxide radical, and has been also linked to some cases of familiar amyotrophic lateral sclerosis. In the present study a set of molecular mechanics parameters for the active site of Cu,Zn-SOD has been derived. Afterward, an extensive molecular dynamics simulation has been carried out in an aqueous environment. The obtained results shed a further light on the structural flexibility of the backbone, where the active site is nested, and the solvation shell occupancy. The relatively small backbone deviation, shown by a root-mean-square deviation below 1.0 A, confirms the accuracy of the parameters. The solvent shell analysis has shown that the first solvation shell is located at about 5 A from the copper ion, generating an empty cavity with enough space to accommodate the superoxide radical. The low residence time means that a high permutation rate of water molecules in both solvation shells is consistent with the efficiency of this catalytic mechanism. Hybrid studies using ONIOM methodologies can now be done to evaluate the mechanistic implications of the explicit inclusion of the whole system.

The dynamics of ligand barrier crossing inside the acetylcholinesterase gorge

The dynamics of ligand movement through the constricted region of the acetylcholinesterase gorge is important in understanding how the ligand gains access to and is released from the active site of the enzyme. Molecular dynamics simulations of the simple ligand, tetramethylammonium, crossing this bottleneck region are conducted using umbrella potential sampling and activated flux techniques. The low potential of mean force obtained is consistent with the fast reaction rate of acetylcholinesterase observed experimentally. From the results of the activated dynamics simulations, local conformational fluctuations of the gorge residues and larger scale collective motions of the protein are found to correlate highly with the ligand crossing.

Molecular dynamics simulations of Cu,Zn superoxide dismutase: effect of temperature on dimer asymmetry

Molecular dynamics simulations of solvated dimeric Cu,Zn superoxide dismutase have been carried out at four temperatures, namely 200, 225, 250 and 300 K. Analysis of the backbone-to-backbone hydrogen bonds number indicates that the symmetry observed in the two subunits at 200 K is gradually lost by heating the system. The C(alpha) atoms displacement cross-correlation maps confirm that the asymmetric behaviour of the two subunits increases as a function of temperature. The dynamic cross-correlation of the subunits volumes indicates a fast correlation between the two subunits at 300 K, which is delayed upon lowering the simulation temperature. These results indicate that temperature plays an essential role in injecting such an asymmetry; the two subunits being asymmetric and in rapid communication at 300 K, and almost symmetric and in slow communication at lower temperatures.

CuZn Superoxide Dismutase Geometry Optimization, Energetics, and Redox Potential Calculations by Density Functional and Electrostatic Methods

The structures, energetics, and orbital- and charge-dependent properties of copper zinc superoxide dismutase (CuZnSOD) have been studied using density functional and electrostatic methods. The CuZnSOD was represented with a model consisting of copper and zinc sites connected by a bridging histidine ligand. In addition to the bridge, three histidine ligands and one water molecule were bonded to the Cu ion in the copper site as first-shell ligands. Two histidine ligands and an aspartate were coordinated to the zinc ion in the zinc site. Full optimization of the model was performed using different functionals, both local and nonlocal. Geometrical parameters calculated with the nonlocal functionals agree well with the experimental X-ray data. In our calculated results, the His61 Nepsilon-Cu bond in the active site breaks during the reduction and protonation, consistent with a number of X-ray structures and with EXAFS and NMR evidence. The reduction potential and pK(a) of the coupled electron/proton reaction catalyzed by CuZnSOD were determined using different models for the extended environment-from an electrostatic representation of continuum solvent, to the full protein/solvent environment using a Poisson-Boltzmann method. The predicted redox potential and pK(a) values determined using the model with the full protein/solvent environment are in excellent agreement with experiment. Inclusion of the full protein environment is essential for an accurate description of the redox process. Although the zinc ion does not play a direct redox role in the dismutation, its electronic contribution is very important for the catalytic mechanism.

A spectroscopic and molecular dynamics study of native and of a mutant of Xenopus laevis Cu,Zn superoxide dismutase: mechanistic consequences of replacing four charged amino acids on the 'electrostatic' loop

Neutralisation by site-directed mutagenesis of four charged and highly conserved residues of the electrostatic loop of Cu,Zn superoxide dismutase from Xenopus laevis, involved in the electrostatic attraction of the substrate: Lys120-->Leu, Asp130-->Gln, Glu131-->Gln and Lys134-->Thr, gives rise to a mutant enzyme which displays an affinity for monovalent inhibitor anions, such as N3-, higher than that of the wild type. Analysis of 300 ps of molecular dynamics simulation carried out on the wild type and on the Xenopus laevis Cu,Zn superoxide dismutase mutant indicates that the two proteins display a distinct dynamical behaviour. In particular the root mean square deviation from the starting structure, the number of residues in random coil conformations, the number of residues in unfavourable regions of the Ramachandran plot indicate that the mutant displays a rigidity higher than the native enzyme. This is also evidenced by the loss of dynamical cross correlations in the simulation of the mutant, which on the other hand are present in the wild type. Moreover the mutant protein shows a different organisation of the backbone-to-backbone hydrogen bonds network that generates a rigid structure leading to an increase of the active site accessibility when compared to the native enzyme. It is suggested that the rigid state in which the mutant is confined, accompanied by the increase of the solvent accessible surface of the active site may explain the difference in reactivity toward the inhibitor anion.

Dimer asymmetry in superoxide dismutase studied by molecular dynamics simulation

Molecular dynamics (MD) simulations of 100 ps have been carried out to study the active-site behaviour of the Cu,Zn superoxide dismutase dimer (SOD) in water. The active site of each subunit was monitored during the whole simulation by calculating the distances between functional residues and the catalytic copper. The results indicate that charge orientation is maintained at each active site but the solvent accessibility varies. Analysis of the MD simulation, carried out by using the atomic displacement covariance matrix, has shown a different intra-subunit correlation pattern for the two monomers and the presence of inter-subunit correlations. The MD simulation presented here indicates an asymmetry in the two active sites and different dynamic behaviour of the two SOD subunits.

Molecular dynamics studies on mutants of Cu,Zn superoxide dismutase: the functional role of charged residues in the electrostatic loop VII

Molecular dynamics (MD) calculations have been performed on mutants of superoxide dismutase (SOD) on some residues present in the electrostatic loop. These calculations have provided the solution structures for the mutants Thr-137-->Ile and Arg; Lys-136-->Ala; Glu-132-->Gln; Glu-133-->Gln; Glu-132, Glu-133-->Gln-132, Gln-133 and-->Gln-132, Lys-133. The structural and dynamic properties of these mutants have been correlated with the catalytic properties and available spectroscopic data. The water molecule present in the active site close to the copper ion in wild type (WT) SOD is missing in the MD average structure of the Thr-137-->Ile mutant, while this molecule is present in the MD average structures of all the other mutants and of WT SOD. This agrees with the experimental data. This is an important result that shows the validity of our calculations and their ability to reproduce even subtle structural features. Addition of one or more positive charges on the 132 and/or 133 positions does not sizably perturb the structure of the active site channel, while the introduction of a positively charged residue (Arg) on position 137 has a large effect on the structure of the electrostatic loop. Analysis of the MD average structures of these mutants has pointed out that the simple electrostatic effects of charged residues in the channel are not the only factor relevant for enzymatic behavior but that the structure of the electrostatic loop and the location of the charged residues also contribute to the catalytic properties of SOD.

Computational chemistry and molecular modeling of electron-transfer proteins

The methods of computational chemistry and molecular modeling are becoming more and more accessible to biochemists with the advent of fast, inexpensive graphics workstations and well-tested computer programs. The state of the art in small molecules allows chemists to use these programs as "black boxes" and be confident of the results at an amazingly high level of precision. This is not the case, however, for biological macro-molecules at this time. Therefore, it is necessary before using the programs listed in Section I that we familiarize ourselves with their theoretical basis and limitations. It is also important that they be used in the context of the accumulated literature on their use. The survey given in this chapter is intended as an introduction to these tools and as a source for initiating the discovery process with the literature cited.

Molecular modeling methods. Basic techniques and challenging problems

An overview is presented of computer modeling and simulation methods that play an increasing role in drug design: quantum chemical methods, molecular mechanics, molecular dynamics and Brownian dynamics. The application of molecular dynamics for the prediction of thermodynamic properties like free energy differences and binding constants is discussed. The Brownian dynamics method is presented in connection with the calculation of effective electrostatic forces using the Poisson-Boltzmann equation, which allows one to sample ligand-binding geometries and to predict the kinetics of diffusion-limited enzyme reactions. New techniques that have recently been extensively developed, such as the global energy minimization and quantum-classical dynamics methods, are also introduced. The molecular modeling methods are illustrated with selected examples.