A comparative analysis of the equilibrium dynamics of a designed protein inferred from NMR,X‐ray,and computations

A detailed analysis of high‐resolution structural data and computationally predicted dynamics was carried out for a designed sugar‐binding protein. The mean‐square deviations in the positions of residues derived from nuclear magnetic resonance (NMR) models and those inferred from X‐ray crystallographic B‐factors for two different crystal forms were compared with the predictions based on the Gaussian Network Model (GNM) and the results from molecular dynamics (MD) simulations. GNM systematically yielded a higher correlation than MD, with experimental data, suggesting that the lack of atomistic details in the coarse‐grained GNM is more than compensated for by the mathematically exact evaluation of fluctuations using the native contacts topology. Evidence is provided that particular loop motions are curtailed by intermolecular contacts in the crystal environment causing a discrepancy between theory and experiments. Interestingly, the information conveyed by X‐ray crystallography becomes more consistent with NMR models and computational predictions when ensembles of X‐ray models are considered. Less precise (broadly distributed) ensembles indeed appear to describe the accessible conformational space under native state conditions better than B‐factors. Our results highlight the importance of using multiple conformations obtained by alternative experimental methods, and analyzing results from both coarse‐grained models and atomic simulations, for accurate assessment of motions accessible to proteins under native state conditions. Proteins 2009. © 2009 Wiley‐Liss, Inc.     

Measurement of chemical exchange between RNA conformers by F NMR

Many noncoding RNA molecules adopt alternative secondary and tertiary conformations that are critical for their roles in gene expression. Although many of these rearrangements are mediated by other biomolecular components, it is important to evaluate the equilibrium relationship of the conformers. To measure the spontaneous interconversion in a bi-stable RNA stem loop sequence into which a single ¹⁹F-uridine label was incorporated, a ¹⁹F–¹⁹F EXSY experiment was employed. The kinetic exchange rate measured from EXSY experiments for this system was 37.3 ± 2.8 s⁻¹. The advantage of this approach is that exchange kinetics can be monitored in any RNA sequence into which a single ¹⁹F nucleotide is incorporated by commercial synthesis. This method is therefore suitable for application to biologically significant systems in which dynamic conformational rearrangement is important for function and may therefore facilitate studies of RNA structure–function relationships.     

Protein loop modeling by using fragment assembly and analytical loop closure

Protein loops are often involved in important biological functions such as molecular recognition, signal transduction, or enzymatic action. The three dimensional structures of loops can provide essential information for understanding molecular mechanisms behind protein functions. In this article, we develop a novel method for protein loop modeling, where the loop conformations are generated by fragment assembly and analytical loop closure. The fragment assembly method reduces the conformational space drastically, and the analytical loop closure method finds the geometrically consistent loop conformations efficiently. We also derive an analytic formula for the gradient of any analytical function of dihedral angles in the space of closed loops. The gradient can be used to optimize various restraints derived from experiments or databases, for example restraints for preferential interactions between specific residues or for preferred backbone angles. We demonstrate that the current loop modeling method outperforms previous methods that employ residue‐based torsion angle maps or different loop closure strategies when tested on two sets of loop targets of lengths ranging from 4 to 12. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

Comparative 5-doxylstearoyllecithin and 3-doxylcholestane EPR spin labeling study of phospholipid bilayer perturbation by different oxidized lecithin species

A 3-doxylcholestane spin label was employed in addition to 5-doxylstearoyl lecithin for a more detailed study of the different effects exerted by variously oxidized lecithins on fatty acid alignment in phospholipid planar bilayers. Either spin label was enclosed in oriented PLPC planar samples also containing in turn a variety of conjugated-dienes lecithins and cleaved chain lecithins, in order to monitor EPR spectral angular dependence loss. Data obtained by use of arachidonoyl-hydroxystearoyl-PC and palmitoyl-hydroxylinoleoyl-PC confirm that the 5-DSPC nitroxide ring almost completely retains its orientation in CD-PCs-containing planar membranes, in contrast with angular dependence loss observed in the presence of the CC-PC molecular species palmitoyl-oxononanoyl-PC and palmitoyl-oxovaleroyl-PC, already seen with cleaved-chain palmitoyl-glutaroyl-PC and palmitoyl-azelaoyl-PC. However, the 3-DC nitroxide ring also loses its orientation with CD-PCs, in addition to being disoriented by cleaved chain-lecithins, similarly to 5DSPC. Joint information from the two spin labels will help to clarify whether OXPC-related disordering involved the whole bilayer structure or only the hydrophobic core. In addition, the propensity of different OXPCs to form bilayer vesicles in water suspension was also determined by Sepharose 4B gel-chromatography. The results suggest that CD-PCs might yield SPB bilayer structures with a disordered hydrophobic core, while pure CC-PC more probably forms non-bilayer disordered structures, possibly micelles or mixed micelle/bilayers.     

Trapping of the quenched conformation associated with non-photochemical quenching of chlorophyll fluorescence at low temperature

The kinetics of non-photochemical quenching (NPQ) of chlorophyll fluorescence was studied in pea leaves at different temperatures between 5 and 25°C and during rapid jumps of the leaf temperature. At 5°C, NPQ relaxed very slowly in the dark and was sustained for up to 30 min. This was independent of the temperature at which quenching was induced. Upon raising the temperature to 25°C, the quenched state relaxed within 1 min, characteristic for qE, the energy-dependent component of NPQ. Measurements of the membrane permeability (ΔA₅₁₅) in dark-adapted and preilluminated leaves and NPQ in the presence of dithiothreitol strongly suggest that the effect of low temperature on NPQ was not because of limitation by the lumenal pH or the de-epoxidation state of the xanthophylls. These data are consistent with the notion that the transition from the quenched to the unquenched state and vice versa involves a structural reorganization in the photosynthetic apparatus. An eight-state reaction scheme for NPQ is proposed, extending the model of Horton and co-workers (FEBS Lett 579:4201–4206, 2005), and a hypothesis is put forward concerning the nature of conformational changes associated with qE. Electronic supplementary material The online version of this article (doi: ) contains supplementary material, which is available to authorized users.     

On a probabilistic model for the numerical estimation of nocturnal migration of birds

The study of nocturnal bird migration by cone methods of observation has a century-long history but has continued to be used up to the present. To describe the flux and estimate the number of passing birds a probabilistic model is proposed. This model is based on the concept of dynamic Poisson ensemble of points in appropriate phase space and has two parameters. One is scalar and the other one is functional. We constructed consistent estimations of these parameters and discuss their use for the numerical estimation of the flux of birds observed in a narrow light cone generated by the bright lunar disk and formed by an open angle of telescope. Selection on the same type of birds was suggested as the necessary condition for the model application. Ground speed of each bird was introduced into the model as a new but obligatory value determining the quantification of the flux of bird.     

New steroid-fused P-heterocycles. Part II. Synthesis and conformational study of oxazaphosphorino[16,17-e]estrone derivatives

16Beta-aminomethyl-17beta-hydroxyestrone 3-methyl ether 6 and its N-propyl (17), N-benzyl (18) and N-arylmethyl derivatives (19-22) were subjected to ring closure reactions with phenylphosphonic dichloride in order to synthetize P-epimeric oxazaphosphorinanes 23a, 24-29 in which the hetero ring is condensed to ring D of the sterane skeleton. The stereostructures of the products were evaluated by 1H, 13C and 31P NMR spectroscopy. The geometry was optimized by utilizing the B3LYP DFT method. The NMR spectral data and the results of the ab initio calculations demonstrated that the stereostructure of the hetero ring was strongly affected by the rigid sterane framework condensed to it, and the phosphoramidate ring proved to adopt predominantly a distorted-boat conformation, regardless of the P-configuration.     

Dynamics of electron transfer from high-potential cytochrome <Emphasis Type="Italic">c</Emphasis> to bacteriochlorophyll dimer in photosynthetic reaction centers as probed using laser-induced temperature jump

Laser-induced temperature jump experiments were used for testing the rates of thermoinduced conformational transitions of reaction center (RC) complexes in chromatophores of Chromatium minutissimum. The thermoinduced transition of the macromolecular RC complex to a state providing effective electron transport from the multiheme cytochrome c to the photoactive bacteriochlorophyll dimer within the temperature range 220–280 K accounts for tens of seconds with activation energy 0.166 eV/molecule. The rate of the thermoinduced transition in the cytochrome–RC complex was found to be three orders of magnitude slower than the rate of similar thermoinduced transition of the electron transfer reaction from the primary to secondary quinone acceptors studied in the preceding work (Chamorovsky et al. in Eur Biophys J 32:537–543, 2003). Parameters of thermoinduced activation of the electron transfer from the multiheme cytochrome c to the photoactive bacteriochlorophyll dimer are discussed in terms of cytochrome c docking onto the RC.     

Allosteric communication in dihydrofolate reductase: signaling network and pathways for closed to occluded transition and back

Escherichia coli dihydrofolate reductase (DHFR) catalyzes the reduction of dihydrofolate to tetrahydrofolate. During the catalytic cycle, DHFR undergoes conformational transitions between the closed (CS) and occluded (OS) states that, respectively, describe whether the active site is closed or occluded by the Met20 loop. The CS→OS and the reverse transition may be viewed as allosteric transitions. Using a sequence-based approach, we identify a network of residues that represents the allostery wiring diagram. Many of the residues in the allostery wiring diagram, which are dispersed throughout the adenosine-binding domain as well as the loop domain, are not conserved. Several of the residues in the network have been previously shown by NMR experiments, mutational studies, and molecular dynamics simulations to be linked to equilibration conformational fluctuations of DHFR. To further probe the nature of events that occur during conformational fluctuations, we use a self-organized polymer model to monitor the kinetics of the CS→OS and the reverse transitions. During the CS→OS transition, coordinated changes in a number of residues in the loop domain enable the Met20 loop to slide along the α-helix in the adenosine-binding domain. Sliding is triggered by pulling of the Met20 loop by the βG-βH loop and the pushing action of the βG-βH loop. The residues that facilitate the Met20 loop motion are part of the network of residues that transmit allosteric signals during the CS→OS transition. Replacement of M16 and G121, whose C^α atoms are about 4.3 Å in the CS, by a disulfide cross-link impedes that CS→OS transition. The order of events in the OS→CS transition is not the reverse of the forward transition. The contact Glu18-Ser49 in the OS persists until the sliding of the Met20 loop is nearly complete. The ensemble of structures in the transition state in both the allosteric transitions is heterogeneous. The most probable transition-state structure resembles the OS (CS) in the CS→OS (OS→CS) transition, which is in accord with the Hammond postulate. Structures resembling the OS (CS) are present as minor (∼ 1-3%) components in equilibrated CS (OS) structures.