RESUMO
Old Yellow Enzymes (OYEs) are flavin-dependent redox enzymes that promote the asymmetric reduction of activated alkenes. Due to the high importance of flavoenzymes in the metabolism of organisms, the interaction between OYEs from the parasites Trypanosoma cruzi and Leishmania braziliensis and three diterpene icetexanes (brussonol and two analogs), were evaluated in the present study, and differences in the binding mechanism and inhibition capacity of these molecules were examined. Although the aforementioned compounds showed poor and negligible activities against T. cruzi and L. braziliensis cells, respectively, the experiments with the purified enzymes indicated that the interaction occurs by divergent mechanisms. Overall, the ligands' inhibitory effect depends on their accessibility to the N5 position of the flavin's isoalloxazine ring. The results also indicated that the OYEs found in both parasites share structural similarities and showed affinities for the diterpene icetexanes in the same range. Nevertheless, the interaction between OYEs and ligands is directed by enthalpy and/or entropy in distinct ways. In conclusion, the binding site of both OYEs exhibits remarkable plasticity, and a large range of different molecules, including that can be substrates and inhibitors, can bind this site. This plasticity should be considered in drug design using OYE as a target.
Assuntos
Doença de Chagas , Leishmania braziliensis , Trypanosoma cruzi , Humanos , NADPH Desidrogenase/química , NADPH Desidrogenase/farmacologia , Doença de Chagas/parasitologia , Flavinas/farmacologiaRESUMO
Molecular crowding is a ubiquitous phenomenon in biological systems, with significant consequences on protein folding and stability. Small compounds, such as the osmolyte trimethylamine N-oxide (TMAO), can also present similar effects. To analyze the effects arising from these solute-like molecules, we performed a series of crowded coarse-grained folding simulations. Two well-known proteins were chosen, CI2 and SH3, modeled by the alpha-carbon-structure-based model. In the simulations, the crowding molecules were represented by low-sized neutral atom beads in different concentrations. The results show that a low level of the volume fraction occupied by neutral agents can change protein stability and folding kinetics for the two systems. However, the kinetics were shown to be unaffected in their respective folding temperatures. The faster kinetics correlates with changes in the folding route for systems at the same temperature, where non-native contacts were shown to be relevant. The transition states of the two systems with and without crowders are similar. In the case of SH3, there are differences in the structuring of two strands, which may be associated with the increase in its folding rate, in addition to the destabilization of the denatured ensemble. The present study also detected a crossover in the thermodynamic stability behavior, previously observed experimentally and theoretically. As the temperature increases, crowders change from destabilizing to stabilizing agents.
Assuntos
Dobramento de Proteína , Cinética , Estabilidade Proteica , Soluções , TermodinâmicaRESUMO
This work reports a simple strategy for Candida auris genomic DNA (gDNA) detection, a multi-resistant fungus associated with nosocomial outbreaks in healthcare settings, presenting high mortality and morbidity rates. The platform was developed using gold electrode sensitized with specific DNA capture probe and ninhydrin as a novel DNA hybridization indicator. The genosensor was able to detect C. auris in urine sample by differential pulse voltammetry and electrochemical impedance spectroscopy. The biosensor's analytical performance was evaluated by differential pulse voltammetry, detecting up to 4.5 pg µL-1 of C. auris gDNA in urine (1:10, V/V). Moreover, the genosensor was reused eight times with no loss in the current signal response. The genosensor showed selectivity and stability, maintaining 100% of its response up to 80 days of storage. In order to analyze interactions of single and double-stranded DNA with ninhydrin, SEM, AFM and molecular dynamics studies followed by docking simulations were performed. Theoretical calculations showed ninhydrin interactions more favorably with dsDNA in an A-T rich binding pocket rather than with the ssDNA. Therefore, the proposed system is a promising electrochemical detection device towards a more accurate detection of C. auris gDNA in biological samples.
Assuntos
Técnicas Biossensoriais , Técnicas Eletroquímicas , Candida/genética , DNA , NinidrinaRESUMO
Chagas disease remains a serious public health concern with unsatisfactory treatment outcomes due to strain-specific drug resistance and various side effects. To identify new therapeutic drugs against Trypanosoma cruzi, we evaluated both the in vitro and in vivo activity of the organometallic gold(III) complex [Au(III)(Hdamp)(L14)]Cl (L1 = SNS-donating thiosemicarbazone), henceforth denoted 4-Cl. Our results demonstrated that 4-Cl was more effective than benznidazole (Bz) in eliminating both the extracellular trypomastigote and intracellular amastigote forms of the parasite without cytotoxic effects on mammalian cells. In in vivo assays, 4-Cl in PBS solution loses the protonation and becomes the 4-neutral. 4-Neutral reduced parasitaemia and tissue parasitism in addition to protecting the liver and heart from tissue damage at 2.8 mg/kg/day. All these changes resulted in the survival of 100% of the mice treated with the gold complex during the acute phase. Analyzing the surviving animals of the acute infection, the parasite load after 150 days of infection was equivalent to those treated with the standard dose of Bz without demonstrating the hepatotoxicity of the latter. In addition, we identified a modulation of interferon gamma (IFN-γ) levels that may be targeting the disease's positive outcome. To the best of our knowledge, this is the first gold organometallic study that shows promise in an in vivo experimental model against Chagas disease.
Assuntos
Doença de Chagas/tratamento farmacológico , Ouro/química , Tiossemicarbazonas/química , Tiossemicarbazonas/farmacologia , Tripanossomicidas/química , Tripanossomicidas/farmacologia , Trypanosoma cruzi/efeitos dos fármacos , Animais , Linhagem Celular , Doença de Chagas/patologia , Cisteína Endopeptidases , Modelos Animais de Doenças , Resistência a Medicamentos/efeitos dos fármacos , Feminino , Coração , Humanos , Interferon gama/metabolismo , Fígado/patologia , Camundongos , Camundongos Endogâmicos BALB C , Simulação de Acoplamento Molecular , Músculo Esquelético/parasitologia , Músculo Esquelético/patologia , Nitroimidazóis , Compostos Organometálicos/química , Compostos Organometálicos/farmacologia , Parasitemia , Proteínas de Protozoários , Análise de SobrevidaRESUMO
Toxoplasma gondii is an intracellular protozoan parasite responsible for toxoplasmosis, which affects humans and animals. Serologic detection of anti-T. gondii immunoglobulins plays a crucial role in the clinical diagnosis of toxoplasmosis. In this work, a novel electrochemical immunosensor for detecting anti-T. gondii immunoglobulins is reported, based on immobilization of an in silico predicted peptide (PepB3), obtained from membrane protein of T. gondii, on the graphite electrode modified with poly(3-hydroxybenzoic acid). Indirect ELISA confirmed infection and binding specificity of peptide PepB3. Molecular modelling and simulations show this peptide binds to the T. gondii human Fab antibody in the surface antigen 1 (SAG1) binding site, remaining a stable complex during the molecular dynamic simulations, especially by hydrogen bonds and hydrophobic interactions. This electrochemical immunosensor was able to discriminate different periods of infection, using infected mouse serum samples, showing selectivity and discriminating infected and uninfected mouse serum.
Assuntos
Anticorpos Antiprotozoários/imunologia , Imunoglobulinas/imunologia , Peptídeos/imunologia , Toxoplasma/imunologia , Toxoplasmose/imunologia , Animais , Antígenos de Protozoários/imunologia , Ensaio de Imunoadsorção Enzimática/métodos , Feminino , Camundongos , Proteínas de Protozoários/imunologia , Sensibilidade e EspecificidadeRESUMO
Reactions of Ni(II) and Pd(II) precursors with S-benzyl-N-(ferrocenyl)methylenedithiocarbazate (HFedtc) led to the formation of heterobimetallic complexes of the type [MII(Fedtc)2] (Mâ¯=â¯Ni and Pd). The characterization of the compounds involved the determination of melting point, FTIR, UV-Vis, 1H NMR, elemental analysis and electrochemical experiments. Furthermore, the crystalline structures of HFedtc and [NiII(Fedtc)2] were determined by single crystal X-ray diffraction. The compounds were evaluated against the intracellular form of Trypanosoma cruzi (Tulahuen Lac-Z strain) and the cytotoxicity assays were assessed using LLC-MK2 cells. The results showed that the coordination of HFedtc to Ni(II) or Pd(II) decreases the in vitro trypanocidal activity while the cytotoxicity against LLC-MK2 cells does not change significantly. [PdII(Fedtc)2] showed the greater potential between the two complexes studied, showing an SI value of 8.9. However, this value is not better than that of the free ligand with an SI of 40, a similar value to that of the standard drug benznidazole (SIâ¯=â¯48). Additionally, molecular docking simulations were performed with Trypanosoma cruzi Old Yellow Enzyme (TcOYE), which predicted that HFedtc binds to the protein, almost parallel to the flavin mononucleotide (FMN) prosthetic group, while the [NiII(Fedtc)2] complex was docked into the enzyme binding site in a significantly different manner. In order to confirm the hypothetical interaction, in vitro experiments of fluorescence quenching and enzymatic activity were performed which indicated that, although HFedtc was not processed by the enzyme, it was able to act as a competitive inhibitor, blocking the hydride transfer from the FMN prosthetic group of the enzyme to the menadione substrate.
Assuntos
Compostos de Benzil/farmacologia , Complexos de Coordenação/farmacologia , Inibidores Enzimáticos/farmacologia , Hidrazinas/farmacologia , Metalocenos/farmacologia , NADPH Desidrogenase/antagonistas & inibidores , Níquel/farmacologia , Paládio/farmacologia , Tripanossomicidas/farmacologia , Trypanosoma cruzi/efeitos dos fármacos , Animais , Sobrevivência Celular/efeitos dos fármacos , Complexos de Coordenação/química , Complexos de Coordenação/metabolismo , Relação Dose-Resposta a Droga , Inibidores Enzimáticos/química , Inibidores Enzimáticos/metabolismo , Humanos , Hidrazinas/química , Macaca mulatta , Metalocenos/química , Simulação de Acoplamento Molecular , Estrutura Molecular , NADPH Desidrogenase/química , NADPH Desidrogenase/metabolismo , Níquel/química , Níquel/metabolismo , Paládio/química , Paládio/metabolismo , Relação Estrutura-Atividade , Tripanossomicidas/química , Tripanossomicidas/metabolismo , Trypanosoma cruzi/metabolismoRESUMO
The class D ß-lactamase OXA-143 has been described as an efficient penicillinase, oxacillinase, and carbapenemase. The D224A variant, known as OXA-231, was described in 2012 as exhibiting less activity toward imipenem and increased oxacillinase activity. Additionally, the P227S mutation was reported as a case of convergent evolution for homologous enzymes. To investigate the impact of both mutations (D224A and P227S), we describe in this paper a deep investigation of the enzymatic activities of these three homologues. OXA-143(P227S) presented enhanced catalytic activity against ampicillin, oxacillins, aztreonam, and carbapenems. In addition, OXA-143(P227S) was the only member capable of hydrolyzing ceftazidime. These enhanced activities were due to a combination of a higher affinity (lower Km) and a higher turnover number (higher kcat). We also determined the crystal structure of apo OXA-231. As expected, the structure of this variant is very similar to the published OXA-143 structure, except for the two M223 conformations and the absence of electron density for three solvent-exposed loop segments. Molecular dynamics calculations showed that both mutants experience higher flexibility compared to that of the wild-type form. Therefore, our results illustrate that D224A and P227S act as deleterious and positive mutations, respectively, within the evolutionary path of the OXA-143 subfamily toward a more efficient carbapenemase.
Assuntos
Acinetobacter baumannii/enzimologia , Carbapenêmicos/metabolismo , Modelos Moleculares , Mutação de Sentido Incorreto , beta-Lactamases/metabolismo , Ampicilina/metabolismo , Aztreonam/metabolismo , Ceftazidima , Hidrólise , Cinética , Simulação de Dinâmica Molecular , Oxacilina/metabolismo , Conformação Proteica em Folha beta , Estabilidade Proteica , Especificidade por Substrato , beta-Lactamases/genéticaRESUMO
Thermostable proteins are advantageous in industrial applications, as pharmaceuticals or biosensors, and as templates for directed evolution. As protein-design methodologies improve, bioengineers are able to design proteins to perform a desired function. Although many rationally designed proteins end up being thermostable, how to intentionally design de novo, thermostable proteins is less clear. UVF is a de novo-designed protein based on the backbone structure of the Engrailed homeodomain (EnHD) and is highly thermostable (Tm > 99°C vs. 52°C for EnHD). Although most proteins generally have polar amino acids on their surfaces and hydrophobic amino acids buried in their cores, protein engineers followed this rule exactly when designing UVF. To investigate the contributions of the fully hydrophobic core versus the fully polar surface to UVF's thermostability, we built two hybrid, chimeric proteins combining the sets of buried and surface residues from UVF and EnHD. Here, we determined a structural, dynamic, and thermodynamic explanation for UVF's thermostability by performing 4 µs of all-atom, explicit-solvent molecular dynamics simulations at 25 and 100°C, Tanford-Kirkwood solvent accessibility Monte Carlo electrostatic calculations, and a thermodynamic analysis of 40 temperature runs by the weighted-histogram analysis method of heavy-atom, structure-based models of UVF, EnHD, and both chimeric proteins. Our models showed that UVF was highly dynamic because of its fully hydrophobic core, leading to a smaller loss of entropy upon folding. The charged residues on its surface made favorable electrostatic interactions that contributed enthalpically to its thermostability. In the chimeric proteins, both the hydrophobic core and charged surface independently imparted thermostability.
Assuntos
Interações Hidrofóbicas e Hidrofílicas , Proteínas/química , Temperatura , Sequência de Aminoácidos , Entropia , Ligação de Hidrogênio , Simulação de Dinâmica Molecular , Método de Monte Carlo , Movimento , Conformação Proteica em alfa-Hélice , Estabilidade Proteica , Proteínas/metabolismo , Eletricidade EstáticaRESUMO
In order to improve the previously observed antichagasic activity of Cu(II) complexes containing 2-chlorobenzhydrazide (2-CH), we report herein the synthesis and anti-Trypanosoma cruzi activity of novel copper complexes containing 2-methoxybenzhydrazide (2-MH), 4-methoxybenzhydrazide (4-MH) and three α-diimine ligands, namely, 1,10-phenanthroline (phen), 2,2-bipyridine (bipy) and 4-4'-dimethoxy-2-2'-bipyridine (dmb). Two of these complexes showed higher in vitro anti-Trypanosoma cruzi activity when compared to benznidazole, the main drug used in Chagas disease treatment. One of them, the copper complex with 4-MH and dmb, [Cu(4-MH)(dmb)(ClO4)2], exhibited a higher selectivity index than that recommended for preclinical studies. Considering this observation, complex [Cu(4-MH)(dmb)(ClO4)2] was selected for preliminary in vivo assays, which verified that this compound was able to reduce parasitemia by 64% at the peak of infection. Further investigations were performed on all compounds. The Cu(II) complexes bind to ct-DNA with Kb values in the range of 103-104 M-1, with [Cu(4-MH)(dmb)(ClO4)2] showing the highest Kb value (1.45 × 104 M-1). Molecular docking simulations predicted that [Cu(4-MH)(dmb)(ClO4)2] binds in the minor groove of the double helix of ct-DNA and forms one hydrogen bond.
Assuntos
Doença de Chagas/tratamento farmacológico , Cobre/química , Tripanossomicidas/farmacologia , Trypanosoma cruzi/efeitos dos fármacos , Animais , Doença de Chagas/parasitologia , Complexos de Coordenação , Feminino , Hidrazinas/síntese química , Hidrazinas/farmacologia , Ligação de Hidrogênio , Ligantes , Camundongos , Camundongos Endogâmicos BALB C , Simulação de Acoplamento Molecular , Nitroimidazóis/farmacologia , Piridinas/síntese química , Piridinas/farmacologia , Tripanossomicidas/síntese químicaRESUMO
New complexes of composition [MX(HL1)] (M = PtII, PdII, X = Cl- or I-) and [MX(L1)] (M = AuIII, X = Cl-; M = PtII, PdII, X = PPh3) have been synthesized using a potentially tridentate thiosemicarbazone (H2L1) containing an additional oxime binding site. Among other analytical methods, all the seven complexes have been structurally characterized by single crystal X-ray diffractometry. Interesting structural features such as the influence of the halide ligands on hydrogen bonds and the formation of supramolecular structures for the phosphine derivatives are discussed. The in vitro trypanocidal activity of the free ligand H2L1 and its derivatives against both extracellular trypomastigote and intracellular amastigote (IC50try/ama) forms of Trypanosoma cruzi (Tulahuen Lac-Z strain) and the cytotoxicity was assessed on LLC-MK2 cell line. The results showed that complexation of the thiosemicarbazone ligand H2L1 to PtII, PdII and AuIII metal centers enhances the in vitro trypanocidal activity and that the cytotoxicity is dependent on both the metal center and coligands. Within the studied series, the AuIII complex showed the greatest potential, being not the most active but the most selective compound with a similar selectivity index to that of the standard drug benznidazole. In order to get a preliminary insight into the mechanism of action of these compounds, in vitro experiments of fluorescence quenching and enzymatic activity were performed using the AuIII complex and Trypanosoma cruzi Old Yellow Enzyme (TcOYE) which indicated that the gold derivative was capable of abstracting the hydride from the prosthetic FMN group of the enzyme. Additionally, molecular docking studies followed by semiempirical simulations showed that the [AuCl(L1)] binds to the binary complex TcOYE/FMN, almost parallel to the FMN prosthetic group, in a close distance that an electron/proton transfer might occur among them.
Assuntos
Compostos Organometálicos/farmacologia , Oximas/farmacologia , Tripanossomicidas/farmacologia , Trypanosoma cruzi/efeitos dos fármacos , Relação Dose-Resposta a Droga , Ouro/química , Ouro/farmacologia , Modelos Moleculares , Estrutura Molecular , Compostos Organometálicos/síntese química , Compostos Organometálicos/química , Oximas/química , Paládio/química , Paládio/farmacologia , Testes de Sensibilidade Parasitária , Platina/química , Platina/farmacologia , Relação Estrutura-Atividade , Tiossemicarbazonas/química , Tiossemicarbazonas/farmacologia , Tripanossomicidas/síntese química , Tripanossomicidas/químicaRESUMO
Tridentate thiosemicarbazone ligands with an ONS donor set, H2L(R) (R = Me and Et) were prepared by reactions of 1-phenyl-1,3-butanedione with 4-R-3-thiosemicarbazides. H2L(R) reacts with Na[AuCl4]·2H2O in MeOH in a 1:1 M ratio under formation of green gold(III) complexes of composition [AuCl(L(R))]. These compounds represent the first examples of gold(III) complexes with ONS chelate-bonded thiosemicarbazones. The in vitro anti-Trypanosoma cruzi activity against both trypomastigote and amastigote forms (IC50try/ama) of CL Brener strains as well as the cytotoxicity against LLC-MK2 cells of the free ligands and complexes was evaluated. The complex [AuCl(L(Me))] was found to be more active and more selective than its precursor ligand and the standard drug benznidazole with a SItry/ama value higher than 200, being considered as a lead candidate for Chagas disease treatment. Moreover the in vitro activity against the replicative amastigote form (IC50ama) of T. cruzi was additionally investigated revealing that [AuCl(L(Me))] was also more potent than benznidazole still with a similar selectivity index. Finally, docking studies showed that free ligands and complexes interact with the same residues of the parasite protease cruzain but with different intensities, suggesting that this protease could be a possible target for the trypanocidal action of the obtained compounds.
Assuntos
Complexos de Coordenação/farmacologia , Ouro/química , Tiossemicarbazonas/farmacologia , Tripanossomicidas/química , Animais , Doença de Chagas/tratamento farmacológico , Complexos de Coordenação/química , Simulação de Acoplamento Molecular , Relação Estrutura-Atividade , Tiossemicarbazonas/química , Tripanossomicidas/farmacologia , Trypanosoma cruzi/efeitos dos fármacosRESUMO
The energy landscape theory has been an invaluable theoretical framework in the understanding of biological processes such as protein folding, oligomerization, and functional transitions. According to the theory, the energy landscape of protein folding is funneled toward the native state, a conformational state that is consistent with the principle of minimal frustration. It has been accepted that real proteins are selected through natural evolution, satisfying the minimum frustration criterion. However, there is evidence that a low degree of frustration accelerates folding. We examined the interplay between topological and energetic protein frustration. We employed a Cα structure-based model for simulations with a controlled nonspecific energetic frustration added to the potential energy function. Thermodynamics and kinetics of a group of 19 proteins are completely characterized as a function of increasing level of energetic frustration. We observed two well-separated groups of proteins: one group where a little frustration enhances folding rates to an optimal value and another where any energetic frustration slows down folding. Protein energetic frustration regimes and their mechanisms are explained by the role of non-native contact interactions in different folding scenarios. These findings strongly correlate with the protein free-energy folding barrier and the absolute contact order parameters. These computational results are corroborated by principal component analysis and partial least square techniques. One simple theoretical model is proposed as a useful tool for experimentalists to predict the limits of improvements in real proteins.
Assuntos
Simulação de Dinâmica Molecular , Dobramento de Proteína , Proteínas , Análise dos Mínimos Quadrados , Análise Multivariada , Análise de Componente Principal , Proteínas/química , Proteínas/metabolismo , TermodinâmicaRESUMO
The energy landscape approach has played a fundamental role in advancing our understanding of protein folding. Here, we quantify protein folding energy landscapes by exploring the underlying density of states. We identify three quantities essential for characterizing landscape topography: the stabilizing energy gap between the native and nonnative ensembles δE, the energetic roughness ΔE, and the scale of landscape measured by the entropy S. We show that the dimensionless ratio between the gap, roughness, and entropy of the system Λ=δE/(ΔEâ(2S)) accurately predicts the thermodynamics, as well as the kinetics of folding. Large Λ implies that the energy gap (or landscape slope towards the native state) is dominant, leading to more funneled landscapes. We investigate the role of topological and energetic roughness for proteins of different sizes and for proteins of the same size, but with different structural topologies. The landscape topography ratio Λ is shown to be monotonically correlated with the thermodynamic stability against trapping, as characterized by the ratio of folding temperature versus trapping temperature. Furthermore, Λ also monotonically correlates with the folding kinetic rates. These results provide the quantitative bridge between the landscape topography and experimental folding measurements.
Assuntos
Modelos Químicos , Dobramento de Proteína , Proteínas/química , EntropiaRESUMO
Configuration-dependent diffusion (CDD) is important for protein folding kinetics with small thermodynamic barriers. CDD can be even more crucial in downhill folding without thermodynamic barriers. We explored the CDD of a downhill protein (BBL), and a two-state protein (CI2). The hidden kinetic barriers due to CDD were revealed. The increased ~1 k(B)T kinetic barrier is in line with experimental value based on other fast folding proteins. Compared to that of CI2, the effective free-energy profile of BBL is found to be significantly influenced by CDD, and the kinetics are totally determined by diffusion. These findings are consistent with both earlier bulk and single-molecule fluorescence measurements. In addition, we found the temperature dependence of CDD. We also found that the ratio of folding transition temperature against optimal kinetic folding temperature can provide both a quantitative measure for the underlying landscape topography and an indicator for the possible appearance of downhill folding. Our study can help for a better understanding of the role of diffusion in protein folding dynamics.
Assuntos
Proteínas/química , Difusão , Cinética , Dobramento de Proteína , Proteínas/metabolismo , Temperatura , TermodinâmicaRESUMO
We present a method for calculating the configurational-dependent diffusion coefficient of a globular protein as a function of the global folding process. Using a coarse-grained structure-based model, we determined the diffusion coefficient, in reaction coordinate space, as a function of the fraction of native contacts formed Q for the cold shock protein (TmCSP). We find nonmonotonic behavior for the diffusion coefficient, with high values for the folded and unfolded ensembles and a lower range of values in the transition state ensemble. We also characterized the folding landscape associated with an energetically frustrated variant of the model. We find that a low-level of frustration can actually stabilize the native ensemble and increase the associated diffusion coefficient. These findings can be understood from a mechanistic standpoint, in that the transition state ensemble has a more homogeneous structural content when frustration is present. Additionally, these findings are consistent with earlier calculations based on lattice models of protein folding and more recent single-molecule fluorescence measurements.
Assuntos
Proteínas de Bactérias/química , Proteínas de Bactérias/metabolismo , Dobramento de Proteína , Thermotoga maritima/metabolismo , Aminoácidos , Simulação por Computador , Difusão , Entropia , Modelos Moleculares , TemperaturaRESUMO
We developed both analytical and simulation methods to explore the diffusion dynamics in protein folding. We found the diffusion as a quantitative measure of escape from local traps along the protein folding funnel with chosen reaction coordinates has two remarkable effects on kinetics. At a fixed coordinate, local escape time depends on the distribution of barriers around it, therefore the diffusion is often time distributed. On the other hand, the environments (local escape barriers) change along the coordinates, therefore diffusion is coordinate dependent. The effects of time-dependent diffusion on folding can lead to non-exponential kinetics and non-Poisson statistics of folding time distribution. The effects of coordinate dependent diffusion on folding can lead to the change of the kinetic barrier height as well as the position of the corresponding transition state and therefore modify the folding kinetic rates as well as the kinetic routes. Our analytical models for folding are based on a generalized Fokker-Planck diffusion equation with diffusion coefficient both dependent on coordinate and time. Our simulation for folding are based on structure-based folding models with a specific fast folding protein CspTm studied experimentally on diffusion and folding with single molecules. The coordinate and time-dependent diffusion are especially important to be considered in fast folding and single molecule studies, when there is a small or no free energy barrier and kinetics is controlled by diffusion while underlying statistics of kinetics become important. Including the coordinate dependence of diffusion 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 will also have to be quantitatively modified. Complex kinetics with multiple time scales may allow us not only to explore the folding kinetics but also probe the local landscape and barrier height distribution with single-molecule experiments.
Assuntos
Dobramento de Proteína , Proteínas e Peptídeos de Choque Frio/química , Simulação por Computador , Difusão , Cinética , Termodinâmica , Thermotoga maritima/químicaRESUMO
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.