PubMed Health⌕ Search

Biomedical subjects

D Thirumalai

Publications and source records attributed to D Thirumalai.

At least 19 recordsLinked to original sources
PubMed2001

Role of counterion condensation in folding of the Tetrahymena ribozyme. I. Equilibrium stabilization by cations.

Folding of RNA into an ordered, compact structure requires substantial neutralization of the negatively charged backbone by positively charged counterions. Using a native gel electrophoresis assay, we have examined the effects of counterion condensation upon the equilibrium folding of the Tetrahymena ribozyme. Incubation of the ribozyme in the presence of mono-, di- and trivalent ions induces a conformational state that is capable of rapidly forming the native structure upon brief exposure to Mg2+. The cation concentration dependence of this transition is directly correlated with the charge of the counterion used to induce folding. Substrate cleavage assays confirm the rapid onset of catalytic activity under these conditions. These results are discussed in terms of classical counterion condensation theory. A model for folding is proposed which predicts effects of charge, ionic radius and temperature on counterion-induced RNA folding transitions.

Animals↗
PubMed2000

Native topology determines force-induced unfolding pathways in globular proteins.

Single-molecule manipulation techniques reveal that stretching unravels individually folded domains in the muscle protein titin and the extracellular matrix protein tenascin. These elastic proteins contain tandem repeats of folded domains with beta-sandwich architecture. Herein, we propose by stretching two model sequences (S1 and S2) with four-stranded beta-barrel topology that unfolding forces and pathways in folded domains can be predicted by using only the structure of the native state. Thermal refolding of S1 and S2 in the absence of force proceeds in an all-or-none fashion. In contrast, phase diagrams in the force-temperature (f,T) plane and steered Langevin dynamics studies of these sequences, which differ in the native registry of the strands, show that S1 unfolds in an allor-none fashion, whereas unfolding of S2 occurs via an obligatory intermediate. Force-induced unfolding is determined by the native topology. After proving that the simulation results for S1 and S2 can be calculated by using native topology alone, we predict the order of unfolding events in Ig domain (Ig27) and two fibronectin III type domains ((9)FnIII and (10)FnIII). The calculated unfolding pathways for these proteins, the location of the transition states, and the pulling speed dependence of the unfolding forces reflect the differences in the way the strands are arranged in the native states. We also predict the mechanisms of force-induced unfolding of the coiled-coil spectrin (a three-helix bundle protein) for all 20 structures deposited in the Protein Data Bank. Our approach suggests a natural way to measure the phase diagram in the (f,C) plane, where C is the concentration of denaturants.

Animals↗
PubMed2000

Mechanisms and kinetics of beta-hairpin formation.

Thermodynamics and kinetics of off-lattice models with side chains for the beta-hairpin fragment of immunoglobulin-binding protein and its variants are reported. For all properties (except refolding time tau(F)) there are no qualitative differences between the full model and the Go version. The validity of the models is established by comparison of the calculated native structure with the Protein Data Bank coordinates and by reproducing the experimental results for the degree of cooperativity and tau(F). For the full model tau(F) approximately 2 micros at the folding temperature (experimental value is 6 micros); the Go model folds 50 times faster. Upon refolding, structural changes take place over three time scales. On the collapse time scale compact structures with intact hydrophobic cluster form. Subsequently, hydrogen bonds form, predominantly originating from the turn by a kinetic zipping mechanism. The assembly of the hairpin is complete when most of the interstrand contacts (the rate-limiting step) is formed. The dominant transition state structure (located by using cluster analysis) is compact and structured. We predict that when hydrophobic cluster is moved to the loop tau(F) marginally increases, whereas moving the hydrophobic cluster closer to the termini results in significant decrease in tau(F) relative to wild type. The mechanism of hairpin formation is predicted to depend on turn stiffness.

Kinetics↗
PubMed2000

Maximizing RNA folding rates: a balancing act.

Large ribozymes typically require very long times to refold into their active conformation in vitro, because the RNA is easily trapped in metastable misfolded structures. Theoretical models show that the probability of misfolding is reduced when local and long-range interactions in the RNA are balanced. Using the folding kinetics of the Tetrahymena ribozyme as an example, we propose that folding rates are maximized when the free energies of forming independent domains are similar to each other. A prediction is that the folding pathway of the ribozyme can be reversed by inverting the relative stability of the tertiary domains. This result suggests strategies for optimizing ribozyme sequences for therapeutics and structural studies.

Animals↗
PubMed1999

Magnesium-dependent folding of self-splicing RNA: exploring the link between cooperativity, thermodynamics, and kinetics.

Folding of the Tetrahymena self-splicing RNA into its active conformation involves a set of discrete intermediate states. The Mg2+-dependent equilibrium transition from the intermediates to the native structure is more cooperative than the formation of the intermediates from the unfolded states. We show that the degree of cooperativity is linked to the free energy of each transition and that the rate of the slow transition from the intermediates to the native state decreases exponentially with increasing Mg2+ concentration. Monovalent salts, which stabilize the folded RNA nonspecifically, induce states that fold in less than 30 s after Mg2+ is added to the RNA. A simple model is proposed that predicts the folding kinetics from the Mg2+-dependent change in the relative stabilities of the intermediate and native states.

Ammonium Sulfate↗
PubMed1999

Stretching single-domain proteins: phase diagram and kinetics of force-induced unfolding.

Single-molecule force spectroscopy reveals unfolding of domains in titin on stretching. We provide a theoretical framework for these experiments by computing the phase diagrams for force-induced unfolding of single-domain proteins using lattice models. The results show that two-state folders (at zero force) unravel cooperatively, whereas stretching of non-two-state folders occurs through intermediates. The stretching rates of individual molecules show great variations reflecting the heterogeneity of force-induced unfolding pathways. The approach to the stretched state occurs in a stepwise "quantized" manner. Unfolding dynamics and forces required to stretch proteins depend sensitively on topology. The unfolding rates increase exponentially with force f till an optimum value, which is determined by the barrier to unfolding when f = 0. A mapping of these results to proteins shows qualitative agreement with force-induced unfolding of Ig-like domains in titin. We show that single-molecule force spectroscopy can be used to map the folding free energy landscape of proteins in the absence of denaturants.

Amino Acid Sequence↗
PubMed1999

Exploring the kinetic requirements for enhancement of protein folding rates in the GroEL cavity.

The chaperonin system, GroEL and GroES of Escherichia coli enable certain proteins to fold under conditions when spontaneous folding is prohibitively slow as to compete with other non-productive channels such as aggregation. We investigated the plausible mechanisms of GroEL-mediated folding using simple lattice models. In particular, we have investigated protein folding in a confined environment, such as those offered by the GroEL, to decipher whether rate and yield enhancement can occur when the substrate protein is allowed to fold within the cavity of the chaperonins. The GroEL cavity is modeled as a cubic box and a simple bead model is used to represent the substrate chain. We consider three distinct characteristic of the confining environment. First, the cavity is taken to be a passive Anfinsen cage in which the walls merely reduce the available conformation space. We find that at temperatures when the native conformation is stable, the folding rate is retarded in the Anfinsen cage. We then assumed that the interior of the wall is hydrophobic. In this case the folding times exhibit a complex behavior. When the strength of the interaction between the polypeptide chain and the cavity is too strong or too weak we find that the rates of folding are retarded compared to spontaneous folding. There is an optimum range of the interaction strength that enhances the rates. Thus, above this value there is an inverse correlation between the folding rates and the strength of the substrate-cavity interactions. The optimal hydrophobic walls essentially pull the kinetically trapped states which leads to a smoother the energy landscape. It is known that upon addition of ATP and GroES the interior cavity of GroEL offers a hydrophilic-like environment to the substrate protein. In order to mimic this within the context of the dynamic Anfinsen cage model, we allow for changes in the hydrophobicity of the walls of the cavity. The duration for which the walls remain hydrophobic during one cycle of ATP hydrolysis is allowed to vary. These calculations show that frequent cycling of the wall hydrophobicity can dramatically reduce the folding times and increase the yield as well under non-permissive conditions. Examination of the structures of the substrate proteins before and after the change in hydrophobicity indicates that there is global unfolding involved. In addition, it is found that a fraction of the molecules kinetically partition to the native state in accordabce with the iterative annealing mechanism. Thus, frequent "unfoldase" activity of chaperonins leading to global unfolding of the polypeptide chain results in enhancement of the folding rates and yield of the folded protein. We suggest that chaperonin efficiency can be greatly enhanced if the cycling time is reduced. The calculations are used to interpret a few experiments on chaperonin-mediated protein folding.

Amino Acid Sequence↗
PubMed1999

Deciphering the timescales and mechanisms of protein folding using minimal off-lattice models.

Considerable insights into the mechanisms and timescales of protein folding have been obtained from detailed studies of minimal off-lattice models. These models are coarse-grained representations of polypeptide chains. Many novel predictions of the mechanisms and timescales of the folding of proteins have been made using simulations of off-lattice models. The concepts derived from these simulations have been used to analyze the recent experiments and simulations of proteins and peptides.

Kinetics↗
PubMed1999

Pair potentials for protein folding: choice of reference states and sensitivity of predicted native states to variations in the interaction schemes.

We examine the similarities and differences between two widely used knowledge-based potentials, which are expressed as contact matrices (consisting of 210 elements) that gives a scale for interaction energies between the naturally occurring amino acid residues. These are the Miyazawa-Jernigan contact interaction matrix M and the potential matrix S derived by Skolnick J et al., 1997, Protein Sci 6:676-688. Although the correlation between the two matrices is good, there is a relatively large dispersion between the elements. We show that when Thr is chosen as a reference solvent within the Miyazawa and Jernigan scheme, the dispersion between the M and S matrices is reduced. The resulting interaction matrix B gives hydrophobicities that are in very good agreement with experiment. The small dispersion between the S and B matrices, which arises due to differing reference states, is shown to have dramatic effect on the predicted native states of lattice models of proteins. These findings and other arguments are used to suggest that for reliable predictions of protein structures, pairwise additive potentials are not sufficient. We also establish that optimized protein sequences can tolerate relatively large random errors in the pair potentials. We conjecture that three body interaction may be needed to predict the folds of proteins in a reliable manner.

Models, Molecular↗
PubMed1998

Lattice models for proteins reveal multiple folding nuclei for nucleation-collapse mechanism.

The nature of the nucleation-collapse mechanism in protein folding is probed using 27-mer and 36-mer lattice models. Three different forms for the interaction potentials are used. Three of the four 27-mer sequences have maximally compact and identical native state while the other has a non-compact native conformation. All the sequences fold thermodynamically and kinetically by a two-state process. Analysis of individual trajectories for each sequence using a self-organizing neural net algorithm shows that upon formation of a critical set of contacts the polypeptide chain rapidly reaches the native conformation which is consistent with a nucleation-collapse mechanism. The algorithm, which reduces the identification of the folding nucleus for each trajectory to one of pattern recognition, is used to show that there are multiple folding nuclei. There is a distribution of nucleation contacts in the transition states with some of them occurring with more probability (when averaged over the denatured ensemble) than others. We also show that there is a distribution in the size of the nuclei with the average number of residues in the folding nuclei being less than about one-third of the chain size. The fluctuations in the sizes of the nuclei are large, suggestive of a broad transition region. The folding nuclei, the structures of each are the corresponding transition states, have varying degree of overlap with the native conformation. The distribution of the radius of gyration of the transition states shows that these structures are an expanded form (by about 25% in the radius of gyration) of the native conformation. Local contacts are most dominant in the folding nuclei while a certain fraction of non-local contacts is necessary to stabilize the transition states. The search for the critical nuclei initially involves the formation of local contacts, while non-local contacts are formed later. The fractional values of PhiF for the two 27-mer mutants found by using the protein engineering protocol are consistent with the microscopic picture of partial formation of structures involving these residues in the transition state. These observations lead to a multiple folding nuclei (MFN) model for nucleation-collapse mechanism in protein folding. The major implication of the MFN model is that, even if the residues whose tertiary interactions are formed nearly completely in the transition state are mutated, it does not disrupt the nature of the nucleation-collapse mechanism. We analyze the experiments on chymotrypsin inhibitor 2 and alpha-spectrin SH3 domain and two circular permutants in light of the MFN model. It is shown that the PhiF-value analysis for these proteins gives considerable support to the MFN model. The theoretical and experimental studies give a coherent picture of the nucleation-collapse mechanism in which there is a distribution of folding nuclei with some more probable than others. The formation of any specific nucleus is not necessary for efficient two-state folding.

Alanine↗
PubMed1998

Hydration for a series of hydrocarbons.

The hydrophobic hydration in a series of hydrocarbons is probed by using molecular dynamics simulations. The solutes considered range from methane to octane. Examination of the shapes of the hydration shell suggests that there is no single stable structure surrounding these solutes. The structure of the water molecules around the solute is not significantly perturbed, even for octane, and the hydrogen bond network is essentially preserved. The solutes are accommodated in the voids of the tetrahedral network of water in such a way as to leave the local environment almost intact. The hydrophobic hydration arises primarily because of the plasticity of the hydrogen bond network. Even for octane we find very little evidence for water-mediated interactions between nonbonded carbon atoms, leading us to suggest that the transition to globular conformations can only occur for very long, linear hydrocarbon chains.

Hydrocarbons↗
PubMed1998

Cooperativity in protein folding: from lattice models with sidechains to real proteins.

BACKGROUND: Over the past few years novel folding mechanisms of globular proteins have been proposed using minimal lattice and off-lattice models. The factors determining the cooperativity of folding in these models and especially their explicit relation to experiments have not been fully established, however. RESULTS: We consider equilibrium folding transitions in lattice models with and without sidechains. A dimensionless measure, omega c, is introduced to quantitatively assess the degree of cooperativity in lattice models and in real proteins. We show that larger values of omega c resembling the values seen in proteins are obtained in lattice models with sidechains. The enhanced cooperativity of such models results from possible denser packing of sidechains in the interior of the model polypeptide chain. We also establish that omega c correlates extremely well with sigma T = (T o - T f) /T o, where T o and T f are collapse and folding transition temperatures, respectively. These theoretical ideas are used to analyze folding transitions in two-state folders (RNase A, chymotrypsin inhibitor 2, fibronectin type III modules and tendamistat) and three-state folders (apomyoglobin and lysozyme). The values of omega c extracted from experiments show a correlation with sigma T (suitably generalized when folding is induced by denaturants or acid). CONCLUSIONS: A quantitative description of the cooperative transition of real proteins can be made by lattice models with sidechains. The degree of cooperativity in minimal models and real proteins can be expressed in terms of the single parameter sigma, which can be estimated from experimental data.

Fibronectins↗
PubMed1998

Virtual atom representation of hydrogen bonds in minimal off-lattice models of alpha helices: effect on stability, cooperativity and kinetics.

BACKGROUND: The most conspicuous feature of a right-handed alpha helix is the presence of hydrogen bonds between the backbone carbonyl oxygen and NH groups along the chain. A simple off-lattice model that includes hydrogen bond interactions using virtual atoms is used to examine the stability, cooperativity and kinetics of the helix-coil transition. RESULTS: We have studied the thermodynamics (using multiple histogram method) and kinetics (by Brownian dynamics simulations) of 16-mer minimal off-lattice models of four-turn alpha-helix sequences. The carbonyl and NH groups are represented as virtual moieties located between two alpha-carbon atoms along the polypeptide chain. The characteristics of the native conformations of the model helices, such as the helical pitch and angular correlations, coincide with those found in real proteins. The transition from coil to helix is quite broad, which is typical of these finite-sized systems. The cooperativity, as measured by a dimensionless parameter, omegac, that takes into account the width and the slope of the transition curves, is enhanced when hydrogen bonds are taken into account. The value of omegac for our model is consistent with that inferred from experiment for an alanine-based helix-forming peptide. The folding time tauF ranges from 6 to 1000 ns in the temperature range 0.7-1.9 T(F), where T(F) is the helix-coil transition temperature. These values are in excellent agreement with the results from recent fast folding experiments. The temperature dependence of tauF exhibits a nearly Arrhenius behavior. Thermally induced unfolding occurs on a time scale that is less than 40-170 ps depending on the final temperature. Our calculations also predict that, although tauF can be altered by changes in the sequence, the dynamic range over which such changes take place is not as large as that predicted for beta-turn formation. CONCLUSIONS: Hydrogen bonds not only affect the stability of alpha-helix formation but also have profound influence on the kinetics. The excellent agreement between our calculations and experiments suggests that these models can be used to investigate the effects of sequence, temperature and viscosity on the helix-coil transition.

Animals↗
PubMed1997

Folding of RNA involves parallel pathways.

Folding kinetics of large RNAs are just beginning to be investigated. We show that the Tetrahymena self-splicing RNA partitions into a population that rapidly reaches the native state, and a slowly folding population that is trapped in metastable misfolded structures. Transitions from the misfolded structures to the native state involve partial unfolding. The total yield of native RNA is increased by iterative annealing of the inactive population, and mildly denaturing conditions increase the rate of folding at physiological temperatures. These results provide the first evidence that an RNA can fold by multiple parallel paths.

Animals↗
PubMed1997

Kinetics of peptide folding: computer simulations of SYPFDV and peptide variants in water.

The folding of Ser-Tyr-Pro-Phe-Asp-Val (SYPFDV), and sequence variants of this peptide (SYPYD and SYPFD) are studied computationally in an explicit water environment. An atomically detailed model of the peptide is embedded in a sphere of TIP3P water molecules and its optimal structure is computed by simulated annealing. At distances from the peptide that are beyond a few solvation shells, a continuum solvent model is employed. The simulations are performed using a mean field approach that enhances the efficiency of sampling peptide conformations. The computations predict a small number of conformations as plausible folded structures. All have a type VI turn conformation for the peptide backbone, similar to that found using NMR. However, some of the structures differ from the experimentally proposed ones in the packing of the proline ring with the aromatic residues. The second most populated structure has, in addition to a correctly folded backbone, the same hydrophobic packing as the conformation measured by NMR. Our simulations suggest a kinetic mechanism that consists of three separate stages. The time-scales associated with these stages are distinct and depend differently on temperature. Electrostatic interactions play an initial role in guiding the peptide chain to a roughly correct structure as measured by the end-to-end distance. At the same time or later the backbone torsions rearrange due to local tendency of the proline ring to form a turn: this step depends on solvation forces and is helped by loose hydrophobic interactions. In the final step, hydrophobic residues pack against each other. We also show the existence of an off the pathway intermediate, suggesting that even in the folding of a small peptide "misfolded" structures can form. The simulations clearly show that parallel folding paths are involved. Our findings suggest that the process of peptide folding shares many of the features expected for the significantly larger protein molecules.

Computer Simulation↗
PubMed1997

Protein folding kinetics: timescales, pathways and energy landscapes in terms of sequence-dependent properties.

BACKGROUND: Recent experimental and theoretical studies have revealed that protein folding kinetics can be quite complex and diverse depending on various factors such as size of the protein sequence and external conditions. For example, some proteins fold apparently in a kinetically two-state manner, whereas others follow complex routes to the native state. We have set out to provide a theoretical basis for understanding the diverse behavior seen in the refolding kinetics of proteins in terms of properties that are intrinsic to the sequence. RESULTS: The folding kinetics of a number of sequences for off-lattice continuum models of proteins is studied using Langevin simulations at two different values of the friction coefficient. We show for these models that there is a remarkable correlation between folding time, tau F, and sigma = (T theta - TF)/T theta, where T theta and TF are the equilibrium collapse and folding transition temperatures, respectively. The microscopic dynamics reveals that several scenarios for the kinetics of refolding arise depending on the range of values of sigma. For relatively small sigma, the chain reaches the native conformation by a direct native conformation nucleation collapse (NCNC) mechanism without being trapped in any detectable intermediates. For moderate and large values of sigma, the kinetics is described by the kinetic partitioning mechanism, according to which a fraction of molecules phi (kinetic partition factor) reach the native conformation via the NCNC mechanism. The remaining fraction attains the native state by off-pathway processes that involve trapping in several misfolded structures. The rate-determining step in the off-pathway processes is the transition from the misfolded structures to the native state. The partition factor phi is also determined by sigma: the smaller the value of sigma, the larger is phi. The qualitative aspects of our results are found to be independent of the friction coefficient. The simulation results and theoretical arguments are used to obtain estimates for timescales for folding via the NCNC mechanism in small proteins, those with less than about 70 amino acid residues. CONCLUSIONS: We have shown that the various scenarios for folding of proteins, and possibly other biomolecules, can be classified solely in terms of sigma. Proteins with small values of sigma reach the native conformation via a nucleation collapse mechanism and their energy landscape is characterized by having one dominant native basin of attraction (NBA). On the other hand, proteins with large sigma get trapped in competing basins of attraction (CBAs) in which they adopt misfolded structures. Only a small fraction of molecules access the native state rapidly when sigma is large. For these sequences, the majority of the molecules approach the native state by a three-stage multipathway mechanism in which the rate-determining step involves a transition from one of the CBAs to the NBA.

Kinetics↗