Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “MD simulation”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 253 records · Page 14Linked to original sources
PubMed2003

A combined atomic force microscopy and molecular dynamics simulation study on a plastocyanin mutant chemisorbed on a gold surface.

A mutant of copper plastocyanin, covalently bound to an Au (111) surface through an engineered disulfide bridge, was investigated in aqueous medium by atomic force microscopy (AFM) and molecular dynamics (MD) simulations. Tapping-mode AFM images revealed adsorption of single molecules which are homogeneously distributed over the substrate and strongly bound to gold and display uniform lateral size. A statistical analysis of the height of the macromolecules on the gold substrate evidenced a distribution around a mean value consistent with that expected from the crystallographic data and with a relatively large standard deviation. A 10-ns classical MD simulation of mutated plastocyanin, hydrated by a layer of water, covalently bound to a gold surface by one or two sulfur atoms, was performed. The simulations indicate that the bound protein retains, in both cases, its overall tertiary structure during the dynamic evolution. Moreover, the macro-molecule can assume different orientations with respect to the gold substrate, which give rise to a distribution of heights on the gold substrate. Experimental and MD simulation results are compared and discussed in connection with the topological and dynamical properties of the protein system.

Adsorption↗
PubMed1999

'Wave-type' structure of a synthetic hexaglycosylated decapeptide: a part of the extracellular domain of human glycophorin A.

The three-dimensional structure of a glycopeptide, His-Thr*-Ser*-Thr*-Ser*-Ser*-Ser*-Val-Thr-Lys, with 2-acetamido-2-deoxy-alpha-D-galactose (GalNAc) residues linked to six adjacent amino acids from Thr-10 to Ser-15, was studied by NMR spectroscopy and molecular dynamics (MD) simulations. The hexaglycosylated decapeptide is part of the extracellular domain of human glycophorin A and shows an extended structure of the peptide backbone due to O-glycosylation. Furthermore, each GalNAc residue exhibits one and only one NOE contact from the NHAc proton to the backbone amide proton of the amino acid that the sugar is directly bound to. This indicates a strong preference for the orientation of all GalNAc residues towards the N-terminus. NOE build-up curves were used to determine 42 inter-proton distances that, in connection with phi angles of the peptide backbone obtained from 3J-coupling constants, resulted in constraints for a MD simulation in water. The NMR data and the MD simulations show a preference for an extended backbone structure. The GalNAc residues are located alternatingly on opposite sides of the backbone and reduce the flexibility of the peptide backbone. The conformation of the molecule is relatively rigid and shows a 'wave-type' 3D structure of the peptide backbone within the glycosylation cluster. This new structural element is also supported by the unusual CD spectrum of the glycopeptide.

Acetylgalactosamine↗
PubMed2003

Putative dynamics of vasopressin in its V1a receptor binding site.

The molecular architecture of the GPCRs, including the dynamic set of interactions between the receptor and the ligand, is one of the key structural questions of biophysical approaches. In the present study, molecular dynamics (MD) simulations were performed on the well-validated molecular model of the vasopressin V1a receptor applying different parameters (i.e., force fields, time variation, use of constraints) in order to sample the conformational space of the endogenous ligand arginine vasopressin (AVP), to explore different putative binding modes, and to analyze the simulation results with respect to experimental data. Noteworthy, it is to mention that for the first time a model of the vasopressin receptor remained stable in a 500 ps MD simulation run under vacuo boundary conditions using the Kollman all-atom FF even though no constraints were imposed. Conclusively, we determined an optimized experimental procedure for studying the dynamics and structure-functionship of this highly important family of GPCRs: the use of MD simulations with the Kollman all-atom force-field parameters on a constrained receptor. Our simplified model may be used as a basis for structure based design of new GPCR ligands and for in silico screening of virtual combinatorial chemistry libraries.

Arginine↗
PubMed2000

Glycosaminoglycan conformation: do aqueous molecular dynamics simulations agree with x-ray fiber diffraction?

Glycosaminoglycan-protein interactions are biologically important and require an appreciation of glycan molecular shape in solution, which is presently unavailable. In previous studies we found strong similarity between aqueous molecular dynamics (MD) simulations and published x-ray diffraction refinements of hyaluronan. We have applied a similar approach here to chondroitin and dermatan, attempting to clarify some of the issues raised by the x-ray diffraction literature relating to chondroitin and dermatan sulfate. We predict that chondroitin has the same beta(1-->4) linkage conformation as hyaluronan, and that their average beta(1-->3) conformations differ. This is explained by changes in hydrogen-bonding across this linkage, resulting from its axial hydroxyl, causing a different sampling of left-handed helices in chondroitin (2.5- to 3.5-fold) as compared with hyaluronan (3.0- to 4.0-fold). Few right-handed helices, which lack intramolecular hydrogen-bonds, were sampled during our MD simulations. Thus, we propose that the 8-fold helix observed in chondroitin-6-sulfate, represented in the literature as an 8(3) helix (right-handed), though it has never been refined, is more likely to be 8(5) (left-handed) helix. Molecular dynamics simulations implied that (4)C(1) and (2)S(O), but not (1)C(4), forms of iduronate could be used in refinements of dermatan x-ray fiber diffraction patterns. Current models of 8-fold dermatan sulfate chains containing (4)C(1) iduronate refine to right-handed helices, which possess no intramolecular hydrogen-bonds. However, MD simulations predict that models containing (2)S(O) iduronate could provide better (8(5) helix) starting structures for refinement. Thus, the 8-fold dermatan sulfate refinement (8(3) helix) could be in error.

Carbohydrate Conformation↗
PubMed2002

Rotational dynamics of propane in Na-Y zeolite: a molecular dynamics and quasielastic neutron-scattering study.

We report results from molecular dynamics (MD) simulations and quasielastic neutron-scattering (QENS) measurements on the rotational dynamics of propane in Na-Y zeolite at room temperature with a loading of four molecules per alpha cage. Rotational part of the intermediate scattering function F(Q,t) obtained from the MD simulation suggests that rotational motion is faster relative to the translational motion. Various rotational models fitted to the MD data suggest that rotation is isotropic. It is found that the hydrogen atoms lie, on the average, on a sphere of radius 1.88+/-0.05 A, which is also the average distance of the hydrogen atoms from the center of mass of the propane molecule. Results from QENS measurements are in excellent agreement with those obtained from MD, suggesting that the intermolecular potential employed in the MD simulation provides a realistic description of propane motion within faujasite. The rotational diffusion constant D(R) is 1.05+/-0.09 x 10(12) sec(-1) from the QENS data, which may be compared with that obtained from the MD data (0.82+/-0.05 x 10(12) sec(-1)).

Journal Article↗
PubMed2004

Molecular dynamics studies of ground state and intermediate of the hyperthermophilic indole-3-glycerol phosphate synthase.

Indole-3-glycerol phosphate synthase catalyzes the terminal ring closure step in tryptophan biosynthesis. In this paper, we compare the results from molecular dynamics (MD) simulations of enzyme-bound substrate at 298, 333, 363, and 385 K and the enzyme-bound intermediate at 385 K, solvated in TIP3P water box with a CHARMM force field. Results from MD simulations agree with experimental studies supporting the observation that Lys-110 is the general acid. Based on its location in the active site during the MD simulations, Glu-210 warrants classification as the general base instead of the previously proposed Glu-159. We find that the relative population of the reactive enzyme-substrate Michaelis conformers [near attack conformers (NACs)] with temperature correlates well (correlation coefficient of 0.96) with the relative activity of this thermophilic enzyme. At higher temperature, the enzyme-substrate electrostatic interaction favors the binding of the substrate in NAC conformation, whereas, at lower temperature, the substrate is distorted and bound in a nonreactive conformation. This change is reflected in the approximately 1,100-fold increase in population of NACs at 385 K relative to 298 K. The easily determined population of NACs at given temperature tells much about the thermophilic property of the enzyme. Thus, the hyperthermophilic enzyme has evolved to have optimum activity at high temperatures, and, with lowering of the temperature, the electrostatic interaction at the active site is enhanced and the structure is deformed. This model can be regarded as a general explanation for the activity of hyperthermophilic enzymes.

Indole-3-Glycerol-Phosphate Synthase↗
PubMed2000

Molecular dynamics simulations of bovine cathepsin B and its complex with CA074.

To promote our better understanding of the dynamic stability of the bovine cathepsin B structure, which is characterized by an extra disulfide bond at Cys148-Cys252 from the other species, and of the binding stability of CA074 (a cathepsin B-specific inhibitor), molecular dynamics (MD) simulations were performed for the enzyme and its CA074 complex, assuming a system in aqueous solution at 300 K. The MD simulation covering 400 ps indicated that the existence of a Cys148-Cys252 disulfide bond increases the conformational flexibility of the occluding loop, although the conformational stability of the overall structure is little affected. The structural characteristics of the complex elucidated by X-ray analysis were suggested to be also intrinsic and stable in the dynamic state; the hydrogen bonding/electrostatic interactions between the main and side chains of CA074 and the Sn and Sn' subsites of the enzyme were maintained throughout the MD simulation. Furthermore, the simulation made clear that the binding of CA074 significantly restricted the conformational flexibility of the substrate binding region, especially the occluding loop, of cathepsin B. Statistical analyses during the simulation suggest that the selectivity of CA074 for cathepsin B stems from the tight P1'-S1' and P2'-S2' interactions, assisted in particular by double hydrogen bonds between the carboxyl two oxygens of the CA074 C-terminus and the imidazole NH groups of His110 and His111 residues.

Animals↗
PubMed2006

Influence of isotherm inflection on the loading dependence of the diffusivities of n-hexane and n-heptane in MFI zeolite. Quasi-elastic neutron scattering experiments supplemented by molecular simulations.

Quasi-Elastic Neutron Scattering (QENS) experiments were carried out to determine (a) Fick diffusivity, D (b) self-diffusivity, Dself, and (c) 1/Gamma, the inverse of the thermodynamic correction factor, for n-hexane (nC6) and n-heptane (nC7) in MFI zeolite (all silica silicalite-1) at 300 K for a variety of loadings. These experimental results are compared with configurational-bias Monte Carlo (CBMC) and molecular dynamics (MD) simulations of, respectively, the adsorption isotherms and diffusivities. For n-hexane, the CBMC simulated isotherm shows a slight inflection at a loading=4 molecules per unit cell; this inflection manifests, also, in the loading dependence of 1/Gamma, obtained from QENS. The trend in the loading dependence of the Fick D and Dself of nC6 obtained from QENS matches the MD simulation results. For nC7 the CBMC simulated isotherm shows a strong inflection at a loading=4 molecules per unit cell. At this loading=4, 1/Gamma tends to zero and there is a very good match between QENS and molecular simulations for the loading dependence of 1/Gamma. Both MD simulations and QENS data on the Fick diffusivity shows a sharp maximum at a loading in the region of=4. For both nC6 and nC7 the simulated values of diffusivity are about an order of magnitude higher than those determined from QENS.

Journal Article↗
PubMed1994

X-ray crystal structure and molecular dynamics simulation of bovine pancreas phospholipase A2-n-dodecylphosphorylcholine complex.

The crystal structure of n-dodecylphosphorylcholine (n-C12PC)-bovine pancreas phospholipase A2 (PLA2) complex provided the following structural characteristics: (1) the dodecyl chain of n-C12PC was located at the PLA2 N-terminal helical region by hydrophobic interactions, which corresponds to the binding pocket of 2-acyl fatty acid chain (beta-chain) of the substrate phospholipid, (2) the region from Lys-53 to Lys-56 creates a choline-receiving pocket of n-C12PC and (3) the N-terminal group of Ala-1 shifts significantly toward the Tyr-52 OH group by the binding of the n-C12PC inhibitor. Since the accuracy of the X-ray analysis (R = 0.275 at 2.3 A resolution) was insufficient to establish these important X-ray insights, the complex structure was further investigated through the molecular dynamics (MD) simulation, assuming a system in aqueous solution at 310K. The MD simulation covering 176 ps showed that the structural characteristics observed by X-ray analysis are intrinsic and also stable in the dynamic state. Furthermore, the MD simulation made clear that the PLA2 binding pocket is large enough to permit the conformational fluctuation of the n-C12PC hydrocarbon chain.

Animals↗
PubMed2006

Li+ cation environment, transport, and mechanical properties of the LiTFSI doped N-methyl-N-alkylpyrrolidinium+TFSI- ionic liquids.

Molecular dynamics (MD) simulations have been performed on N-methyl-N-propylpyrrolidinium bis(trifluoromethanesulfonyl)imide (mppy(+)TFSI(-)) and N,N-dimethyl- pyrrolidinium bis(trifluoromethanesulfonyl)imide (mmpy(+)TFSI(+)) ionic liquids (ILs) doped with 0.25 mol fraction LiTFSI salt at 303-500 K. The liquid density, ion self-diffusion coefficients, and conductivity predicted by MD simulations were found to be in good agreement with experimental data, where available. MD simulations reveal that the Li(+) environment is similar in mppy(+)TFSI(-) and mmpy(+)TFSI(+) ILs doped with LiTFSI. The Li(+) cations were found to be coordinated on average by slightly less than four oxygen atoms with each oxygen atom being contributed by a different TFSI(-) anion. Significant lithium aggregation by sharing up to three TFSI(-) anions bridging two lithiums was observed, particularly at lower temperatures where the lithium aggregates were found to be stable for tens of nanoseconds. Polarization of TFSI(-) anions is largely responsible for the formation of such lithium aggregates. Li(+) transport was found to occur primarily by exchange of TFSI(-) anions in the first coordination shell with a smaller (approximately 30%) contribution also due to Li(+) cations diffusing together with their first coordination shell. In both ILs, ion self-diffusion coefficients followed the order Li(+) < TFSI(-) < mmpy(+) or mppy(+) with all ion diffusion in mmpy(+)TFSI(-) being systematically slower than that in mppy(+)TFSI(-). Conductivity due to the Li(+) cation in LiTFSI doped mppy(+)TFSI(-) IL was found to be greater than that for a model poly(ethylene oxide)(PEO)/LiTFSI polymer electrolyte but significantly lower than that for an ethylene carbonate/LiTFSI liquid electrolyte. Finally, the time-dependent shear modulus for the LiTFSI doped ILs was found to be similar to that for a model poly(ethylene oxide)(PEO)/LiTFSI polymer electrolyte on the subnanosecond time scale.

Journal Article↗
PubMed1998

Time resolved fluorescence properties of phenylalanine in different environments. Comparison with molecular dynamics simulation.

Time resolved fluorescence of the phenylalanine residue (Phe) alone and included in the transmembrane domain (TMD) sequences of the epidermal growth factor receptor (EGFR) and ErbB-2 was studied using the synchrotron radiation source of light, and compared to molecular dynamics (MD) simulations. The fluorescence intensity decay is strongly sensitive to the environment. A mono-exponential decay was obtained for Phe amino acid alone in two different solvents and for Phe included in EGFR transmembrane sequence, with fluorescence lifetime values varying from 1.7 ns (EGFR) to 7.4 ns (Phe dissolved in water). In ErbB-2 transmembrane sequence three lifetimes were detected. The relative amplitude of the shortest one (0.14 ns) is smaller than 10%, whereas the others (0.6 and 2.2 ns) are almost equally represented. They have been attributed to different rotamers exchanging slowly. This interpretation is supported by MD simulations which evidence transitions in time series of the chi 1 dihedral angle of Phe observed in the case of ErbB-2. The anisotropy decays are similar for both peptides and indicate the presence of a correlation time in the nanosecond range (1-4 ns) and the probable existence of a very fast one (< 0.05 ns). Autocorrelation functions computed from MD simulations corroborate these results.

Amino Acid Sequence↗
PubMed2006

Evaluation of efficiency of reference potential spatial warping algorithm in conformational sampling of peptides.

Molecular dynamics (MD) simulation using the reference potential spatial warping algorithm (REPSWA) is a promising method for computing the conformational equilibrium of a system with a rugged energy surface. Its effectiveness has previously been demonstrated using only a simplified model system of a hydrocarbon chain omitting nonbonded interaction terms from the potential energy function. To evaluate the applicability of REPSWA MD simulation to more realistic problems, we applied it to small peptide systems in an aqueous environment. Difficulties were encountered, however, forcing us to devise several modifications. We investigated their effectiveness in comparison to conventional constant-temperature and multicanonical MD simulations. We found that the sampling efficiency of the modified REPSWA MD after a careful optimization of its parameters was better than that of the constant-temperature MD and comparable to that of the multicanonical MD in several cases.

Alanine↗
PubMed2006

Statistical mechanically averaged molecular properties of liquid water calculated using the combined coupled cluster/molecular dynamics method.

Liquid water is investigated theoretically using combined molecular dynamics (MD) simulations and accurate electronic structure methods. The statistical mechanically averaged molecular properties of liquid water are calculated using the combined coupled cluster/molecular mechanics (CC/MM) method for a large number of configurations generated from MD simulations. The method includes electron correlation effects at the coupled cluster singles and doubles level and the use of a large correlation consistent basis set. A polarizable force field has been used for the molecular dynamics part in both the CC/MM method and in the MD simulation. We describe how the methodology can be optimized with respect to computational costs while maintaining the quality of the results. Using the optimized method we study the energetic properties including the heat of vaporization and electronic excitation energies as well as electric dipole and quadrupole moments, the frequency dependent electric (dipole) polarizability, and electric-field-induced second harmonic generation first and second hyperpolarizabilities. Comparisons with experiments are performed where reliable data are available. Furthermore, we discuss the important issue on how to compare the calculated microscopic nonlocal properties to the experimental macroscopic measurements.

Journal Article↗
PubMed2002

Computational study of the catalytic domain of human neutrophil collagenase. specific role of the S3 and S'3 subsites in the interaction with a phosphonate inhibitor.

Human neutrophil collagenase (HNC, MMP-8) is one of the target enzymes for drug treatment of pathologic extracellular matrix degradation. Peptidomimetic inhibitors bind in the S'-side of the enzyme active site occupying the S'1 primary specificity pocket by their large hydrophobic side-chains. The crystal structure of the complex between the catalytic domain of MMP-8 and Pro-Leu-L-TrpP(OH)2 (PLTP) showed that this phosphonate inhibitor binds in the S side of the active site. This finding was unexpected since it represents the first example of accommodation of the bulky Trp indolyl chain in the S1 rather than in the S'1 subsite. Dynamical and structural factors favouring this uncommon mode of binding were therefore investigated. MD simulations performed on the uncomplexed enzyme show that its structure in aqueous solution is only slightly different from the crystal structure found in the complex with PLTP. ED analysis of the MD simulations, performed on PLTP alternatively interacting with the S- or S'-side of the active site, shows that the enzyme fluctuation increases in both cases. The main contribution to the overall enzyme fluctuation is given by the loop 164-173. The fluctuation of this loop is spread over more degrees of freedom when PLTP interacts with the S-side. This dynamical factor can enhance the preference of PLTP for the S subsites of MMP-8. MD simulations also show that ligation of PLTP in the S subsites is further favoured by better zinc chelation, a cation-pi interaction at the S3 subsite and unstrained binding conformations. The role of the S3, S'3 and S'1 subsites in determining the inhibitor binding is discussed.

Catalytic Domain↗
PubMed2001

Structure, dynamics and electrostatics of the active site of glutaredoxin 3 from Escherichia coli: comparison with functionally related proteins.

The chemistry of active-site cysteine residues is central to the activity of thiol-disulfide oxidoreductases of the thioredoxin superfamily. In these reactions, a nucleophilic thiolate is required, but the associated pK(a) values differ vastly in the superfamily, from less than 4 in DsbA to greater than 7 in Trx. The factors that stabilize this thiolate are, however, not clearly established. The glutaredoxins (Grxs), which are members of this superfamily, contain a Cys-Pro-Tyr-Cys motif in their active site. In reduced Grxs, the pK(a) of the N-terminal active-site nucleophilic cysteine residue is lowered significantly, and the stabilization of the corresponding thiolate is expected to influence the redox potential of these enzymes. Here, we use a combination of long molecular dynamics (MD) simulations, pK(a) calculations, and experimental investigations to derive the structure and dynamics of the reduced active site from Escherichia coli Grx3, and investigate the factors that stabilize the thiolate. Several different MD simulations converged toward a consensus conformation for the active-site cysteine residues (Cys11 and Cys14), after a number of local conformational changes. Key features of the model were tested experimentally by measurement of NMR scalar coupling constants, and determination of pK(a) values of selected residues. The pK(a) values of the Grx3 active-site residues were calculated during the MD simulations, and support the underlying structural model. The structure of Grx3, in combination with the pK(a) calculations, indicate that the pK(a) of the N-terminal active-site cysteine residue in Grx3 is intermediate between that of its counterpart in DsbA and Trx. The pK(a) values in best agreement with experiment are obtained with a low (<4) protein dielectric constant. The calculated pK(a) values fluctuate significantly in response to protein dynamics, which underscores the importance of the details of the underlying structures when calculating pK(a) values. The thiolate of Cys11 is stabilized primarily by direct hydrogen bonding with the amide protons of Tyr13 and Cys14 and the thiol proton of Cys14, rather than by long-range interactions from charged groups or from a helix macrodipole. From the comparison of reduced Grx3 with other members of the thioredoxin superfamily, a unifying theme for the structural basis of thiol pK(a) differences in this superfamily begins to emerge.

Bacterial Proteins↗
PubMed2006

RNA unrestrained molecular dynamics ensemble improves agreement with experimental NMR data compared to single static structure: a test case.

Nuclear magnetic resonance (NMR) provides structural and dynamic information reflecting an average, often non-linear, of multiple solution-state conformations. Therefore, a single optimized structure derived from NMR refinement may be misleading if the NMR data actually result from averaging of distinct conformers. It is hypothesized that a conformational ensemble generated by a valid molecular dynamics (MD) simulation should be able to improve agreement with the NMR data set compared with the single optimized starting structure. Using a model system consisting of two sequence-related self-complementary ribonucleotide octamers for which NMR data was available, 0.3 ns particle mesh Ewald MD simulations were performed in the AMBER force field in the presence of explicit water and counterions. Agreement of the averaged properties of the molecular dynamics ensembles with NMR data such as homonuclear proton nuclear Overhauser effect (NOE)-based distance constraints, homonuclear proton and heteronuclear (1)H-(31)P coupling constant (J) data, and qualitative NMR information on hydrogen bond occupancy, was systematically assessed. Despite the short length of the simulation, the ensemble generated from it agreed with the NMR experimental constraints more completely than the single optimized NMR structure. This suggests that short unrestrained MD simulations may be of utility in interpreting NMR results. As expected, a 0.5 ns simulation utilizing a distance dependent dielectric did not improve agreement with the NMR data, consistent with its inferior exploration of conformational space as assessed by 2-D RMSD plots. Thus, ability to rapidly improve agreement with NMR constraints may be a sensitive diagnostic of the MD methods themselves.

Hydrogen Bonding↗
PubMed2005

The influence of amino acid protonation states on molecular dynamics simulations of the bacterial porin OmpF.

Several groups, including our own, have found molecular dynamics (MD) calculations to result in the size of the pore of an outer membrane bacterial porin, OmpF, to be reduced relative to its size in the x-ray crystal structure. At the narrowest portion of its pore, loop L3 was found to move toward the opposite face of the pore, resulting in decreasing the cross-section area by a factor of approximately 2. In an earlier work, we computed the protonation states of titratable residues for this system and obtained values different from those that had been used in previous MD simulations. Here, we show that MD simulations carried out with these recently computed protonation states accurately reproduce the cross-sectional area profile of the channel lumen in agreement with the x-ray structure. Our calculations include the investigation of the effect of assigning different protonation states to the one residue, D(127), whose protonation state could not be modeled in our earlier calculation. We found that both assumptions of charge states for D(127) reproduced the lumen size profile of the x-ray structure. We also found that the charged state of D(127) had a higher degree of hydration and it induced greater mobility of polar side chains in its vicinity, indicating that the apparent polarizability of the D(127) microenvironment is a function of the D(127) protonation state.

Algorithms↗
PubMed2001

Probing the differences between rat liver outer mitochondrial membrane cytochrome b5 and microsomal cytochromes b5.

Two distinct forms of cytochrome b5 exist in the rat hepatocyte. One is associated with the membrane of the endoplasmic reticulum (microsomal, or Mc, cyt b5) while the other is associated with the outer membrane of liver mitochondria (OM cyt b5). Rat OM cyt b5, the only OM cyt b5 identified so far, has a significantly more negative reduction potential and is substantially more stable toward chemical and thermal denaturation than Mc cytochromes b5. In addition, hemin is kinetically trapped in rat OM cyt b5 but not in the Mc proteins. As a result, no transfer of hemin from rat OM cyt b5 to apomyoglobin is observed at pH values as low as 5.2, nor can the thermodyamically favored ratio of hemin orientational isomers be achieved under physiologically relevant conditions. These differences are striking given the similarity of the respective protein folds. A combined theoretical and experimental study has been conducted in order to probe the structural basis behind the remarkably different properties of rat OM and Mc cytochromes b5. Molecular dynamics (MD) simulations starting from the crystal structure of bovine Mc cyt b5 revealed a conformational change that exposes several internal residues to the aqueous environment. The new conformation is equivalent to the "cleft-opened" intermediate observed in a previously reported MD simulation of bovine Mc cyt b5 [Storch, E. M., and Daggett, V. (1995) Biochemistry 34, 9682-9693]. The rat OM protein does not adopt a comparable conformation in MD simulations, thus restricting access of water to the protein interior. Subsequent comparisons of the protein sequences and structures suggested that an extended hydrophobic network encompassing the side chains of Ala-18, Ile-32, Leu-36, and Leu-47 might contribute to the inability of rat OM cyt b5 to adopt the cleft-opened conformation and, hence, stabilize its fold relative to the Mc isoforms. A corresponding network is not present in bovine Mc cyt b5 because positions 18, 32, and 47, are occupied by Ser, Leu, and Arg, respectively. To probe the roles played by Ala-18, Ile-32, and Leu-47 in endowing rat OM cyt b5 with its unusual structural properties, we have replaced them with the corresponding residues in bovine Mc cyt b5. Hence, the I32L (single), A18S/L47R (double), and A18S/L47R/I32L (triple) mutants of rat OM cyt b5 were prepared. The stability of these proteins was found to decrease in the following order: WT rat OM > rat OM I32L > rat OM A18S/L47R > rat OM A18S/L47R/I32L > bovine Mc cyt b5. The decrease in stability of the rat OM protein correlates with the extent to which the hydrophobic cluster involving the side chains of residues 18, 32, 36, and 47 has been disrupted. Complete disruption of the hydrophobic network in the triple mutant is confirmed in a 2.0 A resolution crystal structure of the protein. Disruption of the hydrophobic network also facilitates hemin loss at pH 5.2 for the double and triple mutants, with the less stable triple mutant exhibiting the greater rate of hemin transfer to apomyoglobin. Finally, 1H NMR spectroscopy and side-by-side comparisons of the crystal structures of bovine Mc, rat OM, and rat OM A18S/L47R/I32L cyt b5 allowed us to conclude that the nature of residue 32 plays a key role in controlling the relative stability of hemin orientational isomers A and B in rat OM cyt b5. A similar analysis led to the conclusion that Leu-70 and Ser-71 play a pivotal role in stabilizing isomer A relative to isomer B in Mc cytochromes b5.

Amino Acid Sequence↗