Molecular dynamics and free-energy calculations applied to affinity maturation in antibody 48G7

We investigated the relative free energies of hapten binding to the germ line and mature forms of the 48G7 antibody Fab fragments by applying a continuum model to structures sampled from molecular dynamics simulations in explicit solvent. Reasonable absolute and very good relative free energies were obtained. As a result of nine somatic mutations that do not contact the hapten, the affinity-matured antibody binds the hapten >104 tighter than the germ line antibody. Energetic analysis reveals that van der Waals interactions and nonpolar contributions to solvation are similar and drive the formations of both the germ line and mature antibody–hapten complexes. Affinity maturation of the 48G7 antibody therefore appears to occur through reorganization of the combining site geometry in a manner that optimizes the balance of gaining favorable electrostatic interactions with the hapten and losing those with solvent during the binding process. As reflected by lower rms fluctuations in the antibody–hapten complex, the mature complex undergoes more restricted fluctuations than the germ line complex. The dramatically increased affinity of the 48G7 antibody over its germ line precursor is thus made possible by electrostatic optimization.

Physiological response to long-term peripheral and central leptin infusion in lean and obese mice

Recent data have identified leptin as an afferent signal in a negative-feedback loop regulating the mass of the adipose tissue. High leptin levels are observed in obese humans and rodents, suggesting that, in some cases, obesity is the result of leptin insensitivity. This hypothesis was tested by comparing the response to peripherally and centrally administered leptin among lean and three obese strains of mice: diet-induced obese AKR/J, New Zealand Obese (NZO), and Ay. Subcutaneous leptin infusion to lean mice resulted in a dose-dependent loss of body weight at physiologic plasma levels. Chronic infusions of leptin intracerebroventricularly (i.c.v.) at doses of 3 ng/hr or greater resulted in complete depletion of visible adipose tissue, which was maintained throughout 30 days of continuous i.c.v. infusion. Direct measurement of energy balance indicated that leptin treatment did not increase total energy expenditure but prevented the decrease that follows reduced food intake. Diet-induced obese mice lost weight in response to peripheral leptin but were less sensitive than lean mice. NZO mice were unresponsive to peripheral leptin but were responsive to i.c.v. leptin. Ay mice did not respond to subcutaneous leptin and were 1/100 as sensitive to i.c.v. leptin. The decreased response to leptin in diet-induced obese, NZO, and Ay mice suggests that obesity in these strains is the result of leptin resistance. In NZO mice, leptin resistance may be the result of decreased transport of leptin into the cerebrospinal fluid, whereas in Ay mice, leptin resistance probably results from defects downstream of the leptin receptor in the hypothalamus.

A streptavidin mutant with altered ligand-binding specificity

The biotin-binding site of streptavidin was modified to alter its ligand-binding specificity. In natural streptavidin, the side chains of N23 and S27 make two of the three hydrogen bonds with the ureido oxygen of biotin. These two residues were mutated to severely weaken biotin binding while attempting to maintain the affinity for two biotin analogs, 2-iminobiotin and diaminobiotin. Redesigning of the biotin-binding site used the difference in local electrostatic charge distribution between biotin and these biotin analogs. Free energy calculations predicted that the introduction of a negative charge at the position of S27 plus the mutation N23A should disrupt two of the three hydrogen bonds between natural streptavidin and the ureido oxygen of biotin. In contrast, the imino hydrogen of 2-iminobiotin should form a hydrogen bond with the side chain of an acidic amino acid at position 27. This should reduce the biotin-binding affinity by approximately eight orders of magnitude, while leaving the affinities for these biotin analogs virtually unaffected. In good agreement with these predictions, a streptavidin mutant with the N23A and S27D substitutions binds 2-iminobiotin with an affinity (Ka) of 1 × 106 M−1, two orders of magnitude higher than that for biotin (1 × 104 M−1). In contrast, the binding affinity of this streptavidin mutant for diaminobiotin (2.7 × 104 M−1) was lower than predicted (2.9 × 105 M−1), suggesting the position of the diaminobiotin in the biotin-binding site was not accurately determined by modeling.

Physical reasons for the unusual α-helix stabilization afforded by charged or neutral polar residues in alanine-rich peptides

We have carried out conformational energy calculations on alanine-based copolymers with the sequence Ac-AAAAAXAAAA-NH2 in water, where X stands for lysine or glutamine, to identify the underlying source of stability of alanine-based polypeptides containing charged or highly soluble polar residues in the absence of charge–charge interactions. The results indicate that ionizable or neutral polar residues introduced into the sequence to make them soluble sequester the water away from the CO and NH groups of the backbone, thereby enabling them to form internal hydrogen bonds. This solvation effect dictates the conformational preference and, hence, modifies the conformational propensity of alanine residues. Even though we carried out simulations for specific amino acid sequences, our results provide an understanding of some of the basic principles that govern the process of folding of these short sequences independently of the kind of residues introduced to make them soluble. In addition, we have investigated through simulations the effect of the bulk dielectric constant on the conformational preferences of these peptides. Extensive conformational Monte Carlo searches on terminally blocked 10-mer and 16-mer homopolymers of alanine in the absence of salt were carried out assuming values for the dielectric constant of the solvent ɛ of 80, 40, and 2. Our simulations show a clear tendency of these oligopeptides to augment the α-helix content as the bulk dielectric constant of the solvent is lowered. This behavior is due mainly to a loss of exposure of the CO and NH groups to the aqueous solvent. Experimental evidence indicates that the helical propensity of the amino acids in water shows a dramatic increase on addition of certain alcohols, such us trifluoroethanol. Our results provide a possible explanation of the mechanism by which alcohol/water mixtures affect the free energy of helical alanine oligopeptides relative to nonhelical ones.

Predicting crossover generation in DNA shuffling

We introduce a quantitative framework for assessing the generation of crossovers in DNA shuffling experiments. The approach uses free energy calculations and complete sequence information to model the annealing process. Statistics obtained for the annealing events then are combined with a reassembly algorithm to infer crossover allocation in the reassembled sequences. The fraction of reassembled sequences containing zero, one, two, or more crossovers and the probability that a given nucleotide position in a reassembled sequence is the site of a crossover event are estimated. Comparisons of the predictions against experimental data for five example systems demonstrate good agreement despite the fact that no adjustable parameters are used. An in silico case study of a set of 12 subtilases examines the effect of fragmentation length, annealing temperature, sequence identity and number of shuffled sequences on the number, type, and distribution of crossovers. A computational verification of crossover aggregation in regions of near-perfect sequence identity and the presence of synergistic reassembly in family DNA shuffling is obtained.

Dynamical hierarchy in transition states: Why and how does a system climb over the mountain?

How a reacting system climbs through a transition state during the course of a reaction has been an intriguing subject for decades. Here we present and quantify a technique to identify and characterize local invariances about the transition state of an N-particle Hamiltonian system, using Lie canonical perturbation theory combined with microcanonical molecular dynamics simulation. We show that at least three distinct energy regimes of dynamical behavior occur in the region of the transition state, distinguished by the extent of their local dynamical invariance and regularity. Isomerization of a six-atom Lennard–Jones cluster illustrates this: up to energies high enough to make the system manifestly chaotic, approximate invariants of motion associated with a reaction coordinate in phase space imply a many-body dividing hypersurface in phase space that is free of recrossings even in a sea of chaos. The method makes it possible to visualize the stable and unstable invariant manifolds leading to and from the transition state, i.e., the reaction path in phase space, and how this regularity turns to chaos with increasing total energy of the system. This, in turn, illuminates a new type of phase space bottleneck in the region of a transition state that emerges as the total energy and mode coupling increase, which keeps a reacting system increasingly trapped in that region.

Stability of Coulomb systems in a magnetic field.

I study N electrons and M protons in a magnetic field. It is shown that the total energy per particle is bounded below by a constant independent of M and N, provided the fine structure constant is small. Here, the total energy includes the energy of the magnetic field.

High resolution ultrastructural mapping of total calcium: electron spectroscopic imaging/electron energy loss spectroscopy analysis of a physically/chemically processed nerve-muscle preparation.

We report on a procedure for tissue preparation that combines thoroughly controlled physical and chemical treatments: quick-freezing and freeze-drying followed by fixation with OsO4 vapors and embedding by direct resin infiltration. Specimens of frog cutaneous pectoris muscle thus prepared were analyzed for total calcium using electron spectroscopic imaging/electron energy loss spectroscopy (ESI/EELS) approach. The preservation of the ultrastructure was excellent, with positive K/Na ratios revealed in the fibers by x-ray microanalysis. Clear, high-resolution EELS/ESI calcium signals were recorded from the lumen of terminal cisternae of the sarcoplasmic reticulum but not from longitudinal cisternae, as expected from previous studies carried out with different techniques. In many mitochondria, calcium was below detection whereas in others it was appreciable although at variable level. Within the motor nerve terminals, synaptic vesicles as well as some cisternae of the smooth endoplasmic reticulum yielded positive signals at variance with mitochondria, that were most often below detection. Taken as a whole, the present study reveals the potential of our experimental approach to map with high spatial resolution the total calcium within individual intracellular organelles identified by their established ultrastructure, but only where the element is present at high levels.