Anti-insulin antibody structure and conformation. I. Molecular modeling and mechanics of an insulin antibody.

A knowledge-based three-dimensional model of an anti-insulin antibody, 125, was constructed using the structures of conserved residues found in other known crystallographic immunoglobulins. Molecular modeling and mechanics were done with the 125 amino acid sequences using QUANTA and CHARMm on a Silicon Graphics 4D70GT workstation. A minimal model was made by scaffolding using crystallography coordinates of the antibody HyHEL-5, because it had the highest amino acid sequence homology with 125 (84% light chain, 65% heavy chain). The three hypervariable loop turns that are longer in 125 than in HyHEL-5 (L1, L3, and H3) were modeled separately and incorporated into the HyHEL-5 structure; then other amino acid substitutions were made and torsions optimized. The 125 model maintains all the structural attributes of an antibody and the structures conserved in known antibodies. Although there are many polar amino acids (especially serines) in this site, the overall van der Waals surface shape is determined by positions of aromatic side chains. Based on this model, it is suggested that hydrogen bonding may be key in the interaction between the human insulin A chain loop antigenic epitope and 125.     

Water and polypeptide conformations in the gramicidin channel. A molecular dynamics study.

Theoretical studies of ion channels address several important questions. The mechanism of ion transport, the role of water structure, the fluctuations of the protein channel itself, and the influence of structural changes are accessible from these studies. In this paper, we have carried out a 70-ps molecular dynamics simulation on a model structure of gramicidin A with channel waters. The backbone of the protein has been analyzed with respect to the orientation of the carbonyl and the amide groups. The results are in conformity with the experimental NMR data. The structure of water and the hydrogen bonding network are also investigated. It is found that the water molecules inside the channel act as a collective chain; whereas the conformation in which all the waters are oriented with the dipoles pointing along the axis of the channel is a preferred one, others are also accessed during the dynamics simulation. A collective coordinate involving the channel waters and some of the hydrogen bonding peptide partners is required to describe the transition of waters from one configuration to the other.     

Metabolism of [3-13C]pyruvate in TCA cycle mutants of yeast.

The utilization of pyruvate and acetate by Saccharomyces cerevisiae was examined using 13C and 1H NMR methodology in intact wild-type yeast cells and mutant yeast cells lacking Krebs tricarboxylic acid (TCA) cycle enzymes. These mutant cells lacked either mitochondrial (NAD) isocitrate dehydrogenase (NAD-ICDH1),alpha-ketoglutarate dehydrogenase complex (alpha KGDC), or mitochondrial malate dehydrogenase (MDH1). These mutant strains have the common phenotype of being unable to grow on acetate. [3-13C]-Pyruvate was utilized efficiently by wild-type yeast with the major intermediates being [13C]glutamate, [13C]acetate, and [13C]alanine. Deletion of any one of these Krebs TCA cycle enzymes changed the metabolic pattern such that the major synthetic product was [13C]galactose instead of [13C]glutamate, with some formation of [13C]acetate and [13C]alanine. The fact that glutamate formation did not occur readily in these mutants despite the metabolic capacity to synthesize glutamate from pyruvate is difficult to explain. We discuss the possibility that these data support the metabolon hypothesis of Krebs TCA cycle enzyme organization.     

Gating of the active site of triose phosphate isomerase: Brownian dynamics simulations of flexible peptide loops in the enzyme.

The enzyme triose phosphate isomerase has flexible peptide loops at its active sites. The loops close over these sites upon substrate binding, suggesting that the dynamics of the loops could be of mechanistic and kinetic importance. To investigate these issues, the loop motions in the dimeric enzyme were simulated by Brownian dynamics. The two loops, one on each monomer, were represented by linear chains of appropriately parameterized spheres, each sphere corresponding to an amino acid residue. The loops moved in the electrostatic field of the rest of the enzyme, which was held rigid in its crystallographically observed conformation. In the absence of substrate, the loops exhibited gating of the active site with a period of about 1 ns and occupied "closed" conformations for about half of the time. As the period of gating is much shorter than the enzyme-substrate relaxation time, the motion of the loops does not reduce the rate constant for the approach of substrate from its simple diffusion-controlled value. This suggests that the flexible loops may have evolved to create the appropriate environment for catalysis while, at the same time, minimizing the kinetic penalty for gating the active site.     

Acetylcholinesterase: diffusional encounter rate constants for dumbbell models of ligand.

For some enzymes, virtually every substrate molecule that encounters the entrance to the active site proceeds to reaction, at low substrate concentrations. Such diffusion-limited enzymes display high apparent bimolecular rate constants ((kcat/KM)), which depend strongly upon solvent viscosity. Some experimental studies provide evidence that acetylcholinesterase falls into this category. Interestingly, the asymmetric charge distribution of acetylcholinesterase, apparent from the crystallographic structure, suggests that its electrostatic field accelerates the encounter of its cationic substrate, acetylcholine, with the entrance to the active site. Here we report simulations of the diffusion of substrate in the electrostatic field of acetylcholinesterase. We find that the field indeed guides the substrate to the mouth of the active site. The computed encounter rate constants depend upon the particular relative geometries of substrate and enzyme that are considered to represent successful encounters. With loose reaction criteria, the computed rates exceed those measured experimentally, but the rate constants vary appropriately with ionic strength. Although more restrictive reaction criteria lower the computed rates, they also lead to unrealistic variation of the rate constants with ionic strength. That these simulations do not agree well with experiment suggests that the simple diffusion model is incomplete. Structural fluctuations in the enzyme or events after the encounter may well contribute to rate limitation.     

Prediction of titration properties of structures of a protein derived from molecular dynamics trajectories.

This paper explores the dependence of the molecular dynamics (MD) trajectory of a protein molecule on the titration state assigned to the molecule. Four 100-ps MD trajectories of bovine pancreatic trypsin inhibitor (BPTI) were generated, starting from two different structures, each of which was held in two different charge states. The two starting structures were the X-ray crystal structure and one of the solution structures determined by NMR, and the charge states differed only in the ionization state of N terminus. Although it is evident that the MD simulations were too short to sample fully the equilibrium distribution of structures in each case, standard Poisson-Boltzmann titration state analysis of the resulting configurations shows general agreement between the overall titration behavior of the protein and the charge state assumed during MD simulation: at pH 7, the total net charge of the protein resulting from the titration analysis is consistently lower for the protein with the N terminus assumed to be neutral than for the protein with the N terminus assumed to be charged. For most of the ionizable residues, the differences in the calculated pKaS among the four trajectories are statistically negligible and remain in good agreement with the data obtained by crystal structure titration and by experiment. The exceptions include the N terminus, which responds directly to the change of its imposed charge; the C terminus, which in the NMR structure interacts strongly with the former; and a few other residues (Arg 1, Glu 7, Tyr 35, and Arg 42) whose pKaS reflect the initial structure and the limited trajectory lengths. This study illustrates the importance of the careful assignment of protonation states at the start of MD simulations and points to the need for simulation methods that allow for the variation of the protonation state in the calculation of equilibrium properties.     

Genetics of the “umbrella” branching habit in Capsicum annuum L.

Summary Three inbred lines, MSU 78-101, MSU 79-221 and MSU 74-230 were used to determine the inheritance of the umbrella branching habit in peppers. MSU 79-221, with the umbrella phenotype, was crossed with MSU 78-191 (dwarf) and MSU 74-230 (indeterminate growth habit). Segregating populations were separated on the basis of plant growth habit and fruit bearing habit. Genetic analyses suggested that the umbrella phenotype was controlled by three major recessive genes, ct and dt determining plant habit, and fa determining fruit bearing habit. When the dominant alleles Dt and Ct were in the dominant homozygous or heterozygous condition an indeterminate phenotype was produced. Su, a dominant suppressor gene, apparently acts to suppress the epistatic action of the Ct gene. Modifiers were involved in the control of branching in the umbrella plants. Linkage also was noted between the genes for indeterminate plant habit and non-clustered bearing habit. The information derived from this study will allow the plant breeder to design an efficient breeding program for the development of pepper cultivars suitable for mechanical harvesting systems.Michigan Agricultural Experiment Station Journal Article No. 11073     

Molecular dynamics computations and solid state nuclear magnetic resonance of the gramicidin cation channel.

This paper reports on a coupled approach to determining the structure of the gramicidin A ion channel, utilizing solid state nuclear magnetic resonance (NMR) of isotopically labeled gramicidin channels aligned parallel to the magnetic field direction, and molecular dynamics (MD). MD computations using an idealized right-handed beta-helix as a starting point produce a refined molecular structure that is in excellent agreement with atomic resolution solid state NMR data. The data provided by NMR and MD are complementary to each other. When applied in a coordinated manner they provide a powerful approach to structure determination in molecular systems not readily amenable to x-ray diffraction.     

A molecular dynamics study of thermodynamic and structural aspects of the hydration of cavities in proteins.

The structure and activity of a protein molecule are strongly influenced by the extent of hydration of its cavities. This is, in turn, related to the free energy change on transfer of a water molecule from bulk solvent into a cavity. Such free energy changes have been calculated for two cavities in a sulfate-binding protein. One of these cavities contains a crystallographically observed water molecule while the other does not. Thermodynamic integration and perturbation methods were used to calculate free energies of hydration for each of the cavities from molecular dynamics simulations of two separate events: the removal of a water molecule from pure water, and the introduction of a water molecule into each protein cavity. From the simulations for the pure water system, the excess chemical potential of water was computed to be -6.4 +/- 0.4 kcal/mol, in accord with experiment and with other recent theoretical calculations. For the protein cavity containing an experimentally observed water molecule, the free energy change on hydrating it with one water molecule was calculated as -10.0 +/- 1.3 kcal/mol, indicating the high probability that this cavity is occupied by a water molecule. By contrast, for the cavity in which no water molecules were experimentally observed, the free energy change on hydrating it with one water molecule was calculated as 0.2 +/- 1.5 kcal/mol, indicating its low occupancy by water. The agreement of these results with experiment suggests that thermodynamic simulation methods may become useful for the prediction and analysis of internal hydration in proteins.     

pH‐dependent conformational dynamics of beta‐secretase 1: A molecular dynamics study

Beta‐secretase 1 (BACE‐1) is an aspartyl protease implicated in the overproduction of β‐amyloid fibrils responsible for Alzheimer disease. The process of β‐amyloid genesis is known to be pH dependent, with an activity peak between solution pH of 3.5 and 5.5. We have studied the pH‐dependent dynamics of BACE‐1 to better understand the pH dependent mechanism. We have implemented support for graphics processor unit (GPU) accelerated constant pH molecular dynamics within the AMBER molecular dynamics software package and employed this to determine the relative population of different aspartyl dyad protonation states in the pH range of greatest β‐amyloid production, followed by conventional molecular dynamics to explore the differences among the various aspartyl dyad protonation states. We observed a difference in dynamics between double‐protonated, mono‐protonated, and double‐deprotonated states over the known pH range of higher activity. These differences include Tyr 71‐aspartyl dyad proximity and active water lifetime. This work indicates that Tyr 71 stabilizes catalytic water in the aspartyl dyad active site, enabling BACE‐1 activity.