Investigating Models of Protein Function and Allostery With a Widespread Mutational Analysis of a Light-Activated Protein

To investigate the relationship between a protein’s sequence and its biophysical properties, we studied the effects of more than 100 mutations in Avena sativa light-oxygen-voltage domain 2, a model protein of the Per-Arnt-Sim family. The A. sativa light–oxygen–voltage domain 2 undergoes a photocycle with a conformational change involving the unfolding of the terminal helices. Whereas selection studies typically search for winners in a large population and fail to characterize many sites, we characterized the biophysical consequences of mutations throughout the protein using NMR, circular dichroism, and ultraviolet/visible spectroscopy. Despite our intention to introduce highly disruptive substitutions, most had modest or no effect on function, and many could even be considered to be more photoactive. Substitutions at evolutionarily conserved sites can have minimal effect, whereas those at nonconserved positions can have large effects, contrary to the view that the effects of mutations, especially at conserved positions, are predictable. Using predictive models, we found that the effects of mutations on biophysical function and allostery reflect a complex mixture of multiple characteristics including location, character, electrostatics, and chemistry.     

Hydrogen exchange: the modern legacy of Linderstrøm-Lang.

This discussion, prepared for the Protein Society's symposium honoring the 100th anniversary of Kaj Linderstrøm-Lang, shows how hydrogen exchange approaches initially conceived and implemented by Lang and his colleagues some 50 years ago are contributing to current progress in structural biology. Examples are chosen from the active protein folding field. Hydrogen exchange methods now make it possible to define the structure of protein folding intermediates in various contexts: as tenuous molten globule forms at equilibrium under destabilizing conditions, in kinetic intermediates that exist for less than one second, and as infinitesimally populated excited state forms under native conditions. More generally, similar methods now find broad application in studies of protein structure, energetics, and interactions. This article considers the rise of these capabilities from their inception at the Carlsberg Labs to their contemporary role as a significant tool of modern structural biology.     

Orienting rigid and flexible biological assemblies in ferrofluids for small-angle neutron scattering studies

Small-angle scattering from macromolecules in solution is widely used to study their structures, but the information content is limited because the molecules are generally randomly oriented and hence the data are spherically averaged. The use of oriented rodlike structures for scattering, as in fiber diffraction, greatly increases the amount of structural detail that can be obtained. A new technique using a ferromagnetic fluid has been developed to align elongated structures independent of their intrinsic magnetic properties. This technique is ideal for small-angle neutron scattering because the scattering from the ferrofluid particles can be reduced significantly by matching the neutron scattering length density of the particles to a D₂O solvent (“contrast matching”). The net result is scattering primarily from the ordered biological assembly in a solution environment that can be adjusted to physiological pH and ionic strength. Scattering results from ordered tobacco mosaic virus, tobacco rattle virus, and chromain fibers are presented.     

A Novel Implicit Solvent Model for Simulating the Molecular Dynamics of RNA

Although molecular dynamics simulations can be accelerated by more than an order of magnitude by implicitly describing the influence of the solvent with a continuum model, most currently available implicit solvent simulations cannot robustly simulate the structure and dynamics of nucleic acids. The difficulties become exacerbated especially for RNAs, suggesting the presence of serious physical flaws in the prior continuum models for the influence of the solvent and counter ions on the nucleic acids. We present a novel, to our knowledge, implicit solvent model for simulating nucleic acids by combining the Langevin–Debye model and the Poisson–Boltzmann equation to provide a better estimate of the electrostatic screening of both the water and counter ions. Tests of the model involve comparisons of implicit and explicit solvent simulations for three RNA targets with 20, 29, and 75 nucleotides. The model provides reasonable agreement with explicit solvent simulations, and directions for future improvement are noted.     

The anabolic action of intermittent parathyroid hormone on cortical bone depends partly on its ability to induce nitric oxide‐mediated vasorelaxation in BALB/c mice

There is strong evidence that vasodilatory nitric oxide (NO) donors have anabolic effects on bone in humans. Parathyroid hormone (PTH), the only osteoanabolic drug currently approved, is also a vasodilator. We investigated whether the NO synthase inhibitor L‐NAME might alter the effect of PTH on bone by blocking its vasodilatory effect. BALB/c mice received 28 daily injections of PTH[1–34] (80 µg/kg/day) or L‐NAME (30 mg/kg/day), alone or in combination. Hindlimb blood perfusion was measured by laser Doppler imaging. Bone architecture, turnover and mechanical properties in the femur were analysed respectively by micro‐CT, histomorphometry and three‐point bending. PTH increased hindlimb blood flow by >30% within 10 min of injection (P < 0.001). Co‐treatment with L‐NAME blocked the action of PTH on blood flow, whereas L‐NAME alone had no effect. PTH treatment increased femoral cortical bone volume and formation rate by 20% and 110%, respectively (P < 0.001). PTH had no effect on trabecular bone volume in the femoral metaphysis although trabecular thickness and number were increased and decreased by 25%, respectively. Co‐treatment with L‐NAME restricted the PTH‐stimulated increase in cortical bone formation but had no clear‐cut effects in trabecular bone. Co‐treatment with L‐NAME did not affect the mechanical strength in femurs induced by iPTH. These results suggest that NO‐mediated vasorelaxation plays partly a role in the anabolic action of PTH on cortical bone. © 2016 The Authors. Cell Biochemistry and Function published by John Wiley & Sons, Ltd.     

Investigations into sequence and conformational dependence of backbone entropy,inter-basin dynamics and the Flory isolated-pair hypothesis for peptides

The populations and transitions between Ramachandran basins are studied for combinations of the standard 20 amino acids in monomers, dimers and trimers using an implicit solvent Langevin dynamics algorithm and employing seven commonly used force-fields. Both the basin populations and inter-conversion rates are influenced by the nearest neighbor's conformation and identity, contrary to the Flory isolated-pair hypothesis. This conclusion is robust to the choice of force-field, even though the use of different force-fields produces large variations in the populations and inter-conversion rates between the dominant helical, extended beta, and polyproline II basins. The computed variation of conformational and dynamical properties with different force-fields exceeds the difference between explicit and implicit solvent calculations using the same force-field. For all force-fields, the inter-basin transitions exhibit a directional dependence, with most transitions going through extended beta conformation, even when it is the least populated basin. The implications of these results are discussed in the context of estimates for the backbone entropy of single residues, and for the ability of all-atom simulations to reproduce experimental protein folding data.     

Understanding protein hydrogen bond formation with kinetic H/D amide isotope effects

Through the development of a procedure to measure when hydrogen bonds form under two-state folding conditions, alpha-helices have been determined to form proportionally to denaturant-sensitive surface area buried in the transition state. Previous experiments assessing H/D isotope effects are applied to various model proteins, including lambda and Arc repressor variants, a coiled coil domain, cytochrome c, colicin immunity protein 7, proteins L and G, acylphosphatase, chymotrypsin inhibitor II and a Src SH3 domain. The change in free energy accompanied by backbone deuteration is highly correlated to secondary structure composition when hydrogen bonds are divided into two classes. The number of helical hydrogen bonds correlates with an average equilibrium isotope effect of 8.6 +/- 0.9 cal x mol(-1) x site(-1). However, beta-sheet and long-range hydrogen bonds have little isotope effect. The kinetic isotope effects support our hypothesis that, for helical proteins, hydrophobic association cannot be separated from helix formation in the transition state. Therefore, folding models that describe an incremental build-up of structure in which hydrophobic burial and hydrogen bond formation occur commensurately are more consistent with the data than are models that posit the extensive formation of one quantity before the other.