A statistical mechanical model for β-hairpin kinetics

Understanding the mechanism of protein secondary structure formation is an essential part of the protein-folding puzzle. Here, we describe a simple statistical mechanical model for the formation of a β-hairpin, the minimal structural element of the antiparallel β-pleated sheet. The model accurately describes the thermodynamic and kinetic behavior of a 16-residue, β-hairpin-forming peptide, successfully explaining its two-state behavior and apparent negative activation energy for folding. The model classifies structures according to their backbone conformation, defined by 15 pairs of dihedral angles, and is further simplified by considering only the 120 structures with contiguous stretches of native pairs of backbone dihedral angles. This single sequence approximation is tested by comparison with a more complete model that includes the 215 possible conformations and 15 × 215 possible kinetic transitions. Finally, we use the model to predict the equilibrium unfolding curves and kinetics for several variants of the β-hairpin peptide.

Folding funnels and frustration in off-lattice minimalist protein landscapes

A full quantitative understanding of the protein folding problem is now becoming possible with the help of the energy landscape theory and the protein folding funnel concept. Good folding sequences have a landscape that resembles a rough funnel where the energy bias towards the native state is larger than its ruggedness. Such a landscape leads not only to fast folding and stable native conformations but, more importantly, to sequences that are robust to variations in the protein environment and to sequence mutations. In this paper, an off-lattice model of sequences that fold into a β-barrel native structure is used to describe a framework that can quantitatively distinguish good and bad folders. The two sequences analyzed have the same native structure, but one of them is minimally frustrated whereas the other one exhibits a high degree of frustration.

Cytochrome c′ folding triggered by electron transfer: Fast and slow formation of four-helix bundles

Reduced (FeII) Rhodopseudomonas palustris cytochrome c′ (Cyt c′) is more stable toward unfolding ([GuHCl]1/2 = 2.9(1) M) than the oxidized (FeIII) protein ([GuHCl]1/2 = 1.9(1) M). The difference in folding free energies (ΔΔGf° = 70 meV) is less than half of the difference in reduction potentials of the folded protein (100 mV vs. NHE) and a free heme in aqueous solution (≈−150 mV). The spectroscopic features of unfolded FeII–Cyt c′ indicate a low-spin heme that is axially coordinated to methionine sulfur (Met-15 or Met-25). Time-resolved absorption measurements after CO photodissociation from unfolded FeII(CO)–Cyt c′ confirm that methionine can bind to the ferroheme on the microsecond time scale [kobs = 5(2) × 104 s−1]. Protein folding was initiated by photoreduction (two-photon laser excitation of NADH) of unfolded FeIII–Cyt c′ ([GuHCl] = 2.02–2.54 M). Folding kinetics monitored by heme absorption span a wide time range and are highly heterogeneous; there are fast-folding (≈103 s−1), intermediate-folding (102–101 s−1), and slow-folding (10−1 s−1) populations, with the last two likely containing methionine-ligated (Met-15 or Met-25) ferrohemes. Kinetics after photoreduction of unfolded FeIII–Cyt c′ in the presence of CO are attributable to CO binding [1.4(6) × 103 s−1] and FeII(CO)–Cyt c′ folding [2.8(9) s−1] processes; stopped-flow triggered folding of FeIII–Cyt c′ (which does not contain a protein-derived sixth ligand) is adequately described by a single kinetics phase with an estimated folding time constant of ≈4 ms [ΔGf° = −33(3) kJ mol−1] at zero denaturant.

Configurational diffusion down a folding funnel describes the dynamics of DNA hairpins

Elucidating the mechanism of folding of polynucleotides depends on accurate estimates of free energy surfaces and a quantitative description of the kinetics of structure formation. Here, the kinetics of hairpin formation in single-stranded DNA are measured after a laser temperature jump. The kinetics are modeled as configurational diffusion on a free energy surface obtained from a statistical mechanical description of equilibrium melting profiles. The effective diffusion coefficient is found to be strongly temperature-dependent in the nucleation step as a result of formation of misfolded loops that do not lead to subsequent zipping. This simple system exhibits many of the features predicted from theoretical studies of protein folding, including a funnel-like energy surface with many folding pathways, trapping in misfolded conformations, and non-Arrhenius folding rates.

Lattice model for rapidly folding protein-like heteropolymers.

Protein folding is a relatively fast process considering the astronomical number of conformations in which a protein could find itself. Within the framework of a lattice model, we show that one can design rapidly folding sequences by assigning the strongest attractive couplings to the contacts present in a target native state. Our protein design can be extended to situations with both attractive and repulsive contacts. Frustration is minimized by ensuring that all the native contacts are again strongly attractive. Strikingly, this ensures the inevitability of folding and accelerates the folding process by an order of magnitude. The evolutionary implications of our findings are discussed.

Interaction of GroEL with a highly structured folding intermediate: iterative binding cycles do not involve unfolding.

The GroE proteins are molecular chaperones involved in protein folding. The general mechanism by which they facilitate folding is still enigmatic. One of the central open questions is the conformation of the GroEL-bound nonnative protein. Several suggestions have been made concerning the folding stage at which a protein can interact with GroEL. Furthermore, the possibility exists that binding of the nonnative protein to GroEL results in its unfolding. We have addressed these issues that are basic for understanding the GroE-mediated folding cycle by using folding intermediates of an Fab antibody fragment as molecular probes to define the binding properties of GroEL. We show that, in addition to binding to an early folding intermediate, GroEL is able to recognize and interact with a late quaternary-structured folding intermediate (Dc) without measurably unfolding it. Thus, the prerequisite for binding is not a certain folding stage of a nonnative protein. In contrast, general surface properties of nonnative proteins seem to be crucial for binding. Furthermore, unfolding of a highly structured intermediate does not necessarily occur upon binding to GroEL. Folding of Dc in the presence of GroEL and ATP involves cycles of binding and release. Because in this system no off-pathway reactions or kinetic traps are involved, a quantitative analysis of the reactivation kinetics observed is possible. Our results indicate that the association reaction of Dc and GroEL in the presence of ATP is rather slow, whereas in the absence of ATP association is several orders of magnitude more efficient. Therefore, it seems that ATP functions by inhibiting reassociation rather than promoting release of the bound substrate.

Submillisecond folding of monomeric lambda repressor.

The folding kinetics of a truncated form of the N-terminal domain of phage lambda repressor [lambda 6-85] has been investigated by using the technique of dynamic NMR. lambda 6-85 has been shown previously to fold in a purely two-state fashion. This allows the determination of folding and unfolding rates from simulation of the exchange-broadened aromatic resonances of Tyr-22. The folding kinetics were determined over a range of 1.35 to 3.14 M urea. The urea dependence of both folding and unfolding rate constants is exponential, suggesting that the rate-determining step is invariant at the urea concentrations studied. The folding and unfolding rates extrapolated to 0 M urea at 37 degrees C are 3600 +/- 400 s-1 and 27 +/- 6 s-1, respectively. The observed lambda 6-85 folding rate constant exceeds that of other fast-folding globular proteins by a factor of 14-54. The urea dependence of the folding and unfolding rate constants suggests that the transition state of the rate-determining step is considerably more exposed to solvent than previously studied protein-folding transition states. The surprising rapidity of lambda 6-85 folding and unfolding may be the consequence of its all-helical secondary structure. These kinetic results clearly demonstrate that all of the fundamental events of protein folding can occur on the submillisecond time scale.

Studies of Protein-Protein and Protein-Water Interactions by Small Angle X-Ray Scattering, Terahertz Spectroscopy, ASMOS, And Computer Simulation

The protein folding problem has been one of the most challenging subjects in biological physics due to its complexity. Energy landscape theory based on statistical mechanics provides a thermodynamic interpretation of the protein folding process. We have been working to answer fundamental questions about protein-protein and protein-water interactions, which are very important for describing the energy landscape surface of proteins correctly. At first, we present a new method for computing protein-protein interaction potentials of solvated proteins directly from SAXS data. An ensemble of proteins was modeled by Metropolis Monte Carlo and Molecular Dynamics simulations, and the global X-ray scattering of the whole model ensemble was computed at each snapshot of the simulation. The interaction potential model was optimized and iterated by a Levenberg-Marquardt algorithm. Secondly, we report that terahertz spectroscopy directly probes hydration dynamics around proteins and determines the size of the dynamical hydration shell. We also present the sequence and pH-dependence of the hydration shell and the effect of the hydrophobicity. On the other hand, kinetic terahertz absorption (KITA) spectroscopy is introduced to study the refolding kinetics of ubiquitin and its mutants. KITA results are compared to small angle X-ray scattering, tryptophan fluorescence, and circular dichroism results. We propose that KITA monitors the rearrangement of hydrogen bonding during secondary structure formation. Finally, we present development of the automated single molecule operating system (ASMOS) for a high throughput single molecule detector, which levitates a single protein molecule in a 10 µm diameter droplet by the laser guidance. I also have performed supporting calculations and simulations with my own program codes.

A biochemical analysis into the co-translational folding of a G protein-coupled receptor; GPR35

The folding and targeting of membrane proteins poses a major challenge to the cell, as they must remain insertion competent while their highly hydrophobic transmembrane (TM) domains are transferred from the ribosome, through the aqueous cytosol and into the lipid bilayer. The biogenesis of a mature membrane protein takes place through the insertion and integration into the lipid bilayer. A number of TM proteins have been shown to gain some degree of secondary structure within the ribosome tunnel and to retain this conformation throughout maturation. Although studies into the folding and targeting of a number of membrane proteins have been carried out to date, there is little information on one of the largest class of eukaryotic membrane proteins; the G-protein-coupled receptors (GPCRs). This project studies the early folding events of the human ortholog of GPR35. To analyse the structure of the 1st TM domain, intermediates were generated and assessed by the biochemical method of pegylation (PEG-MAL). A structurally-similar microbial opsin (Bacterioopsin) was also used to investigate the differences in the early protein folding within eukaryotic and prokaryotic translation systems. Results showed that neither the 1st TM domain of GPR35 nor Bacterioopsin were capable of compacting in the ribosome tunnel before their N-terminus reached the ribosome exit point. The results for this assay remained consistent whether the proteins were translated in a eukaryotic or prokaryotic translation system. To examine the communication mechanism between the ribosome, the nascent chain and the protein targeting pathway, crosslinking experiments were carried out using the homobifunctional lysine cross-linker BS3. Specifically, the data generated here show that the nascent chain of GPR35 reaches the ribosomal protein uL23 in an extended conformation and interacts with the SRP protein as it exits the ribosome tunnel. This confirms the role of SRP in the co-translational targeting of GPR35. Using these methods insights into the early folding of GPCRs has been obtained. Further experiments using site-directed mutagenesis to reduce hydrophobicity in the 1st TM domain of GPR35, highlighted the mechanisms by which GPCRs are targeted to the endoplasmic reticulum. Confirming that hydrophobicity within the signal anchor sequence is essential of SRP-dependent targeting. Following the successful interaction of the nascent GPR35 and SRP, GPR35 is successfully targeted to ER membranes, shown here as dog pancreas microsomes (DPMs). Glycosylation of the GPR35 N-terminus was used to determine nascent chain structure as it is inserted into the ER membrane. These glycosylation experiments confirm that TM1 has obtained its compacted state whilst residing in the translocon. Finally, a site-specific cross-linking approach using the homobifunctional cysteine cross-linker, BMH, was used to study the lateral integration of GPR35 into the ER. Cross-linking of GPR35 TM1 and TM2 could be detected adjacent to a protein of ~45kDa, believed to be Sec61α. The loss of this adduct, as the nascent chain extends, showed the lateral movement of GPR35 TM1 from the translocon was dependent on the subsequent synthesis of TM2.

Distinct conformational properties determined by implicit and explicit representation of protein-solvent interactions. An analytical and computer simulation study

In the protein folding problem, solvent-mediated forces are commonly represented by intra-chain pairwise contact energy. Although this approximation has proven to be useful in several circumstances, it is limited in some other aspects of the problem. Here we show that it is possible to achieve two models to represent the chain-solvent system. one of them with implicit and other with explicit solvent, such that both reproduce the same thermodynamic results. Firstly, lattice models treated by analytical methods, were used to show that the implicit and explicitly representation of solvent effects can be energetically equivalent only if local solvent properties are time and spatially invariant. Following, applying the same reasoning Used for the lattice models, two inter-consistent Monte Carlo off-lattice models for implicit and explicit solvent are constructed, being that now in the latter the solvent properties are allowed to fluctuate. Then, it is shown that the chain configurational evolution as well as the globule equilibrium conformation are significantly distinct for implicit and explicit solvent systems. Actually, strongly contrasting with the implicit solvent version, the explicit solvent model predicts: (i) a malleable globule, in agreement with the estimated large protein-volume fluctuations; (ii) thermal conformational stability, resembling the conformational hear resistance of globular proteins, in which radii of gyration are practically insensitive to thermal effects over a relatively wide range of temperatures; and (iii) smaller radii of gyration at higher temperatures, indicating that the chain conformational entropy in the unfolded state is significantly smaller than that estimated from random coil configurations. Finally, we comment on the meaning of these results with respect to the understanding of the folding process. (C) 2009 Elsevier B.V. All rights reserved.

Protein sequence threading, the alignment problem, and a two-step strategy

Conventionally, protein structure prediction via threading relies on some nonoptimal method to align a protein sequence to each member of a library of known structures. We show how a score function (force field) can be modified so as to allow the direct application of a dynamic programming algorithm to the problem. This involves an approximation whose damage can be minimized by an optimization process during score function parameter determination. The method is compared to sequence to structure alignments using a more conventional pair-wise score function and the frozen approximation. The new method produces results comparable to the frozen approximation, but is faster and has fewer adjustable parameters. It is also free of memory of the template's original amino acid sequence, and does not suffer from a problem of nonconvergence, which can be shown to occur with the frozen approximation. Alignments generated by the simplified score function can then be ranked using a second score function with the approximations removed. (C) 1999 John Wiley & Sons, Inc.

Immunoelectron microscopy provides evidence for the presence of mitochondrial heat shock 10-kDa protein (chaperonin 10) in red blood cells and a variety of secretory granules

Hsp10 (10-kDa heat shock protein, also known as chaperonin 10 or Cpn10) is a co-chaperone for Hsp60 in the protein folding process. This protein has also been shown to be identical to the early pregnancy factor, which is an immunosuppressive growth factor found in maternal serum. In this study we have used immunogold electron microscopy to study the subcellular localization of Hsp10 in rat tissues sections embedded in LR Gold resin employing polyclonal antibodies raised against different regions of human Hsp10. In all rat tissues examined including liver, heart, pancreas, kidney, anterior pituitary, salivary gland, thyroid, and adrenal gland, antibodies to Hsp10 showed strong labeling of mitochondria. However, in a number of tissues, in addition to the mitochondrial labeling, strong and highly specific labeling with the Hsp10 antibodies was also observed in several extramitochondrial compartments. These sites included zymogen granules in pancreatic acinar cells, growth hormone granules in anterior pituitary, and secretory granules in PP pancreatic islet cells. Additionally, the mature red blood cells which lack mitochondria, also showed strong reactivity with the Hsp10 antibodies. The observed labeling with the Hsp10 antibodies, both within mitochondria as well as in other compartments/cells, was abolished upon omission of the primary antibodies or upon preadsorption of the primary antibodies with the purified recombinant human Hsp10. These results provide evidence that similar to a number of other recently described mitochondrial proteins (viz., Hsp60, tumor necrosis factor receptor-associated protein- 1, P32 (gC1q-R) protein, and cytochrome c), Hsp10 is also found at a variety of specific extramitochondrial sites in normal rat tissue. These results raise important questions as to how these mitochondrial proteins are translocated to other compartments and their possible function(s) at these sites. The presence of these proteins at extramitochondrial sites in normal tissues has important implications concerning the role of mitochondria in apoptosis and genetic diseases.

Integrating protein structural information

Dissertação apresentada para obtenção de Grau de Doutor em Bioquímica,Bioquímica Estrutural, pela Universidade Nova de Lisboa, Faculdade de Ciências e Tecnologia

The role of the efflux mechanisms in multidrug resistance in Mycobacterium tuberculosis

Abstract The emergence of multi and extensively drug resistant tuberculosis (MDRTB and XDRTB) has increased the concern of public health authorities around the world. The World Health Organization has defined MDRTB as tuberculosis (TB) caused by organisms resistant to at least isoniazid and rifampicin, the main first-line drugs used in TB therapy, whereas XDRTB refers to TB resistant not only to isoniazid and rifampicin, but also to a fluoroquinolone and to at least one of the three injectable second-line drugs, kanamycin, amikacin and capreomycin. Resistance in Mycobacterium tuberculosis is mainly due to the occurrence of spontaneous mutations and followed by selection of mutants by subsequent treatment. However, some resistant clinical isolates do not present mutations in any genes associated with resistance to a given antibiotic, which suggests that other mechanism(s) are involved in the development of drug resistance, namely the presence of efflux pump systems that extrude the drug to the exterior of the cell, preventing access to its target. Increased efflux activity can occur in response to prolonged exposure to subinhibitory concentrations of anti-TB drugs, a situation that may result from inadequate TB therapy. The inhibition of efflux activity with a non-antibiotic inhibitor may restore activity of an antibiotic subject to efflux and thus provide a way to enhance the activity of current anti-TB drugs. The work described in this thesis foccus on the study of efflux mechanisms in the development of multidrug resistance in M. tuberculosis and how phenotypic resistance, mediated by efflux pumps, correlates with genetic resistance. In order to accomplish this goal, several experimental protocols were developed using biological models such as Escherichia coli, the fast growing mycobacteria Mycobacterium smegmatis, and Mycobacterium avium, before their application to M. tuberculosis. This approach allowed the study of the mechanisms that result in the physiological adaptation of E. coli to subinhibitory concentrations of tetracycline (Chapter II), the development of a fluorometric method that allows the detection and quantification of efflux of ethidium bromide (Chapter III), the characterization of the ethidium bromide transport in M. smegmatis (Chapter IV) and the contribution of efflux activity to macrolide resistance in Mycobacterium avium complex (Chapter V). Finally, the methods developed allowed the study of the role of efflux pumps in M. tuberculosis strains induced to isoniazid resistance (Chapter VI). By this manner, in Chapter II it was possible to observe that the physiological adaptation of E. coli to tetracycline results from an interplay between events at the genetic level and protein folding that decrease permeability of the cell envelope and increase efflux pump activity. Furthermore, Chapter III describes the development of a semi-automated fluorometric method that allowed the correlation of this efflux activity with the transport kinetics of ethidium bromide (a known efflux pump substrate) in E. coli and the identification of efflux inhibitors. Concerning M. smegmatis, we have compared the wild-type M. smegmatis mc2155 with knockout mutants for LfrA and MspA for their ability to transport ethidium bromide. The results presented in Chapter IV showed that MspA, the major porin in M. smegmatis, plays an important role in the entrance of ethidium bromide and antibiotics into the cell and that efflux via the LfrA pump is involved in low-level resistance to these compounds in M. smegmatis. Chapter V describes the study of the contribution of efflux pumps to macrolide resistance in clinical M. avium complex isolates. It was demonstrated that resistance to clarithromycin was significantly reduced in the presence of efflux inhibitors such as thioridazine, chlorpromazine and verapamil. These same inhibitors decreased efflux of ethidium bromide and increased the retention of [14C]-erythromycin in these isolates. Finaly, the methods developed with the experimental models mentioned above allowed the study of the role of efflux pumps on M. tuberculosis strains induced to isoniazid resistance. This is described in Chapter VI of this Thesis, where it is demonstrated that induced resistance to isoniazid does not involve mutations in any of the genes known to be associated with isoniazid resistance, but an efflux system that is sensitive to efflux inhibitors. These inhibitors decreased the efflux of ethidium bromide and also reduced the minimum inhibitory concentration of isoniazid in these strains. Moreover, expression analysis showed overexpression of genes that code for efflux pumps in the induced strains relatively to the non-induced parental strains. In conclusion, the work described in this thesis demonstrates that efflux pumps play an important role in the development of drug resistance, namely in mycobacteria. A strategy to overcome efflux-mediated resistance may consist on the use of compounds that inhibit efflux activity, restoring the activity of antimicrobials that are efflux pump substrates, a useful approach particularly in TB where the most effective treatment regimens are becoming uneffective due to the increase of MDRTB/XDRTB.

Estudis funcionals i de plegament dels dominis d'activació de les procarboxipeptidases

Els dominis d’activació (ADs) de les procarboxipeptidases de la subfamília A/B sempre han sorprès ja que representen una quarta part del proenzim. S’han realitzat alguns estudis per intentar descobrir-ne alguna possible funció alternativa, però no han estat fructífers. El descobriment de l’elevada velocitat de plegament del domini d’activació de la procarboxipeptidasa A2 humana, (ADA2h), emperò, va portar a proposar la possibilitat de que realitzessin una funció d’assistència al plegament del domini enzimàtic. Posteriorment, l’anàlisi del plegament d’ADA2h a pH baix va revelar la capacitat d’aquest domini per formar fibres amiloides, a més de demostrar que un increment de l’estabilitat proteica podia prevenir la formació d’aquests agregats. La profunda caracterització del plegament d’ADA2h va fer que aquesta proteïna fos un bon model amiloidogènic, de manera que es van proposar un seguit d’experiments que s’han desenvolupat en el present treball per tal de conèixer millor aquest procés. S’han dut a terme estudis cinètics d’agregació per tal de valorar la contribució dels diferents aminoàcids de la seqüència polipeptídica, utilitzant 29 variants puntuals d’ADA2h. Es va eliminar la contribució de l’estabilitat mitjançant la utilització d’urea, i per dicroïsme circular conjuntament amb un aparell de flux detingut, es van obtenir dues velocitats diferents, v1 i v2, que corresponen a la formació d’un intermediari i a la seva reorganització, respectivament. Experiments complementaris utilitzant espectroscòpia d’infraroig (IR) revelaren la reorganització de l’estat natiu (en aquest cas) per a donar la forma agregada. Les cinètiques d’IR van mostrar que ADA2h forma l’estructura _ típica de les fibres amiloides, previ desplegament les seves hèlixs-_. Finalment, s’han realitzat estudis de biocomputació per tal d’esbrinar possibles funcions alternatives dels ADs. Les superposicions estructurals semblen mostrar similaritat dels ADs amb dominis de reconeixement d’RNA (RRM). Aquesta hipòtesi s’ha comprovat experimentalment amb ADA4h, mostrant una dèbil, però existent, unió a RNA.