The effects of substrate and vertebral number on locomotion in the garter snake Thamnophis elegans

1. Locomotor performance of limbless vertebrates depends on the substrate through which individuals move and may result in selection on vertebral number in different habitats. To evaluate the effect of push-point density on snake locomotion, the density of vegetation and other potential push-points was quantified at two sites in California (coastal and inland), where conspecific snakes differed greatly in vertebral number (230 and 256 average total vertebrae, respectively; Arnold 1988). The coastal site had significantly higher push-point densities than the inland site.
2. Five experimental push-point densities that fell within the natural range of push-point densities were employed in laboratory trials of juvenile snake locomotion. Density of push-points significantly affected both crawling speed and head-to-tail distance (HTD), an indirect measure of lateral bending. The fastest speed was achieved at an intermediate push-point density. The shortest HTD occurred when snakes moved through the lowest push-point density.
3. Sex, total number of vertebrae and total length significantly affected HTD, regardless of push-point density. Snakes with relatively more vertebrae had a shorter HTD, suggesting they were able to achieve greater lateral bending than snakes with fewer vertebrae. Coastal and inland populations did not differ in HTD during locomotion.
4. Numbers of body and tail vertebrae significantly influenced speed at different push-point densities. In general, snakes with more body vertebrae were slower than those with fewer, while snakes with more tail vertebrae were faster than those with fewer. Snakes of greater total length were faster at all densities. Coastal snakes crawled faster than inland snakes at all push-point densities.     

Functional conformational changes of endo-1,4-xylanase II from Trichoderma reesei: A molecular dynamics study

Recent crystallographic studies have revealed a range of structural changes in the three-dimensional structure of endo-1,4-xylanase (XYNII) from Trichoderma reesei. The observed conformational changes can be described as snapshots of an open-close movement of the active site of XYNII. These structures were further analyzed in this study. In addition, a total of four 1 ns molecular dynamics (MD) simulations were performed representing different states of the enzyme. A comparison of the global and local changes found in the X-ray structures and the MD runs suggested that the simulations reproduced a similar kind of active site opening and closing as predicted by the crystal structures. The open-close movement was characterized by the use of distance difference matrixes and the Hingefind program (Wriggers and Schulten, Proteins 29:1–14, 1997) to be a ‘hinge-bending’ motion involving two large rigidly-moving regions and an extended hinge. This conformational feature is probably inherent to this molecular architecture and probably plays a role in the function of XYNII. Proteins 31:434–444, 1998. © 1998 Wiley-Liss, Inc.     

Systematic analysis of domain motions in proteins from conformational change: New results on citrate synthase and T4 lysozyme

Methods developed originally to analyze domain motions from simulation [Proteins 27:425–437, 1997] are adapted and extended for the analysis of X-ray conformers and for proteins with more than two domains. The method can be applied as an automatic procedure to any case where more than one conformation is available. The basis of the methodology is that domains can be recognized from the difference in the parameters governing their quasi-rigid body motion, and in particular their rotation vectors. A clustering algorithm is used to determine clusters of rotation vectors corresponding to main-chain segments that form possible dynamic domains. Domains are accepted for further analysis on the basis of a ratio of interdomain to intradomain fluctuation, and Chasles' theorem is used to determine interdomain screw axes. Finally residues involved in the interdomain motion are identified. The methodology is tested on citrate synthase and the M6I mutant of T4 lysozyme. In both cases new aspects to their conformational change are revealed, as are individual residues intimately involved in their dynamics. For citrate synthase the beta sheet is identified to be part of the hinging mechanism. In the case of T4 lysozyme, one of the four transitions in the pathway from the closed to the open conformation, furnished four dynamic domains rather than the expected two. This result indicates that the number of dynamic domains a protein possesses may not be a constant of the motion. Proteins 30:144–154, 1998. © 1998 Wiley-Liss, Inc.     

Membrane-bound structure and energetics of alpha-synuclein

Aggregation and fibrillation of alpha-synuclein bound to membranes are believed to be involved in Parkinson's and other neurodegenerative diseases. On SDS micelles, the N-terminus of alpha-synuclein forms two curved helices linked by a short loop. However, its structure on lipid bilayers has not been experimentally resolved. Using MD simulations with an implicit membrane model we show here that, on a planar mixed membrane, the truncated alpha-synuclein (residues 1-95) forms a bent helix. Bending of the helix is not due to the protein sequence or membrane binding, but to collective motions of the long helix. The backbone of the helix is approximately 2.5 A above the membrane surface, with some residues partially inserted in the membrane core. The helix periodicity is 11/3 (11 residues complete three full turns) as opposed to 18/5 periodicity of an ideal alpha-helix, with hydrophobic residues towards the membrane, negatively charged residues towards the solvent and lysines on the polar/nonpolar interface. A series of threonines, which are characteristic for alpha-synuclein and perhaps a phosphorylation site, is also located at the hydrophobic/hydrophilic interface with their side chain often hydrogen bonded to the main-chain atom. The calculations show that the energy penalty for change in periodicity from the 18/5 to 11/3 on the anionic membrane is overcome by favorable solvation energy. The binding of truncated alpha-synuclein to membranes is weak. It prefers anionic membranes but it also binds marginally to a neutral membrane, via its C-terminus. Dimerization of helical monomers on the mixed membrane is energetically favorable. However, it slightly interferes with membrane binding. This might promote lateral diffusion of the protein on the membrane surface and facilitate assembly of oligomers that precede fibrillation.     

HingeMaster: normal mode hinge prediction approach and integration of complementary predictors

Protein motion is often the link between structure and function and a substantial fraction of proteins move through a domain hinge bending mechanism. Predicting the location of the hinge from a single structure is thus a logical first step towards predicting motion. Here, we describe ways to predict the hinge location by grouping residues with correlated normal-mode motions. We benchmarked our normal-mode based predictor against a gold standard set of carefully annotated hinge locations taken from the Database of Macromolecular Motions. We then compared it with three existing structure-based hinge predictors (TLSMD, StoneHinge, and FlexOracle), plus HingeSeq, a sequence-based hinge predictor. Each of these methods predicts hinges using very different sources of information-normal modes, experimental thermal factors, bond constraint networks, energetics, and sequence, respectively. Thus it is logical that using these algorithms together would improve predictions. We integrated all the methods into a combined predictor using a weighted voting scheme. Finally, we encapsulated all our results in a web tool which can be used to run all the predictors on submitted proteins and visualize the results.     

Chromophore Structure of Cyanobacterial Phytochrome Cph1 in the Pr State: Reconciling Structural and Spectroscopic Data by QM/MM Calculations

A quantum mechanics (QM)/molecular mechanics (MM) hybrid method was applied to the Pr state of the cyanobacterial phytochrome Cph1 to calculate the Raman spectra of the bound PCB cofactor. Two QM/MM models were derived from the atomic coordinates of the crystal structure. The models differed in the protonation site of His²⁶⁰ in the chromophore-binding pocket such that either the δ-nitrogen (M-HSD) or the ɛ-nitrogen (M-HSE) carried a hydrogen. The optimized structures of the two models display small differences specifically in the orientation of His²⁶⁰ with respect to the PCB cofactor and the hydrogen bond network at the cofactor-binding site. For both models, the calculated Raman spectra of the cofactor reveal a good overall agreement with the experimental resonance Raman (RR) spectra obtained from Cph1 in the crystalline state and in solution, including Cph1 adducts with isotopically labeled PCB. However, a distinctly better reproduction of important details in the experimental spectra is provided by the M-HSD model, which therefore may represent an improved structure of the cofactor site. Thus, QM/MM calculations of chromoproteins may allow for refining crystal structure models in the chromophore-binding pocket guided by the comparison with experimental RR spectra. Analysis of the calculated and experimental spectra also allowed us to identify and assign the modes that sensitively respond to chromophore-protein interactions. The most pronounced effect was noted for the stretching mode of the methine bridge A-B adjacent to the covalent attachment site of PCB. Due a distinct narrowing of the A-B methine bridge bond angle, this mode undergoes a large frequency upshift as compared with the spectrum obtained by QM calculations for the chromophore in vacuo. This protein-induced distortion of the PCB geometry is the main origin of a previous erroneous interpretation of the RR spectra based on QM calculations of the isolated cofactor.Abbreviations: Agp1, phytochrome from Agrobacterium tumefaciens; α-CPC, α-subunit of C-phycocyanin; BV, biliverdin IXα; B3LYP, three-parameter exchange functional according to Becke, Lee, Yang, and Parr; DFT, density functional theory; DrBphP, phytochrome from Deinococcus radiodurans; GAF, domain found in cGMP-specific phosphodiesterases; MM, molecular mechanics; MD, molecular dynamics; N-H ip, N-H in-plane bending; PCB, phycocyanobilin; PED, potential energy distribution; phyA, plant phytochrome; Pr, Pfr, red- and far-red absorbing parent states of phytochrome; PΦB, phytochromobilin; QM, quantum mechanics; RMSD, root mean-square deviation; RR, resonance Raman     

Spalt major controls the development of the notum and of wing hinge primordia of the Drosophila melanogaster wing imaginal disc

The Drosophila wing and the dorsal thorax develop from primordia within the wing imaginal disc. Here we show that spalt major (salm) is expressed within the presumptive dorsal body wall primordium early in wing disc development to specify notum and wing hinge tissue. Upon ectopic salm expression, dorsally located second leg disc cells develop notum and wing hinge tissue instead of sternopleural tissue. Similarly, by salm over-expression within the wing disc, wing blade formation is suppressed and a mirror-image duplication of the notum and wing hinge is formed. In large dorsal clones, which lack salm and its neighboring paralogue spalt related (salr), the cells of the notum primordium do not grow; these dorsal cells are not specified as notum, hence no notum outgrowth develops. These results suggest that the zinc finger factors encoded by the salm/salr complex play important roles in defining cells of the early wing disc as dorsal body wall cells, which develop into a large dorsal body wall territory and form mesonotum and some wing hinge tissue, and in delimiting the wing primordium. We also find that salm activity is down-regulated by its own product and by that of the Pax gene eyegone.