Biomedical subjects

Neelan J Marianayagam

Publications and source records attributed to Neelan J Marianayagam.

8 recordsLinked to original sources

PubMed2007

Dimerization of CtIP may stabilize in vivo interactions with the Retinoblastoma-pocket domain.

CtIP is a tumor suppressor that interacts with Retinoblastoma protein (Rb) to regulate the G1/S-phase transition of the cell cycle. Despite its large size (897 residues) CtIP has few known structured regions. Rather it contains several linear motifs that interact with known binding partners, including an LXCXE motif that binds the pocket domain of Rb-family proteins. This LXCXE motif lies at the C-terminus of the only known structured domain, an N-terminal coiled-coil dimerization domain (DD; residues 45-160). Yeast two-hybrid (Y2H) and GST-pulldown analyses showed that CtIP requires the LXCXE motif to bind the Rb-pocket. Although isothermal titration calorimetry data indicates that the LXCXE motif is the sole determinant of binding affinity for the Rb-pocket domain (K(A) approximately 10(6)M(-1)), Y2H data indicates that the DD is required to stabilize the interaction in vivo. Thus dimerization may increase the apparent stability of the proteins and/or the lifetime of the complexes.

Amino Acid Motifs↗

PubMed2006

Identification of the key LMO2-binding determinants on Ldb1.

The overexpression of LIM-only protein 2 (LMO2) in T-cells, as a result of chromosomal translocations, retroviral insertion during gene therapy, or in transgenic mice models, leads to the onset of T-cell leukemias. LMO2 comprises two protein-binding LIM domains that allow LMO2 to interact with multiple protein partners, including LIM domain-binding protein 1 (Ldb1, also known as CLIM2 and NLI), an essential cofactor for LMO proteins. Sequestration of Ldb1 by LMO2 in T-cells may prevent it binding other key partners, such as LMO4. Here, we show using protein engineering and enzyme-linked immunosorbent assay (ELISA) methodologies that LMO2 binds Ldb1 with a twofold lower affinity than does LMO4. Thus, excess LMO2 rather than an intrinsically higher binding affinity would lead to sequestration of Ldb1. Both LIM domains of LMO2 are required for high-affinity binding to Ldb1 (K(D) = 2.0 x 10(-8) M). However, the first LIM domain of LMO2 is primarily responsible for binding to Ldb1 (K(D) = 2.3 x 10(-7) M), whereas the second LIM domain increases binding by an order of magnitude. We used mutagenesis in combination with yeast two-hybrid analysis, and phage display selection to identify LMO2-binding "hot spots" within Ldb1 that locate to the LIM1-binding region. The delineation of this region reveals some specific differences when compared to the equivalent LMO4:Ldb1 interaction that hold promise for the development of reagents to specifically bind LMO2 in the treatment of leukemia.

Adaptor Proteins, Signal Transducing↗

PubMed2005

Protein folding by distributed computing and the denatured state ensemble.

The distributed computing (DC) paradigm in conjunction with the folding@home (FH) client server has been used to study the folding kinetics of small peptides and proteins, giving excellent agreement with experimentally measured folding rates, although pathways sampled in these simulations are not always consistent with the folding mechanism. In this study, we use a coarse-grain model of protein L, whose two-state kinetics have been characterized in detail by using long-time equilibrium simulations, to rigorously test a FH protocol using approximately 10,000 short-time, uncoupled folding simulations starting from an extended state of the protein. We show that the FH results give non-Poisson distributions and early folding events that are unphysical, whereas longer folding events experience a correct barrier to folding but are not representative of the equilibrium folding ensemble. Using short-time, uncoupled folding simulations started from an equilibrated denatured state ensemble (DSE), we also do not get agreement with the equilibrium two-state kinetics because of overrepresented folding events arising from higher energy subpopulations in the DSE. The DC approach using uncoupled short trajectories can make contact with traditionally measured experimental rates and folding mechanism when starting from an equilibrated DSE, when the simulation time is long enough to sample the lowest energy states of the unfolded basin and the simulated free-energy surface is correct. However, the DC paradigm, together with faster time-resolved and single-molecule experiments, can also reveal the breakdown in the two-state approximation due to observation of folding events from higher energy subpopulations in the DSE.

Kinetics↗

PubMed2005

Native-state dynamics of the ubiquitin family: implications for function and evolution.

Protein dynamics are integral to protein function. In recent years, the use of computer simulation to understand the molecular motions of proteins has become widespread. However, there are few such studies which compare the dynamics of proteins that are structurally and functionally related. In this study, we present native-state molecular dynamic simulations of four proteins which possess a ubiquitin-like fold. Three of these proteins are thought to have evolved from a common ancestral ubiquitin-like protein and have similarities in their function. A fourth protein, which is structurally homologous but which appears to have a different function, is also studied. Local fluctuations in the native state simulations are analysed, and conserved motions of the C-alpha backbone atoms are identified in residues which are important for function. In addition, the global dynamics of the proteins are analysed using the essential-dynamics method. This analysis reveals a slightly higher degree of conservation in dynamics for the three proteins which are functionally related. Both the global and local analyses illustrate how nature has optimized and conserved protein motions for specific biological activity within the ubiquitin family.

Biological Evolution↗

PubMed2005

An independent method for the analysis of protein folding kinetics from all-atom molecular dynamics simulations.

We propose a method for extracting useful kinetic information from all-atom molecular dynamics simulations of protein folding. By calculating the time correlation functions between the evolution of different structural properties during the course of the simulation we can determine the endpoint of the reaction and the mechanism by which it occurs. As a test of our method we use thermal denaturation simulations on a 76 residue protein, ubiquitin. The method we present should be used in combination with current techniques for analyzing molecular dynamics trajectories.

Computer Simulation↗

PubMed2004

The folding pathway of ubiquitin from all-atom molecular dynamics simulations.

The folding (unfolding) pathway of ubiquitin is probed using all-atom molecular dynamics simulations. We dissect the folding pathway using two techniques: first, we probe the folding pathway of ubiquitin by calculating the evolution of structural properties over time and second, we identify the rate determining transition state for folding. The structural properties that we look at are hydrophobic solvent accessible surface area (SASA) and Calpha-root-mean-square deviation (rmsd). These properties on their own tell us relatively little about the folding pathway of ubiquitin; however, when plotted against each other, they become powerful tools for dissecting ubiquitin's folding mechanism. Plots of Calpha-rmsd against SASA serve as a phase space trajectories for the folding of ubiquitin. In this study, these plots show that ubiquitin folds to the native state via the population of an intermediate state. This is shown by an initial hydrophobic collapse phase followed by a second phase of secondary structure arrangement. Analysis of the structure of the intermediate state shows that it is a collapsed species with very little secondary structure. In reconciling these observations with recent experimental data, the transition that we observe in our simulations from the unfolded state (U) to the intermediate state (I) most likely occurs in the dead-time of the stopped flow instrument. The folding pathway of ubiquitin is probed further by identification of the rate-determining transition state for folding. The method used for this is essential dynamics, which utilizes a principal component analysis (PCA) on the atomic fluctuations throughout the simulation. The five transition state structures identified in silico are in good agreement with the experimentally determined transition state. The calculation of phi-values from the structures generated in the simulations is also carried out and it shows a good correlation with the experimentally measured values. An initial analysis of the denatured state shows that it is compact with fluctuating regions of nonnative secondary structure. It is found that the compactness in the denatured state is due to the burial of some hydrophobic residues. We conclude by looking at a correlation between folding kinetics and residual structure in the denatured state. A hierarchical folding mechanism is then proposed for ubiquitin.

Computer Simulation↗

PubMed2004

The power of two: protein dimerization in biology.

The self-association of proteins to form dimers and higher-order oligomers is a very common phenomenon. Recent structural and biophysical studies show that protein dimerization or oligomerization is a key factor in the regulation of proteins such as enzymes, ion channels, receptors and transcription factors. In addition, self-association can help to minimize genome size, while maintaining the advantages of modular complex formation. Oligomerization, however, can also have deleterious consequences when nonnative oligomers associated with pathogenic states are generated. Specific protein dimerization is integral to biological function, structure and control, and must be under substantial selection pressure to be maintained with such frequency throughout biology.

Allosteric Regulation↗

PubMed2002

Fast folding of a four-helical bundle protein.

The FK506-FKBP12 binding-domain of the kinase FRAP (FRB) forms a classic up-down four-helical bundle. The folding pathway of this protein has been investigated using a combination of equilibrium and kinetic studies. The native state of the protein is stable with respect to the unfolded state by some 7 kcal mol(-1) at pH 6.0, 10 degrees C. A kinetic analysis of unfolding and refolding rate constants as a function of chemical denaturant concentration suggests that an intermediate state may be populated during folding at low concentrations of denaturant. The presence of this intermediate state is confirmed by refolding experiments performed in the presence of the hydrophobic dye 8-anilinonaphthalene-1 sulfonate (ANS). ANS binds to the partially folded intermediate state populated during the folding of FRB and undergoes a large change in fluorescence that can be detected using stopped-flow techniques. Analysis of the kinetic data suggests that the intermediate state is compact and it may even be a misfolded species that has to partially unfold before it can reach the transition state. Folding and unfolding rate constants in water are approximately 150-200 s(-1) and 0.005-0.06 s(-1), respectively, at neutral pH and 10 degrees C. The folding of FRB is somewhat slower than for other all-helical proteins, probably as a consequence of the formation of a metastable intermediate state. The folding rate constant in the absence of any populated intermediate can be estimated to be 8800 s(-1). Despite the presence of an intermediate state, which effectively slows folding, the protein still folds rapidly with a half-life of 5 ms at 10 degrees C. The dependence of the rate constants on denaturant concentration indicates that the transition state for folding is compact with some 80% of the surface area exposed in the unfolded state buried in the transition state. Data presented for FRB is compared with kinetic data obtained for other all-helical proteins.

Anilino Naphthalenesulfonates↗