Structural stability of oligomeric chaperonin 10: the role of two beta-strands at the N and C termini in structural stabilization

Chaperonin 10 (cpn10) is a well-conserved subgroup of the molecular chaperone family. GroES, the cpn10 from Escherichia coli, is composed of seven 10kDa subunits, which form a dome-like oligomeric ring structure. From our previous studies, it was found that GroES unfolded completely through a three-state unfolding mechanism involving a partly folded monomer and that this reaction was reversible. In order to study whether these unfolding-refolding characteristics were conserved in other cpn10 proteins, we have examined the structural stabilities of cpn10s from rat mitochondria (RatES) and from hyperthermophilic eubacteria Thermotoga maritima (TmaES), and compared the values to those of GroES. From size-exclusion chromatography experiments in the presence of various concentrations of Gdn-HCl at 25 degrees C, both cpn10s showed unfolding-refolding characteristics similar to those of GroES, i.e. two-stage unfolding reactions that include formation of a partially folded monomer. Although the partially folded monomer of TmaES was considerably more stable compared to GroES and RatES, it was found that the overall stabilities of all three cpn10s were achieved significantly by inter-subunit interactions. We studied this contribution of inter-subunit interactions to overall stability in the GroES heptamer by introducing a mutation that perturbed subunit association, specifically the interaction between the two anti-parallel beta-strands at the N and C termini of this protein. From analyses of the mutants' stabilities, it was revealed that the anti-parallel beta-strands at the subunit interface are crucial for subunit association and stabilization of the heptameric GroES protein.     

Structural stability of covalently linked GroES heptamer: advantages in the formation of oligomeric structure

In order to understand how inter-subunit association stabilizes oligomeric proteins, a single polypeptide chain variant of heptameric co-chaperonin GroES (tandem GroES) was constructed from Escherichia coli heptameric GroES by linking consecutively the C-terminal of one subunit to the N-terminal of the adjacent subunit with a small linker peptide. The tandem GroES (ESC7) showed properties similar to wild-type GroES in structural aspects and co-chaperonin activity. In unfolding and refolding equilibrium experiments using guanidine hydrochloride (Gdn-HCl) as a denaturant at a low protein concentration (50 microg ml(-1)), ESC7 showed a two-state transition with a greater resistance toward Gdn-HCl denaturation (Cm=1.95 M) compared to wild-type GroES (Cm=1.1 M). ESC7 was found to be about 10 kcal mol(-1) more stable than the wild-type GroES heptamer at 50 microg ml(-1). Kinetic unfolding and refolding experiments of ESC7 revealed that the increased stability was mainly attributed to a slower unfolding rate. Also a transient intermediate was detected in the refolding reaction. Interestingly, at the physiological GroES concentration (>1 mg ml(-1)), the free energy of unfolding for GroES heptamer exceeded that for ESC7. These results showed that at low protein concentrations (<1 mg ml(-1)), the covalent linking of subunits contributes to the stability but also complicates the refolding kinetics. At physiological concentrations of GroES, however, the oligomeric state is energetically preferred and the advantages of covalent linkage are lost. This finding highlights a possible advantage in transitioning from multi-domain proteins to oligomeric proteins with small subunits in order to improve structural and kinetic stabilities.     

Potentially amyloidogenic conformational intermediates populate the unfolding landscape of transthyretin: Insights from molecular dynamics simulations

Protein aggregation into insoluble fibrillar structures known as amyloid characterizes several neurodegenerative diseases, including Alzheimer's, Huntington's and Creutzfeldt‐Jakob. Transthyretin (TTR), a homotetrameric plasma protein, is known to be the causative agent of amyloid pathologies such as FAP (familial amyloid polyneuropathy), FAC (familial amyloid cardiomiopathy) and SSA (senile systemic amyloidosis). It is generally accepted that TTR tetramer dissociation and monomer partial unfolding precedes amyloid fibril formation. To explore the TTR unfolding landscape and to identify potential intermediate conformations with high tendency for amyloid formation, we have performed molecular dynamics unfolding simulations of WT‐TTR and L55P‐TTR, a highly amyloidogenic TTR variant. Our simulations in explicit water allow the identification of events that clearly discriminate the unfolding behavior of WT and L55P‐TTR. Analysis of the simulation trajectories show that (i) the L55P monomers unfold earlier and to a larger extent than the WT; (ii) the single α‐helix in the TTR monomer completely unfolds in most of the L55P simulations while remain folded in WT simulations; (iii) L55P forms, early in the simulations, aggregation‐prone conformations characterized by full displacement of strands C and D from the main β‐sandwich core of the monomer; (iv) L55P shows, late in the simulations, severe loss of the H‐bond network and consequent destabilization of the CBEF β‐sheet of the β‐sandwich; (v) WT forms aggregation‐compatible conformations only late in the simulations and upon extensive unfolding of the monomer. These results clearly show that, in comparison with WT, L55P‐TTR does present a much higher probability of forming transient conformations compatible with aggregation and amyloid formation.     

Mechanical unfolding pathway and origin of mechanical stability of proteins of ubiquitin family: An investigation by steered molecular dynamics simulation

Like the muscle protein Titin, proteins of the ubiquitin family exhibit a parallel strand arrangement, but otherwise having a distinctly different fold and not involved in an obvious load‐bearing function, exhibit high resistance to mechanical unfolding. We have applied all‐atom molecular dynamics simulation technique in implicit solvent to present a deep insight into the force‐induced unfolding pathway of three proteins—ubiquitin, NEDD8, and SUMO‐2—all having almost similar structural features. Two intermediates evolve in the unfolding pathway of each of the three proteins. The first intermediate, which has already been identified in case of ubiquitin by earlier simulation results, is similar for ubiquitin and NEDD8, but different in SUMO‐2. We have found a new intermediate with β3–β4 hairpin and some residual α‐helical character; and this intermediate is common for all the three proteins. Thus, proteins of the ubiquitin family pass through a well‐defined conformation in their force‐induced unfolding pathway. Reason behind the higher mechanical stability of the proteins with parallel strand structures like Titin has also been identified. Proteins 2009. © 2008 Wiley‐Liss, Inc.     

Recombination of protein fragments: a promising approach toward engineering proteins with novel nanomechanical properties

Combining single molecule atomic force microscopy (AFM) and protein engineering techniques, here we demonstrate that we can use recombination-based techniques to engineer novel elastomeric proteins by recombining protein fragments from structurally homologous parent proteins. Using I27 and I32 domains from the muscle protein titin as parent template proteins, we systematically shuffled the secondary structural elements of the two parent proteins and engineered 13 hybrid daughter proteins. Although I27 and I32 are highly homologous, and homology modeling predicted that the hybrid daughter proteins fold into structures that are similar to that of parent protein, we found that only eight of the 13 daughter proteins showed beta-sheet dominated structures that are similar to parent proteins, and the other five recombined proteins showed signatures of the formation of significant alpha-helical or random coil-like structure. Single molecule AFM revealed that six recombined daughter proteins are mechanically stable and exhibit mechanical properties that are different from the parent proteins. In contrast, another four of the hybrid proteins were found to be mechanically labile and unfold at forces that are lower than the approximately 20 pN, as we could not detect any unfolding force peaks. The last three hybrid proteins showed interesting duality in their mechanical unfolding behaviors. These results demonstrate the great potential of using recombination-based approaches to engineer novel elastomeric protein domains of diverse mechanical properties. Moreover, our results also revealed the challenges and complexity of developing a recombination-based approach into a laboratory-based directed evolution approach to engineer novel elastomeric proteins.     

Thermal unfolding of dodecameric glutamine synthetase: inhibition of aggregation by urea.

Thermal unfolding of dodecameric manganese glutamine synthetase (622,000 M(r)) at pH 7 and approximately 0.02 ionic strength occurs in two observable steps: a small reversible transition (Tm approximately 42 degrees C; delta H approximately equal to 0.9 J/g) followed by a large irreversible transition (Tm approximately 81 degrees C; delta H approximately equal to 23.4 J/g) in which secondary structure is lost and soluble aggregates form. Secondary structure, hydrophobicity, and oligomeric structure of the equilibrium intermediate are the same as for the native protein, whereas some aromatic residues are more exposed. Urea (3 M) destabilizes the dodecamer (with a tertiary structure similar to that without urea at 55 degrees C) and inhibits aggregation accompanying unfolding at < or = 0.2 mg protein/mL. With increasing temperature (30-70 degrees C) or incubation times at 25 degrees C (5-35 h) in 3 M urea, only dodecamer and unfolded monomer are detected. In addition, the loss in enzyme secondary structure is pseudo-first-order (t1/2 = 1,030 s at 20.0 degrees C in 4.5 M urea). Differential scanning calorimetry of the enzyme in 3 M urea shows one endotherm (Tmax approximately 64 degrees C; delta H = 17 +/- 2 J/g). The enthalpy change for dissociation and unfolding agrees with that determined by urea titrations by isothermal calorimetry (delta H = 57 +/- 15 J/g; Zolkiewski M, Nosworthy NJ, Ginsburg A, 1995, Protein Sci 4: 1544-1552), after correcting for the binding of urea to protein sites exposed during unfolding (-42 J/g). Refolding and assembly to active enzyme occurs upon dilution of urea after thermal unfolding.     

Urea-induced dissociation and unfolding of dodecameric glutamine synthetase from Escherichia coli: calorimetric and spectral studies.

Urea-induced dissociation and unfolding of manganese.glutamine synthetase (Mn.GS) have been studied at 37 degrees C (pH 7) by spectroscopic and calorimetric methods. In 0 to approximately 2 M urea, Mn.GS retains its dodecameric structure and full catalytic activity. Mn.GS is dissociated into subunits in 6 M urea, as evidenced by a 12-fold decrease in 90 degrees light scattering and a monomer molecular weight of 51,800 in sedimentation equilibrium studies. The light scattering decrease in 4 M urea parallels the time course of Trp exposure but occurs more rapidly than changes in secondary structure and Tyr exposure. Early and late kinetic steps appear to involve predominantly disruption of intra-ring and inter-ring subunit contacts, respectively, in the layered hexagonal structure of Mn.GS. The enthalpies for transferring Mn.GS into urea solutions have been measured by titration calorimetry. After correcting for the enthalpy of binding urea to the protein, the enthalpy of dissociation and unfolding of Mn.GS is 14 +/- 4 cal/g. A net proton uptake of approximately 50 H+/dodecamer accompanies unfolding reactions. The calorimetric data are consistent with urea binding to multiple, independent sites in Mn.GS and the number of binding sites increasing approximately 9-fold during the protein unfolding.     

Thermal unfolding molecular dynamics simulation of Escherichia coli dihydrofolate reductase: thermal stability of protein domains and unfolding pathway

Temperature induced unfolding of Escherichia coli dihydrofolate reductase was carried out by using molecular dynamic simulations. The simulations show that the unfolding generally involves an initial end-to-end collapse of the adenine binding domain into partially extended loops, followed by a gradual breakdown of the remaining beta sheet core structure. The core, which consists of beta strands 5-7, was observed to be the most resistant to thermal unfolding. This region, which is made up of part of the N terminus domain and part of the large domain of the E. coli dihydrofolate reductase, may constitute the nucleation site for protein folding and may be important for the eventual formation of both domains. The unfolding of different domains at different stages of the unfolding process suggests that protein domains vary in stability and that the rate at which they unfold can affect the overall outcome of the unfolding pathway. This observation is compared with the recently proposed hierarchical folding model. Finally, the results of the simulation were found to be consistent with a previous experimental study (Frieden, Proc Natl Acad Sci USA 1990;87:4413-4416) which showed that the folding process of E. coli dihydrofolate reductase involves sequential formation of the substrate binding sites.     

pH-induced equilibrium unfolding of apomyoglobin: Substitutions at conserved Trp14 and Met131 and non-conserved Val17 positions

A number of residues in globins family are well conserved but are not directly involved in the primary oxygen-carrying function of these proteins. A possible role for these conserved, non-functional residues has been suggested in promoting a rapid and correct folding process to the native tertiary structure. To test this hypothesis, we have studied pH-induced equilibrium unfolding of mutant apomyoglobins with substitutions of the conserved residues Trp14 and Met131, which are not involved in the function of myoglobin, by various amino acids. This allowed estimating their impact on the stability of various conformational states of the proteins and selecting conditions for a folding kinetics study. The results obtained from circular dichroism, tryptophan fluorescence, and differential scanning microcalorimetry for these mutant proteins were compared with those for the wild type protein and for a mutant with the non-conserved Val17 substituted by Ala. In the native folded state, all of the mutant apoproteins have a compact globular structure, but are destabilized in comparison to the wild type protein. The pH-induced denaturation of the mutant proteins occurs through the formation of a molten globule-like intermediate similar to that of the wild type protein. Thermodynamic parameters for all of the proteins were calculated using the three state model. Stability of equilibrium intermediates at pH ~4.0 was shown to be slightly affected by the mutations. Thus, all of the above substitutions influence the stability of the native state of these proteins. The cooperativity of conformational transitions and the exposed to solvent protein surface were also changed, but not for the substitution at Val17.     

Structure of stable protein folding intermediates by equilibrium phi-analysis: the apoflavodoxin thermal intermediate

Protein intermediates in equilibrium with native states may play important roles in protein dynamics but, in cases, can initiate harmful aggregation events. Investigating equilibrium protein intermediates is thus important for understanding protein behaviour (useful or pernicious) but it is hampered by difficulties in gathering structural information. We show here that the phi-analysis techniques developed to investigate transition states of protein folding can be extended to determine low-resolution three-dimensional structures of protein equilibrium intermediates. The analysis proposed is based solely on equilibrium data and is illustrated by determination of the structure of the apoflavodoxin thermal unfolding intermediate. In this conformation, a large part of the protein remains close to natively folded, but a 40 residue region is clearly unfolded. This structure is fully consistent with the NMR data gathered on an apoflavodoxin mutant designed specifically to stabilise the intermediate. The structure shows that the folded region of the intermediate is much larger than the proton slow-exchange core at 25 degrees C. It also reveals that the unfolded region is made of elements whose packing surface is more polar than average. In addition, it constitutes a useful guide to rationally stabilise the native state relative to the intermediate state, a far from trivial task.     

Pressure versus temperature unfolding of ribonuclease A: an FTIR spectroscopic characterization of 10 variants at the carboxy-terminal site

FTIR spectroscopy was used to characterize and compare the temperature- and pressure-induced unfolding of ribonuclease A and a set of its variants engineered in a hydrophobic region of the C-terminal part of the molecule postulated as a CFIS. The results show for all the ribonucleases investigated, a cooperative, two-state, reversible unfolding transition using both pressure and temperature. The relative stabilities, among the different sites and different variants at the same site, monitored either through the changes in the position of the maximum of the amide I' band and the tyrosine band, or the maximum of the band assigned to the beta-sheet structure, corroborate the results of a previous study using fourth-derivative UV absorbance spectroscopy. In addition, variants at position 108 are the most critical for ribonuclease structure and stability. The V108G variant seems to present a greater conformational flexibility than the other variants. The pressure- and temperature-denaturated states of all the ribonucleases characterized retained some secondary structure. However, their spectral maxima were centered at different wavenumbers, which suggests that pressure- and temperature-denaturated states do not have the same structural characteristics. Nevertheless, there was close correlation between the pressure and temperature midpoint transition values for the whole series of protein variants, which indicated a common tendency of stability toward pressure and heat.     

Studies of the unfolding of an unstable mutant of staphylococcal nuclease: Evidence for low temperature unfolding and compactness of the high temperature unfolded state

Fluorescence and circular dichroism data as a function of temperature were obtained to characterize the unfolding of nuclease A and two of its less stable mutants. These spectroscopic data were obtained with a modified instrument that enables the nearly simultaneous detection of both fluorescence and CD data on the same sample. A global analysis of these multiple datasets yielded an excellent fit of a model that includes a change in the heat capacity change, ΔC_p, between the unfolded and native states. This analysis gives a ΔC_p of 2.2 kcal/mol/·K for thermal unfolding of the WT protein and 1.3 and 1.8 kcal/mol/K for the two mutants. These ΔC_p values are consistent with significant population of the cold unfolded state at ∼0°C. Independent evidence for the existence of a cold unfolded state is the observation of a separately migrating peak in size exclusion chromatography. The new chromatographic peak is seen near 0°C, has a partition coefficient corresponding to a larger hydrodynamic radius, and shows a red-shifted fluorescence spectrum, as compared to the native protein. Data also indicate that the high-temperature unfolded form of mutant nuclease is relatively compact. Size exclusion chromatography shows the high temperature unfolded form to have a hydrodynamic radius that is larger than that for the native form, but smaller than that for the urea or pH-induced unfolded forms. Addition of chemical denaturants to the high-temperature unfolded form causes a further unfolding of the protein, as indicated by an increase in the apparent hydrodynamic radius and a decrease in the rotational correlation time for Trp140 (as determined by fluorescence anisotropy decay measurements). Proteins 28:227–240, 1997 © 1997 Wiley-Liss Inc.     

The unfolding of oxidized c-type cytochromes: the instructive case of Bacillus pasteurii

The reversible unfolding of oxidized Bacillus pasteurii cytochrome c(553) by guanidinium chloride under equilibrium conditions has been monitored by NMR and optical spectroscopy. The results obtained indicate that unfolding takes place through a mechanism involving the detachment from heme iron coordination of the sulfur of the Met71 axial ligand and yielding either a high spin (HS) or a low spin (LS(1)) species, depending on the pH value. In the LS(1) form the Met71 is replaced by another protein ligand, possibly Lys. The ligand exchange reaction does not reach completion until the protein backbone reaches a largely unfolded state, as monitored through 1H-15N NMR experiments, thus demonstrating that there is a significant correlation between formation of the Fe-S bond and native structure stability. 1H/2H exchange data, however, show that helix alpha(3), the C-terminal region of helix alpha(4), and helix alpha(5) maintain low exchangeability of the amide protons in the LS(1) form. This finding most likely implies that these regions maintain some ordered non-covalent structure, in which the amide moieties are involved in H-bonds. Finally, a folding mechanism is proposed and discussed in terms of analogies and differences with the larger mitochondrial cytochrome c proteins. It is concluded that the thermodynamic stability of the region around the metal cofactor is determined by the chemical nature of the residues around the axial methionine residue.     

Computational diagnosis of protein conformational diseases: short molecular dynamics simulations reveal a fast unfolding of r-LDL mutants that cause familial hypercholesterolemia

The molecular basis of conformational diseases frequently resides in mutant proteins constituting a subset of the vast mutational space. While the subtleties of protein structure point to molecular dynamics (MD) techniques as promising tools for an efficient exploration of such a space, the average size of proteins and the time scale of unfolding events make this goal difficult with present computational capabilities. We show here, nevertheless, that an efficient approach is already feasible for modular proteins. Familial hypercholesterolemia (FH) is a conformational disease linked to mutations in the gene encoding the low density lipoprotein receptor. A high percentage of these mutations has been found in the seven small modular binding repeats of the receptor. Taking advantage of its small size, we have performed an in depth MD study of the fifth binding repeat. Fast unfolding dynamics have been observed in the absence of a structural bound calcium ion, which agrees with its reported essential role in the stability of the module. In addition, several mutations detected in FH patients have been analyzed, starting from the native conformation. Our results indicate that in contrast with the wild type protein and an innocuous control mutant, disease-related mutants experience, in short simulation times (2-8 ns), gross departures from the native state that lead to unfolded conformations and, in some cases, to binding site desorganization deriving in calcium release. Computational diagnosis of mutations leading to conformational diseases seems thus feasible, at least for small or modular pathogenic proteins.     

Measuring unfolding of proteins in the presence of denaturant using fluorescence correlation spectroscopy

IFABP is a small (15 kDa) protein consisting mostly of antiparallel beta-strands that surround a large cavity into which ligands bind. We have previously used FCS to show that the native protein, labeled with fluorescein, exhibits dynamic fluctuation with a relaxation time of 35 micros. Here we report the use of FCS to study the unfolding of the protein induced by guanidine hydrochloride. Although the application of this technique to measure diffusion coefficients and molecular dynamics is straightforward, the FCS results need to be corrected for both viscosity and refractive index changes as the guanidine hydrochloride concentration increases. We present here a detailed study of the effects of viscosity and refractive index of guanidine hydrochloride solutions to calibrate FCS data. After correction, the increase in the diffusion time of IFABP corresponds well with the unfolding transition monitored by far ultraviolet circular dichroism. We also show that the magnitude of the 35 micros phase, reflecting the conformational fluctuation in the native state, decreases sharply as the concentration of denaturant increases and the protein unfolds. Although FCS experiments indicate that the unfolded state at pH 2 is rather compact and native-like, the radius in the presence of guanidine hydrochloride falls well within the range expected for a random coil.     

Mechanical unfolding of a titin Ig domain: structure of unfolding intermediate revealed by combining AFM,molecular dynamics simulations,NMR and protein engineering

The mechanical unfolding of an immunoglobulin domain from the human muscle protein titin (TI I27) has been shown to proceed via a metastable intermediate in which the A-strand is detached. The structure and properties of this intermediate are characterised in this study. A conservative destabilising mutation in the A-strand has no effect on the unfolding force, nor the dependence of the unfolding force on the pulling speed, indicating that the unfolding forces measured in an AFM experiment are those required for the unfolding of the intermediate and not the native state. A mutant of TI I27 with the A-strand deleted (TI I27-A) is studied by NMR and standard biophysical techniques, combined with protein engineering. Molecular dynamics simulations show TI I27-A to be a good model for the intermediate. It has a structure very similar to the native state, and is surprisingly stable. Comparison with a Phi-value analysis of the unfolding pathway clearly shows that the protein unfolds by a different pathway under an applied force than on addition of denaturant.     

Structural impact of cations on lipid bilayer models: Nanomechanical properties by AFM-force spectroscopy

Atomic Force Microscopy (AFM) has become an invaluable tool for studying the micro- and nanoworlds. As a stand-alone, high-resolution imaging technique and force transducer, it defies most other surface instrumentation in ease of use, sensitivity and versatility. The main strength of AFM relies on the possibility to operate in an aqueous environment on a wide variety of biological samples, from single molecules – DNA or proteins – to macromolecular assemblies like biological membranes. Understanding the effect of mechanical stress on membranes is of primary importance in biophysics, since cells are known to perform their function under a complex combination of forces. In the later years, AFM-based Force-Spectroscopy (AFM-FS) has provided a new vista on membrane mechanics in a confined area within the nanometer realm, where most of the specific molecular interactions take place. Lipid membranes are electrostatically charged entities that physiologically coexist with electrolyte solutions. Thus, specific interactions with ions are a matter of considerable interest. The distribution of ions in the solution and their interaction with the membranes are factors that substantially modify the structure and dynamics of the cell membranes. Furthermore, signaling processes are modified by the membrane capability of retaining ions. Supported Lipid Bilayers (SLBs) are a versatile tool to investigate phospholipid membranes mimicking biological surfaces. In the present contribution, we review selected experiments on the mechanical stability of SLBs as models of lipid membranes by means of AFM-FS, with special focus on the effect of cations and ionic strength in the overall nanomechanical stability.     

Temperature and pressure denaturation of chignolin: Folding and unfolding simulation by multibaric-multithermal molecular dynamics method

A multibaric‐multithermal molecular dynamics (MD) simulation of a 10‐residue protein, chignolin, was performed. All‐atom model with the Amber parm99SB force field was used for the protein and the TIP3P model was used for the explicit water molecules. This MD simulation covered wide ranges of temperature between 260 and 560 K and pressure between 0.1 and 600 MPa and sampled many conformations without getting trapped in local‐minimum free‐energy states. Folding events to the native β‐hairpin structure occurred five times and unfolding events were observed four times. As the temperature and/or pressure increases, fraction of folded chignolin decreases. The partial molar enthalpy change ΔH and partial molar volume change ΔV of unfolding were calculated as ΔH = 24.1 ± 4.9 kJ/mol and ΔV = ?5.6 ± 1.5 cm³/mol, respectively. These values agree well with recent experimental results. Illustrating typical local‐minimum free‐energy conformations, folding and unfolding pathways were revealed. When chignolin unfolds from the β‐hairpin structure, only the C terminus or both C and N termini open first. It may undergo an α‐helix or 3₁₀‐helix structure and finally unfolds to the extended structure. Difference of the mechanism between temperature denaturation and pressure denaturation is also discussed. Temperature denaturation is caused by making the protein transferred to a higher entropy state and making it move around more with larger space. The reason for pressure denaturation is that water molecules approach the hydrophobic residues, which are not well hydrated at the folded state, and some hydrophobic contacts are broken. Proteins 2012;. © 2012 Wiley Periodicals, Inc.     

Mechanical properties of a mature biofilm from a wastewater system: from microscale to macroscale level

A fundamental understanding of biofilm mechanical stability is critical in order to describe detachment and develop biofouling control strategies. It is thus important to characterise the elastic deformation and flow behaviour of the biofilm under different modes of applied force. In this study, the mechanical properties of a mature wastewater biofilm were investigated with methods including macroscale compression and microscale indentation using atomic force microscopy (AFM). The mature biofilm was found to be mechanically isotropic at the macroscale level as its mechanical properties did not depend on the scales and modes of loading. However, the biofilm showed a tendency for mechanical inhomogeneity at the microscale level as indentation progressed deeper into the matrix. Moreover, it was observed that the adhesion force had a significant influence on the elastic properties of the biofilm at the surface, subjected to microscale tensile loading. These results are expected to inform a damage-based model for biofilm detachment.     

Mechanical measurement of single-molecule binding rates: kinetics of DNA helix-destabilization by T4 gene 32 protein

Bacteriophage T4 gene 32 protein (gp32) is a single-stranded DNA (ssDNA) binding protein, and is essential for DNA replication, recombination and repair. While gp32 binds preferentially and cooperatively to ssDNA, it has not been observed to lower the thermal melting temperature of natural double-stranded DNA (dsDNA). However, in single-molecule stretching experiments, gp32 significantly destabilizes lambda DNA. In this study, we develop a theory of the effect of the protein on single dsDNA stretching curves, and apply it to the measured dependence of the DNA overstretching force on pulling rate in the presence of the full-length and two truncated forms of the protein. This allows us to calculate the rate of cooperative growth of single clusters of protein along ssDNA that are formed as the dsDNA molecule is stretched, as well as determine the site size of the protein binding to ssDNA. The rate of cooperative binding (ka) of both gp32 and of its proteolytic fragment *I (which lacks 48 residues from the C terminus) varies non-linearly with protein concentration, and appears to exceed the diffusion limit. We develop a model of protein association with the ends of growing clusters of cooperatively bound protein enhanced by 1-D diffusion along dsDNA, under the condition of protein excess. Upon globally fitting ka versus protein concentration, we determine the binding site size and the non-cooperative binding constants to dsDNA for gp32 and I. Our experiment mimics the growth of clusters of gp32 that likely exist at the DNA replication fork in vivo, and explains the origin of the "kinetic block" to dsDNA melting by gene 32 protein observed in thermal melting experiments.