Spectral densities of nitrogen nuclei in Escherichia coli ribonuclease HI obtained by 15N NMR relaxation and molecular dynamics

Summary Spectral densities of the ¹⁵N amide in Escherichia coli ribonuclease HI, obtained from NMR relaxation experiments, were compared with those calculated using a molecular dynamics (MD) simulation. All calculations and comparisons assumed that the auto-correlation function describing the internal motions of the molecule was independent of the auto-correlation function associated with overall rotational diffusion. Comparisons were limited to those residues for which the auto-correlation function of internal motions rapidly relaxed and reached a steady state within 205 ps. The results show the importance of frequency components as well as amplitudes of internal motions in order to obtain a meaningful comparison of MD simulations with NMR data.     

Fundamental and biotechnological applications of neutron scattering measurements for macromolecular dynamics

To explore macromolecular dynamics on the picosecond timescale, we used neutron spectroscopy. First, molecular dynamics were analyzed for the hyperthermophile malate dehydrogenase from Methanococcus jannaschii and a mesophilic homologue, the lactate dehydrogenase from Oryctolagus cunniculus muscle. Hyperthermophiles have elaborate molecular mechanisms of adaptation to extremely high temperature. Using a novel elastic neutron scattering approach that provides independent measurements of the global flexibility and of the structural resilience (rigidity), we have demonstrated that macromolecular dynamics represents one of these molecular mechanisms of thermoadaptation. The flexibilities were found to be similar for both enzymes at their optimal activity temperature and the resilience is higher for the hyperthermophilic protein. Secondly, macromolecular motions were examined in a native and immobilized dihydrofolate reductase (DHFR) from Escherichia coli. The immobilized mesophilic enzyme has increased stability and decreased activity, so that its properties are changed to resemble those of the thermophilic enzyme. Are these changes reflected in dynamical behavior? For this study, we performed quasielastic neutron scattering measurements to probe the protein motions. The residence time is 7.95 ps for the native DHFR and 20.36 ps for the immobilized DHFR. The average height of the potential barrier to local motions is therefore increased in the immobilized DHFR, with a difference in activation energy equal to 0.54 kcal/mol, which is, using the theoretical rate equation, of the same order than expected from calculation.     

Hemoglobin dynamics in red blood cells: correlation to body temperature

A transition in hemoglobin behavior at close to body temperature has been discovered recently by micropipette aspiration experiments on single red blood cells (RBCs) and circular dichroism spectroscopy on hemoglobin solutions. The transition temperature was directly correlated to the body temperatures of a variety of species. In an exploration of the molecular basis for the transition, we present neutron scattering measurements of the temperature dependence of hemoglobin dynamics in whole human RBCs in vivo. The data reveal a change in the geometry of internal protein motions at 36.9°C, at human body temperature. Above that temperature, amino acid side-chain motions occupy larger volumes than expected from normal temperature dependence, indicating partial unfolding of the protein. Global protein diffusion in RBCs was also measured and the findings compared favorably with theoretical predictions for short-time self-diffusion of noncharged hard-sphere colloids. The results demonstrated that changes in molecular dynamics in the picosecond time range and angstrom length scale might well be connected to a macroscopic effect on whole RBCs that occurs at body temperature.     

Location and dynamics of acyl chain NBD-labeled phosphatidylcholine (NBD-PC) in DPPC bilayers. A molecular dynamics and time-resolved fluorescence anisotropy study

100-ns molecular dynamics simulations of fluid 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) bilayers, both pure and containing 7-nitrobenz-2-oxa-1,3-diazol-4-yl (NBD) acyl-chain labeled fluorescent analogs (C6-NBD-PC and C12-NBD-PC), are described. These molecules are widely used as probes for lipid structure and dynamics. The results obtained here for pure DPPC agree with both experimental and theoretical published works. We verified that the NBD fluorophore of both derivatives loops to a transverse location closer to the interface than to the center of the bilayer. Whereas this was observed previously in experimental literature works, conflicting transverse locations were proposed for the NBD group. According to our results, the maximum of the transverse distribution of NBD is located around the glycerol backbone/carbonyl region, and the nitro group is the most external part of the fluorophore. Hydrogen bonds from the NH group of NBD (mostly to glycerol backbone lipid O atoms) and to the nitro O atoms of NBD (from water OH groups) are continuously observed. Rotation of NBD occurs with ∼ 2.5-5 ns average correlation time for these probes, but very fast, unresolved reorientation motions occur in < 20 ps, in agreement with time-resolved fluorescence anisotropy measurements. Finally, within the uncertainty of the analysis, both probes show lateral diffusion dynamics identical to DPPC.     

Energy transfer processes in Gloeobacter violaceus PCC 7421 that possesses phycobilisomes with a unique morphology

We examined energy transfer dynamics in phycobilisomes (PBSs) of cyanobacteria in relation to the morphology and pigment compositions of PBSs. We used Gloeobacter violaceus PCC 7421 and measured time-resolved fluorescence spectra in three types of samples, i.e., intact cells, PBSs, and rod assemblies separated from cores. Fremyella diplosiphon, a cyanobacterial species well known for its complementary chromatic adaptation, was used for comparison after growing under red or green light. Spectral data were analyzed by the fluorescence decay-associated spectra with components common in lifetimes with a time resolution of 3 ps/channel and a spectral resolution of 2 nm/channel. This ensured a higher resolution of the energy transfer kinetics than those obtained by global analysis with fewer sampling intervals. We resolved four spectral components in phycoerythrin (PE), three in phycocyanin (PC), two in allophycocyanin, and two in photosystem II. The bundle-like PBSs of G. violaceus showed multiple energy transfer pathways; fast (≈ 10 ps) and slow (≈ 100 ps and ≈ 500 ps) pathways were found in rods consisting of PE and PC. Energy transfer time from PE to PC was two times slower in G. violaceus than in F. diplosiphon grown under green light.     

Lactose Binding to Galectin-1 Modulates Structural Dynamics, Increases Conformational Entropy, and Occurs with Apparent Negative Cooperativity

Galectins are a family of lectins with a conserved carbohydrate recognition domain that interacts with β-galactosides. By binding cell surface glycoconjugates, galectin-1 (gal-1) is involved in cell adhesion and migration processes and is an important regulator of tumor angiogenesis. Here, we used heteronuclear NMR spectroscopy and molecular modeling to investigate lactose binding to gal-1 and to derive solution NMR structures of gal-1 in the lactose-bound and unbound states. Structure analysis shows that the β-strands and loops around the lactose binding site, which are more open and dynamic in the unbound state, fold in around the bound lactose molecule, dampening internal motions at that site and increasing motions elsewhere throughout the protein to contribute entropically to the binding free energy. CD data support the view of an overall more open structure in the lactose-bound state. Analysis of heteronuclear single quantum coherence titration binding data indicates that lactose binds the two carbohydrate recognition domains of the gal-1 dimer with negative cooperativity, in that the first lactose molecule binds more strongly (K₁ = 21 ± 6 × 10³ M^− 1) than the second (K₂ = 4 ± 2 × 10³ M^− 1). Isothermal calorimetry data fit using a sequential binding model present a similar picture, yielding K₁ = 20 ± 10 × 10³ M^− 1 and K₂ = 1.67 ± 0.07 × 10³ M^− 1. Molecular dynamics simulations provide insight into structural dynamics of the half-loaded lactose state and, together with NMR data, suggest that lactose binding at one site transmits a signal through the β-sandwich and loops to the second binding site. Overall, our results provide new insight into gal-1 structure-function relationships and to protein-carbohydrate interactions in general.     

Subpicosecond midinfrared spectroscopy of the Pfr reaction of phytochrome Agp1 from Agrobacterium tumefaciens

Phytochromes are light-sensing pigments found in plants and bacteria. For the first time, the P_fr photoreaction of a phytochrome has been subject to ultrafast infrared vibrational spectroscopy. Three time constants of 0.3 ps, 1.3 ps, and 4.0 ps were derived from the kinetics of structurally specific marker bands of the biliverdin chromophore of Agp1-BV from Agrobacterium tumefaciens after excitation at 765 nm. VIS-pump-VIS-probe experiments yield time constants of 0.44 ps and 3.3 ps for the underlying electronic-state dynamics. A reaction scheme is proposed including two kinetic steps on the S₁ excited-state surface and the cooling of a vibrationally hot P_fr ground state. It is concluded that the upper limit of the E-Z isomerization of the C₁₅ = C₁₆ methine bridge is given by the intermediate time constant of 1.3 ps. The reaction scheme is reminiscent of that of the corresponding P_r reaction of Agp1-BV as published earlier.     

Solvent effects in the slow dynamics of proteins

The influence of solvent on the slow internal dynamics of proteins is studied by comparing molecular dynamics simulations of solvated and unsolvated lysozyme. The dynamical trajectories are projected onto the protein's normal modes in order to obtain a separate analysis for each of the associated time scales. The results show that solvent effects are important for the slowest motions (below approximately 1 ps(-1)) but negligible for faster motions. The damping effects seen in the latter show that the principal source of friction in protein dynamics is not the solvent, but the protein itself.     

A Direct Coupling between Global and Internal Motions in a Single Domain Protein? MD Investigation of Extreme Scenarios

Proteins are not rigid molecules, but exhibit internal motions on timescales ranging from femto- to milliseconds and beyond. In solution, proteins also experience global translational and rotational motions, sometimes on timescales comparable to those of the internal fluctuations. The possibility that internal and global motions may be directly coupled has intriguing implications, given that enzymes and cell signaling proteins typically associate with binding partners and cellular scaffolds. Such processes alter their global motion and may affect protein function. Here, we present molecular dynamics simulations of extreme case scenarios to examine whether a possible relationship exists. In our model protein, a ubiquitin-like RhoGTPase binding domain of plexin-B1, we removed either internal or global motions. Comparisons with unrestrained simulations show that internal and global motions are not appreciably coupled in this single-domain protein. This lack of coupling is consistent with the observation that the dynamics of water around the protein, which is thought to permit, if not stimulate, internal dynamics, is also largely independent of global motion. We discuss implications of these results for the structure and function of proteins.     

OmpA: Gating and dynamics via molecular dynamics simulations

Outer membrane proteins (OMPs) of Gram-negative bacteria have a variety of functions including passive transport, active transport, catalysis, pathogenesis and signal transduction. Whilst the structures of ∼ 25 OMPs are currently known, there is relatively little known about their dynamics in different environments. The outer membrane protein, OmpA from Escherichia coli has been studied extensively in different environments both experimentally and computationally, and thus provides an ideal test case for the study of the dynamics and environmental interactions of outer membrane proteins. We review molecular dynamics simulations of OmpA and its homologues in a variety of different environments and discuss possible mechanisms of pore gating. The transmembrane domain of E. coli OmpA shows subtle differences in dynamics and interactions between a detergent micelle and a lipid bilayer environment. Simulations of the crystallographic unit cell reveal a micelle-like network of detergent molecules interacting with the protein monomers. Simulation and modelling studies emphasise the role of an electrostatic-switch mechanism in the pore-gating mechanism. Simulation studies have been extended to comparative models of OmpA homologues from Pseudomonas aeruginosa (OprF) and Pasteurella multocida (PmOmpA), the latter model including the periplasmic C-terminal domain.     

Hydration effect on low-frequency protein dynamics observed in simulated neutron scattering spectra

Hydration effects on protein dynamics were investigated by comparing the frequency dependence of the calculated neutron scattering spectra between full and minimal hydration states at temperatures between 100 and 300 K. The protein boson peak is observed in the frequency range 1-4 meV at 100 K in both states. The peak frequency in the minimal hydration state shifts to lower than that in the full hydration state. Protein motions with a frequency higher than 4 meV were shown to undergo almost harmonic motion in both states at all temperatures simulated, whereas those with a frequency lower than 1 meV dominate the total fluctuations above 220 K and contribute to the origin of the glass-like transition. At 300 K, the boson peak becomes buried in the quasielastic contributions in the full hydration state but is still observed in the minimal hydration state. The boson peak is observed when protein dynamics are trapped within a local minimum of its energy surface. Protein motions, which contribute to the boson peak, are distributed throughout the whole protein. The fine structure of the dynamics structure factor is expected to be detected by the experiment if a high resolution instrument (<∼20 μeV) is developed in the near future.     

Photonic characteristics and ex vivo imaging of Escherichia coli-Xen14 within the bovine reproductive tract

The objectives of this study were to (1) characterize the photonic properties of Escherichia coli-Xen14 and (2) conduct photonic imaging of E. coli-Xen14 within bovine reproductive tract segments (RTS) ex vivo (Bos indicus). E. coli-Xen14 was grown for 24 h in Luria Bertani medium (LB), with or without kanamycin (KAN). Every 24 h, for an 8-d interval, inoculums were imaged and photonic emissions (PE) collected. Inoculums were subcultured and plated daily to determine the colony forming units (CFU) and ratio of photon emitters to nonemitters. In the second objective, abattoir-derived bovine reproductive tracts (n = 9) were separated into posterior and anterior vagina, cervix, uterine body, and uterine horns. Two concentrations (3.2 × 10⁸ and 3.2 × 10⁶ CFU/200 μL for relative [High] and [Low], respectively) of E. coli-Xen14 were placed in translucent tubes for detection of PE through RTS. The CFU did not differ (P = 0.31) over time with or without KAN presence; they remained stable with 99.93% and 99.98% photon emitters, respectively. However, PE were lower (P < 0.0001) in cultures containing KAN than in those containing no KAN (629.8 ± 117.7 vs. 3012.0 ± 423.5 relative lights units per second [RLU/sec], respectively). On average, the percentage of PE between RTS, for both concentrations, was higher (P < 0.05) in the uterine body. In summary, E. coli-Xen14 remained stable with respect to the proportions of photon emitters with or without KAN (used to selectively culture E. coli-Xen14). However, KAN presence suppressed photonic activity. The ability to detect PE through various segments of the reproductive tract demonstrated the feasibility of monitoring the presence of E. coli-Xen14 in the bovine reproductive tract ex vivo.     

How long does DNA keep the memory of its conformation? A time-dependent canonical correlation analysis of molecular dynamics simulation

The time dependence of the correlation between motions of different parts of DNA is analyzed from a 200 ps molecular dynamics simulation of the double-stranded self-complementary d(CTGATCAG) in the B form. Each nucleotide is decomposed into three subunits corresponding to the furanose ring (SU), the base (BA), and the backbone (SK). The motion of each subunit is considered as the superimposition of rigid body translation, rigid body rotation, and internal deformation. Canonical time-dependent correlation functions calculated with coordinates describing the different components of the subunits motion are defined and computed. This allows us to probe how long a particular type of motion of one subunit influences the other types of motions of other subunits (cross correlation functions) or how long a particular subunit keeps the memory of its own conformation or location (autocorrelation functions). From auto-correlation analysis it is found that deformation decorrelates within a few tenths of picoseconds, rotational correlation times are on the order of 8 ps, while translational motions are long-time correlated. The deformation of a subunit is not correlated to the deformation of another one (at the 200 ps time scale of our simulation), but influences slightly their translation and orientation as time increases. © 1996 John Wiley & Sons, Inc.     

Studies on chromophore coupling in isolated phycobiliproteins: III. Picosecond excited state kinetics and time-resolved fluorescence spectra of different allophycocyanins from Mastigocladus laminosus

The excited state kinetics of three different allophycocyanin (AP) complexes has been studied by picosecond fluorescence spectroscopy. Both the fluorescence kinetics and the decay-associated fluorescence spectra of the different complexes can be understood on the basis of a structural model for AP which uses (a) an analogy to the known x-ray determined structure of C-phycocyanin, (b) the biochemical analogies of AP and C-phycocyanin, and (c) the biochemical composition of AP-B (AP-681). A model is developed that describes the excited state kinetics as a mixture of internal conversion processes within a coupled exciton pair and energy transfer processes between exciton pairs. We found excited state relaxation times in the range of 13 ps (AP with linker peptide) up to 66 ps (AP-B). The trimeric aggregates AP 660 and AP 665 show one fast relaxation component each, as was expected on the basis of their symmetry properties. The lower symmetry of AP-B (AP-681) gives rise to two fast lifetime components (τ₁ = 23 ps and τ₂ = 66 ps) which are attributed to internal conversion and/or energy transfer between excitonic states formed by the coupling of symmetrically and spectrally nonequivalent chromophores. It is proposed that the internal conversion between exciton states of strongly coupled chromophores fulfills the requirements of the small energy gap limit. Thus, internal conversion rates in the order of tens of picoseconds are feasible. The influence of the interaction of the linker peptide on the properties of the AP trimer are manifested in the fluorescence kinetics. Lack of the linker peptide in AP 660 gives rise to a heterogeneity in the chromophore conformations and chromophore-chromophore interactions.     

Backbone dynamics of escherichia coli adenylate kinase at the extreme stages of the catalytic cycle studied by (15)N NMR relaxation

Adenylate kinase from Escherichia coli (AKeco), consisting of a single 23.6 kDa polypeptide chain folded into domains CORE, AMPbd, and LID, catalyzes the reaction AMP + ATP --> 2ADP. Domains LID and AMPbd execute large-scale movements during catalysis. Backbone dynamics of ligand-free and AP(5)A-inhibitor-bound AKeco were studied comparatively with (15)N NMR relaxation methods. Overall diffusion with correlation times of 15.05 (11.42) ns and anisotropy D(parallel)/D(perp) = 1.25 (1.10), and fast internal motions with correlation times up to 100 ps (50 ps), were determined for AKeco (AKecoAP(5)A). Fast internal motions affect 93% of the AKeco sites, with pronounced preference for domains AMPbd and LID, and 47% of the AKecoAP(5)A sites, with limited variability along the chain. The mean squared generalized order parameters, , of secondary structure elements and loops are affected by ligand binding differentially and in a domain-specific manner. Nanosecond motions predominate within AMPbd. Prominent exchange contributions, associated in particular with residue G10 of the nucleotide-binding P-loop motif, are interpreted to reflect hydrogen-bond dynamics at the inhibitor-binding site. The hypothesis of energetic counter balancing of substrate binding based on crystallographic data is strongly supported by the solution NMR results. Correlations between backbone dynamics and domain displacement are established.     

Photoactivation Reduces Side-Chain Dynamics of a LOV Photoreceptor

We used neutron-scattering experiments to probe the conformational dynamics of the light, oxygen, voltage (LOV) photoreceptor PpSB1-LOV from Pseudomonas putida in both the dark and light states. Global protein diffusion and internal macromolecular dynamics were measured using incoherent neutron time-of-flight and backscattering spectroscopy on the picosecond to nanosecond timescales. Global protein diffusion of PpSB1-LOV is not influenced by photoactivation. Observation-time-dependent global diffusion coefficients were found, which converge on the nanosecond timescale toward diffusion coefficients determined by dynamic light scattering. Mean-square displacements of localized internal motions and effective force constants, <k′>, describing the resilience of the proteins were determined on the respective timescales. Photoactivation significantly modifies the flexibility and the resilience of PpSB1-LOV. On the fast, picosecond timescale, small changes in the mean-square displacement and <k′> are observed, which are enhanced on the slower, nanosecond timescale. Photoactivation results in a slightly larger resilience of the photoreceptor on the fast, picosecond timescale, whereas in the nanosecond range, a significantly less resilient structure of the light-state protein is observed. For a residue-resolved interpretation of the experimental neutron-scattering data, we analyzed molecular dynamics simulations of the PpSB1-LOV X-ray structure. Based on these data, it is tempting to speculate that light-induced changes in the protein result in altered side-chain mobility mostly for residues on the protruding Jα helix and on the LOV-LOV dimer interface. Our results provide strong experimental evidence that side-chain dynamics play a crucial role in photoactivation and signaling of PpSB1-LOV via modulation of conformational entropy.     

Excitation dynamics in Photosystem I from Chlamydomonas reinhardtii. Comparative studies of isolated complexes and whole cells

Identical time-resolved fluorescence measurements with ~ 3.5-ps resolution were performed for three types of PSI preparations from the green alga, Chlamydomonas reinhardtii: isolated PSI cores, isolated PSI–LHCI complexes and PSI–LHCI complexes in whole living cells. Fluorescence decay in these types of PSI preparations has been previously investigated but never under the same experimental conditions. As a result we present consistent picture of excitation dynamics in algal PSI. Temporal evolution of fluorescence spectra can be generally described by three decay components with similar lifetimes in all samples (6–8 ps, 25–30 ps, 166–314 ps). In the PSI cores, the fluorescence decay is dominated by the two fastest components (~ 90%), which can be assigned to excitation energy trapping in the reaction center by reversible primary charge separation. Excitation dynamics in the PSI–LHCI preparations is more complex because of the energy transfer between the LHCI antenna system and the core. The average trapping time of excitations created in the well coupled LHCI antenna system is about 12–15 ps longer than excitations formed in the PSI core antenna. Excitation dynamics in PSI–LHCI complexes in whole living cells is very similar to that observed in isolated complexes. Our data support the view that chlorophylls responsible for the long-wavelength emission are located mostly in LHCI. We also compared in detail our results with the literature data obtained for plant PSI.     

Energy transfer in the peridinin-chlorophyll protein complex reconstituted with mixed chlorophyll sites

We use femtosecond transient absorption spectroscopy to study chlorophyll (Chl)-Chl energy transfer in the peridinin-chlorophyll protein (PCP) reconstituted with mixtures of either chlorophyll b (Chlb) and Chld or Chla and bacteriochlorophyll a (BChla). Analysis of absorption and transient absorption spectra demonstrated that reconstitution with chlorophyll mixtures produces a significant fraction of PCP complexes that contains a different Chl in each domain of the PCP monomer. The data also suggest that binding affinity of Chla is less than that of the other three Chl species. By exciting the Chl species lying at higher energy, we obtained energy transfer times of 40 ± 5 ps (Chlb-Chld) and 59 ± 3 ps (Chla-BChla). The experimental values match those obtained from the Förster equation, 36 and 50 ps, respectively, showing that energy transfer proceeds via the Förster mechanism. Excitation of peridinin in the PCP complex reconstituted with Chla/BChla mixture provided time constants of 2.6 and 0.4 ps for the peridinin-Chla and peridinin-BChla energy transfer, matching those obtained from studies of PCP complexes reconstituted with single chlorophyll species.