Cooperative ligand reorientations in cytochrome c3: a molecular dynamics simulation.

Molecular dynamics simulations of a tetraheme cytochrome c3 were performed to investigate dynamic aspects of the motion of the axial heme iron ligands. It was found that persistent transitions between alternate axial imidazole orientations of the histidine incorporated in the CXXCH heme binding sequence occurred via correlated motions. The correlated motions involved virtually all of the atoms comprising the polypeptide backbone of the heme binding sequence as well as the histidine imidazole side-chain.     

Molecular dynamics simulations and free energy calculations of base flipping in dsRNA

The family of adenosine deaminases acting on RNA (ADARs) targets adenosines in RNA that is mainly double stranded. Some substrates are promiscuously deaminated whereas others, such as the mammalian glutamate receptor B (gluR-B) pre-mRNA, are more selectively deaminated. Many DNA/RNA-base modification enzymes use a base flipping mechanism to be able to reach their target base and it is believed that ADARs function in a similar way. In this study we used molecular dynamics (MD) simulations to describe two sites on the gluR-B pre-mRNA, the selectively targeted R/G site and the nontargeted 46 site, in an attempt to explain the substrate specificity. We used regular MD and also a forced base flipping method with umbrella sampling to calculate the free energy of base opening. Spontaneous opening of the mismatched adenosine was observed for the R/G site but not for the 46 site.     

Molecular dynamics simulations of the cytochrome c3-rubredoxin complex from Desulfovibrio vulgaris.

Molecular dynamics simulations have been carried out on the complex formed between the tetraheme cytochrome c3 and the iron protein rubredoxin from the sulfate-reducing bacterium Desulfovibrio vulgaris. These simulations were performed both with explicit solvent water molecules included, and without solvent molecules using a distance-dependent dielectric constant to approximate the screening effects of solvent. The results of both simulations are strikingly different, indicating that the representation of environmental effects is important in such simulations. For example, a striking adaptation of the two proteins seen in the nonsolvated simulation is not seen when explicit solvent water is included; in fact, the complex appears to become weaker in the solvated simulation. Nonetheless, the iron-iron distance decreases more significantly in the solvated simulation than in the nonsolvated simulation. It was found that in both cases molecular dynamics optimized the structures further than energy minimization alone.     

Molecular dynamics simulation of protein folding by essential dynamics sampling: folding landscape of horse heart cytochrome c

A new method for simulating the folding process of a protein is reported. The method is based on the essential dynamics sampling technique. In essential dynamics sampling, a usual molecular dynamics simulation is performed, but only those steps, not increasing the distance from a target structure, are accepted. The distance is calculated in a configurational subspace defined by a set of generalized coordinates obtained by an essential dynamics analysis of an equilibrated trajectory. The method was applied to the folding process of horse heart cytochrome c, a protein with approximately 3000 degrees of freedom. Starting from structures, with a root-mean-square deviation of approximately 20 A from the crystal structure, the correct folding was obtained, by utilizing only 106 generalized degrees of freedom, chosen among those accounting for the backbone carbon atoms motions, hence not containing any information on the side chains. The folding pathways found are in agreement with experimental data on the same molecule.     

Absolute binding free energy calculations using molecular dynamics simulations with restraining potentials

The absolute (standard) binding free energy of eight FK506-related ligands to FKBP12 is calculated using free energy perturbation molecular dynamics (FEP/MD) simulations with explicit solvent. A number of features are implemented to improve the accuracy and enhance the convergence of the calculations. First, the absolute binding free energy is decomposed into sequential steps during which the ligand-surrounding interactions as well as various biasing potentials restraining the translation, orientation, and conformation of the ligand are turned "on" and "off." Second, sampling of the ligand conformation is enforced by a restraining potential based on the root mean-square deviation relative to the bound state conformation. The effect of all the restraining potentials is rigorously unbiased, and it is shown explicitly that the final results are independent of all artificial restraints. Third, the repulsive and dispersive free energy contribution arising from the Lennard-Jones interactions of the ligand with its surrounding (protein and solvent) is calculated using the Weeks-Chandler-Andersen separation. This separation also improves convergence of the FEP/MD calculations. Fourth, to decrease the computational cost, only a small number of atoms in the vicinity of the binding site are simulated explicitly, while all the influence of the remaining atoms is incorporated implicitly using the generalized solvent boundary potential (GSBP) method. With GSBP, the size of the simulated FKBP12/ligand systems is significantly reduced, from approximately 25,000 to 2500. The computations are very efficient and the statistical error is small ( approximately 1 kcal/mol). The calculated binding free energies are generally in good agreement with available experimental data and previous calculations (within approximately 2 kcal/mol). The present results indicate that a strategy based on FEP/MD simulations of a reduced GSBP atomic model sampled with conformational, translational, and orientational restraining potentials can be computationally inexpensive and accurate.     

Molecular aspects of cytochrome c oxidase: Structure and dynamics

Summary In the last few years much attention has been dedicated to the elucidation of some of the molecular aspects of cytochrome c oxidase. It has been shown conclusively that the enzyme from several sources (yeast, Neurospora, heart, liver) contains seven different subunits, which are asymmetrically inserted in the membrane. All of these are in contact with the lipid bilayer (except subunits V and VI) and to a greater or lesser extent with the water phase as well (except for subunit I). Subunit II of the enzyme appears to be involved in the formation of the binding site of cytochrome c. The location of the redox groups of the enzyme is still a matter of controversy. Their distance from the cytochrome c heme group is approximately 35 Å such that electron tunneling appears to be the only possible mechanism for transporting electrons across such a distance.A proton pump appears to be associated with electron transport and approximately one proton is extruded per electron equivalent reducing oxygen via the enzyme. N,N, dicyclohexylcarbodiimide a well-established inhibitor of H⁺-translocating ATPases inhibits the proton pump and labels specifically subunit III of the enzyme.Abbreviations CCCP carbonyl cyanide mchlorophenylhydrazone - DADH₂ diamino-durene (reduced form) - DCCD N,N-dicyclohexylcarbodiimide - FCCP carbonyl cyanide p-trifluoromethoxyphenylhydrazone - Hepes 4-(2-hydroxyethyl)-1-piperazine-ethanesulphonic acid - NEM N-ethyl-maleimide - pH transmembrane chemical-potential gradient of H⁺ - transmembrane electricalpotential gradient - transmembrane electrochemical-potential gradient of H⁺ - Q/QH₂ oxidised/reduced form of ubiquinone - TMPD/TMPDH₂ non-protonated/protonated form of tetramethyl-p-phenylenediamine - SDS sodium dodecyl sulphate     

Insights into scFv:drug binding using the molecular dynamics simulation and free energy calculation

Molecular dynamics simulations and free energy calculation have been performed to study how the single-chain variable fragment (scFv) binds methamphetamine (METH) and amphetamine (AMP). The structures of the scFv:METH and the scFv:AMP complexes are analyzed by examining the time-dependence of their RMSDs, by analyzing the distance between some key atoms of the selected residues, and by comparing the averaged structures with their corresponding crystallographic structures. It is observed that binding an AMP to the scFv does not cause significant changes to the binding pocket of the scFv:ligand complex. The binding free energy of scFv:AMP without introducing an extra water into the binding pocket is much stronger than scFv:METH. This is against the first of the two scenarios postulated in the experimental work of Celikel et al. (Protein Science 18, 2336 (2009)). However, adding a water to the AMP (at the position of the methyl group of METH), the binding free energy of the scFv:AMP-H2O complex, is found to be significantly weaker than scFv:METH. This is consistent with the second of the two scenarios given by Celikel et al. Decomposition of the binding energy into ligand-residue pair interactions shows that two residues (Tyr175 and Tyr177) have nearly-zero interactions with AMP in the scFv:AMP-H2O complex, whereas their interactions with METH in the scFv:METH complex are as large as -0.8 and -0.74 kcal mol^-1. The insights gained from this study may be helpful in designing more potent antibodies in treating METH abuse.     

3D-QSAR-aided design of potent c-Met inhibitors using molecular dynamics simulation and binding free energy calculation

Abstract

Mesenchymal-epithelial transition factor (c-Met) is a member of receptor tyrosine kinase. It involves in various cellular signaling pathways which includes proliferation, motility, migration, and invasion. Over-expression of c-Met has been reported in various cancers. Hence, it is an ideal therapeutic target for cancer. The main objective of the study is to identify crucial residues involved in the inhibition of c-Met kinase and to design a series of potent imidazo [4,5-b] pyrazine derivatives as c-Met inhibitors. Docking was used to identify important active site residues involved in the inhibition of c-Met kinase which was further validated by 100 ns of molecular dynamics simulation and free energy calculation using molecular mechanics generalized born surface area. Furthermore, binding energy decomposition identified that residues Tyr1230, Met1211, Asp1222, Tyr1159, Met1160, Val1092, Ala1108, and Leu1157 contributed favorably to the binding stability of compound 32. Receptor-guided Comparative Molecular Field Analysis (CoMFA) (q² = 0.751, NOC = 6, r² = 0.933) and Comparative Molecular Similarity Indices Analysis (COMSIA) (q² = 0.744, NOC = 6, r² = 0.950) models were generated based on the docked conformation of the most active compound 32. The robustness of these models was tested using various validation techniques and found to be predictive. The results of CoMFA and CoMSIA contour maps exposed the regions favorable to enhance the activity. Based on this information, 27 novel c-Met inhibitors were designed. These designed compounds exhibited potent activity than the most active compound of the existing dataset.

Communicated by Ramaswamy H. Sarma     

Measurement of free radical oxygen generation by cytochrome c reduction requirement for cytochrome c oxidase blockade

Data is presented showing that one commercial preparation of cytochrome c, used to trap and measure free radical superoxide anion, can be contaminated with cytochrome c oxidase activity. This activity can vary from lot to lot, can introduce variability into the measurement of superoxide anion and can result in falsely low estimations of free radical formation. This cytochrome c oxidase activity can be inhibited by low (0.2 mM) concentrations of KCN. Blockade of the cytochrome c oxidase activity allows reproducible measurement of superoxide anion formation at low levels by red cells.     

Multiwavelength analysis of the kinetics of reduction of cytochrome aa3 by cytochrome c.

Some new approaches to the kinetic study of the reduction of cytochrome aa3 by cytochrome c are presented. The primary innovations are the use of a spectrometer which can acquire multiwavelength data as fast as every 10 microseconds, and the application of a variety of analytical methods which can utilize simultaneously all of the time-resolved spectral data. These techniques include singular value decomposition (SVD), deconvolutions based on pure Gaussian models for absorption peaks, deconvolutions based on isolated absorption spectra for the pure components, and simulations of SVD-deduced and actual experimental difference spectra. The reduction characteristics of the anaerobic resting enzyme can be distinguished from those of pulsed forms. In the former case, only two electrons can be bound by cytochrome aa3, whereas in the latter case complete reduction of the enzyme is achieved.     

Molecular dynamics simulation studies on structural and conformational changes in tyrosine-67 nitrated cytochrome c

Protein tyrosine nitration is well-established post-translational modification occurring in a number of diseases, viz. neurodegenerative, cardiovascular diseases, ageing, etc. Tyrosine-67 (Tyr-67) nitration of cytochrome c (cyt c) was observed under oxidative stress affecting its structure and electron transfer properties. Hence, in this study, molecular dynamics (MD) simulations were carried out at room temperature to investigate the structural and conformational changes in the nitrated cyt c's. MD results revealed that the bond between FE (Heme-105) and S (Met-80) considerably weakened, radius of gyration, backbone and Cα root-mean-square deviations decreased and hydrogen bonding increased in the nitrated cyt c's relative to wild type (WT) cyt c. Ramachandran plot analysis revealed that N- and C-terminal helices also affected by nitration at CE2 carbon atom. Furthermore, essential dynamics analysis showed that amplitude of concerted motion decreased in the nitrated cyt c's, perhaps due to the increase in the hydrogen bonding interaction. Taken together, the structural and conformational changes in the active site Tyr-67 nitrated cyt c may have implications in the loss of electron/proton transfer and gain of apoptotic properties.     

Molecular dynamics simulation of interleukin-2 and its complex and determination of the binding free energy

Interleukin-2 (IL-2) protein belongs to the signal modulator cytokine's family and therefore it is prevalent for immunological responses. It has been identified as a centrally important potential drug target for the inhibition of protein-protein interactions; so as to suppress the immunological responses associated with autoimmune, inflammatory and immunological diseases, and cancer. In the present work, we have performed two independent 100?ns of molecular dynamics (MD) simulations on the apo IL-2 protein and its ligand-bound complex (with a potent inhibitor FRG), to study the effect of inhibitor binding on the dynamics and stability of the protein. The calculation of binding free energy via post-processing end state method of Molecular Mechanics Poisson Boltzmann Surface Area (MM-PBSA) and Molecular Mechanics Generalised Born Surface Area (MM-GBSA) has inferred a good correlation in accordance with the already reported experimental data, demonstrating that the free energy of binding calculated by the two methods has no significant difference. The investigation of individual components of free energy revealed that the association of IL-2 protein with FRG ligand is primarily driven by the van der Waals energy contribution that represents the non-polar/hydrophobic energy contribution as dominant in this case of ligand binding.     

Molecular dynamics and free energy calculations of the B and Z forms of C8-arylguanine modified oligonucleotides

Arylhydrazines found in the mushroom Agaricus bisporus have been shown to be carcinogenic. Upon metabolic activation, arylhydrazines are transformed into aryl radicals, forming 8-arylpurines, which may play a role in arylhydrazine carcinogenesis. These adducts are poorly read and inhibit chain extension but do alter the conformational preferences of oligonucleotides. We have shown that C8-phenylguanine modification of d(CGCGCG*CGCG) (G*= 8-phenylguanine) stabilizes it in the Z-DNA conformation (B/Z-DNA=1:1, 200 mM NaCl, pH 7.4). Here we have conducted molecular dynamics and free energy calculations to determine the sources(s) of these conformational affects and to predict the affect of the related C8-tolyl and C8-hydroxymethylphenyl guanine adducts on B/Z-DNA equilibrium. Force field parameters for the modified guanines were first developed using Guassian98 employing the B3LYP method and the standard 6-31G* basis set and fit to the Cornell 94 force field with RESP. Molecular dynamics simulations and free energy calculations, using the suite of programs contained in Amber 6 and 7 with the Cornell 94 force field, were used to determine the structural and thermodynamic properties of the DNA. The principal factors that drive conformation are stacking of the aryl group over the 5'-cytosine in the phenyl and tolyl modified oligonucleotides while hydrogen bonding opposes stacking in the hydroxymethylphenyl derivative. The phenyl and tolyl-modified DNA's favored the Z-DNA form as did the hydroxymethylphenyl derivative when hydrogen bonding was not present. The B-DNA conformation was preferred by the unmodified oligonucleotide and by the hydroxymethylphenyl-modified oligonucleotide when hydrogen bonding was considered. Z-DNA stability was not found to directly correlated with carcinogenicity and additional biological factors, such as recognition and repair, may also need to be considered in addition to Z-DNA formation.     

Molecular dynamics study of Desulfovibrio africanus cytochrome c3 in oxidized and reduced forms

A 5-ns molecular dynamics study of a tetraheme cytochrome in fully oxidized and reduced forms was performed using the CHARMM molecular modeling program, with explicit water molecules, Langevin dynamics thermalization, Particle Mesh Ewald long-range electrostatics, and quantum mechanical determination of heme partial charges. The simulations used, as starting points, crystallographic structures of the oxidized and reduced forms of the acidic cytochrome c(3) from Desulfovibrio africanus obtained at pH 5.6. In this paper we also report structures for the two forms obtained at pH 8. In contrast to previous cytochrome c(3) dynamics simulations, our model is stable. The simulation structures agree reasonably well with the crystallographic ones, but our models show higher flexibility and the water molecules are more labile. We have compared in detail the differences between the simulated and experimental structures of the two redox states and observe that the hydration structure is highly dependent on the redox state. We have also analyzed the interaction energy terms between the hemes, the protein residues, and water. The direct electrostatic interaction between hemes is weak and nearly insensitive to the redox state, but the remaining terms are large and contribute in a complex way to the overall potential energy differences that we see between the redox states.     

Computational investigation of the HIV-1 Rev multimerization using molecular dynamics simulations and binding free energy calculations

The HIV Rev protein mediates the nuclear export of viral mRNA, and is thereby essential for the production of late viral proteins in the replication cycle. Rev forms a large organized multimeric protein-protein complex for proper functioning. Recently, the three-dimensional structures of a Rev dimer and tetramer have been resolved and provide the basis for a thorough structural analysis of the binding interaction. Here, molecular dynamics (MD) and binding free energy calculations were performed to elucidate the forces thriving dimerization and higher order multimerization of the Rev protein. It is found that despite the structural differences between each crystal structure, both display a similar behavior according to our calculations. Our analysis based on a molecular mechanics-generalized Born surface area (MM/GBSA) and a configurational entropy approach demonstrates that the higher order multimerization site is much weaker than the dimerization site. In addition, a quantitative hot spot analysis combined with a mutational analysis reveals the most contributing amino acid residues for protein interactions in agreement with experimental results. Additional residues were found in each interface, which are important for the protein interaction. The investigation of the thermodynamics of the Rev multimerization interactions performed here could be a further step in the development of novel antiretrovirals using structure based drug design. Moreover, the variability of the angle between each Rev monomer as measured during the MD simulations suggests a role of the Rev protein in allowing flexibility of the arginine rich domain (ARM) to accommodate RNA binding.     

Kinetic study of the reduction mechanism for Desulfovibrio gigas cytochrome c3.

The kinetic aspects of the reduction process in cytochrome c3 from Desulfovibrio gigas have been investigated over a wide range of pH values ranging between pH 5.8 and pH 9.8. The data have been analyzed in the framework of an I2H4 interaction network coupled to a proton-linked equilibrium between two tertiary structures (Cornish-Bowden, A. & Koshland, D.E. Jr (1970) J. Biol. Chem. 245, 6241-6250). The kinetic rate constants for the reduction of the four hemes for the two tertiary conformations have been characterized in the framework of the thermodynamic network obtained from the equilibrium analysis (Coletta, M., Catarino, T., LeGall, J.J. & Xavier, A.V. (1991) Eur. J. Biochem. 202, 1101-1106). The intrinsic reduction rate constants determined by reaction with sodium dithionite for two hemes (namely heme 4 and heme 1) are significantly faster than those for the other two heme residues. In view of the equilibrium redox properties, heme 4 (with the fastest reduction rate) may then work as the kinetic electron-capturing site for the electrons from sodium dithionite. The transfer to hemes 2 and 3 then occurs by virtue of their free-energy levels at equilibrium. At our experimental conditions, there is also transfer of electrons to hemes 2 and 3 from heme 1, which is reduced at a slower rate than heme 4, thus contributing to the biphasic kinetics observed for the overall process. The kinetic parameters obtained are discussed in terms of the mechanism proposed for the coupling between the electron and proton transfer, as induced by the heme/heme cooperativity network.     

Molecular dynamics simulations of the resting and hydrogen peroxide-bound states of cytochrome c peroxidase.

The current hypothesis for the formation of the catalytically active compound I of peroxidases from the resting state and peroxide involves formation of a reversible "inner-sphere" complex in which the peroxide is bound to the heme iron. It is this precursor that is postulated to then form compound I. However, this crucial putative transient intermediate has not yet been definitively detected or characterized by experimental methods. We report here the use of energy minimization and molecular dynamics simulation together with the known X-ray structure of cytochrome c peroxidase to investigate the nature of this complex and comparisons of it with the resting state in which a water is bound as a ligand. Among the properties monitored in these simulations are the mode of binding of the peroxide to the heme iron, its interactions with neighboring amino acid residues, and the extent to which the binding of the peroxide perturbs both the local environment around the heme unit and more distant regions. The results of this study indicate that solvated, full protein dynamics is required to obtain reliable results for the known resting-state complex and hence for the uncharacterized peroxide complex. In this complex, the peroxide binds to the heme iron in a dynamically averaged end-on fashion, rather than a bridged structure, with approximately equal probability of each oxygen serving as the ligand to the iron. Binding of the peroxide as a ligand disrupts the H-bonded network of waters in the distal binding pocket which are present in the resting state, but there is no dramatic perturbation of the nearby amino acid residues.(ABSTRACT TRUNCATED AT 250 WORDS)     

Probing the interaction mechanism of small molecule inhibitors with matriptase based on molecular dynamics simulation and free energy calculations

Matriptase is a serine protease associated with a wide variety of human tumors and carcinoma progression. Up to now, many promising anti-cancer drugs have been developed. However, the detailed structure–function relationship between inhibitors and matriptase remains elusive. In this work, molecular dynamics simulation and binding free energy studies were performed to investigate the biochemistry behaviors of two class inhibitors binding to matriptase. The binding free energies predicted by MM/GBSA methods are in good agreement with the experimental bioactivities, and the analysis of the individual energy terms suggests that the van der Waals interaction is the major driving force for ligand binding. The key residues responsible for achieving strong binding have been identified by the MM/GBSA free energy decomposition analysis. Especially, Trp215 and Phe99 had an important impact on active site architecture and ligand binding. The results clearly identify the two class inhibitors exist different binding modes. Through summarizing the two different modes, we have mastered some important and favorable interaction patterns between matriptase and inhibitors. Our findings would be helpful for understanding the interaction mechanism between the inhibitor and matriptase and afford important guidance for the rational design of potent matriptase inhibitors.