Communication between noncontacting macromolecules

We present a quantitative experimental demonstration of solvent-mediated communication between noncontacting biopolymers. We show that changes in the activity of a solvent component brought about by a conformational change in one biopolymer can result in changes in the physical properties of a second noncontacting biopolymer present in solution. Specifically, we show that the release of protons on denaturation of a donor polymer (in this case, a four-stranded DNA tetraplex, iDNA) modulates the melting temperature of a noncontacting, acceptor polymer [in this case poly(A)]. In addition to such proton-mediated cross talk, we also demonstrate counterion-mediated cross talk between noncontacting biopolymers. Specifically, we show that counterion association/release on denaturation of native salmon sperm DNA (the donor polymer) can modulate the melting temperature of poly(dA)⋅poly(dT) (the acceptor polymer). Taken together, these two examples demonstrate how poly(A) and poly(dA)⋅poly(dT) can serve as molecular probes that report the pH and free salt concentrations in solution, respectively. Further, we demonstrate how such through-solvent dialogue between biopolymers that do not directly interact can be used to evaluate (in a model-free manner) association/dissociation reactions of solvent components (e.g., protons, sodium cations) with one of the two biopolymers. We propose that such through-solution dialogue is a general property of all biopolymers. As a result, such solvent-mediated cross talk should be considered when assessing reactions of multicomponent systems such as those that exist in essentially all biological processes.

Simple model of protein folding kinetics.

A simple model of the kinetics of protein folding is presented. The reaction coordinate is the "correctness" of a configuration compared with the native state. The model has a gap in the energy spectrum, a large configurational entropy, a free energy barrier between folded and partially folded states, and a good thermodynamic folding transition. Folding kinetics is described by a master equation. The folding time is estimated by means of a local thermodynamic equilibrium assumption and then is calculated both numerically and analytically by solving the master equation. The folding time has a maximum near the folding transition temperature and can have a minimum at a lower temperature.

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.

Non-Arrhenius kinetics for the loop closure of a DNA hairpin

Intramolecular chain diffusion is an elementary process in the conformational fluctuations of the DNA hairpin-loop. We have studied the temperature and viscosity dependence of a model DNA hairpin-loop by FRET (fluorescence resonance energy transfer) fluctuation spectroscopy (FRETfs). Apparent thermodynamic parameters were obtained by analyzing the correlation amplitude through a two-state model and are consistent with steady-state fluorescence measurements. The kinetics of closing the loop show non-Arrhenius behavior, in agreement with theoretical prediction and other experimental measurements on peptide folding. The fluctuation rates show a fractional power dependence (β = 0.83) on the solution viscosity. A much slower intrachain diffusion coefficient in comparison to that of polypeptides was derived based on the first passage time theory of SSS [Szabo, A., Schulten, K. & Schulten, Z. (1980) J. Chem. Phys. 72, 4350–4357], suggesting that intrachain interactions, especially stacking interaction in the loop, might increase the roughness of the free energy surface of the DNA hairpin-loop.

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.

Species specificity in the cell-free conversion of prion protein to protease-resistant forms: a model for the scrapie species barrier.

Scrapie is a transmissible neurodegenerative disease that appears to result from an accumulation in the brain of an abnormal protease-resistant isoform of prion protein (PrP) called PrPsc. Conversion of the normal, protease-sensitive form of PrP (PrPc) to protease-resistant forms like PrPsc has been demonstrated in a cell-free reaction composed largely of hamster PrPc and PrPsc. We now report studies of the species specificity of this cell-free reaction using mouse, hamster, and chimeric PrP molecules. Combinations of hamster PrPc with hamster PrPsc and mouse PrPc with mouse PrPsc resulted in the conversion of PrPc to protease-resistant forms. Protease-resistant PrP species were also generated in the nonhomologous reaction of hamster PrPc with mouse PrPsc, but little conversion was observed in the reciprocal reaction. Glycosylation of the PrPc precursors was not required for species specificity in the conversion reaction. The relative conversion efficiencies correlated with the relative transmissibilities of these strains of scrapie between mice and hamsters. Conversion experiments performed with chimeric mouse/hamster PrPc precursors indicated that differences between PrPc and PrPsc at residues 139, 155, and 170 affected the conversion efficiency and the size of the resultant protease-resistant PrP species. We conclude that there is species specificity in the cell-free interactions that lead to the conversion of PrPc to protease-resistant forms. This specificity may be the molecular basis for the barriers to interspecies transmission of scrapie and other transmissible spongiform encephalopathies in vivo.

Rapid treadmilling of brain microtubules free of microtubule-associated proteins in vitro and its suppression by tau

We have determined the treadmilling rate of brain microtubules (MTs) free of MT-associated proteins (MAPs) at polymer mass steady state in vitro by using [3H]GTP-exchange. We developed buffer conditions that suppressed dynamic instability behavior by ≈10-fold to minimize the contribution of dynamic instability to total tubulin-GTP exchange. The MTs treadmilled rapidly under the suppressed dynamic instability conditions, at a minimum rate of 0.2 μm/min. Thus, rapid treadmilling is an intrinsic property of MAP-free MTs. Further, we show that tau, an axonal stabilizing MAP involved in Alzheimer’s disease, strongly suppresses the treadmilling rate. These results indicate that tau’s function in axons might involve suppression of axonal MT treadmilling. We describe mathematically how treadmilling and dynamic instability are mechanistically distinct MT behaviors. Finally, we present a model that explains how small changes in the critical tubulin subunit concentration at MT minus ends, caused by intrinsic differences in rate constants or regulatory proteins, could produce large changes in the treadmilling rate.

Exploring the origins of topological frustration: Design of a minimally frustrated model of fragment B of protein A

Topological frustration in an energetically unfrustrated off-lattice model of the helical protein fragment B of protein A from Staphylococcus aureus was investigated. This Gō-type model exhibited thermodynamic and kinetic signatures of a well-designed two-state folder with concurrent collapse and folding transitions and single exponential kinetics at the transition temperature. Topological frustration is determined in the absence of energetic frustration by the distribution of Fersht φ values. Topologically unfrustrated systems present a unimodal distribution sharply peaked at intermediate φ, whereas highly frustrated systems display a bimodal distribution peaked at low and high φ values. The distribution of φ values in protein A was determined both thermodynamically and kinetically. Both methods yielded a unimodal distribution centered at φ = 0.3 with tails extending to low and high φ values, indicating the presence of a small amount of topological frustration. The contacts with high φ values were located in the turn regions between helices I and II and II and III, intimating that these hairpins are in large part required in the transition state. Our results are in good agreement with all-atom simulations of protein A, as well as lattice simulations of a three- letter code 27-mer (which can be compared with a 60-residue helical protein). The relatively broad unimodal distribution of φ values obtained from the all-atom simulations and that from the minimalist model for the same native fold suggest that the structure of the transition state ensemble is determined mostly by the protein topology and not energetic frustration.

RNA folding kinetics regulates translation of phage MS2 maturation gene

The gene for the maturation protein of the single-stranded RNA coliphage MS2 is preceded by an untranslated leader of 130 nt, which folds into a cloverleaf, i.e., three stem–loop structures enclosed by a long distance interaction (LDI). This LDI prevents translation because its 3′ moiety contains the Shine–Dalgarno sequence of the maturation gene. Previously, several observations suggested that folding of the cloverleaf is kinetically delayed, providing a time window for ribosomes to access the RNA. Here we present direct evidence for this model. In vitro experiments show that ribosome binding to the maturation gene is faster than refolding of the denatured cloverleaf. This folding delay appears related to special properties of the leader sequence. We have replaced the three stem–loop structures by a single five nt loop. This change does not affect the equilibrium structure of the LDI. Nevertheless, in this construct, the folding delay has virtually disappeared, suggesting that now the RNA folds faster than ribosomes can bind. Perturbation of the cloverleaf by an insertion makes the maturation start permanently accessible. A pseudorevertant that evolved from an infectious clone carrying the insertion had overcome this defect. It showed a wild-type folding delay before closing down the maturation gene. This experiment reveals the biological significance of retarded cloverleaf formation.

Exploring the folding free energy surface of a three-helix bundle protein

The multidimensional free energy surface for a small fast folding helical protein is explored based on first-principle calculations. The model represents the 46-residue segment from fragment B of staphylococcal protein A. The relationship between collapse and tertiary structure formation, and the order of collapse and secondary structure formation, are investigated. We find that the initial collapse process gives rise to a transition state with about 30% of the native tertiary structure and 50–70% of the native helix content. We also observe two distinct distributions of native helix in this collapsed state (Rg ≈ 12 Å), one with about 20% of the native helical hydrogen bonds, the other with near 70%. The former corresponds to a local minimum. The barrier from this metastable state to the native state is about 2 kBT. In the latter case, folding is essentially a downhill process involving topological assembly. In addition, the order of formation of secondary structure among the three helices is examined. We observe cooperative formation of the secondary structure in helix I and helix II. Secondary structure in helix III starts to form following the formation of certain secondary structure in both helix I and helix II. Comparisons of our results with those from theory and experiment are made.

Identification of the prooxidant site of human ceruloplasmin: A model for oxidative damage by copper bound to protein surfaces

Free transition metal ions oxidize lipids and lipoproteins in vitro; however, recent evidence suggests that free metal ion-independent mechanisms are more likely in vivo. We have shown previously that human ceruloplasmin (Cp), a serum protein containing seven Cu atoms, induces low density lipoprotein oxidation in vitro and that the activity depends on the presence of a single, chelatable Cu atom. We here use biochemical and molecular approaches to determine the site responsible for Cp prooxidant activity. Experiments with the His-specific reagent diethylpyrocarbonate (DEPC) showed that one or more His residues was specifically required. Quantitative [14C]DEPC binding studies indicated the importance of a single His residue because only one was exposed upon removal of the prooxidant Cu. Plasmin digestion of [14C]DEPC-treated Cp (and N-terminal sequence analysis of the fragments) showed that the critical His was in a 17-kDa region containing four His residues in the second major sequence homology domain of Cp. A full length human Cp cDNA was modified by site-directed mutagenesis to give His-to-Ala substitutions at each of the four positions and was transfected into COS-7 cells, and low density lipoprotein oxidation was measured. The prooxidant site was localized to a region containing His426 because CpH426A almost completely lacked prooxidant activity whereas the other mutants expressed normal activity. These observations support the hypothesis that Cu bound at specific sites on protein surfaces can cause oxidative damage to macromolecules in their environment. Cp may serve as a model protein for understanding mechanisms of oxidant damage by copper-containing (or -binding) proteins such as Cu, Zn superoxide dismutase, and amyloid precursor protein.

From residue matching patterns to protein folding topographies: General model and bovine pancreatic trypsin inhibitor

A coarse-grained model for protein-folding dynamics is introduced based on a discretized representation of torsional modes. The model, based on the Ramachandran map of the local torsional potential surface and the class (hydrophobic/polar/neutral) of each residue, recognizes patterns of both torsional conformations and hydrophobic-polar contacts, with tolerance for imperfect patterns. It incorporates empirical rates for formation of secondary and tertiary structure. The method yields a topological representation of the evolving local torsional configuration of the folding protein, modulo the basins of the Ramachandran map. The folding process is modeled as a sequence of transitions from one contact pattern to another, as the torsional patterns evolve. We test the model by applying it to the folding process of bovine pancreatic trypsin inhibitor, obtaining a kinetic description of the transitions between the contact patterns visited by the protein along the dominant folding pathway. The kinetics and detailed balance make it possible to invert the result to obtain a coarse topographic description of the potential energy surface along the dominant folding pathway, in effect to go backward or forward between a topological representation of the chain conformation and a topographical description of the potential energy surface governing the folding process. As a result, the strong structure-seeking character of bovine pancreatic trypsin inhibitor and the principal features of its folding pathway are reproduced in a reasonably quantitative way.

Molecular dynamics and free-energy calculations applied to affinity maturation in antibody 48G7

We investigated the relative free energies of hapten binding to the germ line and mature forms of the 48G7 antibody Fab fragments by applying a continuum model to structures sampled from molecular dynamics simulations in explicit solvent. Reasonable absolute and very good relative free energies were obtained. As a result of nine somatic mutations that do not contact the hapten, the affinity-matured antibody binds the hapten >104 tighter than the germ line antibody. Energetic analysis reveals that van der Waals interactions and nonpolar contributions to solvation are similar and drive the formations of both the germ line and mature antibody–hapten complexes. Affinity maturation of the 48G7 antibody therefore appears to occur through reorganization of the combining site geometry in a manner that optimizes the balance of gaining favorable electrostatic interactions with the hapten and losing those with solvent during the binding process. As reflected by lower rms fluctuations in the antibody–hapten complex, the mature complex undergoes more restricted fluctuations than the germ line complex. The dramatically increased affinity of the 48G7 antibody over its germ line precursor is thus made possible by electrostatic optimization.

Pressure-induced protein-folding/unfolding kinetics

We use an off-lattice minimalist model to describe the effects of pressure in slowing down the folding/unfolding kinetics of proteins when subjected to increasingly larger pressures. The potential energy function used to describe the interactions between beads in the model includes the effects of pressure on the pairwise interaction of hydrophobic groups in water. We show that pressure affects the participation of contacts in the transition state. More significantly, pressure exponentially decreases the chain reconfigurational diffusion coefficient. These results are consistent with experimental results on the kinetics of pressure-denaturation of staphylococcal nuclease.

Flickering fusion pores comparable with initial exocytotic pores occur in protein-free phospholipid bilayers

For the act of membrane fusion, there are two competing, mutually exclusive molecular models that differ in the structure of the initial pore, the pathway for ionic continuity between formerly separated volumes. Because biological “fusion pores” can be as small as ionic channels or gap junctions, one model posits a proteinaceous initial fusion pore. Because biological fusion pore conductance varies widely, another model proposes a lipidic initial pore. We have found pore opening and flickering during the fusion of protein-free phospholipid vesicles with planar phospholipid bilayers. Fusion pore formation appears to follow the coalescence of contacting monolayers to create a zone of hemifusion where continuity between the two adherent membranes is lipidic, but not aqueous. Hypotonic stress, causing tension in the vesicle membrane, promotes complete fusion. Pores closed soon after opening (flickering), and the distribution of fusion pore conductance appears similar to the distribution of initial fusion pores in biological fusion. Because small flickering pores can form in the absence of protein, the existence of small pores in biological fusion cannot be an argument in support of models based on proteinaceous pores. Rather, these results support the model of a lipidic fusion pore developing within a hemifused contact site.