Lorna J. Smith

Understanding protein folding via free-energy surfaces from theory and experiment.

The ability of protein molecules to fold into their highly structured functional states is one of the most remarkable evolutionary achievements of biology. In recent years, our understanding of the way in which this complex self-assembly process takes place has increased dramatically. Much of the reason for this advance has been the development of energy surfaces (landscapes), which allow the folding reaction to be described and visualized in a meaningful manner. Analysis of these surfaces, derived from the constructive interplay between theory and experiment, has led to the development of a unified mechanism for folding and a recognition of the underlying factors that control the rates and products of the folding process.

Charge-state dependent compaction and dissociation of protein complexes: Insights from ion mobility and molecular dynamics.

Collapse to compact states in the gas phase, with smaller collision cross sections than calculated for their native-like structure, has been reported previously for some protein complexes although not rationalized. Here we combine experimental and theoretical studies to investigate the gas-phase structures of four multimeric protein complexes during collisional activation. Importantly, using ion mobility-mass spectrometry (IM-MS), we find that all four macromolecular complexes retain their native-like topologies at low energy. Upon increasing the collision energy, two of the four complexes adopt a more compact state. This collapse was most noticeable for pentameric serum amyloid P (SAP) which contains a large central cavity. The extent of collapse was found to be highly correlated with charge state, with the surprising observation that the lowest charge states were those which experience the greatest degree of compaction. We compared these experimental results with in vacuo molecular dynamics (MD) simulations of SAP, during which the temperature was increased. Simulations showed that low charge states of SAP exhibited compact states, corresponding to collapse of the ring, while intermediate and high charge states unfolded to more extended structures, maintaining their ring-like topology, as observed experimentally. To simulate the collision-induced dissociation (CID) of different charge states of SAP, we used MS to measure the charge state of the ejected monomer and assigned this charge to one subunit, distributing the residual charges evenly among the remaining four subunits. Under these conditions, MD simulations captured the unfolding and ejection of a single subunit for intermediate charge states of SAP. The highest charge states recapitulated the ejection of compact monomers and dimers, which we observed in CID experiments of high charge states of SAP, accessed by supercharging. This strong correlation between theory and experiment has implications for further studies as well as for understanding the process of CID and for applications to gas-phase structural biology more generally.

Defect‐Rich Ultrathin ZnAl‐Layered Double Hydroxide Nanosheets for Efficient Photoreduction of CO2 to CO with Water

Defect-rich ultrathin ZnAl-layered double hydroxide nanosheets are successfully prepared. Under UV-vis irradiation, these nanosheets are superior efficient catalysts for the photoreduction of CO2 to CO with water. The formed oxygen vacancies lead to the formation of coordinatively unsaturated Zn(+) centers within the nanosheets, responsible for the very high photocatalytic activities.

A refined solution structure of hen lysozyme determined using residual dipolar coupling data.

A high resolution NMR structure of hen lysozyme has been determined using 209 residual 1H–15N dipolar coupling restraints from measurements made in two different dilute liquid crystalline phases (bicelles) in conjunction with a data set of 1632 NOE distance restraints, 110 torsion angle restraints, and 60 hydrogen bond restraints. The ensemble of 50 low‐energy calculated structures has an average backbone RMSD of 0.50±0.13Å to the mean structure and of 1.49±0.10Å to the crystal structure of hen lysozyme. To assess the importance of the dipolar coupling data in the structure determination, the final structures are compared with an ensemble calculated using an identical protocol but excluding the dipolar coupling restraints. The comparison shows that structures calculated with the dipolar coupling data are more similar to the crystal structure than those calculated without, and have better stereochemical quality. The structures also show improved quality factors when compared with additional dipolar coupling data that were not included in the structure calculations, with orientation‐dependent 15N chemical shift changes measured in the bicelle solutions, and with T1/T2 values obtained from 15N relaxation measurements. Analysis of the ensemble of NMR structures and comparisons with crystal structures, 15N relaxation data, and molecular dynamics simulations of hen lysozyme provides a detailed description of the solution structure of this protein and insights into its dynamical behavior.

Characterisation of protein unfolding by NMR diffusion measurements

The characterisation of non-native states of proteins is a key problem instudies of protein folding. Complete characterisation of these states requiresa description of both local and global properties, including moleculardimensions. Here we present results from pulsed field gradient experimentsdesigned to compare the effective hydrodynamic radii of a protein in nativeand non-native states. Measurements performed on lysozyme indicate that theeffective hydrodynamic radius increases by 38±1% on unfolding in urea,a result completely consistent with a recent study by small-angle X-rayscattering.

Amyloid fibril formation by a helical cytochrome

The substitution of alanines for the two cysteines which form thioether linkages to the haem group in cytochrome c 552 from Hydogenobacter thermophilus destabilises the native protein fold. The holo form of this variant slowly converts into a partially folded apo state that over prolonged periods of time aggregates into fibrillar structures. Characterisation of these structures by electron microscopy and thioflavin‐T binding assays shows that they are amyloid fibrils. The data demonstrate that when the native state of this cytochrome is destabilised by loss of haem, even this highly α‐helical protein can form β‐sheet structures of the type most commonly associated with protein deposition diseases.

Heme proteins—Diversity in structural characteristics, function, and folding

The characteristics of heme prosthetic groups and their binding sites have been analyzed in detail in a data set of nonhomologous heme proteins. Variations in the shape, volume, and chemical composition of the binding site, in the mode of heme binding and in the number and nature of heme–protein interactions are found to result in significantly different heme environments in proteins with different functions in biology. Differences are also seen in the properties of the apo states of the proteins. The apo states of proteins that bind heme permanently in their functional form show some disorder, ranging from local unfolding in the heme binding pocket to complete unfolding to give a random coil. In contrast, proteins that bind heme transiently are fully folded in their apo and holo states, presumably allowing both apo and holo forms to remain biologically active resisting aggregation or proteolysis. The principles identified here provide a framework for the design of de novo proteins that will exhibit tight heme ligand binding and for the identification of the function of structural genomic target proteins with heme ligands. Proteins 2010.

Assessing equilibration and convergence in biomolecular simulations

If molecular dynamics simulations are used to characterize the folding of peptides or proteins, a wide range of conformational states needs to be sampled. This study reports an analysis of peptide simulations to identify the best methods for assessing equilibration and sampling in these systems where there is significant conformational disorder. Four trajectories of a β peptide in methanol and four trajectories of an α peptide in water, each of 5 ns in length, have been studied. Comparisons have also been made with two 50‐ns trajectories of the β peptide in methanol. The convergence rates of quantities that probe both the extent of conformational sampling and the local dynamical properties have been characterized. These include the numbers of hydrogen bonds populated, clusters identified, and main chain torsion angle transitions in the trajectories. The relative equilibrium rates of different quantities are found to vary significantly between the two systems studied reflecting both the differences in peptide primary structure and the different solvents used. A cluster analysis of the simulation trajectories is identified as a very effective method for judging the convergence of the simulations. This is particularly the case if the analysis includes a comparison of multiple trajectories calculated for the same system from different starting structures. Proteins 2002;48:487–496.

Stimulation and inhibition of fibril formation by a peptide in the presence of different concentrations of SDS

Sodium dodecyl sulphate (SDS), a detergent that mimics some characteristics of biological membranes, has been found to affect significantly fibril formation by a peptide from human complement receptor 1. In aqueous solution the peptide is unfolded but slowly aggregates to form fibrils. In sub‐micellar concentrations of SDS the peptide is initially α‐helical but converts rapidly to a β‐sheet structure and large quantities of fibrils form. In SDS above the critical micellar concentration the peptide adopts a stable α‐helical structure and no fibrils are observed. These findings demonstrate the sensitivity of fibril formation to solution conditions and suggest a possible role for membrane components in amyloid fibril formation in living systems.

Entropy calculations on the molten globule state of a protein: Side-chain entropies of alpha-lactalbumin

We present entropy estimates based on molecular dynamics simulations of models of the molten globule state of the protein α‐lactalbumin at low pH. The entropy calculations use the covariance matrix of atom‐positional fluctuations and yield the complete configurational entropy. The configurational entropy of the entire protein and of each of its side chains is calculated. Exposed side chains show a larger entropy compared to buried side chains. A comparison to data from rotamer counting is made and significant differences are found. Proteins 2002;46:215–224.