Radhakrishnan, Ravi
Email Address
ORCID
Disciplines
relationships.isProjectOf
relationships.isOrgUnitOf
Position
Introduction
Research Interests
Collection
12 results
Search Results
Now showing 1 - 10 of 12
Publication Effect of oxidatively damaged DNA on the active site preorganization during nucleotide incorporation in a high fidelity polymerase from Bacillus stearothermophilus(2007-12-07) Venkatramani, Ravindra; Radhakrishnan, RaviWe study the effect of the oxidative lesion 8-oxoguanine (8oxoG) on the preorganization of the active site for DNA replication in the closed (active) state of the Bacillus fragment (BF), a Klenow analog from Bacillus stearothermophilus. Our molecular dynamics and free energy simulations of explicitly solvated model ternary complexes of BF bound to correct dCTP/incorrect dATP opposite guanine (G) and 8oxoG bases in DNA suggest that the lesion introduces structural and energetic changes at the catalytic site to favor dATP insertion. Despite the formation of a stable Watson- Crick pairing in the 8oxoG:dCTP system, the catalytic geometry is severely distorted to possibly slow down catalysis. Indeed, our calculated free energy landscapes associated with active site preorganization suggest additional barriers to assemble an efficient catalytic site, which need to be overcome during dCTP incorporation opposite 8oxoG relative to that opposite undamaged G. In contrast, the catalytic geometry for the Hoogsteen pairing in the 8oxoG:dATP system is highly organized and poised for efficient nucleotide incorporation via the "two-metal-ion" catalyzed phosphoryl transfer mechanism. However, the free energy calculations suggest that the catalytic geometry during dATP incorporation opposite 8oxoG is considerably less plastic than that during dCTP incorporation opposite G despite a very similar, well organized catalytic site for both systems. A correlation analysis of the dynamics trajectories suggests the presence of significant coupling between motions of the polymerase fingers and the primary distance for nucleophilic attack (i.e., between the terminal primer O3′ and the dNTP Pα atoms) during correct dCTP incorporation opposite undamaged G. This coupling is shown to be disrupted during nucleotide incorporation by the polymerase with oxidatively damaged DNA/dNTP substrates. We also suggest that the lesion affects DNA interactions with key polymerase residues, thereby affecting the enzymes ability to discriminate against noncomplementary DNA/dNTP substrates. Taken together, our results provide a unified structural, energetic, and dynamic platform to rationalize experimentally observed relative nucleotide incorporation rates for correct dCTP/incorrect dATP insertion opposite an undamaged/oxidatively damaged template G by BF.Publication Computational Study of the Force Dependence of Phosphoryl Transfer during DNA Synthesis by a High Fidelity Polymerase(2008-02-29) Venkatramani, Ravindra; Radhakrishnan, RaviHigh fidelity polymerases are efficient catalysts of phosphodiester bond formation during DNA replication or repair. We interpret molecular dynamics simulations of a polymerase bound to its substrate DNA and incoming nucleotide using a quasiharmonic model to study the effect of external forces applied to the bound DNA on the kinetics of phosphoryl transfer. The origin of the force dependence is shown to be an intriguing coupling between slow, delocalized polymerase-DNA modes and fast catalytic site motions. Using noncognate DNA substrates we show that the force dependence is context specific.Publication A Multiscale Computational Approach to Dissect Early Events in the Erb Family Receptor Mediated Activation, Differential Signaling, and Relevance to Oncogenic Transformations(2006-12-12) Liu, Yingting; Purvis, Jeremy; Shih, Andrew; Weinstein, Joshua; Agrawal, Neeraj; Radhakrishnan, RaviWe describe a hierarchical multiscale computational approach based on molecular dynamics simulations, free energy-based molecular docking simulations, deterministic network-based kinetic modeling, and hybrid discrete/continuum stochastic dynamics protocols to study the dimermediated receptor activation characteristics of the Erb family receptors, specifically the epidermal growth factor receptor (EGFR). Through these modeling approaches, we are able to extend the prior modeling of EGF-mediated signal transduction by considering specific EGFR tyrosine kinase (EGFRTK) docking interactions mediated by differential binding and phosphorylation of different C-terminal peptide tyrosines on the RTK tail. By modeling signal flows through branching pathways of the EGFRTK resolved on a molecular basis, we are able to transcribe the effects of molecular alterations in the receptor (e.g., mutant forms of the receptor) to differing kinetic behavior and downstream signaling response. Our molecular dynamics simulations show that the drug sensitizing mutation (L834R) of EGFR stabilizes the active conformation to make the system constitutively active. Docking simulations show preferential characteristics (for wildtype vs. mutant receptors) in inhibitor binding as well as preferential enhancement of phosphorylation of particular substrate tyrosines over others. We find that in comparison to the wildtype system, the L834R mutant RTK preferentially binds the inhibitor erlotinib, as well as preferentially phosphorylates the substrate tyrosine Y1068 but not Y1173. We predict that these molecular level changes result in preferential activation of the Akt signaling pathway in comparison to the Erk signaling pathway for cells with normal EGFR expression. For cells with EGFR over expression, the mutant over activates both Erk and Akt pathways, in comparison to wildtype. These results are consistent with qualitative experimental measurements reported in the literature. We discuss these consequences in light of how the network topology and signaling characteristics of altered (mutant) cell lines are shaped differently in relationship to native cell lines.Publication Calculation of free energies in fluid membranes subject to heterogeneous curvature fields(2009-07-30) Agrawal, Neeraj J; Radhakrishnan, RaviWe present a computational methodology for incorporating thermal effects and calculating relative free energies for elastic fluid membranes subject to spatially dependent intrinsic curvature fields using the method of thermodynamic integration. Based on a simple model for the intrinsic curvature imposed only in a localized region of the membrane, we employ thermodynamic integration to calculate the free-energy change as a function of increasing strength of the intrinsic curvature field and a thermodynamic cycle to compute free-energy changes for different sizes of the localized region. By explicitly computing the free-energy changes and by quantifying the loss of entropy accompanied with increasing membrane deformation, we show that the membrane stiffness increases with increasing intrinsic field, thereby, renormalizing the membrane bending rigidity. The second main conclusion of this work is that the entropy of the membrane decreases with increasing size of the localized region subject to the curvature field. Our results help to quantify the free-energy change when a planar membrane deforms under the influence of curvature-inducing proteins at a finite temperature.Publication The effect of oxidatively damaged DNA on the active site pre-organization during nucleotide incorporation in a high fidelity polymerase from Bacillus stearothermophilus(2008-05-15) Venkatramani, Ravindra; Radhakrishnan, RaviWe study the effect of the oxidative lesion 8-oxoguanine (8oxoG) on the pre-organization of the active site for DNA replication in the closed (active) state of the Bacillus Fragment (BF), a Klenow analog from Bacillus stearothermophilus. Our molecular dynamics and free energy simulations of explicitly solvated model ternary complexes of BF bound to correct dCTP/incorrect dATP opposite guanine (G) and 8oxoG bases in DNA suggest that the lesion introduces structural and energetic changes at the catalytic site to favor dATP insertion. Despite the formation of a stable Watson-Crick pairing in the 8oxoG:dCTP system, the catalytic geometry is severely distorted to possibly slow down catalysis. Indeed, our calculated free energy landscapes associated with active site pre-organization suggest additional barriers to assemble an efficient catalytic site, which need to be overcome during dCTP incorporation opposite 8oxoG relative to that opposite undamaged G. In contrast, the catalytic geometry for the Hoogsteen pairing in the 8oxoG:dATP system is highly organized and poised for efficient nucleotide incorporation via the "twometal- ion" catalyzed phosphoryl transfer mechanism. However, the free energy calculations suggest that the catalytic geometry during dATP incorporation opposite 8oxoG is considerably less plastic than that during dCTP incorporation opposite G despite a very similar, well organized catalytic site for both systems. A correlation analysis of the dynamics trajectories suggests the presence of significant coupling between motions of the polymerase fingers and the primary distance for nucleophilic attack (i.e., between the terminal primer O3´ and the dNTP Pα atoms) during correct dCTP incorporation opposite undamaged G. This coupling is shown to be disrupted during nucleotide incorporation by the polymerase with oxidatively damaged DNA/dNTP substrates. We also suggest that the lesion affects DNA interactions with key polymerase residues, thereby affecting the enzymes ability to discriminate against noncomplementary DNA/dNTP substrates. Taken together, our results provide a unified structural, energetic, and dynamic platform to rationalize experimentally observed relative nucleotide incorporation rates for correct dCTP/incorrect dATP insertion opposite an undamaged/oxidatively damaged template G by BF.Publication The ErbB3/HER3 Intracellular Domain is Competent to Bind ATP and Catalyze Autophosphorylation(2010-04-27) Shi, Fumin; Telesco, Shannon; Liu, Yingting; Radhakrishnan, Ravi; Lemmon, MarkErbB3/HER3 is one of four members of the human epidermal growth factor receptor (EGFR/HER) or ErbB receptor tyrosine kinase family. ErbB3 binds neuregulins via its extracellular region and signals primarily by heterodimerizing with ErbB2/HER2/Neu. A recently appreciated role for ErbB3 in resistance of tumor cells to EGFR/ErbB2-targeted therapeutics has made it a focus of attention. However, efforts to inactivate ErbB3 therapeutically in parallel with other ErbB receptors are challenging because its intracellular kinase domain is thought to be an inactive pseudokinase that lacks several key conserved (and catalytically important) residues-including the catalytic base aspartate. We report here that, despite these sequence alterations, ErbB3 retains sufficient kinase activity to robustly trans-autophosphorylate its intracellular region--although it is substantially less active than EGFR and does not phosphorylate exogenous peptides. The ErbB3 kinase domain binds ATP with a K(d) of approximately 1.1 microM. We describe a crystal structure of ErbB3 kinase bound to an ATP analogue, which resembles the inactive EGFR and ErbB4 kinase domains (but with a shortened alphaC-helix). Whereas mutations that destabilize this configuration activate EGFR and ErbB4 (and promote EGFR-dependent lung cancers), a similar mutation conversely inactivates ErbB3. Using quantum mechanics/molecular mechanics simulations, we delineate a reaction pathway for ErbB3-catalyzed phosphoryl transfer that does not require the conserved catalytic base and can be catalyzed by the "inactive-like" configuration observed crystallographically. These findings suggest that ErbB3 kinase activity within receptor dimers may be crucial for signaling and could represent an important therapeutic target.Publication Molecular Systems Biology of ErbB1 Signaling: Bridging the Gap through Multiscale Modeling and High-Performance Computing(2008-12-01) Shih, Andrew; Radhakrishnan, Ravi; Purvis, JeremyThe complexity in intracellular signaling mechanisms relevant for the conquest of many diseases resides at different levels of organization with scales ranging from the subatomic realm relevant to catalytic functions of enzymes to the mesoscopic realm relevant to the cooperative association of molecular assemblies and membrane processes. Consequently, the challenge of representing and quantifying functional or dysfunctional modules within the networks remains due to the current limitations in our understanding of mesoscopic biology, i.e., how the components assemble into functional molecular ensembles. A multiscale approach is necessary to treat a hierarchy of interactions ranging from molecular (nm, ns) to signaling (μm, ms) length and time scales, which necessitates the development and application of specialized modeling tools. Complementary to multiscale experimentation (encompassing structural biology, mechanistic enzymology, cell biology, and single molecule studies) multiscale modeling offers a powerful and quantitative alternative for the study of functional intracellular signaling modules. Here, we describe the application of a multiscale approach to signaling mediated by the ErbB1 receptor which constitutes a network hub for the cell’s proliferative, migratory, and survival programs. Through our multiscale model, we mechanistically describe how point-mutations in the ErbB1 receptor can profoundly alter signaling characteristics leading to the onset of oncogenic transformations. Specifically, we describe how the point mutations induce cascading fragility mechanisms at the molecular scale as well as at the scale of the signaling network to preferentially activate the survival factor Akt. We provide a quantitative explanation for how the hallmark of preferential Akt activation in cell-lines harboring the constitutively active mutant ErbB1 receptors causes these cell-lines to be addicted to ErbB1-mediated generation of survival signals. Consequently, inhibition of ErbB1 activity leads to a remarkable therapeutic response in the addicted cell lines.Publication Geometry of mediating protein affects the probability of loop formation in DNA.(2008-04-15) Agrawal, Neeraj J; Radhakrishnan, Ravi; Purohit, Prashant KRecent single molecule experiments have determined the probability of loop formation in DNA as a function of the DNA contour length for different types of looping proteins. The optimal contour length for loop formation as well as the probability density functions have been found to be strongly dependent on the type of looping protein used. We show, using Monte Carlo simulations and analytical calculations, that these observations can be replicated using the wormlike-chain model for double-stranded DNA if we account for the nonzero size of the looping protein. The simulations have been performed in two dimensions so that bending is the only mode of deformation available to the DNA while the geometry of the looping protein enters through a single variable which is representative of its size. We observe two important effects that seem to directly depend on the size of the enzyme: 1), the overall propensity of loop formation at any given value of the DNA contour length increases with the size of the enzyme; and 2), the contour length corresponding to the first peak as well as the first well in the probability density functions increases with the size of the enzyme. Additionally, the eigenmodes of the fluctuating shape of the looped DNA calculated from simulations and theory are in excellent agreement, and reveal that most of the fluctuations in the DNA occur in regions of low curvature.Publication Coupling of Fast and Slow Modes in the Reaction Pathway of the Minimal Hammerhead Ribozyme Cleavage(2007-10-01) Radhakrishnan, RaviBy employing classical molecular dynamics, correlation analysis of coupling between slow and fast dynamical modes, and free energy (umbrella) sampling using classical as well as mixed quantum mechanics molecular mechanics force fields, we uncover a possible pathway for phosphoryl transfer in the self-cleaving reaction of the minimal hammerhead ribozyme. The significance of this pathway is that it initiates from the minimal hammerhead crystal structure and describes the reaction landscape as a conformational rearrangement followed by a covalent transformation. The delineated mechanism is catalyzed by two metal (Mg2+) ions, proceeds via an in-line-attack by CYT 17 O2′ on the scissile phosphorous (ADE 1.1 P), and is therefore consistent with the experimentally observed inversion configuration. According to the delineated mechanism, the coupling between slow modes involving the hammerhead backbone with fast modes in the cleavage site appears to be crucial for setting up the in-line nucleophilic attack.Publication Steady-State Kinetic Modeling Constrains Cellular Resting States and Dynamic Behavior(2009-03-06) Purvis, Jeremy E; Radhakrishnan, Ravi; Diamond, Scott LA defining characteristic of living cells is the ability to respond dynamically to external stimuli while maintaining homeostasis under resting conditions. Capturing both of these features in a single kinetic model is difficult because the model must be able to reproduce both behaviors using the same set of molecular components. Here, we show how combining small, well-defined steady-state networks provides an efficient means of constructing large-scale kinetic models that exhibit realistic resting and dynamic behaviors. By requiring each kinetic module to be homeostatic (at steady state under resting conditions), the method proceeds by (i) computing steady-state solutions to a system of ordinary differential equations for each module, (ii) applying principal component analysis to each set of solutions to capture the steady-state solution space of each module network, and (iii) combining optimal search directions from all modules to form a global steady-state space that is searched for accurate simulation of the time-dependent behavior of the whole system upon perturbation. Importantly, this stepwise approach retains the nonlinear rate expressions that govern each reaction in the system and enforces constraints on the range of allowable concentration states for the full-scale model. These constraints not only reduce the computational cost of fitting experimental time-series data but can also provide insight into limitations on system concentrations and architecture. To demonstrate application of the method, we show how small kinetic perturbations in a modular model of platelet P2Y(1) signaling can cause widespread compensatory effects on cellular resting states.