Search PubMed⌕ Search

Biomedical subjects

Lorna J Smith

Publications and source records attributed to Lorna J Smith.

12 recordsLinked to original sources

Molecular dynamics simulations of Hydrogenobacter thermophilus cytochrome c552: comparisons of the wild-type protein, a b-type variant, and the apo state.

Molecular dynamic simulations have been performed for wild-type Hydrogenobacter thermophilus cytochrome c(552), a b-type variant of the protein, and the apo state with the heme prosthetic group removed. In the b-type variant, Cys 10 and Cys 13 were mutated to alanine residues, and so the heme group was no longer covalently bound to the protein. Two 8-ns simulations have been performed for each system at 298 and 360 K. The simulations of the wild-type protein at 298 K show a very close agreement with experimental NMR data. A fluxional process involving the side chain of Met 59, which coordinates to the heme iron, is observed in accord with proposals from NMR studies. Overall, the structure and dynamical behavior of the protein during the simulations of the b-type variant is closely similar to that of the wild-type protein. However, side chains in the heme-binding site show larger fluctuations in the b-type variant simulation at 360 K. In addition, structural changes are seen for a number of residues close to the heme group, particularly Gly 22 and Ser 51. The simulations of the apo state show significant conformational changes for residues 50-59. These residues form a loop region, which packs over the heme group in the wild-type protein and hydrogen bonds to the heme propionate groups. In the absence of heme, in the apo state simulations, these residues form short but persistent regions of beta-sheet secondary structure. These could provide nucleation sites for the conversion to amyloid fibrils.

Alanine↗

Oxidative refolding of amyloidogenic variants of human lysozyme.

The oxidative refolding of human lysozyme and its two best characterised amyloidogenic variants, Ile56Thr and Asp67His, has been investigated in vitro by means of the concerted application of a range of biophysical techniques. The results show that in each case the ensemble of reduced denatured conformers initially collapses into a large number of unstructured intermediates with one or two disulphide bonds, the majority of which then fold to form the native-like three-disulphide intermediate, des-[77-95]. The slow step in the overall folding reaction involves the rearrangement of the latter to the fully oxidised native protein containing four disulphide bonds. The Ile56Thr and Asp67His variants were found to fold faster than the wild-type protein by a factor of 2 and 3 respectively, an observation that can be attributed primarily to the reduction in the barriers to conformational rearrangements that results from both the mutations. The efficient folding of these variants despite their enhanced propensities to aggregate when compared to the wild-type protein is consistent with their ability to be secreted in sufficient quantities to give rise to the systemic amyloidoses with which they are associated.

Amyloid↗

Characterization of the denaturation of human alpha-lactalbumin in urea by molecular dynamics simulations.

Molecular dynamics (MD) simulations were used to characterize the non-cooperative denaturation of the molten globule A-state of human alpha-lactalbumin by urea. A solvent of explicit urea and water molecules was used, corresponding to a urea concentration of approximately 6M. Three simulations were performed at temperatures of 293K, 360K and 400K, with lengths of 2 ns, 8 ns and 8 ns respectively. The results of the simulations were compared with experimental data from NMR studies of human alpha-lactalbumin and related peptides. During the simulations, hydrogen bonds were formed from the protein to both urea and water molecules as intra-protein hydrogen bonds were lost. Urea was shown to compete efficiently with water as both a hydrogen bond donor and acceptor. Radial distribution functions of water and urea around hydrophobic side chain atoms showed a significant increase in urea molecules in the solvation shell as the side chains became exposed during denaturation. A considerable portion of the native-like secondary structure persisted throughout the simulations. However, in the simulations at 360K and 400K, there were substantial changes in the packing of aromatic and other hydrophobic side chains in the protein, and many native contacts were lost. The results suggest that during the non-cooperative denaturation of the molten globule, secondary structure elements are stabilized by non-specific, non-native interactions.

Circular Dichroism↗

Photo-CIDNP NMR spectroscopy of a heme-containing protein.

There are relatively few examples of the application of photo-CIDNP NMR spectroscopy to chromophore-containing proteins. The most likely reason for this is that simultaneous absorption of light by the photosensitiser molecule and the protein chromophore reduces the effectiveness of the photochemical reaction that produces the observed nuclear polarisation. We present details of experiments performed on the air-oxidised form of a small cytochrome, from the thermophilic bacterium Hydrogenobacter thermophilus, using both the wild-type protein and apo and holo forms of a double alanine b-type mutant. We show that, along with the apo state, it is possible to generate CIDNP in the air-oxidised form of the b-type mutant, but not in the corresponding c-type cytochrome. This finding is supported by control experiments using horse-heart cytochrome c.

Animals↗

Characterization of the molten globule state of retinol-binding protein using a molecular dynamics simulation approach.

Retinol-binding protein transports retinol, and circulates in the plasma as a macromolecular complex with the protein transthyretin. Under acidic conditions retinol-binding protein undergoes a transition to the molten globule state, and releases the bound retinol ligand. A biased molecular dynamics simulation method has been used to generate models for the ensemble of conformers populated within this molten globule state. Simulation conformers, with a radius of gyration at least 1.1 A greater than that of the native state, contain on average 37%beta-sheet secondary structure. In these conformers the central regions of the two orthogonal beta-sheets that make up the beta-barrel in the native protein are highly persistent. However, there are sizable fluctuations for residues in the outer regions of the beta-sheets, and large variations in side chain packing even in the protein core. Significant conformational changes are seen in the simulation conformers for residues 85-104 (beta-strands E and F and the E-F loop). These changes give an opening of the retinol-binding site. Comparisons with experimental data suggest that the unfolding in this region may provide a mechanism by which the complex of retinol-binding protein and transthyretin dissociates, and retinol is released at the cell surface.

Computer Simulation↗

NMR analysis shows that a b-type variant of Hydrogenobacter thermophilus cytochrome c552 retains its native structure.

Conversion of Hydrogenobacter thermophilus cytochrome c(552) into a b-type cytochrome by mutagenesis of both heme-binding cysteines to alanines significantly reduces the stability of the protein (Tomlinson, E. J., and Ferguson, S. J. (2000) Proc. Natl. Acad. Sci. U. S. A. 97, 5156-5160). To understand the effects of this change on the structure and dynamics of the protein, hetero-nuclear (15)N-edited NMR techniques have been used to characterize this b-type variant. The backbone (15)N, (1)H(N), and (1)H(alpha), and (1)H(beta) resonances of the protein have been assigned. Analysis of (3)J(HN)alpha coupling constants, nuclear Overhauser enhancement intensities, and chemical shift index data demonstrates that the four alpha-helices present in the wild-type protein are retained in the b-type variant. Comparison of the chemical shifts for the b-type and wild-type proteins indicates that the tertiary structures of the two proteins are closely similar. Some subtle differences are, however, observed for residues in the N-terminal region and in the vicinity of the heme-binding pocket. Hydrogen exchange studies show that there are 25 backbone amide protons that exchange very slowly in the b-type variant and confirm that the fluctuations within the b-type protein are of a similar extent to those in the wild-type protein. These data demonstrate the notable retention of the native secondary structure and tertiary fold despite the absence of covalent linkages between the heme group and the protein.

Alanine↗

Asparagine and glutamine side-chain conformation in solution and crystal: a comparison for hen egg-white lysozyme using residual dipolar couplings.

Experimental (15)N-(1)H and (1)H-(1)H residual dipolar couplings (RDCs) for the asparagine (Asn) and glutamine (Gln) side chains of hen egg-white lysozyme are measured and analysed in conjunction with (1)N relaxation data, information about chi(1) torsion angles in solution and molecular dynamics simulations. The RDCs are compared to values predicted from 16 high-resolution crystal structures. Two distinct groups of Asn and Gln side chains are identified. The first contains residues whose side chains show a fixed, relatively rigid, conformation in solution. For these residues there is good agreement between the experimental and predicted RDCs. This agreement improves when the experimental order parameter, S, is included in the calculation of the RDCs from the crystal structures. The comparison of the experimental RDCs with values calculated from the X-ray structures shows that the similarity between the oxygen and nitrogen electron densities is a limitation to the correct assignment of the Asn and Gln side-chain orientation in X-ray structures. In the majority of X-ray structures a 180 degrees rotation about chi(2) or chi(3), leading to the swapping of N(delta/epsilon 2) and O(delta/epsilon 1), is necessary for at least one Asn or Gln residue in order to achieve good agreement between experimental and predicted RDCs. The second group contains residues whose side chains do not adopt a single, well-defined, conformation in solution. These residues do not show a correlation between the experimental and predicted RDCs. In many cases the family of crystal structures shows a range of orientations for these side chains, but in others the crystal structures show a well-defined side-chain position. In the latter case, this is found to arise from crystallographic contacts and does not represent the behaviour of the side chain in solution.

Animals↗

Computational methods for generating models of denatured and partially folded proteins.

Partially folded and denatured proteins can give important insights into protein folding, misfolding, and aggregation. Such non-native states of proteins are however very difficult to characterise in detail as they are dynamic, heterogeneous systems comprising of ensembles of interconverting conformers. This article describes methods that produce models for non-native proteins in atomic detail. A variety of molecular dynamics based protocols are discussed together with some recent procedures that include restraints from experimental data. These models provide an important framework for interpreting experimental data from studies of non-native states using nuclear magnetic resonance spectroscopy, fluorescence, circular dichroism, and small angle scattering techniques.

Animals↗

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 alpha-helical but converts rapidly to a beta-sheet structure and large quantities of fibrils form. In SDS above the critical micellar concentration the peptide adopts a stable alpha-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.

Circular Dichroism↗

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 beta peptide in methanol and four trajectories of an alpha peptide in water, each of 5 ns in length, have been studied. Comparisons have also been made with two 50-ns trajectories of the beta 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.

Amino Acid Sequence↗

Long-range interactions within a nonnative protein.

Protein folding and unfolding are coupled to a range of biological phenomena, from the regulation of cellular activity to the onset of neurodegenerative diseases. Defining the nature of the conformations sampled in nonnative proteins is crucial for understanding the origins of such phenomena. We have used a combination of nuclear magnetic resonance (NMR) spectroscopy and site-directed mutagenesis to study unfolded states of the protein lysozyme. Extensive clusters of hydrophobic structure exist within the wild-type protein even under strongly denaturing conditions. These clusters involve distinct regions of the sequence but are all disrupted by a single point mutation that replaced residue Trp62 with Gly located at the interface of the two major structural domains in the native state. Thus, nativelike structure in the denatured protein is stabilized by the involvement of Trp62 in nonnative and long-range interactions.

Amino Acid Substitution↗

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 alpha-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.

Entropy↗