Computer modeling of actinomycin D interactions with double-helical DNA

We have performed molecular mechanical calculations on intercalation complexes of actinomycin D with a series of base-paired hexanucleoside pentaphosphates; d(GCGCGC)2, d(GCCGGC)2, d(GCATGC)2, d(GCTAGC)2 and d(ATGCAT)2. Our results are in good agreement with previous experimental work on sequence selectivity. The results provide a rationalization for the strong preference of actinomycin D to intercalate on the 3' side of guanine residues, consistent with previously proposed models. Finally, the computed structures for d(ATGCAT)2-actinomycin D complexes have been compared with two-dimensional nuclear magnetic resonance nuclear Overhauser effect experimental results. To our knowledge, this is the first extensive comparison of molecular mechanical model structures for a drug-DNA complex with experimental solution phase data. We find generally good agreement between our computational models and the experimental solution phase structures.     

A highly conserved nuclear gene for low-level phylogenetics: elongation factor-1 alpha recovers morphology-based tree for heliothine moths

Molecular systematists need increased access to nuclear genes. Highly conserved, low copy number protein-encoding nuclear genes have attractive features for phylogenetic inference but have heretofore been applied mostly to very ancient divergences. By virtue of their synonymous substitutions, such genes should contain a wealth of information about lower-level taxonomic relationships as well, with the advantage that amino acid conservatism makes both alignment and primer definition straightforward. We tested this postulate for the elongation factor-1 alpha (EF-1 alpha) gene in the noctuid moth subfamily Heliothinae, which has probably diversified since the middle Tertiary. We sequenced 1,240 bp in 18 taxa representing heliothine groupings strongly supported by previous morphological and allozyme studies. The single most parsimonious gene tree and the neighbor-joining tree for all nucleotides show almost complete concordance with the morphological tree. Homoplasy and pairwise divergence levels are low, transition/transversion ratios are high, and phylogenetic information is spread evenly across gene regions. The EF-1 alpha gene and presumably other highly conserved genes hold much promise for phylogenetics of Tertiary age eukaryote groups.      

Plant virus infection development as affected by heavy metal stress

The work was focused on the investigation of possible dependencies between the development of viral infection in plants and the presence of high heavy metal concentrations in soil. Field experiments have been conducted in order to study the development of systemic tobacco mosaic virus (TMV) infection in Lycopersicon esculentum L. cv. Miliana plants under effect of separate salts of heavy metals Cu, Zn and Pb deposited in soil. As it is shown, simultaneous effect of viral infection and heavy metals in tenfold maximum permissible concentration leads to decrease of total chlorophyll content in experiment plants mainly due to the degradation of chlorophyll a. The reduction of chlorophyll concentration under the combined influence of both stress factors was more serious comparing to the separate effect of every single factor. Plants' treatment with toxic concentrations of lead and zinc leaded to slight delay in the development of systemic TMV infection together with more than twofold increase of virus content in plants that may be an evidence of synergism between these heavy metal's and virus' effects. Contrary, copper although decreased total chlorophyll content but showed protective properties and significantly reduced amount of virus in plants.     

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.     

Simulation of nitroxide electron paramagnetic resonance spectra from brownian trajectories and molecular dynamics simulations

A simulated continuous wave electron paramagnetic resonance spectrum of a nitroxide spin label can be obtained from the Fourier transform of a free induction decay. It has been previously shown that the free induction decay can be calculated by solving the time-dependent stochastic Liouville equation for a set of Brownian trajectories defining the rotational dynamics of the label. In this work, a quaternion-based Monte Carlo algorithm has been developed to generate Brownian trajectories describing the global rotational diffusion of a spin-labeled protein. Also, molecular dynamics simulations of two spin-labeled mutants of T4 lysozyme, T4L F153R1, and T4L K65R1 have been used to generate trajectories describing the internal dynamics of the protein and the local dynamics of the spin-label side chain. Trajectories from the molecular dynamics simulations combined with trajectories describing the global rotational diffusion of the protein are used to account for all of the dynamics of a spin-labeled protein. Spectra calculated from these combined trajectories correspond well to the experimental spectra for the buried site T4L F153R1 and the helix surface site T4L K65R1. This work provides a framework to further explore the modeling of the dynamics of the spin-label side chain in the wide variety of labeling environments encountered in site-directed spin labeling studies.     

Three-dimensional structure for the beta 2 adrenergic receptor protein based on computer modeling studies.

Computer-aided model building techniques have been used to construct three-dimensional model structures for hamster beta 2 adrenergic receptor. Experimental data were used as constraints to guide the model building procedure, and a number of rather strict criteria were applied to assess the physical plausibility of model structures. We present details of our best model structure to date, which is consistent with a large body of experimental data. We also discuss in detail our model building procedures and evaluation criteria, which we believe may be of general utility in modeling projects.     

Molecular modeling of flexible arm‐mediated interactions between bacterial chemoreceptors and their modification enzyme

Sensory adaptation in bacterial chemotaxis is mediated by methylation and demethylation of specific glutamyl residues in the cytoplasmic domain of chemoreceptors. Methylation is catalyzed by methyltransferase CheR. In E. coli and related organisms, methylation sufficiently rapid to be physiologically effective requires a carboxyl terminal pentapeptide sequence on the receptor being modified or, via adaptational assistance, on a neighboring homodimer in a receptor cluster. Pentapeptide‐enhanced methylation is thought to be mediated by a ～30 residue, potentially disordered sequence that serves as a flexible arm connecting the receptor body and pentapeptide‐bound methyltransferase, thus allowing diffusionally restricted enzyme to reach methyl‐accepting sites. However, it was not known how many or which sites on the same or neighboring receptors were accessible to the tethered enzyme. We investigated using molecular modeling and found that, in a hexagonal array of trimers of receptor dimers, CheR tethered to a dimer of chemoreceptor Tar by its native 30‐residue flexible‐arm sequence could reach all methyl‐accepting sites on the dimer to which it was tethered plus 48 methyl‐accepting sites distributed among nine neighboring dimers, equivalent to the total sites carried by six receptors. This modeling‐determined methylation neighborhood of one enzyme‐binding dimer and six neighbors corresponds precisely with the experimentally identified neighborhood of seven. Thus, the experimentally observed adaptational assistance can occur by docking of pentapeptide‐bound, diffusionally restricted enzyme to methyl‐accepting sites on neighboring receptors. Our analysis introduces the notion that physiologically relevant adaptational assistance could occur even if only a subset of sites on a particular receptor are within reach.     

A net ecosystem carbon budget for snow dominated forested headwater catchments: linking water and carbon fluxes to critical zone carbon storage

Climate-driven changes in carbon (C) cycling of forested ecosystems have the potential to alter long-term C sequestration and the global C balance. Prior studies have shown that C uptake and partitioning in response to hydrologic variation are system specific, suggesting that a comprehensive assessment is required for distinct ecosystems. Many sub-humid montane forest ecosystems in the US are projected to experience increased water limitation over the next decades and existing water-limited forests can be used as a model for how changes in the hydrologic cycle will impact such ecosystems more broadly. Toward that goal we monitored precipitation, net ecosystem exchange and lateral soil and stream C fluxes in three semi-arid to sub-humid montane forest catchments for several years (WY 2009–2013) to investigate how the amount and timing of water delivery affect C stores and fluxes. The key control on aqueous and gaseous C fluxes was the distribution of water between winter and summer precipitation, affecting ecosystem C uptake versus heterotrophic respiration. We furthermore assessed C stores in soil and above- and below-ground biomass to assess how spatial patterns in water availability influence C stores. Topographically-driven patterns in catchment wetness correlated with modeled soil C stores, reflecting both long-term trends in local C uptake as well as lateral redistribution of C leached from upslope organic soil horizons to convergent landscape positions. The results suggest that changes in the seasonality of precipitation from winter snow to summer rain will influence both the amount and the spatial distribution of soil C stores.     

Streptavidin-biotin binding energetics

The high affinity energetics in the streptavidin-biotin system provide an excellent model system for studying how proteins balance enthalpic and entropic components to generate an impressive overall free energy for ligand binding. We review here concerted site-directed mutagenesis, biophysical, and computational studies of aromatic and hydrogen bonding interaction energetics between streptavidin and biotin. These results also have provided insight into how streptavidin builds a large activation barrier to dissociation by managing the enthalpic and entropic activation components. Finally, we review recent studies of the biotin dissociation pathway that address the fundamental question of how ligands exit protein binding pockets.     

Ovalbumin(323-339) peptide binds to the major histocompatibility complex class II I-A(d) protein using two functionally distinct registers

Proteins of the class II major histocompatibility complex (MHC) bind antigenic peptides that are subsequently presented to T cells. Previous studies have shown that most of the residues required for binding of the chicken ovalbumin (Ova) 323-339 peptide to the I-A(d) MHC class II protein are contained within the shorter 325-336 peptide. This observation is somewhat inconsistent with the X-ray structure of the Ova peptide covalently attached to I-A(d) ( structure) in which residues 323 and 324 form binding interactions with the protein. A second register for the Ova(325-336) peptide is proposed where residues 326 and 327 occupy positions similar to residues 323 and 324 in the structure. Two Ova peptides that minimally encompass the and alternate registers, Ova(323-335) and Ova(325-336), respectively, were found to dissociate from I-A(d) with distinct kinetics. The dissociation rates for both peptides were enhanced when the His81 residue of the MHC beta-chain was replaced with an asparagine. In the structure the betaH81 residue forms a hydrogen bond to the backbone carbonyl of I323. If the Ova(325-336) peptide were also bound in the register, there would be no comparable hydrogen-bond acceptor for the betaH81 side chain that could explain this peptide's sensitivity to the betaH81 replacement. The Ova(323-335) peptide that binds in the register does not stimulate a T-cell hybridoma that is stimulated by Ova(325-336) bound in the alternate register. These results demonstrate that a single peptide can bind to an MHC peptide in alternate registers producing distinct T-cell responses.