Biomedical subjects

Aaron R Dinner

Publications and source records attributed to Aaron R Dinner.

17 recordsLinked to original sources

PubMed2006

Implications of alternative substrate binding modes for catalysis by uracil-DNA glycosylase: an apparent discrepancy resolved.

A theoretical study showed that the base excision repair enzyme uracil-DNA glycosylase (UDG) exploits electrostatic interactions with backbone phosphate groups in the substrate for catalysis. Although experiments performed to test the calculated results confirmed the predicted importance of the -2, -1, and +1 phosphate groups, there was an apparent disagreement with regard to the +2 phosphate group. The calculations indicated that it made an important contribution, while experimentally, the effect of its deletion or neutralization was small. The +2 phosphate group interacts directly with an active site histidine (H148 in humans) in the crystal structure of UDG in complex with double-stranded (ds) DNA. We previously calculated that H148 has a strong anticatalytic effect due to its protonation, and here we use alchemical free energy simulations to estimate its site-specific pKa. We find that it is positively charged over the entire experimental pH range (4-10), so its deprotonation cannot compensate for deletion or neutralization of the +2 phosphate group. The free energy simulations are facilitated by an efficient charge-scaling procedure that allows quantitative correction for the implicit treatment of solvent far from the active site; improvements are made to that method to account carefully for differences in the truncation of electrostatic interactions in the contributing molecular-mechanical and continuum-electrostatic approaches. Additional simulations are used to demonstrate that the +2 phosphate group is fully solvent exposed in complexes with single-stranded DNA substrates like those used in the experiments. In contrast, it is well-structured and buried in the dsDNA complex used in the original simulations. Differences in solvent shielding thus account for the apparent lack of an effect observed experimentally upon neutralization or deletion of this group.

Catalysis↗

PubMed2006

Bias annealing: a method for obtaining transition paths de novo.

Computational studies of dynamics in complex systems require means for generating reactive trajectories with minimum knowledge about the processes of interest. Here, we introduce a method for generating transition paths when an existing one is not already available. Starting from biased paths obtained from steered molecular dynamics, we use a Monte Carlo procedure in the space of whole trajectories to shift gradually to sampling an ensemble of unbiased paths. Application to basin-to-basin hopping in a two-dimensional model system and nucleotide-flipping by a DNA repair protein demonstrates that the method can efficiently yield unbiased reactive trajectories even when the initial steered dynamics differ significantly. The relation of the method to others and the physical basis for its success are discussed.

Journal Article↗

PubMed2006

A model for TCR gene segment use.

The TCR alpha-chain is assembled by somatic recombination of variable (V) and joining (J) gene segments at the CD4+ CD8+ stage of development. In this study, we present the first analytical model for deletional rearrangement and show that it is consistent with almost all available data on Valpha Jalpha use in mice and humans. A key feature of the model is that both "local" and "express service" models of rearrangement can be obtained by varying a single parameter that describes the number of gene segments accessible at a time. We find that the window is much larger for Valpha segments than Jalpha segments, which reconciles seemingly conflicting data for the former. Implications for the properties of the repertoire as a whole and experiments that seek to probe them are discussed. Special considerations for allelic inclusion are treated in the Appendices.

Alleles↗

PubMed2006

Multilineage transcriptional priming and determination of alternate hematopoietic cell fates.

Hematopoietic stem cells and their progenitors exhibit multilineage patterns of gene expression. Molecular mechanisms underlying the generation and refinement of these patterns during cell fate determination remain unexplored because of the absence of suitable experimental systems. Using PU.1(-/-) progenitors, we demonstrate that at subthreshold levels, this Ets transcription factor regulates a mixed pattern (macrophage/neutrophil) of gene expression within individual myeloid progenitors. Increased PU.1 levels refine the pattern and promote macrophage differentiation by modulating a novel regulatory circuit comprised of counter antagonistic repressors, Egr-1,2/Nab-2 and Gfi-1. Egr-1 and Egr-2 function redundantly to activate macrophage genes and to repress the neutrophil program. These results are used to assemble and mathematically model a gene regulatory network that exhibits both graded and bistable behaviors and accounts for the onset and resolution of mixed lineage patterns during cell fate determination.

Animals↗

PubMed2006

A hybrid deterministic-stochastic algorithm for modeling cell signaling dynamics in spatially inhomogeneous environments and under the influence of external fields.

Cell signaling dynamics mediate myriad processes in biology. It has become increasingly clear that inter- and intracellular signaling reactions often occur in a spatially inhomogeneous environment and that it is important to account for stochastic fluctuations of certain species involved in signaling reactions. The importance of these effects enhances the difficulty of gleaning mechanistic information from observations of a few experimental reporters and highlights the significance of synergistic experimental and computational studies. When both stochastic fluctuations and spatial inhomogeneity must be included in a model simultaneously, however, the resulting computational demands quickly become overwhelming. In many situations the failure of standard coarse-graining methods (i.e., ignoring spatial variation or stochastic fluctuations) when applied to all components of a complex system does not exclude the possibility of successfully applying such coarse-graining to some components of the system. Following this approach alleviates computational cost but requires "hybrid" algorithms where some variables are treated at a coarse-grained level while others are not. We present an efficient algorithm for simulation of stochastic, spatially inhomogeneous reaction-diffusion kinetics coupled to coarse-grained fields described by (stochastic or deterministic) partial differential equations (PDEs). The PDEs could represent mean-field descriptions of reactive species present in large copy numbers or evolution of hydrodynamic variables that influence signaling (e.g., membrane shape or cytoskeletal motion). We discuss the approximations made to derive our algorithm and test its efficacy by applying it to problems that include many features typical of realistic cell signaling processes.

Algorithms↗

PubMed2006

Conformational sampling via a self-regulating effective energy surface.

The difficulty of efficiently sampling the phase space of complex systems with rough energy surfaces is well known. Typical solutions to the problem involve accelerating the crossing of barriers, but such methods often have the secondary problem that the low-energy states of interest are inadequately sampled, unless the parameters of the search algorithm are modified as the system evolves. A method is presented to improve the sampling with particular emphasis on the low-energy conformations, which make the most important contributions to the thermodynamics of the system. The algorithm proposed here samples the details of the minima, while easily surmounting barriers. This is achieved by introducing a self-regulating sampling variable which depends on the current state of the system. Two replicas of the system are introduced and the sampling variable is treated as a particle coupled to the physical system. The method is illustrated with a simple model system and is applied to the realistic example of barrier crossing in a protein-ligand complex.

Algorithms↗

PubMed2006

Dynamic coupling between coordinates in a model for biomolecular isomerization.

To understand a complex reaction, it is necessary to project the dynamics of the system onto a low-dimensional subspace of physically meaningful coordinates. We recently introduced an automatic method for identifying coordinates that relate closely to stable-state commitment probabilities and successfully applied it to a model for biomolecular isomerization, the C(7eq)-->alpha(R) transition of the alanine dipeptide [A. Ma and A. R. Dinner, J. Phys. Chem. B 109, 6769 (2005)]. Here, we explore approximate means for estimating diffusion tensors for systems subject to restraints in one and two dimensions and then use the results together with an extension of Kramers theory for unimolecular reaction rates [A. Berezhkovskii and A. Szabo, J. Chem. Phys. 122, 014503 (2005)] to show explicitly that both the potential of mean force and the diffusion tensor are essential for describing the dynamics of the alanine dipeptide quantitatively. In particular, the signficance of off-diagonal elements of the diffusion tensor suggests that the coordinates of interest are coupled by the hydrodynamic-like response of the bath of remaining degrees of freedom.

Acetamides↗

PubMed2006

Monte Carlo simulations of biomolecules: The MC module in CHARMM.

We describe the implementation of a general and flexible Monte Carlo (MC) module for the program CHARMM, which is used widely for modeling biomolecular systems with empirical energy functions. Construction and use of an almost arbitrary move set with only a few commands is made possible by providing several predefined types of moves that can be combined. Sampling can be enhanced by noncanonical acceptance criteria, automatic optimization of step sizes, and energy minimization. A systematic procedure for improving MC move sets is introduced and applied to simulations of two peptides. The resulting move sets allow MC to sample the configuration spaces of these systems much more rapidly than Langevin dynamics. The rate of convergence of the difference in free energy between ethane and methanol in explicit solvent is also examined, and comparable performances are observed for MC and the Nosé-Hoover algorithm. Its ease of use combined with its sampling efficiency make the MC module in CHARMM an attractive alternative for exploring the behavior of biomolecular systems.

Alanine↗

PubMed2006

B cell antigen receptor signaling and internalization are mutually exclusive events.

Engagement of the B cell antigen receptor initiates two concurrent processes, signaling and receptor internalization. While both are required for normal humoral immune responses, the relationship between these two processes is unknown. Herein, we demonstrate that following receptor ligation, a small subpopulation of B cell antigen receptors are inductively phosphorylated and selectively retained at the cell surface where they can serve as scaffolds for the assembly of signaling molecules. In contrast, the larger population of non-phosphorylated receptors is rapidly endocytosed. Each receptor can undergo only one of two mutually exclusive fates because the tyrosine-based motifs that mediate signaling when phosphorylated mediate internalization when not phosphorylated. Mathematical modeling indicates that the observed competition between receptor phosphorylation and internalization enhances signaling responses to low avidity ligands.

Amino Acid Motifs↗

PubMed2005

Enhancement of diffusion-controlled reaction rates by surface-induced orientational restriction.

We explore the means by which immobilization of a substrate on a surface can increase the rate of a diffusion-controlled enzymatic reaction. A quasichemical approach is developed and compared with Brownian dynamics simulations. We use these methods to show that restricting only the orientation of the enzyme by long-range interactions with the surface is sufficient for enhancing catalysis.

Biophysics↗

PubMed2005

Control of genotypic allelic inclusion through TCR surface expression.

To gain insight into the molecular causes and functional consequences of allelic inclusion of TCR alpha-chains, we develop a computational model for thymocyte selection in which the signal that determines cell fate depends on surface expression. Analysis of receptor pairs on selected dual TCR cells reveals that allelic inclusion permits both autoreactive TCR and receptors not in the single TCR cell repertoire to be selected. However, in comparison with earlier theoretical studies, relatively few dual TCR cells display receptors with high avidity for thymic ligands because their alpha-chains compete aggressively for the beta-chain, which hinders rescue from clonal deletion. This feature of the model makes clear that allelic inclusion does not in itself compromise central tolerance. A specific experiment based on modulation of TCR surface expression levels is proposed to test the model.

Alleles↗

PubMed2005

Automatic method for identifying reaction coordinates in complex systems.

To interpret simulations of a complex system to determine the physical mechanism of a dynamical process, it is necessary to identify the small number of coordinates that distinguish the stable states from the transition states. We develop an automatic method for identifying these degrees of freedom from a database of candidate physical variables. In the method neural networks are used to determine the functional dependence of the probability of committing to a stable state (committor) on a set of coordinates, and a genetic algorithm selects the combination of inputs that yields the best fit. The method enables us to obtain the first set of coordinates that is demonstrably sufficient to specify the transition state of the C(7eq)--> alpha(R) isomerization of the alanine dipeptide in the presence of explicit water molecules. It is revealed that the solute-solvent coupling can be described by a solvent-derived electrostatic torque around one of the main-chain bonds, and the collective, long-ranged nature of this interaction accounts for previous failures to characterize this reaction.

Alanine↗

PubMed2004

Grand canonical Monte Carlo simulations of water in protein environments.

The grand canonical simulation algorithm is considered as a general methodology to sample the configuration of water molecules confined within protein environments. First, the probability distribution of the number of water molecules and their configuration in a region of interest for biochemical simulations, such as the active site of a protein, is derived by considering a finite subvolume in open equilibrium with a large system serving as a bulk reservoir. It is shown that the influence of the bulk reservoir can be represented as a many-body potential of mean force acting on the atoms located inside the subvolume. The grand canonical Monte Carlo (GCMC) algorithm, augmented by a number of technical advances to increase the acceptance of insertion attempts, is implemented, and tested for simple systems. In particular, the method is illustrated in the case of a pure water box with periodic boundary conditions. In addition, finite spherical systems of pure water and containing a dialanine peptide, are simulated with GCMC while the influence of the surrounding infinite bulk is incorporated using the generalized solvent boundary potential [W. Im, S. Berneche, and B. Roux, J. Chem. Phys. 114, 2924 (2001)]. As a last illustration of water confined in the interior of a protein, the hydration of the central cavity of the KcsA potassium channel is simulated.

Algorithms↗

PubMed2004

CD4 enhances T cell sensitivity to antigen by coordinating Lck accumulation at the immunological synapse.

How T cells respond with extraordinary sensitivity to minute amounts of agonist peptide and major histocompatibility complex (pMHC) molecules on the surface of antigen-presenting cells bearing large numbers of endogenous pMHC molecules is not understood. Here we present evidence that CD4 affects the responsiveness of T helper cells by controlling spatial localization of the tyrosine kinase Lck in the synapse. This finding, as well as further in silico and in vitro experiments, led us to develop a molecular model in which endogenous and agonist pMHC molecules act cooperatively to amplify T cell receptor signaling. At the same time, activation due to endogenous pMHC molecules alone is inhibited. A key feature is that the binding of agonist pMHC molecules to the T cell receptor results in CD4-mediated spatial localization of Lck, which in turn enables endogenous pMHC molecules to trigger many T cell receptors. We also discuss broader implications for T cell biology, including thymic selection, diversity of the repertoire of self pMHC molecules and serial triggering.

Animals↗

PubMed2003

Atomistic understanding of kinetic pathways for single base-pair binding and unbinding in DNA.

We combine free-energy calculations and molecular dynamics to elucidate a mechanism for DNA base-pair binding and unbinding in atomic detail. Specifically, transition-path sampling is used to overcome computational limitations associated with conventional techniques to harvest many trajectories for the flipping of a terminal cytosine in a 3-bp oligomer in explicit water. Comparison with free-energy projections obtained with umbrella sampling reveals four coordinates that separate true dynamic transition states from stable reactant and product states. Unbinding proceeds via two qualitatively different pathways: one in which the flipping base breaks its intramolecular hydrogen bonds before it unstacks and another in which it ruptures both sets of interactions simultaneously. Both on- and off-pathway intermediates are observed. The relation of the results to coarse-grained models for DNA-based biosensors is discussed.

Base Pairing↗

PubMed2003

The immunological synapse balances T cell receptor signaling and degradation.

The immunological synapse is a specialized cell-cell junction between T cell and antigen-presenting cell surfaces. It is characterized by a central cluster of antigen receptors, a ring of integrin family adhesion molecules, and temporal stability over hours. The role of this specific organization in signaling for T cell activation has been controversial. We use in vitro and in silico experiments to determine that the immunological synapse acts as a type of adaptive controller that both boosts T cell receptor triggering and attenuates strong signals.

Adaptor Proteins, Signal Transducing↗

PubMed2003

Elastic energy storage in beta-sheets with application to F1-ATPase.

We present a methodology for obtaining the elastic properties of protein motifs. We combine the use of interpolated structures (IS), molecular dynamics (MD) and collective coordinates to deduce the elastic properties of the beta-sheet in F(1) ATPase. We find that about 3.5 kcal/mol (6 k(B) T at room temperature) of elastic energy is stored in the beta-sheet as the beta-subunit undergoes its hinge bending motion, in good agreement with the finite element model of Wang and Oster [Nature (1998) 396:279-282]. The technique should be useful for beta-sheets in other proteins and aid in the construction of phenomenological models for molecular motors that are computationally prohibitive for MD alone.

Kinetics↗