Proline-Rich Salivary Proteins Have Extended Conformations

Three basic proline-rich salivary proteins have been produced through the recombinant route. IB5 is a small basic proline-rich protein that is involved in the binding of plant tannins in the oral cavity. II-1 is a larger protein with a closely related backbone; it is glycosylated, and it is also able to bind plant tannins. II-1ng has the same polypeptidic backbone as II-1, but it is not glycosylated. Small angle x-ray scattering experiments on dilute solutions of these proteins confirm that they are intrinsically disordered. IB5 and II-1ng can be described through a chain model including a persistence length and cross section. The measured radii of gyration (R_g = 27.9 and 41.0 ± 1 Å respectively) and largest distances (r_max = 110 and 155 ± 10 Å respectively) show that their average conformations are rather extended. The length of the statistical segment (twice the persistence length) is b = 30 Å, which is larger than the usual value (18 Å − 20 Å) for unstructured polypeptide chains. These characteristics are presumably related to the presence of polyproline helices within the polypeptidic backbones. For both proteins, the radius of gyration of the chain cross-section is R_c = 2.7 ± 0.2Å. The glycosylated protein II-1 has similar conformations but the presence of large polyoside sidegroups yields the structure of a branched macromolecule with the same hydrophobic backbone and hydrophilic branches. It is proposed that the unusually extended conformations of these proteins in solution facilitate the capture of plant tannins in the oral cavity.     

Investigating Single Molecule Adhesion by Atomic Force Spectroscopy

Atomic force spectroscopy is an ideal tool to study molecules at surfaces and interfaces. An experimental protocol to couple a large variety of single molecules covalently onto an AFM tip is presented. At the same time the AFM tip is passivated to prevent unspecific interactions between the tip and the substrate, which is a prerequisite to study single molecules attached to the AFM tip. Analyses to determine the adhesion force, the adhesion length, and the free energy of these molecules on solid surfaces and bio-interfaces are shortly presented and external references for further reading are provided. Example molecules are the poly(amino acid) polytyrosine, the graft polymer PI-g-PS and the phospholipid POPE (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine). These molecules are desorbed from different surfaces like CH₃-SAMs, hydrogen terminated diamond and supported lipid bilayers under various solvent conditions. Finally, the advantages of force spectroscopic single molecule experiments are discussed including means to decide if truly a single molecule has been studied in the experiment.     

Urea Impedes the Hydrophobic Collapse of Partially Unfolded Proteins

Proteins are denatured in aqueous urea solution. The nature of the molecular driving forces has received substantial attention in the past, whereas the question how urea acts at different phases of unfolding is not yet well understood at the atomic level. In particular, it is unclear whether urea actively attacks folded proteins or instead stabilizes unfolded conformations. Here we investigated the effect of urea at different phases of unfolding by molecular dynamics simulations, and the behavior of partially unfolded states in both aqueous urea solution and in pure water was compared. Whereas the partially unfolded protein in water exhibited hydrophobic collapses as primary refolding events, it remained stable or even underwent further unfolding steps in aqueous urea solution. Further, initial unfolding steps of the folded protein were found not to be triggered by urea, but instead, stabilized. The underlying mechanism of this stabilization is a favorable interaction of urea with transiently exposed, less-polar residues and the protein backbone, thereby impeding back-reactions. Taken together, these results suggest that, quite generally, urea-induced protein unfolding proceeds primarily not by active attack. Rather, thermal fluctuations toward the unfolded state are stabilized and the hydrophobic collapse of partially unfolded proteins toward the native state is impeded. As a result, the equilibrium is shifted toward the unfolded state.     

Forces and Dynamics of Glucose and Inhibitor Binding to Sodium Glucose Co-transporter SGLT1 Studied by Single Molecule Force Spectroscopy

Single molecule force spectroscopy was employed to investigate the dynamics of the sodium glucose co-transporter (SGLT1) upon substrate and inhibitor binding on the single molecule level. CHO cells stably expressing rbSGLT1 were probed by using atomic force microscopy tips carrying either thioglucose, 2′-aminoethyl β-d-glucopyranoside, or aminophlorizin. Poly(ethylene glycol) (PEG) chains of different length and varying end groups were used as tether. Experiments were performed at 10, 25 and 37 °C to address different conformational states of SGLT1. Unbinding forces between ligands and SGLT1 were recorded at different loading rates by changing the retraction velocity, yielding binding probability, width of energy barrier of the binding pocket, and the kinetic off rate constant of the binding reaction. With increasing temperature, width of energy barrier and average life time increased for the interaction of SGLT1 with thioglucose (coupled via acrylamide to a long PEG) but decreased for aminophlorizin binding. The former indicates that in the membrane-bound SGLT1 the pathway to sugar translocation involves several steps with different temperature sensitivity. The latter suggests that also the aglucon binding sites for transport inhibitors have specific, temperature-sensitive conformations.     

Complex Stability of Single Proteins Explored by Forced Unfolding Experiments

In the last decade atomic force microscopy has been used to measure the mechanical stability of single proteins. These force spectroscopy experiments have shown that many water-soluble and membrane proteins unfold via one or more intermediates. Recently, Li and co-workers found a linear correlation between the unfolding force of the native state and the intermediate in fibronectin, which they suggested indicated the presence of a molecular memory or multiple unfolding pathways (1). Here, we apply two independent methods in combination with Monte Carlo simulations to analyze the unfolding of α-helices E and D of bacteriorhodopsin (BR). We show that correlation analysis of unfolding forces is very sensitive to errors in force calibration of the instrument. In contrast, a comparison of relative forces provides a robust measure for the stability of unfolding intermediates. The proposed approach detects three energetically different states of α-helices E and D in trimeric BR. These states are not observed for monomeric BR and indicate that substantial information is hidden in forced unfolding experiments of single proteins.     

Containment of Extended Length Polymorphisms in Silk Proteins

The spider silk gene family to the current date has been developed by gene duplication and homogenization events as well as conservation of crucial sequence parts. These evolutionary processes have created an amazing diversity of silk types each associated with specific properties and functions. In addition, they have led to allelic and gene variants within a species as exemplified by the major ampullate spidroin 1 gene of Nephila clavipes. Due to limited numbers of individuals screened to date little is known about the extent of these heterogeneities and how they are finally manifested in the proteins. Using expanded sample sizes, we show that sequence variations expressed as deletions or insertions of tri-nucleotides lead to different sized and structured repetitive units throughout a silk protein. Moreover, major ampullate spidroins 1 can quite dramatically differ in their overall lengths; however, extreme variants do not spread widely in a spider population. This suggests that a certain size range stabilized by purifying selection is important for spidroin 1 gene integrity and protein function. More than one locus for spidroin 1 genes possibly exist within one individual genome, which are homogenized in size, are differentially expressed and give a spider a certain degree of adaptation on silk’s composition and properties. Such mechanisms are shared to a lesser extent by the second major ampullate spidroin gene.     

An Extended Guinier Analysis for Intrinsically Disordered Proteins

Guinier analysis allows model-free determination of the radius of gyration (R_g) of a biomolecule from X-ray or neutron scattering data, in the limit of very small scattering angles. Its range of validity is well understood for globular proteins, but is known to be more restricted for unfolded or intrinsically disordered proteins (IDPs). We have used ensembles of disordered structures from molecular dynamics simulations to investigate which structural properties cause deviations from the Guinier approximation at small scattering angles. We find that the deviation from the Guinier approximation is correlated with the polymer scaling exponent ν describing the unfolded ensemble. We therefore introduce an empirical, ν-dependent, higher-order correction term, to augment the standard Guinier analysis. We test the new fitting scheme using all-atom simulation data for several IDPs and experimental data for both an IDP and a destabilized mutant of a folded protein. In all cases tested, we achieve an accuracy of the inferred R_g within ～ 3% of the true R_g. The method is straightforward to implement and extends the range of validity to a maximum qR_g of ～ 2 versus ～ 1.1 for Guinier analysis. Compared with the Guinier or Debye approaches, our method allows data from wider angles with lower noise to be used to analyze scattering data accurately. In addition to R_g, our fitting scheme also yields estimates of the scaling exponent ν in excellent agreement with the reference ν determined from the underlying molecular ensemble.     

Effective Binding Force Calculation in Dimeric Proteins

We apply the Blue Moon constrained Molecular Dynamics technique to study a particular case of molecular recognition, one of the main issues of modern molecular biology. We investigate the effects of mutation of interface residues on the binding strength of the dimeric protein superoxide dismutase from Photobacterium leiognathi. With our technique we produce a specific path describing the separation of the dimers and we calculate the effective mean force involved in the process. We apply the method to two mutants and compare the results with those obtained in an earlier calculation on the native enzyme. The method is sensitive to the mutations and allows us to establish a semi-quantitative hierarchy for the association strengths of the three enzymes.     

Prediction of methyl-side Chain Dynamics in Proteins

A simple analytical model is presented for the prediction of methyl-side chain dynamics in comparison with S(2) order parameters obtained by NMR relaxation spectroscopy. The model, which is an extension of the local contact model for backbone order parameter prediction, uses a static 3D protein structure as input. It expresses the methyl-group S(2) order parameters as a function of local contacts of the methyl carbon with respect to the neighboring atoms in combination with the number of consecutive mobile dihedral angles between the methyl group and the protein backbone. For six out of seven proteins the prediction results are good when compared with experimentally determined methyl-group S(2) values with an average correlation coefficient r = 0.65+/-0.14. For the unusually rigid cytochrome c(2) no significant correlation between prediction and experiment is found. The presented model provides independent support for the reliability of current side-chain relaxation methods along with their interpretation by the model-free formalism.     

Energy Landscape and Dynamics of Biomolecules Extended Abstract

Proteins are not isolated homogeneous systems. Each protein can exist in a very large number of conformations (conformational substates) that are characterized by an energy landscape. The main conformational motions, similar to the α and β fluctuations in glasses, are linked to fluctuations in the bulk solvent and the hydration shell.     

Rotamer Dynamics: Analysis of Rotamers in Molecular Dynamics Simulations of Proteins

Given by χ torsional angles, rotamers describe the side-chain conformations of amino acid residues in a protein based on the rotational isomers (hence the word rotamer). Constructed rotamer libraries, based on either protein crystal structures or dynamics studies, are the tools for classifying rotamers (torsional angles) in a way that reflect their frequency in nature. Rotamer libraries are routinely used in structure modeling and evaluation. In this perspective article, we would like to encourage researchers to apply rotamer analyses beyond their traditional use. Molecular dynamics (MD) of proteins highlight the in silico behavior of molecules in solution and thus can identify favorable side-chain conformations. In this article, we used simple computational tools to study rotamer dynamics (RD) in MD simulations. First, we isolated each frame in the MD trajectories in separate Protein Data Bank files via the cpptraj module in AMBER. Then, we extracted torsional angles via the Bio3D module in R language. The classification of torsional angles was also done in R according to the penultimate rotamer library. RD analysis is useful for various applications such as protein folding, study of rotamer-rotamer relationship in protein-protein interaction, real-time correlation between secondary structures and rotamers, study of flexibility of side chains in binding site for molecular docking preparations, use of RD as guide in functional analysis and study of structural changes caused by mutations, providing parameters for improving coarse-grained MD accuracy and speed, and many others. Major challenges facing RD to emerge as a new scientific field involve the validation of results via easy, inexpensive wet-lab methods. This realm is yet to be explored.     

Measuring “Unmeasurable” Folding Kinetics of Proteins by Single-Molecule Force Spectroscopy

Folding and unfolding are fundamental biological processes in cell and are important for the biological functions of proteins. Characterizing the folding and unfolding kinetics of proteins is important for understanding the energetic landscape leading to the active native conformations of these molecules. However, the thermal or chemical-induced unfolding of many proteins is irreversible in vitro, precluding characterization of the folding kinetics of such proteins, just as it is impossible to “un-boil” an egg. Irreversible unfolding often manifests as irreversible aggregation of unfolded polypeptide chains, which typically occurs between denatured protein molecules in response to the exposure of hydrophobic residues to solvent. An example of such a protein where thermal denaturation results in irreversible aggregation is the β-1,4 endoxylanase from Bacillus circulans (BCX). Here, we report the use of single-molecule atomic force microscopy to directly measure the folding kinetics of BCX in vitro. By mechanically unfolding BCX, we essentially allowed only one unfolded molecule to exist in solution at a given time, effectively eliminating the possibility for aggregation. We found that BCX can readily refold back to the native state, allowing us to measure its folding kinetics for the first time. Our results demonstrate that single-molecule force-spectroscopy-based methods can adequately tackle the challenge of “un-boiling eggs”, providing a general methodology to characterize the folding kinetics of many proteins that suffer from irreversible denaturation and thus cannot be characterized using traditional equilibrium methodologies.     

Microextraction of Nuclear Proteins from Single Maize Embryos

The analysis of DNA binding proteins can be difficult when only small quantities of tissue expressing the desired protein are available. We present a protocol for the preparation of nuclear extracts from as little as 100 mg of tissue. This protocol is well suited for extraction of DNA binding proteins from tissues that are difficult to obtain in large quantities such as maize embryos.     

Effective Population Size Dynamics and the Demographic Collapse of Bornean Orang-Utans

Bornean orang-utans experienced a major demographic decline and local extirpations during the Pleistocene and Holocene due to climate change, the arrival of modern humans, of farmers and recent commercially-driven habitat loss and fragmentation. The recent loss of habitat and its dramatic fragmentation has affected the patterns of genetic variability and differentiation among the remaining populations and increased the extinction risk of the most isolated ones. However, the contribution of recent demographic events to such genetic patterns is still not fully clear. Indeed, it can be difficult to separate the effects of recent anthropogenic fragmentation from the genetic signature of prehistoric demographic events. Here, we investigated the genetic structure and population size dynamics of orang-utans from different sites. Altogether 126 individuals were analyzed and a full-likelihood Bayesian approach was applied. All sites exhibited clear signals of population decline. Population structure is known to generate spurious bottleneck signals and we found that it does indeed contribute to the signals observed. However, population structure alone does not easily explain the observed patterns. The dating of the population decline varied across sites but was always within the 200–2000 years period. This suggests that in some sites at least, orang-utan populations were affected by demographic events that started before the recent anthropogenic effects that occurred in Borneo. These results do not mean that the recent forest exploitation did not leave its genetic mark on orang-utans but suggests that the genetic pool of orang-utans is also impacted by more ancient events. While we cannot identify the main cause for this decline, our results suggests that the decline may be related to the arrival of the first farmers or climatic events, and that more theoretical work is needed to understand how multiple demographic events impact the genome of species and how we can assess their relative contributions.     

Fluorescence Spectral Dynamics of Single LHCII Trimers

Single-molecule spectroscopy was employed to elucidate the fluorescence spectral heterogeneity and dynamics of individual, immobilized trimeric complexes of the main light-harvesting complex of plants in solution near room temperature. Rapid reversible spectral shifts between various emitting states, each of which was quasi-stable for seconds to tens of seconds, were observed for a fraction of the complexes. Most deviating states were characterized by the appearance of an additional, red-shifted emission band. Reversible shifts of up to 75 nm were detected. By combining modified Redfield theory with a disordered exciton model, fluorescence spectra with peaks between 670 nm and 705 nm could be explained by changes in the realization of the static disorder of the pigment-site energies. Spectral bands beyond this wavelength window suggest the presence of special protein conformations. We attribute the large red shifts to the mixing of an excitonic state with a charge-transfer state in two or more strongly coupled chlorophylls. Spectral bluing is explained by the formation of an energy trap before excitation energy equilibration is completed.