Electrostatic steering and ionic tethering in enzyme–ligand binding: Insights from simulations

To bind at an enzyme’s active site, a ligand must diffuse or be transported to the enzyme’s surface, and, if the binding site is buried, the ligand must diffuse through the protein to reach it. Although the driving force for ligand binding is often ascribed to the hydrophobic effect, electrostatic interactions also influence the binding process of both charged and nonpolar ligands. First, electrostatic steering of charged substrates into enzyme active sites is discussed. This is of particular relevance for diffusion-influenced enzymes. By comparing the results of Brownian dynamics simulations and electrostatic potential similarity analysis for triose-phosphate isomerases, superoxide dismutases, and β-lactamases from different species, we identify the conserved features responsible for the electrostatic substrate-steering fields. The conserved potentials are localized at the active sites and are the primary determinants of the bimolecular association rates. Then we focus on a more subtle effect, which we will refer to as “ionic tethering.” We explore, by means of molecular and Brownian dynamics simulations and electrostatic continuum calculations, how salt links can act as tethers between structural elements of an enzyme that undergo conformational change upon substrate binding, and thereby regulate or modulate substrate binding. This is illustrated for the lipase and cytochrome P450 enzymes. Ionic tethering can provide a control mechanism for substrate binding that is sensitive to the electrostatic properties of the enzyme’s surroundings even when the substrate is nonpolar.

Regular and stochastic dynamics in the real quadratic family

We prove the Regulat or Stochastic Conjecture for the real quadratic family which asserts that almost every real quadratic map Pc, c ∈ [−2, 1/4], has either an attracting cycle or an absolutely continuous invariant measure.

Extracellular HIV-1 virus protein R causes a large inward current and cell death in cultured hippocampal neurons: Implications for AIDS pathology

The small HIV-1 accessory protein Vpr (virus protein R) is a multifunctional protein that is present in the serum and cerebrospinal fluid of AIDS patients. We previously showed that Vpr can form cation-selective ion channels across planar lipid bilayers, introducing the possibility that, if incorporated into the membranes of living cells, Vpr might form ion channels and consequently perturb the maintained ionic gradient. In this study, we demonstrate, by a variety of approaches, that Vpr added extracellularly to intact cells does indeed form ion channels. We use confocal laser scanning microscopy to examine the subcellular localization of fluorescently labeled Vpr. Plasmalemma depolarization and damage are examined using the anionic potential-sensitive dye bis(1,3-dibutylbarbituric acid) trimethine oxonol and propidium iodide (PI), respectively, and the effect of Vpr on whole-cell current is demonstrated directly by using the patch-clamp technique. We show that recombinant purified extracellular Vpr associates with the plasmalemma of hippocampal neurons to cause a large inward cation current and depolarization of the plasmalemma, eventually resulting in cell death. Thus, we demonstrate a physiological action of extracellular Vpr and present its mechanistic basis. These findings may have important implications for neuropathologies in AIDS patients who possess significant amounts of Vpr in the cerebrospinal fluid.

The use of site-directed fluorophore labeling and donor–donor energy migration to investigate solution structure and dynamics in proteins

The use of molecular genetics for introducing fluorescent molecules enables the use of donor–donor energy migration to determine intramolecular distances in a variety of proteins. This approach can be applied to examine the overall molecular dimensions of proteins and to investigate structural changes upon interactions with specific target molecules. In this report, the donor–donor energy migration method is demonstrated by experiments with the latent form of plasminogen activator inhibitor type 1. Based on the known x-ray structure of plasminogen activator inhibitor type 1, three positions forming the corners of a triangle were chosen. Double Cys substitution mutants (V106C-H185C, H185C-M266C, and M266C-V106C) and corresponding single substitution mutants (V106C, H185C, and M266C) were created and labeled with a sulfhydryl specific derivative of BODIPY (=the D molecule). The side lengths of this triangle were obtained from analyses of the experimental data. The analyses account for the local anisotropic order and rotational motions of the D molecules, as well as for the influence of a partial DD-labeling. The distances, as determined from x-ray diffraction, between the Cα-atoms of the positions V106C–H185C, H185C–M266C, and M266C–V106C were 60.9, 30.8, and 55.1 Å, respectively. These are in good agreement with the distances of 54 ± 4, 38 ± 3, and 55 ± 3 Å, as determined between the BODIPY groups attached via linkers to the same residues. Although the positions of the D-molecules and the Cα-atoms physically cannot coincide, there is a reasonable agreement between the methods.

Stability and dynamics in a hyperthermophilic protein with melting temperature close to 200°C

The rubredoxin protein from the hyperthermophilic archaebacterium Pyrococcus furiosus was examined by a hydrogen exchange method. Even though the protein does not exhibit reversible thermal unfolding, one can determine its stability parameters—free energy, enthalpy, entropy, and melting temperature—and also the distribution of stability throughout the protein, by using hydrogen exchange to measure the reversible cycling of the protein between native and unfolded states that occurs even under native conditions.

Vibrational population relaxation of carbon monoxide in the heme pocket of photolyzed carbonmonoxy myoglobin: Comparison of time-resolved mid-IR absorbance experiments and molecular dynamics simulations

The vibrational energy relaxation of carbon monoxide in the heme pocket of sperm whale myoglobin was studied by using molecular dynamics simulation and normal mode analysis methods. Molecular dynamics trajectories of solvated myoglobin were run at 300 K for both the δ- and ɛ-tautomers of the distal His-64. Vibrational population relaxation times of 335 ± 115 ps for the δ-tautomer and 640 ± 185 ps for the ɛ-tautomer were estimated by using the Landau–Teller model. Normal mode analysis was used to identify those protein residues that act as the primary “doorway” modes in the vibrational relaxation of the oscillator. Although the CO relaxation rates in both the ɛ- and δ-tautomers are similar in magnitude, the simulations predict that the vibrational relaxation of the CO is faster in the δ-tautomer with the distal His playing an important role in the energy relaxation mechanism. Time-resolved mid-IR absorbance measurements were performed on photolyzed carbonmonoxy hemoglobin (Hb13CO). From these measurements, a T1 time of 600 ± 150 ps was determined. The simulation and experimental estimates are compared and discussed.

Induced fit DNA recognition by a minor groove binding analogue of Hoechst 33258: fluctuations in DNA A tract structure investigated by NMR and molecular dynamics simulations

NMR analysis and molecular dynamics simulations of d(GGTAATTACC)2 and its complex with a tetrahydropyrimidinium analogue of Hoechst 33258 suggest that DNA minor groove recognition in solution involves a combination of conformational selection and induced fit, rather than binding to a preorganised site. Analysis of structural fluctuations in the bound and unbound states suggests that the degree of induced fit observed is primarily a consequence of optimising van der Waals contacts with the walls of the minor groove resulting in groove narrowing through: (i) changes in base step parameters, including increased helical twist and propeller twist; (ii) changes to the sugar–phosphate backbone conformation to engulf the bound ligand; (iii) suppression of bending modes at the TpA steps. In contrast, the geometrical arrangement of hydrogen bond acceptors on the groove floor appears to be relatively insensitive to DNA conformation (helical twist and propeller twist). We suggest that effective recognition of DNA sequences (in this case an A tract structure) appears to depend to a significant extent on the sequence being flexible enough to be able to adopt the geometrically optimal conformation compatible with the various binding interactions, rather than involving ‘lock and key’ recognition.

Arp2/3 complex and actin dynamics are required for actin-based mitochondrial motility in yeast

The Arp2/3 complex is implicated in actin polymerization-driven movement of Listeria monocytogenes. Here, we find that Arp2p and Arc15p, two subunits of this complex, show tight, actin-independent association with isolated yeast mitochondria. Arp2p colocalizes with mitochondria. Consistent with this result, we detect Arp2p-dependent formation of actin clouds around mitochondria in intact yeast. Cells bearing mutations in ARP2 or ARC15 genes show decreased velocities of mitochondrial movement, loss of all directed movement and defects in mitochondrial morphology. Finally, we observe a decrease in the velocity and extent of mitochondrial movement in yeast in which actin dynamics are reduced but actin cytoskeletal structure is intact. These results support the idea that the movement of mitochondria in yeast is actin polymerization driven and that this movement requires Arp2/3 complex.

Quantitative ER ↔ Golgi Transport Kinetics and Protein Separation upon Golgi Exit Revealed by Vesicular Integral Membrane Protein 36 Dynamics in Live CellsV⃞

To quantitatively investigate the trafficking of the transmembrane lectin VIP36 and its relation to cargo-containing transport carriers (TCs), we analyzed a C-terminal fluorescent-protein (FP) fusion, VIP36-SP-FP. When expressed at moderate levels, VIP36-SP-FP localized to the endoplasmic reticulum, Golgi apparatus, and intermediate transport structures, and colocalized with epitope-tagged VIP36. Temperature shift and pharmacological experiments indicated VIP36-SP-FP recycled in the early secretory pathway, exhibiting trafficking representative of a class of transmembrane cargo receptors, including the closely related lectin ERGIC53. VIP36-SP-FP trafficking structures comprised tubules and globular elements, which translocated in a saltatory manner. Simultaneous visualization of anterograde secretory cargo and VIP36-SP-FP indicated that the globular structures were pre-Golgi carriers, and that VIP36-SP-FP segregated from cargo within the Golgi and was not included in post-Golgi TCs. Organelle-specific bleach experiments directly measured the exchange of VIP36-SP-FP between the Golgi and endoplasmic reticulum (ER). Fitting a two-compartment model to the recovery data predicted first order rate constants of 1.22 ± 0.44%/min for ER → Golgi, and 7.68 ± 1.94%/min for Golgi → ER transport, revealing a half-time of 113 ± 70 min for leaving the ER and 1.67 ± 0.45 min for leaving the Golgi, and accounting for the measured steady-state distribution of VIP36-SP-FP (13% Golgi/87% ER). Perturbing transport with AlF4− treatment altered VIP36-SP-GFP distribution and changed the rate constants. The parameters of the model suggest that relatively small differences in the first order rate constants, perhaps manifested in subtle differences in the tendency to enter distinct TCs, result in large differences in the steady-state localization of secretory components.

Polygalacturonase-Mediated Solubilization and Depolymerization of Pectic Polymers in Tomato Fruit Cell Walls1 : Regulation by pH and Ionic Conditions

The hydrolysis of cell wall pectins by tomato (Lycopersicon esculentum) polygalacturonase (PG) in vitro is more extensive than the degradation affecting these polymers during ripening. We examined the hydrolysis of polygalacturonic acid and cell walls by PG isozyme 2 (PG2) under conditions widely adopted in the literature (pH 4.5 and containing Na+) and under conditions approximating the apoplastic environment of tomato fruit (pH 6.0 and K+ as the predominate cation). The pH optima for PG2 in the presence of K+ were 1.5 and 0.5 units higher for the hydrolysis of polygalacturonic acid and cell walls, respectively, compared with activity in the presence of Na+. Increasing K+ concentration stimulated pectin solubilization at pH 4.5 but had little influence at pH 6.0. Pectin depolymerization by PG2 was extensive at pH values from 4.0 to 5.0 and was further enhanced at high K+ levels. Oligomers were abundant products in in vitro reactions at pH 4.0 to 5.0, decreased sharply at pH 5.5, and were negligible at pH 6.0. EDTA stimulated PG-mediated pectin solubilization at pH 6.0 but did not promote oligomer production. Ca2+ suppressed PG-mediated pectin release at pH 4.5 yet had minimal influence on the proportional recovery of oligomers. Extensive pectin breakdown in processed tomato might be explained in part by cation- and low-pH-induced stimulation of PG and other wall-associated enzymes.

Domain structure and dynamics in the helical filaments formed by RecA and Rad51 on DNA

Both the bacterial RecA protein and the eukaryotic Rad51 protein form helical nucleoprotein filaments on DNA that catalyze strand transfer between two homologous DNA molecules. However, only the ATP-binding cores of these proteins have been conserved, and this same core is also found within helicases and the F1-ATPase. The C-terminal domain of the RecA protein forms lobes within the helical RecA filament. However, the Rad51 proteins do not have the C-terminal domain found in RecA, but have an N-terminal extension that is absent in the RecA protein. Both the RecA C-terminal domain and the Rad51 N-terminal domain bind DNA. We have used electron microscopy to show that the lobes of the yeast and human Rad51 filaments appear to be formed by N-terminal domains. These lobes are conformationally flexible in both RecA and Rad51. Within RecA filaments, the change between the “active” and “inactive” states appears to mainly involve a large movement of the C-terminal lobe. The N-terminal domain of Rad51 and the C-terminal domain of RecA may have arisen from convergent evolution to play similar roles in the filaments.

Intrinsic compressibility and volume compression in solvated proteins by molecular dynamics simulation at high pressure.

Constant pressure and temperature molecular dynamics techniques have been employed to investigate the changes in structure and volumes of two globular proteins, superoxide dismutase and lysozyme, under pressure. Compression (the relative changes in the proteins' volumes), computed with the Voronoi technique, is closely related with the so-called protein intrinsic compressibility, estimated by sound velocity measurements. In particular, compression computed with Voronoi volumes predicts, in agreement with experimental estimates, a negative bound water contribution to the apparent protein compression. While the use of van der Waals and molecular volumes underestimates the intrinsic compressibilities of proteins, Voronoi volumes produce results closer to experimental estimates. Remarkably, for two globular proteins of very different secondary structures, we compute identical (within statistical error) protein intrinsic compressions, as predicted by recent experimental studies. Changes in the protein interatomic distances under compression are also investigated. It is found that, on average, short distances compress less than longer ones. This nonuniform contraction underlines the peculiar nature of the structural changes due to pressure in contrast with temperature effects, which instead produce spatially uniform changes in proteins. The structural effects observed in the simulations at high pressure can explain protein compressibility measurements carried out by fluorimetric and hole burning techniques. Finally, the calculation of the proteins static structure factor shows significant shifts in the peaks at short wavenumber as pressure changes. These effects might provide an alternative way to obtain information concerning compressibilities of selected protein regions.

Ultra-fast excited state dynamics in green fluorescent protein: multiple states and proton transfer.

The green fluorescent protein (GFP) of the jellyfish Aequorea Victoria has attracted widespread interest since the discovery that its chromophore is generated by the autocatalytic, posttranslational cyclization and oxidation of a hexapeptide unit. This permits fusion of the DNA sequence of GFP with that of any protein whose expression or transport can then be readily monitored by sensitive fluorescence methods without the need to add exogenous fluorescent dyes. The excited state dynamics of GFP were studied following photo-excitation of each of its two strong absorption bands in the visible using fluorescence upconversion spectroscopy (about 100 fs time resolution). It is shown that excitation of the higher energy feature leads very rapidly to a form of the lower energy species, and that the excited state interconversion rate can be markedly slowed by replacing exchangeable protons with deuterons. This observation and others lead to a model in which the two visible absorption bands correspond to GFP in two ground-state conformations. These conformations can be slowly interconverted in the ground state, but the process is much faster in the excited state. The observed isotope effect suggests that the initial excited state process involves a proton transfer reaction that is followed by additional structural changes. These observations may help to rationalize and motivate mutations that alter the absorption properties and improve the photo stability of GFP.