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.

Molecular mechanics calculations of the riboacetal internucleotide linkage in double and triple helices.

Structures of Watson-Crick base paired 15-nucleobase oligomer strands in A-type or B-type conformation in which one strand [a strand of alternating nucleotide and riboacetal thymidine nucleoside (RT) units, RP] is DNA and the other is composed of alternating nucleotides and riboacetal nucleosides have been studied by molecular mechanics. Analogously, oligomer strands of RNA in place of DNA have been modeled. The calculations indicate that the RP strand is more stable when complexed in an A-type duplex relative to a B-type form and that this conformational preference is presumably due to the more uniform nature of the former. Nearly planar ribose rings were more commonly observed in the minimized structures of the B-type DNA.RP duplexes as compared with A-type duplexes, despite the fact that planar ribofuranose rings are known to be energetically unfavorable in oligonucleotides. Computed relative stabilities of all duplexes containing the RP strand suggest that such heteroduplexes are less stable than the corresponding double-stranded DNA and double-stranded RNA species. These findings are in agreement with experimental results which show, when equivalent sequences were compared, that a DNA.RNA control forms a more stable duplex than RP hound to a complementary single-stranded RNA strand. In contrast, molecular mechanics studies of complementary triple-helical (DNA)2.RP, (DNA)2.DNA, and (DNA)2.RNA structures indicate that the binding of RP as a Hoogsteen strand stabilizes the underlying duplex to a greater extent compared with native oligonucleotides. These calculations suggest that puckering of the ribose ring in the riboacetal linkage leads to a more favorable interaction with a complementary nucleic acid target than the proposed planar geometry and that this puckering may account for the enhanced binding of RP to a double-stranded target.

Direct molecular level measurements of the electrostatic properties of a protein surface

In this work, we used direct measurements with the surface force apparatus to determine the pH-dependent electrostatic charge density of a single binding face of streptavidin. Mean field calculations have been used with considerable success to model electrostatic potential fields near protein surfaces, but these models and their inherent assumptions have not been tested directly at the molecular level. Using the force apparatus and immobilized, oriented monolayers of streptavidin, we measured a pI of 5–5.5 for the biotin-binding face of the protein. This differs from the pI of 6.3 for the soluble protein and confirms that we probed the local electrostatic features of the macromolecule. With finite difference solutions of the linearized Poisson–Boltzmann equation, we then calculated the pH-dependent charge densities adjacent to the same face of the protein. These calculated values agreed quantitatively with those obtained by direct force measurements. Although our study focuses on the pH-dependence of surface electrostatics, this direct approach to probing the electrostatic features of proteins is applicable to investigations of any perturbations that alter the charge distribution of the surfaces of immobilized molecules.

Femtosecond observation of benzyne intermediates in a molecular beam: Bergman rearrangement in the isolated molecule

In this communication, we report our femtosecond real-time observation of the dynamics for the three didehydrobenzene molecules (p-, m-, and o-benzyne) generated from 1,4-, 1,3-, and 1,2-dibromobenzene, respectively, in a molecular beam, by using femtosecond time-resolved mass spectrometry. The time required for the first and the second C-Br bond breakage is less than 100 fs; the benzyne molecules are produced within 100 fs and then decay with a lifetime of 400 ps or more. Density functional theory and high-level ab initio calculations are also reported herein to elucidate the energetics along the reaction path. We discuss the dynamics and possible reaction mechanisms for the disappearance of benzyne intermediates. Our effort focuses on the isolated molecule dynamics of the three isomers on the femtosecond time scale.

The importance of reactant positioning in enzyme catalysis: A hybrid quantum mechanics/molecular mechanics study of a haloalkane dehalogenase

Hybrid quantum mechanics/molecular mechanics calculations using Austin Model 1 system-specific parameters were performed to study the SN2 displacement reaction of chloride from 1,2-dichloroethane (DCE) by nucleophilic attack of the carboxylate of acetate in the gas phase and by Asp-124 in the active site of haloalkane dehalogenase from Xanthobacter autotrophicus GJ10. The activation barrier for nucleophilic attack of acetate on DCE depends greatly on the reactants having a geometry resembling that in the enzyme or an optimized gas-phase structure. It was found in the gas-phase calculations that the activation barrier is 9 kcal/mol lower when dihedral constraints are used to restrict the carboxylate nucleophile geometry to that in the enzyme relative to the geometries for the reactants without dihedral constraints. The calculated quantum mechanics/molecular mechanics activation barriers for the enzymatic reaction are 16.2 and 19.4 kcal/mol when the geometry of the reactants is in a near attack conformer from molecular dynamics and in a conformer similar to the crystal structure (DCE is gauche), respectively. This haloalkane dehalogenase lowers the activation barrier for dehalogenation of DCE by 2–4 kcal/mol relative to the single point energies of the enzyme's quantum mechanics atoms in the gas phase. SN2 displacements of this sort in water are infinitely slower than in the gas phase. The modest lowering of the activation barrier by the enzyme relative to the reaction in the gas phase is consistent with mutation experiments.