Second-contact shell mutation diminishes streptavidin-biotin binding affinity through transmitted effects on equilibrium dynamics.

We report a point mutation in the second contact shell of the high-affinity streptavidin-biotin complex that appears to reduce binding affinity through transmitted effects on equilibrium dynamics. The Y54F streptavidin mutation causes a 75-fold loss of binding affinity with 73-fold faster dissociation, a large loss of binding enthalpy (ΔΔH = 3.4 kcal/mol at 37 °C), and a small gain in binding entropy (TΔΔS = 0.7 kcal/mol). The removed Y54 hydroxyl is replaced by a water molecule in the bound structure, but there are no observable changes in structure in the first contact shell and no additional changes surrounding the mutation. Molecular dynamics simulations reveal a large increase in the atomic fluctuation amplitudes for W79, a key biotin contact residue, compared to the fluctuation amplitudes in the wild-type. The increased W79 atomic fluctuation amplitudes are caused by loss of water-mediated hydrogen bonds between the Y54 hydroxyl group and peptide backbone atoms in and near W79. We propose that the increased atomic fluctuation amplitudes diminish the integrity of the W79-biotin interaction and represents a loosening of the "tryptophan collar" that is critical to the slow dissociation and high affinity of streptavidin-biotin binding. These results illustrate how changes in protein dynamics distal to the ligand binding pocket can have a profound impact on ligand binding, even when equilibrium structure is unperturbed.

Theoretical determination of the standard reduction potentials of pheophytin-a in N,N-dimethyl formamide and membrane.

Quantum mechanical/molecular mechanics (QM/MM) calculations were performed on the neutral, anionic, and dianionic forms of Pheophytin-a (Pheo-a) in N,N-dimethyl formamide (DMF) in order to calculate the absolute free energy of reduction of Pheo-a in solution. The geometry of the solvated species was optimized by restricted open-shell density functional treatment (ROB3LYP) using the 6-31G(d) basis set for the molecular species while the primary solvent shell consisting of 45 DMF molecules was treated by the MM method using the universal force field (UFF). Electronic energies of the neutral, anionic, and dianionic species were obtained by carrying out single point density functional theory (DFT) calculations using the 6-311+G(2d,2p) basis set on the respective ONIOM optimized geometries. The CHARMM27 force field was used to account for the dynamical nature of the primary solvation shell of 45 DMF molecules. In the calculations using solvent shells, the required atomic charges for each solvent molecule were obtained from ROB3LYP/6-31G(d) calculation on a single isolated DMF molecule. Randomly sampled configurations generated by the Monte Carlo (MC) technique were used to determine the contribution of the primary shell to the free energy of solvation of the three species. The dynamical nature of the primary shell significantly corrects the free energy of solvation. Frequency calculations at the ROB3LYP/6-31G(d) level were carried out on the optimized geometries of truncated 47-atom models of the neutral and ionic species in vacuum so as to determine the differences in thermal energy and molecular entropy. The Born energy of ion-dielectric interaction, the Onsager energy of dipole-dielectric interaction, and the Debye-Hückel energy of ion-ionic cloud interaction for the pheophytin-solvent aggregate were added as perturbative corrections. The Born interaction also makes a large contribution to the absolute free energy of reduction. An implicit solvation model (DPCM) was also employed for the calculation of standard reduction potentials in DMF. Both the models were successful in reproducing the standard reduction potentials. An explicit solvent treatment(QM/MM/MC + Born + Onsager + Debye corrections) yielded the one electron reduction potential of Pheo-a as -0.92 +/- 0.27 V and the two electron reduction potential as -1.34 +/- 0.25 V at 298.15 K, while the implicit solvent treatment yielded the corresponding values as -1.03 +/- 0.17 and -1.30 +/- 0.17 V, respectively. The calculated values more or less agree with the experimental midpoint potentials of -0.90 and -1.25 V, respectively. Moreover, a numerical finite difference Poisson-Boltzmann solver (FDPB) along with the DPCM methodology was employed to obtain the reduction potential of pheophytin in the thylakoid membrane. The calculated reduction potential value of -0.58 V is in excellent agreement with the reported value -0.61 V.

Unique rate expression for glucose production in C(4) plants.

A physicochemical interpretation of a recently formulated temperature-dependent, steady-state rate expression for the production of glucose equivalent in C(4) plants is given here. We show that the rate equation is applicable to a wide range of C(4) plants.

Integrated kinetics for the production of glucose in plant cells and the effect of temperature.

We prepare a temperature-dependent formulation of the integrated kinetics for the overall process of photosynthesis in eukaryotic cells. To avoid complexity, the C4 plants are chosen because their rate of photosynthesis is independent of the partial pressure of O2. A systematically simplified but comprehensive scheme for both light and dark reactions is considered. The reaction rate per reaction center in the thylakoid membrane is related to the rate of exciton transfer between chlorophyll neighbors. An expression is formulated for the light reaction rate (R1'). The NADPH formation rate is related to R1' and the survival probability of the membrane. Rates of different steps in the simplified scheme can be related to each other by applying a few steady state conditions. The saturation probability of CO2 in a bundle sheath is also considered. The photochemical efficiency (phi) appears in terms of these probabilities. We find the glucose production rate as R(glucose) = (8/3) upsilon L: [corrected] R1'phi g(T)([G3P]/[P(i)]2) exp(-deltaG(E)S/RT), where g(T) is the activation quotient of the involved enzymes, G3P and P(i) represent glyceraldehyde-3-phosphate and inorganic phosphates, and deltaG(E)S is the free energy for the apparent equilibrium between G3P and glucose. This is the first time that such a comprehensive expression for R(glucose) has been derived. The probabilities are generally given by sigmoid curves. The corresponding parameters can be easily determined. The quotient g(T) incorporates a Gaussian distribution for temperature dependence and a sigmoid function describing deactivation. The theoretical plots of photochemical efficiency and glucose production rate versus temperature are in excellent agreement with the experimental ones, thereby validating the formalism.