Hinge-bending motion in citrate synthase arising from normal mode calculations

A normal mode analysis of the closed form of dimeric citrate synthase has been performed. The largest-amplitude collective motion predicted by this method compares well with the crystallographically observed hinge-bending motion. Such a result supports those obtained previously in the case of hinge-bending motions of smaller systems, such as lysozyme or hexokinase. Taken together, all these results suggest that low-frequency normal modes may become useful for determining a first approximation of the conformational path between the closed and open forms of these proteins.     

Dynamic structures of granulocyte colony-stimulating factor proteins studied by normal mode analysis: two domain-type motions in low frequency modes

In this study, granulocyte colony-stimulating factor (GCSF) proteins were chosen as subjects for normal mode analysis. As helical cytokines with a four helix bundled type topology, they were classified into long chain and short chain groups by Sprang and Bazan. Normal mode calculations were also carried out with leukemia inhibitory factor (LIF), interleukin-6 (IL-6), and growth hormone (GH) as members of the long chain group and GCSF and IL-2 and IL-4 as members of the short chain group. For the GCSF families it was found that the fluctuations in the helical region are smaller than in the loop region, and it is clear that on the whole the smaller fluctuation residues belong to a large hydrophobic core region. Thus, it can be imagined how the receptor binding sites approach the receptor within the normal time-scale of pico seconds. In addition, two similar domain-type motions in low frequency modes were found with proteins in the long chain group, although we never observed any sequence similarity in the two separate two-domain regions in each protein of the long chain group. On the other hand, these two domain-type motions were not clear in proteins of the short chain group.     

Dynamic structure of subtilisin-eglin c complex studied by normal mode analysis

Normal mode analysis of subtilisin-eglin c complex was performed to investigate the dynamics at the interface between the enzyme and the inhibitor. The internal motions of the complex calculated from the normal modes were divided into three parts: the internal motions changing the shape of each molecule, the external rigid-body motions changing their mutual dispositions, and the coupling between the internal and external motions. From the results of the analysis, the following characteristic features were found in the dynamics at the interface regions: 1) negative correlation between the internal and external motions within each molecule, and 2) positive correlation between the external motions of the two molecules. The former decreases the apparent amplitudes of motions at the interface. The latter minimizes the interference between individual motions of the two molecules. These dynamic characteristics allow the enzyme and the inhibitor to move as freely as possible. This finding suggests that the experimental evidence of the large entropy gain on binding should be attributed not only to strong hydrophobic interactions, but also to the dynamic structure of the complex, which is found to minimize an unavoidable loss of the conformational entropy on binding.     

Unique mode of GCC box recognition by the DNA-binding domain of ethylene-responsive element-binding factor (ERF domain) in plant

Ethylene-responsive element-binding proteins (EREBPs)have novel DNA-binding domains (ERF domains), which are widely conserved in plants, and interact specifically with sequences containing AGCCGCC motifs (GCC box). Deletion experiments show that some flanking region at the N terminus of the conserved 59-amino acid ERF domain is required for stable binding to the GCC box. Three ERF domain-containing fragments of EREBP2, EREBP4, and AtERF1 from tobacco and Arabidopsis, bind to the sequence containing the GCC box with a high binding affinity in the pM range. The high affinity binding is conferred by a monomeric ERF domain fragment, and DNA truncation experiments show that only 11-base pair DNA containing the GCC box is sufficient for stable ERF domain interaction. Systematic DNA mutation analyses demonstrate that the specific amino acid contacts are confined within the 6-base pair GCCGCC region of the GCC box, and the first G, the fourth G, and the sixth C exhibit highest binding specificity common in all three ERF domain-containing fragments studied. Other bases within the GCC box exhibit modulated binding specificity varying from protein to protein, implying that these positions are important for differential binding by different EREBPs. The conserved N-terminal half is likely responsible for formation of a stable complex with the GCC box and the divergent C-terminal half for modulating the specificity.     

Trunk strength in combined motions of rotation and flexion/extension in normal young adults

Thirty-eight normal healthy young subjects (14 males, 24 females) with mean ages of 23 years (males) and 21 years (females), performed 36 functional rotational tasks of the trunk. The subject's lower extremities were stabilized in a stabilizing platform, allowing the entire motion of flexion-rotation and extension-rotation to take place in the trunk. Of these tasks, 18 were isometric and the other 18 were isokinetic. The isometric tasks consisted of flexion-rotation and extension-rotation from a 20 degrees, 40 degrees and 60 degrees flexed trunk in 20 degrees, 40 degrees and 60 degrees of axial rotation. The isokinetic activity consisted of flexion-rotation and extension-rotation from upright and flexed postures respectively in 20 degrees, 40 degrees and 60 degrees rotation planes at 15 degrees, 30 degrees and 60 degrees/s angular velocities. The results revealed that the males were significantly stronger than females (p < 0.01) and isometric activities produced significantly greater torque compared to isokinetic efforts (p < 0.01). The degree of trunk flexion was not significant; the angle of rotation, although significant, had only a small effect. The 60 degrees trunk rotation was significantly different from 20 degrees and 40 degrees of trunk rotation. The multiple regressions were all significant (p < 0.01); however, they predicted only 40 to 60% of the variations. Based on the results and analysis, it is suggested that the motion involved rather than the torque may have a consequential effect in the precipitation of back injuries.     

Heteronuclear relaxation study of the PH domain of beta-spectrin: restriction of loop motions upon binding inositol trisphosphate

The structural dynamics of protein ligand binding sites is one factor determining the specificity towards related ligands. In this context, the spectrin PH domain, which binds to a number of phosphatidylinositol lipid head groups, was investigated with respect to the dynamics of the binding loops. The latter were found to be of intermediate flexibility on a picosecond to nanosecond time-scale in the free protein and become more rigid upon ligand binding. Significant 15N and proton chemical shift changes occur in the binding loops. The internal correlation time, determined from 15N heteronuclear relaxation data using the standard model-free approach, decreases upon ligand binding. For several residues a concomitant rise in the generalized order parameter is observed. This is interpreted as a dampening effect of the ligand on a slow loop motion, while a fast component is not affected. Molecular dynamics simulations were performed to further investigate this situation. In fact, two time-scales of loop motions in the free state are observed in a 9 ns molecular dynamics trajectory. Agreement with generalized order parameters obtained from the experiment improves when a subtrajectory is analyzed that excludes rare dihedral transitions.     

A novel mode of DNA recognition by a beta-sheet revealed by the solution structure of the GCC-box binding domain in complex with DNA

The 3D solution structure of the GCC-box binding domain of a protein from Arabidopsis thaliana in complex with its target DNA fragment has been determined by heteronuclear multidimensional NMR in combination with simulated annealing and restrained molecular dynamic calculation. The domain consists of a three-stranded anti-parallel beta-sheet and an alpha-helix packed approximately parallel to the beta-sheet. Arginine and tryptophan residues in the beta-sheet are identified to contact eight of the nine consecutive base pairs in the major groove, and at the same time bind to the sugar phosphate backbones. The target DNA bends slightly at the central CG step, thereby allowing the DNA to follow the curvature of the beta-sheet.     

Characterizing semilocal motions in proteins by NMR relaxation studies

The understanding of protein function is incomplete without the study of protein dynamics. NMR spectroscopy is valuable for probing nanosecond and picosecond dynamics via relaxation studies. The use of 15N relaxation to study backbone dynamics has become virtually standard. Here, we propose to measure the relaxation of additional nuclei on each peptide plane allowing for the observation of anisotropic local motions. This allows the nature of local motions to be characterized in proteins. As an example, semilocal rotational motion was detected for part of a helix of the protein Escherichia coli flavodoxin.     

Discrimination of forearm''s motions by surface EMG signals using neural network

The effect of the antiviral, antitumoural xanthate D609 on the activity of phospholipase A2, C (PC- and Pi-specific) and D was investigated. D609 is the first model substance of a new concept of antiviral therapy that interferes with cellular regulation mechanisms, rather than with virus coded enzymes. Exclusively phosphatidylcholine (PC) specific phospholipase C (PC-PLC) was found to be inhibited in a dose-dependent manner. Enzyme activity was determined either as the rate of acid release from PC or as the rate of phosphorylcholine production form 3H labelled PC. Lineweaver-Burk plots revealed D609 as a competitive inhibitor of PC-PLC with a Ki of 6.4 microM. In addition, D609 competitively inhibited PC-PLC mediated cleavage of P-nitrophenylphosphorylcholine (p-NPP), a pseudo-substrate of PC-PLC with a Ki of 8.8 microM. These data suggest that D609 competes with the phosphorylcholine residue of PC for binding to PC-PLC.     

Enzyme specificity under dynamic control: a normal mode analysis of alpha-lytic protease

We have used alpha-lytic protease as a model system for exploring the relationship between the internal dynamics of an enzyme and its substrate specificity. The wild-type enzyme is highly specific for small substrates in its primary specificity pocket, while the M190A mutant has a much broader specificity, efficiently catalyzing cleavage of both large and small substrates. Normal modes have been calculated for both the wild-type and the mutant enzyme to determine how internal vibrations contribute to these contrasting specificity profiles. We find that for the atoms lining the walls of the specificity pocket, the wild-type normal modes have a more symmetric character, with the walls vibrating in phase, and the size of the pocket remaining relatively fixed. This is in agreement with X-ray crystallographic data on conformational substates trapped at 120 K. In contrast, we find that in the mutant, the binding pocket normal modes have a more antisymmetric character, with the walls vibrating out of phase, and the pocket able to expand and contract. These results suggest that the internal vibrations of a molecule may play an important role in determining substrate binding and specificity. A small change in protein structure can have a significant effect on the pattern of molecular vibrations, and thus on enzymatic properties, even if the overall amplitudes of the vibrations, as measured by NMR relaxation or crystallographic B-factors, remain largely unchanged.     

Dynamics of biomolecules: assignment of local motions by fluorescence anisotropy decay

Many biological systems have multiple fluorophores that experience multiple depolarizing motions, requiring multiple lifetimes and correlation times to define the fluorescence intensity and anisotropy decays, respectively. To simplify analyses, an assumption often made is that all fluorophores experience all depolarizing motions. However, this assumption usually is invalid, because each lifetime is not necessarily associated with each correlation time. To help establish the correct associations and recover accurate kinetic parameters, a general kinetic scheme that can examine all possible associations is presented. Using synthetic data sets, the ability of the scheme to discriminate among all nine association models possible for two lifetimes and two correlation times has been evaluated. Correct determination of the association model, and accurate recovery of the decay parameters, required the global analysis of related data sets. This general kinetic scheme was then used for global analyses of liver alcohol dehydrogenase anisotropy data sets. The results indicate that only one of the two tryptophan residues in each subunit is depolarized by process(es) independent of the enzyme's rotations. By applying the proper kinetic scheme and appropriate analysis procedures to time-resolved fluorescence anisotropy data, it is therefore possible to examine the dynamics of specific portions of a macromolecule in solution.