布氏田鼠肥满度分析和小型兽类肥满度指标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.     

Molecular dynamics simulations of phenylimidazole inhibitor complexes of cytochrome P450 cam

Molecular dynamics simulations have been performed on three phenylimidazole inhibitor complexes ofP450 _cam, utilizing the X-ray structures and the AMBER suite of programs. Compared to their corresponding optimized X-ray structures, very similar features were observed for the 1-phenylimidazole (1-PI) and 2-phenylimidazole (2-PI) complexes during a 100 ps MD simulation. The 1-PI inhibitor binds as a Type II complex with the imidazole nitrogen as a ligand of the heme iron. Analysis of the inhibitor-enzyme interctions during the MD simulations reveals that electrostatic interactions of the imidazole with the heme and van der Waals interactions of the phenyl ring with nearby hydrophobic residues are dominant. By contrast, 2-PI binds as a Type I inhibitor in the substrate binding pocket, but not as a ligand of the iron. The interactions of this inhibitor are qualitatively different from that of the Type II 1-PI, being mainly electrostatic/H-bonding interactions with a bound water and polar residues. Although the third compound, 4-PI, in common with 1-PI, also binds as a Type II inhibitor, with one nitrogen of the imidazole as a ligand to the iron, the MD average binding orientation deviates significantly from the X-ray structure. The most important changes observed include: (1) the rotation of the imidazole ring of this inhibitor by about 90° to enhance electrostatic interactions of the imidazole NH group with the carbonyl group of LEU244, and (2) the rotation of the carbonyl group of ASP251 to form a H-bond with VAL254. An analysis of the H-bonding network surrounding this substrate in the optimized crystal structure revealed that there is no H-bonding partner either for the free polar NH group in the imidazole ring of 4-phenylimidazole or for the polar carbonyl group of the nearby ASP251 residue. The deviation of the dynamically averaged inhibitor-enzyme structure of the 4-PI complex from the optimized crystal structure can therefore be rationalized as a consequence of the optimization of the electrostatic interactions among the polar groups.     

β-glucuronidase as a marker for clonal analysis of tomato lateral roots

In higher plants, the root-shoot axis established during embryogenesis is extended and modified by the development of primary and lateral apical meristems. While the structure of several shoot apical meristems has been deduced by combining histological studies with clonal analysis, the application of this approach to root apical meristems has been limited by a lack of visible genetic markers. We have tested the feasibility of using a synthetic gene consisting of the maize transposable elementActivator (Ac) inserted between a 35S CaMV promoter and the coding region of a -glucuronidase (GUS) reporter gene as a means of marking cell lineages in roots. The GUS gene was activated in individual cells byAc excision, and the resulting sectors of GUS-expressing cells were detected with the histochemical stain X-Gluc. Sectors in lateral roots originated from bothAc excision in meristematic cells and from parent root sectors that bisect the founder cell population for the lateral root initial. Analysis of root tip sectors confirmed that the root cap, and root proper have separate initials. Large sectors in the body of the lateral root encompassed both cortex and vascular tissues. The number of primary initial cells predicted from the size and arrangement of the sectors observed ranged from two to four and appeared to vary between roots. We conclude that transposon-based clonal analysis using GUS expression as a genetic marker is an effective approach for deducing the functional organization of root apical meristems.