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.

Ionophores and receptors using cation-π interactions: Collarenes

Cation-π interactions are important forces in molecular recognition by biological receptors, enzyme catalysis, and crystal engineering. We have harnessed these interactions in designing molecular systems with circular arrangement of benzene units that are capable of acting as ionophores and models for biological receptors. [n]Collarenes are promising candidates with high selectivity for a specific cation, depending on n, because of their structural rigidity and well-defined cavity size. The interaction energies of [n]collarenes with cations have been evaluated by using ab initio calculations. The selectivity of these [n]collarenes in aqueous solution was revealed by using statistical perturbation theory in conjunction with Monte Carlo and molecular dynamics simulations. It has been observed that in [n]collarenes the ratio of the interaction energies of a cation with it and the cation with the basic building unit (benzene) can be correlated to its ion selectivity. We find that collarenes are excellent and efficient ionophores that bind cations through cation-π interactions. [6]Collarene is found to be a selective host for Li+ and Mg2+, [8]collarene for K+ and Sr2+, and [10]collarene for Cs+ and Ba2+. This finding indicates that [10]collarene and [8]collarene could be used for effective separation of highly radioactive isotopes, 137Cs and 90Sr, which are major constituents of nuclear wastes. More interestingly, collarenes of larger cavity size can be useful in capturing organic cations. [12]Collarene exhibits a pronounced affinity for tetramethylammonium cation and acetylcholine, which implies that it could serve as a model for acetylcholinestrase. Thus, collarenes can prove to be novel and effective ionophores/model-receptors capable of heralding a new direction in molecular recognition and host-guest chemistry.

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.

Molecular evidence for multiple origins of woodiness and a New World biogeographic connection of the Macaronesian Island endemic Pericallis (Asteraceae: Senecioneae)

The prevalence of woody species in oceanic islands has attracted the attention of evolutionary biologists for more than a century. We used a phylogeny based on sequences of the internal-transcribed spacer region of nuclear ribosomal DNA to trace the evolution of woodiness in Pericallis (Asteraceae: Senecioneae), a genus endemic to the Macaronesian archipelagos of the Azores, Madeira, and Canaries. Our results show that woodiness in Pericallis originated independently at least twice in these islands, further weakening some previous hypotheses concerning the value of this character for tracing the continental ancestry of island endemics. The same data suggest that the origin of woodiness is correlated with ecological shifts from open to species-rich habitats and that the ancestor of Pericallis was an herbaceous species adapted to marginal habitats of the laurel forest. Our results also support Pericallis as closely related to New World genera of the tribe Senecioneae.

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.

Molecular mechanism for silencing virally transduced genes involves histone deacetylation and chromatin condensation

Virally transduced genes are often silenced after integration into the host genome. Chromatin immunoprecipitation and nuclease sensitivity experiments now demonstrate that silencing of the transgene is characterized by deacetylation of histone H4 lysines and chromatin condensation. Trichostatin A treatment results in dramatic reactivation of gene expression that is preceded by histone acetylation and chromatin decondensation. Analysis of individual histone H4 lysines demonstrate that chromatin domain opening is coincident with rapid acetylation of histone H4 K5, K12, and K16 and that maintenance of the open domain is correlated with acetylation of histone H4 K8. Removal of trichostatin A results in rapid deacetylation of histone H4 K8, chromatin condensation, and transcription silencing. The results suggest that deacetylation of histone H4 lysines and coincident chromatin condensation are critically involved in the silencing of virally transduced genes.

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.

Ozone Sensitivity in Hybrid Poplar Is Correlated with a Lack of Defense-Gene Activation1

Ozone is a major gaseous pollutant thought to contribute to forest decline. Although the physiological and morphological responses of forest trees to ozone have been well characterized, little is known about the molecular basis for these responses. Our studies compared the response to ozone of ozone-sensitive and ozone-tolerant clones of hybrid poplar (Populus maximowizii × Populus trichocarpa) at the physiological and molecular levels. Gas-exchange analyses demonstrated clear differences between the ozone-sensitive clone 388 and the ozone-tolerant clone 245. Although ozone induced a decrease in photosynthetic rate and stomatal conductance in both clones, the magnitude of the decrease in stomatal conductance was significantly greater in the ozone-tolerant clone. RNA-blot analysis established that ozone-induced mRNA levels for phenylalanine ammonia-lyase, O-methyltransferase, a pathogenesis-related protein, and a wound-inducible gene were significantly higher in the ozone-tolerant than in the ozone-sensitive plants. Wound- and pathogen-induced levels of these mRNAs were also higher in the ozone-tolerant compared with the ozone-sensitive plants. The different physiological and molecular responses to ozone exposure exhibited by clones 245 and 388 suggest that ozone tolerance involves the activation of salicylic-acid- and jasmonic-acid-mediated signaling pathways, which may be important in triggering defense responses against oxidative stress.

Influence of Water Content and Temperature on Molecular Mobility and Intracellular Glasses in Seeds and Pollen1

Although the occurrence of intracellular glasses in seeds and pollen has been established, physical properties such as rotational correlation times and viscosity have not been studied extensively. Using electron paramagnetic resonance spectroscopy, we examined changes in the molecular mobility of the hydrophilic nitroxide spin probe 3-carboxy-proxyl during melting of intracellular glasses in axes of pea (Pisum sativum L.) seeds and cattail (Typha latifolia L.) pollen. The rotational correlation time of the spin probe in intracellular glasses of both organisms was approximately 10−3 s. Using the distance between the outer extrema of the electron paramagnetic resonance spectrum (2Azz) as a measure of molecular mobility, we found a sharp increase in mobility at a definite temperature during heating. This temperature increased with decreasing water content of the samples. Differential scanning calorimetry data on these samples indicated that this sharp increase corresponded to melting of the glassy matrix. Molecular mobility was found to be inversely correlated with storage stability. With decreasing water content, the molecular mobility reached a minimum, and increased again at very low water content. Minimum mobility and maximum storage stability occurred at a similar water content. This correlation suggests that storage stability might be at least partially controlled by molecular mobility. At low temperatures, when storage longevity cannot be determined on a realistic time scale, 2Azz measurements can provide an estimate of the optimum storage conditions.

Purification and Characterization of Peroxidases Correlated with Lignification in Poplar Xylem1

Lignin is an integral cell wall component of all vascular plants. Peroxidases are widely believed to catalyze the last enzymatic step in the biosynthesis of lignin, the dehydrogenation of the p-coumaryl alcohols. As the first stage in identifying lignin-specific peroxidase isoenzymes, the classical anionic peroxidases found in the xylem of poplar (Populus trichocarpa Trichobel) were purified and characterized. Five different poplar xylem peroxidases (PXP 1, PXP 2, PXP 3–4, PXP 5, and PXP 6) were isolated. All five peroxidases were strongly glycosylated (3.6% to 4.9% N-glucosamine), with apparent molecular masses between 46 and 54 kD and pI values between pH 3.1 and 3.8. Two of the five isolated peroxidases (PXP 3–4 and PXP 5) could oxidize the lignin monomer analog syringaldazine, an activity previously correlated with lignification in poplar. Because these isoenzymes were specifically or preferentially expressed in xylem, PXP 3–4 and PXP 5 are suggested to be involved in lignin polymerization.

Molecular and Physiological Responses to Water Deficit in Drought-Tolerant and Drought-Sensitive Lines of Sunflower1 : Accumulation of Dehydrin Transcripts Correlates with Tolerance

To investigate correlations between phenotypic adaptation to water limitation and drought-induced gene expression, we have studied a model system consisting of a drought-tolerant line (R1) and a drought-sensitive line (S1) of sunflowers (Helianthus annuus L.) subjected to progressive drought. R1 tolerance is characterized by the maintenance of shoot cellular turgor. Drought-induced genes (HaElip1, HaDhn1, and HaDhn2) were previously identified in the tolerant line. The accumulation of the corresponding transcripts was compared as a function of soil and leaf water status in R1 and S1 plants during progressive drought. In leaves of R1 plants the accumulation of HaDhn1 and HaDhn2 transcripts, but not HaElip1 transcripts, was correlated with the drought-adaptive response. Drought-induced abscisic acid (ABA) concentration was not associated with the varietal difference in drought tolerance. Stomata of both lines displayed similar sensitivity to ABA. ABA-induced accumulation of HaDhn2 transcripts was higher in the tolerant than in the sensitive genotype. HaDhn1 transcripts were similarly accumulated in the tolerant and in the sensitive plants in response to ABA, suggesting that additional factors involved in drought regulation of HaDhn1 expression might exist in tolerant plants.

N2 fixation in marine heterotrophic bacteria: dynamics of environmental and molecular regulation.

Molecular and immunological techniques were used to examine N2 fixation in a ubiquitous heterotrophic marine bacterium, the facultative anaerobic Vibrio natriegens. When batch cultures were shifted from aerobic N-replete to anaerobic N-deplete conditions, transcriptional and post-translational regulation of N2 fixation was observed. Levels of nifHDK mRNA encoding the nitrogenase enzyme were highest at 140 min postshift and undetectable between 6 and 9 h later. Immunologically determined levels of nitrogenase enzyme (Fe protein) were highest between 6 and 15 h postshift, and nitrogenase activity peaked between 6 and 9 h postshift, declining by a factor of 2 after 12-15 h. Unlike their regulation in cyanobacteria, Fe protein and nitrogenase activity were present when nifHDK mRNA was absent in V. natriegens, indicating that nitrogenase is stored and stable under anaerobic conditions. Both nifHDK mRNA and Fe protein disappeared within 40 min after cultures were shifted from N2-fixing conditions (anaerobic, N-deplete) to non- N2-fixing conditions (aerobic, N-enriched) but reappeared when shifted to conditions favoring N2 fixation. Thus, unlike other N2-fixing heterotrophic bacteria, nitrogenase must be resynthesized after aerobic exposure in V. natriegens. Immunological detection based on immunoblot (Western) analysis and immunogold labeling correlated positively with nitrogenase activity; no localization of nitrogenase was observed. Because V. natriegens continues to fix N2 for many hours after anaerobic induction, this species may play an important role in providing "new" nitrogen in marine ecosystems.

Molecular cloning and functional characterization of the Xenopus Ca(2+)-binding protein frequenin.

Frequenin was originally identified in Drosophila melanogaster as a Ca(2+)-binding protein facilitating transmitter release at the neuromuscular junction. We have cloned the Xenopus frequenin (Xfreq) by PCR using degenerate primers combined with low-stringency hybridization. The deduced protein has 70% identity with Drosophila frequenin and about 38-58% identity with other Ca(2+)-binding proteins. The most prominent features are the four EF-hands, Ca(2+)-binding motifs. Xfreq mRNA is abundant in the brain and virtually nondetectable from adult muscle. Western blot analysis indicated that Xfreq is highly concentrated in the adult brain and is absent from nonneural tissues such as heart and kidney. During development, the expression of the protein correlated well with the maturation of neuromuscular synapses. To determine the function of Xfreq at the developing neuromuscular junction, the recombinant protein was introduced into Xenopus embryonic spinal neurons by early blastomere injection. Synapses made by spinal neurons containing exogenous Xfreq exhibited a much higher synaptic efficacy. These results provide direct evidence that frequenin enhances transmitter release at the vertebrate neuromuscular synapse and suggest its potential role in synaptic development and plasticity.

Heterogeneity in molecular recognition by restriction endonucleases: osmotic and hydrostatic pressure effects on BamHI, Pvu II, and EcoRV specificity.

The cleavage specificity of the Pvu II and BamHI restriction endonucleases is found to be dramatically reduced at elevated osmotic pressure. Relaxation in specificity of these otherwise highly accurate and specific enzymes, previously termed "star activity," is uniquely correlated with osmotic pressure between 0 and 100 atmospheres. No other colligative solvent property exhibits a uniform correlation with star activity for all of the compounds tested. Application of hydrostatic pressure counteracts the effects of osmotic pressure and restores the natural selectivity of the enzymes for their canonical recognition sequences. These results indicate that water solvation plays an important role in the site-specific recognition of DNA by many restriction enzymes. Osmotic pressure did not induce an analogous effect on the specificity of the EcoRV endonuclease, implying that selective hydration effects do not participate in DNA recognition in this system. Hydrostatic pressure was found to have little effect on the star activity induced by changes in ionic strength, pH, or divalent cation, suggesting that distinct mechanisms may exist for these observed alterations in specificity. Recent evidence has indicated that BamHI and EcoRI share similar structural motifs, while Pvu II and EcoRV belong to a different structural family. Evidently, the use of hydration water to assist in site-specific recognition is a motif neither limited to nor defined by structural families.