The upper and lower lip regions in lysozyme from bacteriophage lambda (λ-lysozyme) are flexible in solution and exhibit two different conformations in crystal structures of the protein. MD simulations have been used to characterize the structure and dynamics of these lip regions, which surround the active site. Ten different simulations have been run including those with restraining to experimental NOE distance and 1H-15N order parameter data. The simulations show that the lower lip region, although undergoing considerable backbone fluctuations, contains two persistent β-strands. In the upper lip region, a wide range of conformations are populated and it is not clear from the available data whether some helical secondary structure is present. The work provides a clear example of the advantages of combining MD simulations with experimental data to obtain a structural interpretation of the latter. In this case, time-averaged order parameter restraining has played an essential role in enabling convergence between two different starting structures and identifying the extent to which flexible regions in solution can contain persistent secondary structure.

Experimental studies of barley and maize lipid transfer proteins (LTPs) show that the two proteins bind the ligand palmitate in opposite orientations in their internal cavities. Moreover, maize LTP is reported to bind the ligand caprate in the internal cavity in a mixture of two orientations with approximately equal occupancy. Six 30 ns molecular dynamics (MD) simulations of maize and barley LTP with ligands bound in two orientations (modes M and B) have been used to understand the different ligand binding preferences. The simulations show that both maize and barley LTP could bind palmitate in the orientation observed experimentally for maize LTP (mode M), with the predominant interaction being a salt bridge between the ligand carboxylate headgroup and a conserved arginine side chain. However, the simulation of barley LTP with palmitate in the mode B orientation shows the most favorable protein–ligand interaction energy. In contrast, the simulations of maize LTP with palmitate and with caprate in the mode B orientation show no persistent ligand binding, the ligands leaving the cavity during the simulations. Sequence differences between maize and barley LTP in the AB loop region, in residues at the base of the hydrophobic cavity, and in the helix A region are identified as contributing to the different behavior. The simulations reproduce well the experimentally observed binding preferences for palmitate and suggest that the experimental data for maize LTP with caprate reflect ligand mobility in binding mode M rather than the population of binding modes M and B.

The structure and stability of the allergenic nonspecific lipid transfer protein (LTP) of peach were compared with the homologous LTP1 of barley and its liganded form LTP1b. All three proteins were resistant to gastric pepsinolysis and were only slowly digested at 1 to 2 out of 14 potential tryptic and chymotryptic cleavage sites under duodenal conditions. Peach LTP was initially cleaved at Tyr79-Lys80 and then at Arg39-Thr40 (a site lost in barley LTP1). Molecular dynamics simulations of the proteins under folded conditions showed that the backbone flexibility is limited, explaining the resistance to duodenal proteolysis. Arg39 and Lys80 side chains were more flexible in simulations of peach compared with barley LTP1. This may explain differences in the rates of cleavage observed experimentally for the two proteins and suggests that the flexibility of individual amino acid side chains could be important in determining preferred proteolytic cleavage sites. In order to understand resistance to pepsinolysis, proteins were characterized by NMR spectroscopy at pH 1.8. This showed that the helical regions of both proteins remain folded at this pH. NMR hydrogen exchange studies confirmed the rigidity of the structures at acidic pH, with barley LTP1 showing some regions with greater protection. Collectively, these data suggest that the rigidity of the LTP scaffold is responsible for their resistance to proteolysis. Gastroduodenal digestion conditions do not disrupt the 3D structure of peach LTP, explaining why LTPs retain their ability to bind IgE after digestion and hence their allergenic potential.

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.

NMR techniques have been used to characterise the effects of a lipid-like post-translational modification on barley lipid transfer protein (LTP1b). NMR chemical shift data indicate that the lipid-like molecule lies in the hydrophobic cavity of LTP1b, with Tyr 79 being displaced to accommodate the ligand in the cavity. The modified protein has a reduced level of backbone amide hydrogen exchange protection, presumably reflecting increased dynamics in the protein. This may result from a loosening of the protein structure and may explain the enhanced surface properties observed for LTP1b.