PubMed Health⌕ Search

Biomedical subjects

Timothy A Cross

Publications and source records attributed to Timothy A Cross.

10 recordsLinked to original sources

Molecular dynamics simulations of Trp side-chain conformational flexibility in the gramicidin A channel.

Gramicidin A (gA) is prototypical peptide antibiotic and a model ion channel former. Configured in the solid-state NMR beta(6.5)-helix channel conformation, gA was subjected to 1-ns molecular dynamics (MD) gas phase simulations using the all-atom charmm22 force field to ascertain the conformational stability of the Trp side chains as governed by backbone and neighboring side-chain contacts. Three microcanonical trajectories were computed using different initial atomic velocities for each of twenty different initial structures. For each set, one of the four Trp side chains in each monomer was initially positioned in one of the five non-native conformations (A. E. Dorigo et al., Biophysical Journal, 1999, Vol. 76, 1897-1908), the other Trps being positioned in the native state, o1. In three additional control simulations, all Trps were initiated in the native conformation. After equilibration, constraints were removed and subsequent conformational changes of the initially constrained Trp were measured. The chi(1) was more flexible than chi(2.1). The energetically optimal orientation, o1 (Dorigo et al., 1999), was the most stable in all four Trp positions (9, 11, 13, 15) and remained unchanged for the entire 1 ns simulation in 19 of 24 trials. Changes in chi(1) from each of the 5 suboptimal states occur readily. Two of the non-native conformations reverted readily to o1, whereas the other three converted to an intermediate state, i2. There were frequent interconversions between i2 and o1. We speculate that experimentally observed Trp stability is caused by interactions with the lipid-water interface, and that stabilization of one of the suboptimal conformations in gA, such as i2, by lipid headgroups could produce a secondary, metastable conformational state. This could explain recent experimental studies of differences in the channel conductance dispersity between gA and a Trp-to-Phe gA analog, gramicidin M (gM, J. C. Markham et al., Biochimica et Biophysica Acta, 2001, Vol. 1513, 185-192).

Dimerization↗

Atomic refinement with correlated solid-state NMR restraints.

The orientation data provided by solid-state NMR can provide a great deal of structural information about membrane proteins. The quality of the information provided is, however, somewhat degraded by sign degeneracies in measurements of the dipolar coupling tensor. This is reflected in the dipolar coupling penalty function used in atomic refinement, which is less capable of properly restraining atoms when dipolar sign degeneracies are present. In this report we generate simulated solid-state NMR data using a variety of procedures, including back-calculation from crystal structures of alpha-helical and beta-sheet membrane proteins. We demonstrate that a large fraction of the dipolar sign degeneracies are resolved if anisotropic dipolar coupling measurements are correlated with anisotropic chemical shift measurements, and that all sign degeneracies can be resolved if three data types are correlated. The advantages of correlating data are demonstrated with atomic refinement of two test membrane proteins. When refinement is performed using correlated dipolar couplings and chemical shifts, perturbed structures converge to conformations with a larger fraction of correct dipolar signs than when data are uncorrelated. In addition, the final structures are closer to the original unperturbed structures when correlated data are used in the refinement. Thus, refinement with correlated data leads to improved atomic structures. The software used to correlate dipolar coupling and chemical shift data and to set up energy functions and their derivatives for refinement, CNS-SS02, is available at our web site.

Algorithms↗

Initial structural and dynamic characterization of the M2 protein transmembrane and amphipathic helices in lipid bilayers.

Amphipathic helices in membrane proteins that interact with the hydrophobic/hydrophilic interface of the lipid bilayer have been difficult to structurally characterize. Here, the backbone structure and orientation of an amphipathic helix in the full-length M2 protein from influenza A virus has been characterized. The protein has been studied in hydrated DMPC/DMPG lipid bilayers above the gel to liquid-crystalline phase transition temperature by solid-state NMR spectroscopy. Characteristic PISA (Polar Index Slant Angle) wheels reflecting helical wheels have been observed in uniformly aligned bilayer preparations of both uniformly 15N labeled and amino acid specific labeled M2 samples. Hydrogen/deuterium exchange studies have shown the very slow exchange of some residues in the amphipathic helix and more rapid exchange for the transmembrane helix. These latter results clearly suggest the presence of an aqueous pore. A variation in exchange rate about the transmembrane helical axis provides additional support for this claim and suggests that motions occur about the helical axes in this tetramer to expose the entire backbone to the pore.

Influenza A virus↗

NMR spin locking of proton magnetization under a frequency-switched Lee-Goldburg pulse sequence.

The spin dynamics of NMR spin locking of proton magnetization under a frequency-switched Lee-Goldburg (FSLG) pulse sequence is investigated for a better understanding of the line-narrowing mechanism in PISEMA experiments. For the sample of oriented 15N(1,3,5,7)-labeled gramicidin A in hydrated DMPC bilayers, it is found that the spin-lattice relaxation time T(1rho)(H) in the tilted rotating frame is about five times shorter when the 1H magnetization is spin locked at the magic angle by the FSLG sequence compared to the simple Lee-Goldburg sequence. It is believed that the rapid phase alternation of the effective fields during the FSLG cycles results in averaging of the spin lock field so that the spin lock becomes less efficient. A FSLG supercycle has been suggested here to slow the phase alternation. It has been demonstrated experimentally that a modified PISEMA pulse sequence with such supercycles gives rise to about 30% line narrowing in the dipolar dimension in the PISEMA spectra compared to a standard PISEMA pulse sequence.

Gramicidin↗

Noncontact dipole effects on channel permeation. VI. 5F- and 6F-Trp gramicidin channel currents.

Fluorination of peptide side chains has been shown to perturb gramicidin channel conductance without significantly changing the average side chain structure, which, it is hoped, will allow detailed analysis of electrostatic modulation of current flow. Here we report a 1312-point potassium current-voltage-concentration data set for homodimeric channels formed from gramicidin A (gA) or any of eight fluorinated Trp analogs in both lecithin and monoglyceride bilayers. We fit the data with a three-barrier, two-site, two-ion (3B2S) kinetic model. The fluorination-induced changes in the rate constants were constrained by the same factor in both lipids. The rate constant changes were converted to transition-state free-energy differences for comparison with previous electrostatic potential energy differences based on an ab initio force field. The model allowed a reasonably good fit (chi = 8.29 with 1271 degrees of freedom). The measured changes were subtle. Nevertheless, the fitted energy perturbations agree well with electrostatic predictions for five of the eight peptides. For the other three analogs, the fitted changes suggested a reduced translocation barrier rather than the reduced exit barrier as predicted by electrostatics.

Algorithms↗

Uniformity, ideality, and hydrogen bonds in transmembrane alpha-helices.

Protein environments substantially influence the balance of molecular interactions that generate structural stability. Transmembrane helices exist in the relatively uniform low dielectric interstices of the lipid bilayer, largely devoid of water and with a very hydrophobic distribution of amino acid residues. Here, through an analysis of bacteriorhodopsin crystal structures and the transmembrane helix structure from M2 protein of influenza A, some helices are shown to be exceptionally uniform in hydrogen bond geometry, peptide plane tilt angle, and backbone torsion angles. Evidence from both the x-ray crystal structures and solid-state NMR structure suggests that the intramolecular backbone hydrogen bonds are shorter than their counterparts in water-soluble proteins. Moreover, the geometry is consistent with a dominance of electrostatic versus covalent contributions to these bonds. A comparison of structure as a function of resolution shows that as the structures become better characterized the helices become much more uniform, suggesting that there is a possibility that many more uniform helices will be observed, even among the moderate resolution membrane protein structures that are currently in the Protein Data Bank that do not show such features.

Bacteriorhodopsins↗

The role of Trp side chains in tuning single proton conduction through gramicidin channels.

We present an extensive set of measurements of proton conduction through gramicidin A (gA), B (gB), and M (gM) homodimer channels which have 4, 3, or 0 Trp residues at each end of the channel, respectively. In gA we find a shoulder separating two domains of conductance increasing with concentration, confirming the results of Eisenman, G., B. Enos, J. Hagglund, and J. Sandblom. 1980. Ann. NY. Acad. Sci. 339:8-20. In gB, the shoulder is shifted by approximately 1/2 pH unit to higher H(+) concentrations and is very sharply defined. No shoulder appears in the gM data, but an associated transition from sublinear to superlinear I-V values occurs at a 100-fold higher [H(+)] in gM than in gA. The data in the low concentration domain are analyzed using a configuration space model of single-proton conduction, assuming that the difference in the proton potential of mean force (PMF) between gA and its analogs is constant, similar to the results of Anderson, D., R. B. Shirts, T. A. Cross, and D. D. Busath. 2001. Biophys. J. 81:1255-1264. Our results suggest that the average amplitudes of the calculated proton PMFs are nearly correct, but that the water reorientation barrier calculated for gA by molecular dynamics using the PM6 water model (Pomès, R., and B. Roux. 1997. Biophys. J. 72:246a) must be reduced in amplitude by 1.5 kcal/mol or more, and is not rate-limiting for gA.

Biophysical Phenomena↗

Solid-state 19F-NMR analysis of 19F-labeled tryptophan in gramicidin A in oriented membranes.

The response of membrane-associated peptides toward the lipid environment or other binding partners can be monitored by solid-state NMR of suitably labeled side chains. Tryptophan is a prominent amino acid in transmembrane helices, and its (19)F-labeled analogues are generally biocompatible and cause little structural perturbation. Hence, we use 5F-Trp as a highly sensitive NMR probe to monitor the conformation and dynamics of the indole ring. To establish this (19)F-NMR strategy, gramicidin A was labeled with 5F-Trp in position 13 or 15, whose chi(1)/chi(2) torsion angles are known from previous (2)H-NMR studies. First, the alignment of the (19)F chemical shift anisotropy tensor within the membrane was deduced by lineshape analysis of oriented samples. Next, the three principal axes of the (19)F chemical shift anisotropy tensor were assigned within the molecular frame of the indole ring. Finally, determination of chi(1)/chi(2) for 5F-Trp in the lipid gel phase showed that the side chain alignment differs by up to 20 degrees from its known conformation in the liquid crystalline state. The sensitivity gain of (19)F-NMR and the reduction in the amount of material was at least 10-fold compared with previous (2)H-NMR studies on the same system and 100-fold compared with (15)N-NMR.

Cell Membrane↗

The effect of Hartmann-Hahn mismatching on polarization inversion spin exchange at the magic angle.

The effect of the Hartmann-Hahn mismatch delta = omega(eff)-omega(1S) during polarization inversion spin exchange at the magic angle (PISEMA) has been investigated, where omega(eff) and omega(1S) represent the amplitudes of the 1H effective spin-locking field at the magic angle and the 15N RF spin-locking field, respectively. During the PISEMA evolution period, the exact Hartmann-Hahn match condition (i.e., delta = 0) yields a maximum dipolar scaling factor of 0.816 for PISEMA experiments, while any mismatch results in two different effective fields for the first and second half of each frequency switched Lee-Goldburg (FSLG) cycle. The mismatch effect on the scaling factor depends strongly on the transition angle from one effective field to the other within each FSLG cycle as well as on the cycle time. At low RF spin-lock amplitudes in which the FSLG cycle time is relatively long, the scaling factor rapidly becomes smaller as omega(1S) becomes greater than omega(eff). On the other hand, when omega(1S) < omega(eff), there is relatively little effect on the scaling factor with variation in delta. As a result, the presence of RF inhomogeneities may significantly broaden the line-width in the dipolar dimension because of the mismatch effect. Higher RF spin-lock amplitudes result in a relatively small variation for the scaling factor. Furthermore, ramped amplitude of the 15N RF spin-lock field in synchronization with the flip-flop of the FSLG sequence minimizes the transition angle between the two effective fields within the FSLG cycle. It is shown experimentally that such a ramped amplitude not only gives rise to the same scaling factor but also results in a narrower dipolar line-width in comparison with the rectangular amplitude.

Gramicidin↗