Monitoring calcium-induced conformational changes in recoverin by electrospray mass spectrometry.

Recoverin is a calcium-binding protein that regulates the vertebrate photoresponse by inhibiting rhodopsin kinase in response to high calcium concentrations. It is heterogeneously N-acylated by myristoyl and related fatty acyl residues that are thought to act as "calcium-myristoyl switches," whereby, in the presence of Ca2+, the N-terminal acyl group is extended away from recoverin and, in the absence of calcium, it is more closely associated with the protein. Here we use electrospray ionization mass spectrometry (ESI/MS) to examine hydrogen isotopic exchange rates for specific regions of both acylated and nonacylated recoverin in the presence and absence of calcium. The deuterium exchange rates of three regions in the hydrophobic myristoyl binding pocket of acylated recoverin decreased in the absence of calcium. This effect is most likely due to the closer association of the acyl group with the protein under these conditions. In contrast, rates of deuterium incorporation increased in the absence of calcium for other regions, including the two functional calcium-binding sites. In addition to supporting the calcium-myristoyl switch hypothesis, a comparison of the behavior of acylated and unacylated recoverin revealed that the N-acyl group (N-lauroyl or N-myristoyl) exerts a significant stabilizing influence on the dynamics of recoverin. We demonstrate that the new technique of monitoring hydrogen isotopic exchange by ESI/MS can be used to obtain useful information concerning protein structures in solution using smaller amounts of protein and under more physiologically relevant conditions than is typically possible with NMR or X-ray crystallography.     

Ditylenchus dipsaci Infestation of Trifolium repens. I. Temperature Effects,Seedling Invasion,and a Field Survey

Rates of development of stem nematode (Ditylenchus dipsaci) in white clover (Trifolium repens) seedlings were found to be linearly related to temperature. Basal developmental temperature (T_b) was 3 °C, and the thermal constant (S) for development of gravid adult females from freshly laid eggs was 270 accumulated day-degrees above the T_b. Only 12% at 20 °C and 4% at 4 °C of the gravid female nematodes inoculated into seedling axils successfully penetrated seedling epidermis. These nematodes slowly migrated within the seedling and after a lag of 5 days at 20 °C started to lay eggs. The maximal rate of egg production was temperature-dependent, being 0.8 and 3.1 eggs female⁻¹ day⁻¹ at 10 and 20 °C, respectively. Nematodes emigrated rapidly from infested stolons when they were immersed in water, with rates being highest at 25 °C and lowest at 4 °C. The sensitivity to temperature of many of the parameters that govern nematode population dynamics indicates that climatic changes will have a marked effect upon this host-parasite system. A study of infested stolons from the field indicated that nematode numbers increased up to 3,000 or more before tissue senesence, triggered by nematode damage, caused a mass emigration of nematodes from the stolon.     

布氏田鼠肥满度分析和小型兽类肥满度指标K与KWL（重长指标）的比较

通过对两个肥满度指标的理论和生物学意义分析,以及对布氏田鼠肥满度的研究和实际应用的讨论,认为描述动物的肥满度时,重长指标ＫＷＬ优于指标Ｋ。两指标的最大差别是成体的ＫＷＬ值大于幼体,而成体的Ｋ值小于幼体。布氏田鼠肥满度没有性别差异;有异著的年龄差异,成体鼠的肥满度高于幼鼠;有显著的季节变化,鼠种群春季肥满度最高,夏季降低,秋季回升;有显著的年际变化,高数量年的肥满度高于低数量年。     

Self-organization of an oscillatory neural system

Hebbian dynamics is used to derive the differential equations for the synaptic strengths in the neural circuitry of the locomotive oscillator. Initially, neural connection are random. Under a specified arborization hypothesis relating to the density of neural connections, the differential equations are shown to model the self-organization and the stability of the oscillator.     

Long timestep dynamics of peptides by the dynamics driver approach

Previous experience with the Langevin/implicit-Euler scheme for dynamics (“LI”) on model systems (butane, water) has shown that LI is numerically stable for timesteps in the 5–20 fs range but quenches high-frequency modes. To explore applications to polypeptides, we apply LI to model systems (several dipeptides, a tetrapeptide, and a 13-residue oligoalanine) and also develop a new dynamics driver approach (“DA”). The DA scheme, based on LI, addresses the important issue of proper sampling, which is unlikely to be solved by small-time step integration methods or implicit methods with intrinsic damping at room temperature, such as LI. Equilibrium averages, time-dependent molecular properties, and sampling trends at room temperature are reported for both LI and DA dynamics simulations, which are then compared to those generated by a standard explicit discretization of the Langevin equation with a 1 fs timestep. We find that LI's quenching effects are severe on both the fast and slow (due to vibrational coupling) frequency modes of all-atom polypeptides and lead to more restricted dynamics at moderate timesteps (40 fs). The DA approach empirically counteracts these damping effects by adding random atomic perturbations to the coordinates at each step (before the minimization of a dynamics function). By restricting the energetic fluctuations and controlling the kinetic energy, we are able with a 60 fs timestep to generate continuous trajectories that sample more of the relevant conformational space and also reproduce reasonably Boltzmann statistics. Although the timescale for transition may be accelerated by the DA approach, the transitional. information obtained for the alanine dipeptide and the tetrapeptide is consistent with that obtained by several other theoretical approaches that focus specifically on the determination of pathways. While the trajectory for oligoalanine by the explicit scheme over the nanosecond timeframe remains in the vicinity of the full α_R-helix starting structure, and a high-temperature (6000°K) MD trajectory departs slowly from the a helical structure, the DA-generated trajectory for the same CPU time exhibits unfolding and refolding and reveals a range of conformations with an intermediate helix content. Significantly, this range of states is more consistent with spectroscopic experiments on small peptides, as well as the cooperative two-state model for helix–coil transition. The good, near-Boltzmann statistics reported for the smaller systems above, in combination with the interesting oligoalanine results, suggest that DA is a promising tool for efficiently exploring conformational spaces of biomolecules and exploring folding/unfolding processes of polypeptides. © 1995 Wiley-Liss, Inc.     

Protein conformational landscapes: Energy minimization and clustering of a long molecular dynamics trajectory

Using energy minimization and cluster analysis, we have analyzed a 1020 ps molecular dynamics trajectory of solvated bovine pancreatic trypsin inhibitor. Elucidation of conformational sub states in this way both illustrates the degree of conformational convergence in the simulation and reduces the structural data to a tractable subset. The relative movement of structures upon energy minimization was used to estimate the sizes of features on the protein potential energy surface. The structures were analyzed using their pairwise root-mean-square C_α deviations, which gave a global measure of conformational changes that would not be apparent by monitoring single degrees of freedom. At time scales of 0.1 ps, energy minimization detected sharp transitions between energy minima separated by 0.1 Å rms deviation. Larger conformational clusters containing these smaller minima and separated by 0.25 Å were seen at 1 ps time scales. Both of these small features of the conformational landscape were characterized by movements in loop regions associated with small, correlated backbone dihedral angle shifts. On a nanosecond time scale, the main features of the protein energy landscape were clusters separated by over 0.7 Å rms deviation, with only seven of these sub states visited over the 1 ns trajectory. These substates, discernible both before and after energy minimization, differ mainly in a monotonic pivot of the loop residues 11–18 over the course of the simulation. This loop contains lysine 17, which specifically binds to trypsin in the active site. The trajectory did not return to previously visited clusters, indicating that this trajectory has not been shown to have completely sampled the conformational substates available to it. Because the apparent convergence to a single region of conformation space depends on both the time scale of observation and the size of the conformational features examined, convergence must be operationally defined within the context of the simulation. © 1995 Wiley-Liss, Inc.     

Dynamic properties of some β-chain mutant hemoglobins

The thermal behavior of the Soret band relative to the carbonmonoxy derivatives of some β-chain mutant hemoglobins is studied in the temperature range 300–10 K and compared to that of wild-type carbonmonoxy hemoglobin. The band profile at various temperatures is modeled as a Voigt function that accounts for homogeneous broadening and for the coupling with high- and low-frequency vibrational modes, while inhomogeneous broadening is taken into account with a gaussian distribution of purely electronic transition frequencies. The various contributions to the overall bandwidth are singled out With this analysis and their temperature dependence, in turn, gives information on structural and dynamic properties of the system studied. In the wildtype and mutant hemoglobins, the values of homogeneous bandwidth and of the coupling constants to high-frequency vibrational modes are not modified with respect to natural human hemoglobin, thus indicating that the local electronic and vibrational properties of the heme–CO complex are not altered by the recombinant procedures. On the contrary, differences in the protein dynamic behavior are observed. The most relevant are those relative to the “polar isosteric” βVal-67(Ell) →Thr substitution, localized in the heme pocket, which results in decreased coupling with low-frequency modes and increased anharmonic motions. Mutations involving residue βLys-144(HC1) at the C-terminal and residue βCys-112(G14) at the α₁β₁ interface have a smaller effect consisting in an increased coupling with low-frequency modes. Mutations at the β-N-terminal and at the α₁β₂ interface have no effect on the dynamic properties of the heme pocket. © 1995 Wiley-Liss, Inc.     

Combined use of 13C chemical shift and 1Hα−13Cα heteronuclear NOE data in monitoring a protein NMR structure refinement

Summary A large portion of the ¹³C resonance assignments for murine epidermal growth factor (mEGF) at pH 3.1 and 28°C has been determined at natural isotope abundance. Sequence-specific ¹³C assignments are reported for 100% of the assignable C, 96% of the C, 86% of the aromatic and 70% of the remaining peripheral aliphatic resonances of mEGF. A good correlation was observed between experimental and back-calculated C chemical shifts for regions of regular -sheet structure. These assignments also provide the basis for interpreting ¹H–¹³C heteronuclear NOE (HNOE) values in mEGF at natural isotope abundance. Some of the backbone polypeptide segments with high internal mobility, indicated by these ¹H–¹³C HNOE measurements, correlate with locations of residues involved in the putative mEGF-receptor binding site. Using four families of mEGF structures obtained over the last few years, we demonstrate that standard deviations between experimental and back-calculated C values can be used to monitor the refinement of this protein's structure, particularly for -sheet regions. Improved agreement between calculated and observed values of C is correlated with other measures of structure quality, including lowered values of residual constraint violations and more negative values of conformational energy. These results support the view that experimental conformation-dependent chemical shifts, C, can provide a reliable source of information for monitoring the process of protein structure refinement and are potentially useful restraints for driving the refinement.Abbreviations HSQC heteronuclear single-quantum coherence spectroscopy - PFG pulsed-field gradient - TOCSY ¹H-¹H total correlation spectroscopy - EGF epidermal growth factor - mEGF murine EGF - hEGF human EGF - hTGF human type- transforming growth factor - DIPSI spm-locking pulse sequence - NOE nuclear Overhauser effect - HNOE heteronuclear Overhauser effect     

Larmor frequency selective model free analysis of protein NMR relaxation

Summary The Lipari-Szabo dynamical formalism is extended by setting the time constants of the Lorentzian terms to and . This analysis is compared to the earlier proposed three-parameter extended model free formalism with regard to the range of equivalence and the advantages of the simplified two-parameter (S _inff ^sup2 ,S _infH ^sup2 ) and (S _inff ^sup2 ,S _infN ^sup2 ) representations. Spectral density components are calculated and compared to those obtained from the spectral density analysis formalism. Protein relaxation data, commonly analyzed in terms of the two-parameter representation, may correspond to a dynamically heterogeneous behaviour that is more appropriately represented in terms of a fast limit order parameter and a second, lower frequency order parameter.     

Dynamic modelling of a helical peptide in solution using NMR data: Multiple conformations and multi-spin effects

Summary Nuclear Overhauser effect (NOE) measurements on molecules in solution provide information about only the ensemble-averaged properties of these molecules. An algorithm is presented that uses a list of NOEs to produce an ensemble of molecules that on average agrees with these NOEs, taking into account the effect of surrounding spins on the buildup of each NOE (spin diffusion). A simplified molecular dynamics simulation on several copies of the molecule in parallel is restrained by forces that are derived directly from differences between calculated and measured NOEs. The algorithm is tested on experimental NOE data of a helical peptide derived from bovine pancreatic trypsin inhibitor.