PubMed Health⌕ Search

Biomedical subjects

S H Northrup

Publications and source records attributed to S H Northrup.

16 recordsLinked to original sources
PubMed1993

Effects of charged amino acid mutations on the bimolecular kinetics of reduction of yeast iso-1-ferricytochrome c by bovine ferrocytochrome b5.

The reduction of wild-type yeast iso-1-ferricytochrome c (ycytc) and several mutants by trypsin-solubilized bovine liver ferrocytochrome b5 (cytb5) has been studied under conditions in which the electron-transfer reaction is bimolecular. The effect of electrostatic charge modifications and steric changes on the kinetics has been determined by experimental and theoretical observations of the electron-transfer rates of ycytc mutants K79A, K'72A, K79A/K'72A, and R38A (K' is used to signify trimethyllysine (Tml)). A structurally robust Brownian dynamics (BD) method simulating diffusional docking and electron transfer was employed to predict the mutation effect on the rate constants. A realistic model of the electron-transfer event embodied in an intrinsic unimolecular rate constant is used which varies exponentially with donor-acceptor distance. The BD method quantitatively predicts rate constants over a considerable range of ionic strengths. Semiquantitative agreement is obtained in predicting the perturbing influence of the mutations on the rate constants. Both the experimentally observed rate constants and those predicted by BD descend in the following order: native ycytc > K79A > K'72A > K79A/K'72A. Variant R38A was studied at a different ionic strength than this series of mutations, and the theory agreed with experiment in predicting a smaller rate constant for the mutant. In all cases the predicted effect of mutation was in the correct direction, but not as large as that observed. The BD simulations predict that the two proteins dock through essentially a single domain, with a distance of closest approach of the two heme groups in rigid body docking typically around 12 A. Two predominant classes of complexes were calculated, the most frequent involving the quartet of cytb5/ycytc interactions, Glu48-Arg13, Glu56-Lys87, Asp60-Lys86, and heme-Tml72, having an average electrostatic energy of -13.0 kcal/mol. The second most important complexes were of the type previously postulated (Salemme, 1976; Mauk et al., 1986; Rodgers et al., 1988) with interactions Glu44-Lys27, Glu48-Arg13, Asp60-Tml72, and heme-Lys79 and having an energy of -6.4 kcal/mol. The ionic strength dependence of the bimolecular reaction rate was well reproduced using a discontinuous dielectric model, but poorly so for a uniform dielectric model.

Amino Acid Sequence↗
PubMed1992

Kinetics of protein-protein association explained by Brownian dynamics computer simulation.

Protein-protein bond formations, such as antibody-antigen complexation or aggregation of protein monomers into dimers and larger aggregates, occur with bimolecular rate constants on the order of 10(6) M-1.s-1, which is only 3 orders of magnitude slower than the diffusion-limited Smoluchowski rate. However, since the protein-protein bond requires rotational alignment to within a few angstroms of tolerance, purely geometric estimates would suggest that the observed rates might be 6 orders of magnitude below the Smoluchowski rate. Previous theoretical treatments have not been solved for the highly specific docking criteria of protein-protein association--the entire subunit interface must be aligned within 2 A of the correct position. Several studies have suggested that diffusion alone could not produce the rapid association kinetics and have postulated "lengthy collisions" and/or the operation of electrostatic or hydrophobic steering forces to accelerate the association. In the present study, the Brownian dynamics simulation method is used to compute the rate of association of neutral spherical model proteins with the stated docking criteria. The Brownian simulation predicts a rate of 2 x 10(6) M-1.s-1 for this generic protein-protein association, a rate that is 2000 times faster than that predicted by the simplest geometric calculation and is essentially equal to the rates observed for protein-protein association in aqueous solution. This high rate is obtained by simple diffusive processes and does not require any attractive or steering forces beyond those achieved for a partially formed bond. The rate enhancement is attributed to a diffusive entrapment effect, in which a protein pair surrounded and trapped by water undergoes multiple collisions with rotational reorientation during each encounter.

Computer Simulation↗
PubMed1991

Reduction of horse heart ferricytochrome c by bovine liver ferrocytochrome b5. Experimental and theoretical analysis.

The reduction of horse heart ferricytochrome c by the tryptic fragment of bovine liver cytochrome b5 and its dimethyl ester heme (DME)-substituted derivative has been studied as a function of ionic strength, pH, and temperature under solution conditions where the reaction is bimolecular. The rate constant for ferricytochrome c reduction by native ferrocytochrome b5 is 1.8 (+/- 0.2) x 10(7) M-1 s-1 (25 degrees C) with delta H++ = 7.5 (+/- 0.2) kcal/mol and delta S++ = -0.3 (+/- 0.6) eu (pH 7.0, I = 0.348 M). Under the same solution conditions, the reduction of ferricytochrome c by DME-ferrocytochrome b5 proceeds with a rate constant of 1.7 (+/- 0.1) x 10(7) M-1 s-1 with delta H++ = 7.9 (+/- 0.4) kcal/mol and delta S++ = 1 (+/- 1) eu. The rate constants for both reactions are strongly dependent on ionic strength. A detailed electrostatic analysis of the proteins has been performed. Two relatively simple Brownian dynamics simulation models predict rate constants for the reaction between the two native proteins that demonstrate a dependence on ionic strength similar to that observed experimentally. In one of these models, the proteins are treated as spheres with reactive surface patches that are defined by a 5 degrees cone generated about the dipole vector calculated for each protein and aligned with the presumed electron-transfer site near the partially exposed heme edge. The second model replaces the reactive patch assumption with an exponential distance dependence for the probability of reaction that permits estimation of a value for the distance-dependence factor alpha. Calculations with this latter model in combination with the aligned dipole assumption provide a reasonable approximation to the observed ionic strength dependence for the reaction and are consistent with a value of alpha = 1.2 A-1.

Animals↗
PubMed1990

Electrostatic field around cytochrome c: theory and energy transfer experiment.

Energy transfer in the "rapid diffusion" limit from electronically excited terbium(III) chelates in three different charge states to horse heart ferricytochrome c was measured as a function of ionic strength. Theoretical rate constants calculated by numerical integration of the Forster integral (containing the Poisson-Boltzmann-generated protein electrostatic potential) were compared with the experimental data to evaluate the accuracy of protein electrostatic field calculations at the protein/solvent interface. Two dielectric formalisms were used: a simple coulombic/Debye-Hückel procedure and a finite difference method [Warwicker, J. & Watson, H. C. (1982) J. Mol. Biol. 157, 671-679] that accounts for the low-dielectric protein interior and the irregular protein/solvent boundary. Good agreement with experiment was obtained and the ionic-strength dependence of the reaction was successfully reproduced. The sensitivity of theoretical rate constants to the choices of effective donor sphere size and the energy transfer distance criterion was analyzed. Electrostatic potential and rate-constant calculations were carried out on sets of structures collected along two molecular dynamics trajectories of cytochrome c. Protein conformational fluctuations were shown to produce large variations in the calculated energy transfer rate constant. We conclude that protein fluctuations and the resulting transient structures can play significant roles in biological or catalytic activities that are not apparent from examination of a static structure. For calculating protein electrostatics, large-scale low-frequency conformational fluctuations, such as charged side-chain reorientation, are established to be as important as the computational method for incorporating dielectric boundary effects.

Animals↗
PubMed1988

Brownian dynamics of cytochrome c and cytochrome c peroxidase association.

Brownian dynamics computer simulations of the diffusional association of electron transport proteins cytochrome c (cyt c) and cytochrome c peroxidase (cyt c per) were performed. A highly detailed and realistic model of the protein structures and their electrostatic interactions was used that was based on an atomic-level spatial description. Several structural features played a role in enhancing and optimizing the electron transfer efficiency of this reaction. Favorable electrostatic interactions facilitated long-lived nonspecific encounters between the proteins that allowed the severe orientational criteria for reaction to be overcome by rotational diffusion during encounters. Thus a "reduction-in-dimensionality" effect operated. The proteins achieved plausible electron transfer orientations in a multitude of electrostatically stable encounter complexes, rather than in a single dominant complex.

Chemical Phenomena↗
PubMed1988

Brownian dynamics simulation of protein association.

The Brownian Dynamics (BD) method is applied to study the diffusive dynamics and interaction of two proteins, cytochrome c (CYTC) and cytochrome c peroxidase (CYP). We examine the role of protein electrostatic charge distribution in the facilitation of protein-protein docking prior to the electron transfer step, assessing the influence of individual charged amino acid residues. Accurate interaction potentials are computed by iterating the linearized Poisson-Boltzmann (PB) equation around the larger protein CYP. The low dielectric constant inside proteins, electrolyte screening effects and irregular protein surface topography are taken into account. We observe a large ensemble of electrostatically stable encounter complexes seemingly with acceptable geometric requirements for electron transfer rather than a single dominant complex. Stabilities of the large variety of docking complexes are rationalized in terms of generalized charged residue complementarities. However, it is found that the electrostatic interactions giving rise to complex stabilities are somewhat nonspecific in nature. A large series of additional simulations are performed in which individual charged residues on CYTC have been chemically modified. Resulting perturbations of the association rate are significant and qualitatively similar to results observed in comparable kinetics experiments. We therefore demonstrate the potential of the Brownian dynamics method to estimate the effects of site-directed mutagenesis on protein-protein and protein-ligand diffusional association rates.

Amino Acids↗
PubMed1986

Simulation of biomolecular diffusion and complex formation.

Diffusion is a phenomenon of very widespread importance in molecular biophysics. Diffusion can determine the rates and character of the assembly of multisubunit structures, the binding of ligands to receptors, and the internal motions of molecules and assemblies that involve solvent surface displacements. Current computer simulation techniques provide much more detailed descriptions of diffusional processes than have been available in the past. Models can be constructed to include such realistic features as structural subunits at the submolecular level (domains, monomers, or atoms); detailed electrostatic charge distributions and corresponding solvent-screened inter- and intramolecular interactions; and hydrodynamic interactions. The trajectories can be analyzed either to provide direct information on biomolecular function (e.g., the bimolecular rate constant for formation of an electron-transfer complex between two proteins), or to provide or test models for the interpretation of experimental data (e.g., the time dependence of fluorescence depolarization for segments of DNA). Here, we first review the theory of diffusional simulations, with special emphasis on new techniques such as those for obtaining transport properties of flexible assemblies and rate constants of diffusion-controlled reactions. Then we survey a variety of recent applications, including studies of large-scale motion in DNA segments and substrate "steering" in enzyme-substrate binding. We conclude with a discussion of current work (e.g., formation of protein complexes) and possible areas for future work.

Biopolymers↗
PubMed1982

Dynamical theory of activated processes in globular proteins.

A method is described for calculating the reaction rate in globular proteins of activated processes such as ligand binding or enzymatic catalysis. The method is based on the determination of the probability that the system is in the transition state and of the magnitude of the reactive flux for transition-state systems. An "umbrella sampling" simulation procedure is outlined for evaluating the transition-state probability. The reactive flux is obtained from an approach described previously for calculating the dynamics of transition-state trajectories. An application to the rotational isomerization of an aromatic ring in the bovine pancreatic trypsin inhibitor is presented. The results demonstrate the feasibility of calculating rate constants for reactions in proteins and point to the importance of solvent effects for reactions that occur near the protein surface.

Chemical Phenomena↗
PubMed1980

Internal mobility of ferrocytochrome c.

In the refinement of the X-ray diffraction structures of molecules, it is conventional to introduce atomic 'temperature factors' of the Debye-Waller form to characterize the widths of the electron density peaks corresponding to the atoms. Although these factors are known to include a variety of contributions other than thermal fluctuations of the atomic positions, recent progress in the refinement of protein structures has led to inferences concerning atomic mobilities from the temperature factor data for several proteins. Atomic position fluctuations can be calculated independently by the molecular dynamics method, in which the classical equations of motion for the atoms of an equilibrated protein are solved on a computer. We now show that the X-ray diffraction and dynamical simulation methods yield similar pictures of the atomic mobility in tuna ferrocytochrome c.

Animals↗