The Origin of Nonmonotonic Complex Behavior and the Effects of Nonnative Interactions on the Diffusive Properties of Protein Folding

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Topography of funneled landscapes determines the thermodynamics and kinetics of protein folding

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

The water factor in the protein-folding problem

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Configuration-dependent diffusion can shift the kinetic transition state and barrier height of protein folding

We show that diffusion can play an important role in protein-folding kinetics. We explicitly calculate the diffusion coefficient of protein folding in a lattice model. We found that diffusion typically is configuration- or reaction coordinate-dependent. The diffusion coefficient is found to be decreasing with respect to the progression of folding toward the native state, which is caused by the collapse to a compact state constraining the configurational space for exploration. The configuration- or position-dependent diffusion coefficient has a significant contribution to the kinetics in addition to the thermodynamic free-energy barrier. It effectively changes (increases in this case) the kinetic barrier height as well as the position of the corresponding transition state and therefore modifies the folding kinetic rates as well as the kinetic routes. The resulting folding time, by considering both kinetic diffusion and the thermodynamic folding free-energy profile, thus is slower than the estimation from the thermodynamic free-energy barrier with constant diffusion but is consistent with the results from kinetic simulations. The configuration- or coordinate-dependent diffusion is especially important with respect to fast folding, when there is a small or no free-energy barrier and kinetics is controlled by diffusion. Including the configurational dependence will challenge the transition state theory of protein folding. The classical transition state theory will have to be modified to be consistent. The more detailed folding mechanistic studies involving phi value analysis based on the classical transition state theory also will have to be modified quantitatively.

Electrochemical genosensors for the detection of Bonamia parasite. Selection of single strand-DNA (ssDNA) probes by simulation of the secondary structure folding

A post-PCR nucleic acid work by comparing experimental data, from electrochemical genosensors, and bioinformatics data, derived from the simulation of the secondary structure folding and prediction of hybridisation reaction, was carried out in order to rationalize the selection of ssDNA probes for the detection of two Bonamia species, B. exitiosa and B. ostreae, parasites of Ostrea edulis.Six ssDNA probes (from 11 to 25 bases in length, 2 thiolated and 4 biotinylated) were selected within different regions of B. ostreae and B. exitiosa PCR amplicons (300 and 304 bases, respectively) with the aim to discriminate between these parasite species. ssDNA amplicons and probes were analyzed separately using the "Mfold Web Server" simulating the secondary structure folding behaviour. The hybridisation of amplicon-probe was predicted by means of "Dinamelt Web Server". The results were evaluated considering the number of hydrogen bonds broken and formed in the simulated folding and hybridisation process, variance in gaps for each sequence and number of available bases. In the experimental part, thermally denatured PCR products were captured at the sensor interface via sandwich hybridisation with surface-tethered probes (thiolated probes) and biotinylated signalling probes. A convergence between analytical signals and simulated results was observed, indicating the possibility to use bioinformatic data for ssDNA probes selection to be incorporated in genosensors. (C) 2011 Elsevier B.V. All rights reserved.

Frustration and hydrophobicity interplay in protein folding and protein evolution

A lattice model is used to study mutations and compacting effects on protein folding rates and folding temperature. In the context of protein evolution, we address the question regarding the best scenario for a polypeptide chain to fold: either a fast nonspecific collapse followed by a slow rearrangement to form the native structure or a specific collapse from the unfolded state with the simultaneous formation of the native state. This question is investigated for optimized sequences, whose native state has no frustrated contacts between monomers, and also for mutated sequences, whose native state has some degree of frustration. It is found that the best scenario for folding may depend on the amount of frustration of the native structure. The implication of this result on protein evolution is discussed. (c) 2006 American Institute of Physics.

Probing the kinetics of single molecule protein folding

We propose an approach to integrate the theory, simulations, and experiments in protein-folding kinetics. This is realized by measuring the mean and high-order moments of the first-passage time and its associated distribution. The full kinetics is revealed in the current theoretical framework through these measurements. In the experiments, information about the statistical properties of first-passage times can be obtained from the kinetic folding trajectories of single molecule experiments ( for example, fluorescence). Theoretical/simulation and experimental approaches can be directly related. We study in particular the temperature-varying kinetics to probe the underlying structure of the folding energy landscape. At high temperatures, exponential kinetics is observed; there are multiple parallel kinetic paths leading to the native state. At intermediate temperatures, nonexponential kinetics appears, revealing the nature of the distribution of local traps on the landscape and, as a result, discrete kinetic paths emerge. At very low temperatures, exponential kinetics is again observed; the dynamics on the underlying landscape is dominated by a single barrier. The ratio between first-passage-time moments is proposed to be a good variable to quantitatively probe these kinetic changes. The temperature-dependent kinetics is consistent with the strange kinetics found in folding dynamics experiments. The potential applications of the current results to single-molecule protein folding are discussed.

Crystallization and preliminary X-ray diffraction analysis of suramin, a highly charged polysulfonated napthylurea, complexed with a myotoxic PLA(2) from Bothrops asper venom

Suramin is a highly charged polysulfonated napthylurea that interferes in a number of physiologically relevant processes such as myotoxicity, blood coagulation and several kinds of cancers. This synthetic compound was complexed with a myotoxic Lys49 PLA(2) from Bothrops asper venom and crystallized by the hanging-drop vapor diffusion method at 18 degreesC. The crystals belong to the orthorhombic space group P2(1)2(1)2(1), with unit cell parameters a=49.05, b=63.84 and c=85.67 Angstrom, Diffraction data was collected to 1.78 Angstrom. (C) 2004 Elsevier B.V. All rights reserved.

Assignment of the disulfide bridges in bothropstoxin-I, a myonecrotic Lys49 PLA(2) homolog from Bothrops jararacussu snake venom

Bothropstoxin-I (BthTX-1), a Lys49 phospholipase A(2) homolog with no apparent catalytic activity, was first isolated from Bothrops jararacussu snake venom and completely sequenced in this laboratory. It is a 121-amino-acid single polypeptide chain, highly myonecrotic, despite its inability to catalyze hydrolysis of egg yolk phospholipids, and has 14 half-cystine residues identified at positions 27, 29, 44, 45, 50, 51, 61, 84, 91, 96, 98, 105, 123, and 131 (numbering according to the conventional alignment including gaps, so that the last residue is Cys 131). In order to access its seven disulfide bridges, two strategies were followed: (1) Sequencing of isolated peptides from (tryptic + SV8) and chymotryptic digests by Edman-dansyl degradation; (2) crystallization of the protein and determination of the crystal structure so that at least two additional disulfide bridges could be identified in the final electron density map. Identification of the disulfide-containing peptides from the enzymatic digests was achieved following the disappearance of the original peptides from the HPLC profile after reduction and carboxymethylation of the digest. Following this procedure, four bridges were initially identified from the tryptic and SV8 digests: Cys50-Cys131, Cys51-Cys98, Cys61-Cys91, and Cys84-Cys96. From the chymotryptic digest other peptides were isolated either containing some of the above bridges, therefore confirming the results from the tryptic digest, or presenting a new bond between Cys27 and Cys123. The two remaining bridges were identified as Cys29-Cys45 and Cys44-Cys105 by determination of the crystal structure, showing that BthTX-1 disulfide bonds follow the normal pattern of group II PLA(2)s.

Crytallization and high-resolution X-ray diffraction data collection of an Asp49 PLA(2) from Bothrops jararacussu venom both in the presence and absence of Ca2+ ions

Snake venom PLA(2)s have been extensively studied due to their role in mediating and disrupting physiological processes such as coagulation, platelet aggregation and myotoxicity. The Ca2+ ion bound to the putative calcium-binding loop is essential for hydrolytic activity. We report the crystallization in the presence and absence of Ca2+ and X-ray diffraction data collection at 1.60 Angstrom (with Ca2+) and 1.36 Angstrom (without Ca2+) of an Asp49 PLA(2) from Bothrops jararacussu venom. The crystals belong to orthorhombic space group C222(1). Initial refinement and electron density analysis indicate significant conformational. changes upon Ca2+ binding. (C) 2004 Elsevier B.V. All fights reserved.

Pauli nonlocality in heavy-ion rainbow scattering: A further test of the folding model

Nonlocal interactions are an intrinsically quantum phenomenon. In this work we point out that, in the context of heavy ions, such interactions can be studied through the refractive elastic scattering of these systems at intermediate energies. We show that most of the observed energy dependence of the local equivalent bare potential arises from the exchange nonlocality. The nonlocality parameter extracted from the data was found to be very close to the one obtained from folding models. The effective mass of the colliding, heavy-ion, system was found to be close to the nucleon effective mass in nuclear matter.