PubMed Health⌕ Search

Biomedical subjects

R Elber

Publications and source records attributed to R Elber.

27 records · Page 2Linked to original sources
PubMed1993

Ligand binding and conformation change in the dimeric hemoglobin of the clam Scapharca inaequivalvis.

The reaction with carbon monoxide of the cooperative dimeric hemoglobin from Scapharca inaequivalvis has been examined by flash photolysis. In the nanosecond time range, geminate rebinding of 5% of dissociated CO occurs with a rate constant of 1.4 x 10(7) s-1. There is a change in absorbance of deoxyhemoglobin following photolysis at a rate of 1.2 x 10(6) s-1, consistent with a shift in the position of the Soret band to longer wavelengths. The amplitude of the change is proportional to the population of deoxydimer. In much of the Soret region this change is greater than the absorbance excursion associated with geminate recombination. There is at least one other slower change associated with the singly liganded species. Geminate rebinding of NO has components of 50, 8, and 0.035 ns-1, accounting for 75%, 25%, and less than 1% of the total reaction observed after a 35-ps photolysis flash. Simulation of diffusion of NO by molecular dynamics shows the ligands moving from the heme pocket to a subsidiary space between the edge of the heme and the surface of the protein.

Animals↗
PubMed1992

Distal pocket residues affect picosecond ligand recombination in myoglobin. An experimental and molecular dynamics study of position 29 mutants.

Time courses for intramolecular NO and O2 recombination to native and three position 29 mutants of sperm whale myoglobins were measured after laser photolysis on picosecond and nanosecond time scales. The rates for the first phase of NO recombination were 1.8, 2.5, 29, and > or = 100 ns-1 for Ala29, Val29, Leu29 (native), and Phe29 myoglobin, respectively, at room temperature. This order is not correlated with the overall association rate constants for NO binding which were all in the range 20-50 x 10(6) M-1 s-1 and is the opposite of that observed for the rate constants for the overall thermal dissociation of NO which were 5.0, 2.8, 0.98, and 0.21 x 10(-4) s-1 for Ala29, Val29, Leu29 (native), and Phe29 myoglobin, respectively, at 20 degrees C. This inverse correlation suggests that photo- and thermally dissociated ligand molecules experience similar kinetic and equilibrium barriers to rebinding. The larger side chains of Leu29 and Phe29 inhibit rapid movement of the ligand away from the iron atom facilitating geminate recombination. The smaller side chains of Val29 and Ala29 increase the space available to the ligand, decreasing the rate of geminate recombination and enhancing complete dissociation. Diffusion of NO in the distal pocket of myoglobin was simulated using a variant of the molecular dynamics program CHARMM that includes the locally enhanced sampling protocol (Elber, R., and Karplus, M. (1991) J. Am. Chem. Soc. 112, 9161-9175; Roitberg, A., and Elber, R. (1991) J. Chem. Phys. 95, 9277-9287) and the x-ray structures of Carver et al. (Carver, T. E., Brantley, R. E., Jr., Singleton, E. W., Arduini, R. M., Quillin, M. L., Phillips, G. N., Jr., and Olson, J. S. (1992) J. Biol. Chem. 267, 14443-14450). Both accelerated (5,000 K) and room temperature ligands were used, and comparisons were made between simulations with a complete hydration shell surrounding the protein and those with only eight water molecules near the distal histidine. Photodissociated ligands initially move away from the heme plane, past Leu29, and toward Leu32, Phe33, Ile107, and Ile111. These theoretical results confirm that a complete description of picosecond ligand recombination must include the dynamics of ligand movement in the distal portion of the heme pocket.

Animals↗
PubMed1991

Computational studies of ligand diffusion in globins: I. Leghemoglobin.

The thermally assisted diffusion of a small ligand (carbon monoxide) through a protein matrix (lupine leghemoglobin) is investigated computationally. The diffusion paths are calculated by a variant of the time-dependent Hartree approximation which we call LES (locally enhanced sampling). The variant which was recently introduced by Elber and Karplus is based on the classical TDSCF approximation of Gerber et al. The simulation enables more significant search for diffusion pathways than was possible before. This is done by increasing the number of ligand trajectories using a single trajectory for the protein. We compare qualitatively diffusion rates in leghemoglobin and in myoglobin. The calculation shows that the diffusion in leghemoglobin is much faster than the diffusion in myoglobin, in agreement with experiment. The gate in leghemoglobin is opened by fluctuations at a close contact between the B/C and the G helices. The most relevant fluctuation is the rigid shift of the C helix with respect to the G helix. This path is not observed in a comparable calculation for myoglobin. This finding is rationalized by the lack of the D helix in leghemoglobin and a significantly more flexible CE loop. Supporting experimental evidence for the importance of the CE loop in leghemoglobin can be found in the kinetics studies of Gibson et al.

Amino Acid Sequence↗
PubMed1990

Molecular dynamics study of secondary structure motions in proteins: application to myohemerythrin.

The concept of secondary structure motions is examined in a molecular dynamics simulation of the protein myohemerythrin. We extracted from the simulation a corresponding trajectory of helices and demonstrated that the fluctuations of the protein are dominated by a rigid shift of these secondary structure elements. The relative motions of the helices are irregular, with no clear periodicity. They are bounded by approximately 2 A for the center of mass motions and by 20 degrees for the relative orientations. The potential of mean force for the interactions of the helices was calculated, and the correlations between the different extended motions were investigated. It is shown that the one-dimensional mean force potentials are close to quadratic for most of the helices coordinates. The anharmonicity is reflected by changes in the direction of the normal modes as a function of the energy and by the existence of multiple free energy minima for the helices packing. The multiple conformations are associated with a single type of secondary structure coordinate: the angle that describes the relative orientation of the helices in a plane perpendicular to the line connecting their center of mass.

Binding Sites↗
PubMed1989

Reaction path study of conformational transitions and helix formation in a tetrapeptide.

Conformational transitions between the 112 stable states of the tetrapeptide isobutyryl-Ala3-NH-methyl (IAN) are studied theoretically. The objective of the investigation is to advance the understanding of helix formation and of conformational transitions in polypeptides. The possible reaction paths between extended chain and helical configurations are examined in detail. The study of the multiple reaction paths in this 48-atom molecule became possible due to development of a new computational algorithm. It is shown that the helix-coil transitions in IAN follow a sequence of local dihedral flips and that the number of the available routes for the transition is significantly lower than in a random search. A quasi-melting point is obtained at 5 +/- 1 kcal (1 cal = 4.18 J)/mol above the lowest energy minimum. Below this point the molecule is trapped in one or very few minima, and above it the molecule hops between a large number of configurations.

Calorimetry↗
PubMed1987

Multiple conformational states of proteins: a molecular dynamics analysis of myoglobin.

A molecular dynamics simulation of myoglobin provides the first direct demonstration that the potential energy surface of a protein is characterized by a large number of thermally accessible minima in the neighborhood of the native structure (for example, approximately 2000 minima were sampled in a 300-picosecond trajectory). This is expected to have important consequences for the interpretation of the activity of transport proteins and enzymes. Different minima correspond to changes in the relative orientation of the helices coupled with side-chain rearrangements that preserve the close packing of the protein interior. The conformational space sampled by the simulation is similar to that found in the evolutionary development of the globins. Glasslike behavior is expected at low temperatures. The minima obtained from the trajectory do not satisfy certain criteria for ultrametricity.

Models, Structural↗