Translocation and interactions of L-arabinose in OmpF porin: A molecular dynamics study

The passage of a natural substrate, L-arabinose (L-ARA) through Escherichia coli porin embedded in an artificial bilayer, is studied by equilibrium molecular dynamics simulations. We investigate the early stage of translocation process of L-ARA from intra-cellular to extra-cellular side (Int-to-Ext) across the bilayer. The average trajectory path over all L-ARA molecules along with quantum-mechanical configuration-optimizations at PM3 level predict the existence of at least three trapping zones. The common feature within all these zones is that L-ARA remains perpendicular to the channel axis. It is remarkable how the orientation and translational-rotational motion of L-ARA molecule play a role in its transport through OmpF channel. These simulations are important for better understanding of permeation process in OmpF channel. They also provide an insight into the chiral recognition of translocation process in protein nanochannels from substrate and protein prospects and help interpret experiments on permeation process of small dipolar molecules across biological membranes.     

Brownian dynamics of interactions between aldolase mutants and F-actin

Previous Brownian dynamics (BD) simulations (Ouporov IG, Knull HR and Thomasson KA 1999. Biophys. J. 76: 17-27) of complex formation between rabbit aldolase and F-actin have identified three lysine residues (K288, K293 and K341) on aldolase and acidic residues (DEDE) at the N-terminus of actin as important to binding. BD simulations of computer models of aldolase mutants with any of these lysine residues replaced by alanine show reduced binding energy; the greatest effect of a single substitution is for K341A, and replacement of all three lysines greatly reduces binding. BD simulations of wild-type rabbit aldolase vs altered F-actin show that binding is decreased if any one of the four N-terminal acidic residues is replaced by alanine and binding is greatly reduced if three or more of the N-terminal acidic residues are replaced; none of the four actin residues appear more critical for binding than the others.     

Energetics of charge-charge interactions between residues adjacent in sequence

Statistical analysis of the residue separation between a pair of ionizable side chains within 4 ? of each other was performed on a set of 1560 non-homologous PDB structures. We found that the frequency of pairs of like charges (i.e., pairs consisting of acidic residues Asp and Glu or pairs consisting of basic residues Arg and Lys) is two orders of magnitude lower than the pairs of oppositely charged residues (salt-bridges). We also found that for pairs of like charges the distribution is skewed dramatically towards short residue separation (<3). On the basis of these observations, we hypothesize that at short residue separation the repulsion between charges does not contribute much to the protein stability and the effects are largely dominated by the long range charge-charge interactions with other ionizable groups in the protein molecule. To test this hypothesis, we incorporated various pairs of charged residues at position 63 and 64 of ubiquitin and compared the stabilities of these variants. We also performed calculations of the expected changes in the charge-charge interactions. A very good correlation between experimental changes in the stability of ubiquitin variants, and changes in the energy of charge-charge interactions provides support for the hypothesis that a pair of ionizable residues next to each other in sequence modulates protein stability via long range charge-charge interactions with the rest of the protein.     

Intragenic suppressors of an OmpF assembly mutant and assessment of the roles of various OmpF residues in assembly through informational suppressors

We employed two separate genetic approaches to examine the roles of various OmpF residues in assembly. In one approach, intragenic suppressors of a temperature-sensitive OmpF assembly mutant carrying a W214E substitution were sought at 42 degrees C, or at 37 degrees C in a genetic background lacking the periplasmic folding factor SurA. In the majority of cases (58 out of 61 revertants), the suppressors mapped either at the original site (position 214) or two residues downstream from it. In the remaining three revertants that were obtained in a surA background, an alteration of N230Y was located 16 residues away from the original site. The N230Y suppressor also corrected OmpF315 assembly at 42 degrees C in a surA(+) background, indicating that the two different physiological environments imposed similar assembly constraints. The specificity of N230Y was tested against five different residues at position 214 of mature OmpF. Clear specificity was displayed, with maximum suppression observed for the original substitution at position 214 (E214) against which the N230Y suppressor was isolated, and no negative effect on OmpF assembly was noted when the wild-type W214 residue was present. The mechanism of suppression may involve compensation for a specific conformational defect. The second approach involved the application of informational suppressors (Su-tRNA) in combination with ompF amber mutations to generate variant OmpF proteins. In this approach we targeted the Y40, Q66, W214, and Y231 residues of mature OmpF and replaced them with S, Q, L, and Y through the action of Su-tRNAs. Thus, a total of 16 variant OmpF proteins were generated, of which three were identical to the parental protein, and two variants carrying W214Q and Y231Q substitutions were similar to assembly-defective proteins isolated previously (R. Misra, J. Bacteriol. 175:5049-5056, 1993). The results obtained from these analyses provided useful information regarding the compatibility of various alterations in OmpF assembly.     

Conductance and selectivity fluctuations in D127 mutants of the bacterial porin OmpF

A recent molecular dynamics study questioned the protonation state and physiological role of aspartate 127 (D127) of E. coli porin OmpF. To address that question we isolated two OmpF mutants with D127 either neutralized (D127N) or replaced by a positively charged lysine (D127K). The charge state of the residue at position 127 has clear effects on both conductance and selectivity. The D127K but not the D127N mutant expresses resilient conductance and selectivity fluctuations. These fluctuations reflect, we think, either changes in the ionization state of K127 and/or transitions between unstable subconformations as induced by the electrostatic repulsion between two positively charged residues, K127 and the nearby R167. Our results slightly favor the view that in WT OmpF residue D127 is deprotonated. As for the role of D127 in OmpF functionality, the gating of both mutants shows very similar sensitivity toward voltage as WT OmpF. Moreover, the current fluctuations of the D127K mutant were observed also in the absence of an applied electric field. We therefore dismiss D127 as a key residue in the control mechanism of the voltage-dependent gating of OmpF.     

NMR and computational studies of interactions between remote residues in gangliosides

Conformational preferences of the gangliosides GM1, GM1b, and GD1a have been investigated by using a systematic combination of NMR distance constraints and molecular mechanics calculations. These gangliosides share a common four-sugar core but differ in the number or placement of sialic acid residues attached to the core. Placement of the sialic acid residues is shown to influence the preferred core conformation. The origin of these effects is postulated to be intramolecular interactions of the sialic acid residues with other remote residues. In the case of GM1, hydrogen bonds between the internal sialic acid and an N-acetyl group on GalNAc are suggested. In the case of GD1a, a hydrogen-bonding network between the terminal and internal sialic acids is suggested to play a role.     

The geometry of interactions between catalytic residues and their substrates

The process of deducing the catalytic mechanism of an enzyme from its structure is highly complex and requires extensive experimental work to validate a proposed mechanism. As one step towards improving the reliability of this process, we have gathered statistics describing the typical geometry of catalytic residues with regard to the substrate and one another. In order to analyse residue-substrate interactions, we have assembled a dataset of structures of enzymes of known mechanism bound to substrate, product, or a substrate analogue. Despite the challenges presented in obtaining such experimental data, we were able to include 42 enzyme structures. We have also assembled a separate dataset of catalytic residues which act upon other catalytic residues, using a set of 60 enzyme structures. For both datasets, we have extracted the distances between residues with a given catalytic function and their target moieties. The geometry of residues whose function involves the transfer or sharing of hydrogens (either with substrate or another residue) was analysed more closely. The results showed that the geometry for such productive interactions (prior to the transition state) closely resembles that seen in non-catalytic hydrogen bonds, with distances and angles in the normal expected range. Such statistics provide limits on "expected geometries" for catalytic residues, which will help to identify these residues and elucidate enzyme mechanisms.     

Spontaneous degradation of polypeptides at aspartyl and asparaginyl residues: effects of the solvent dielectric.

We have investigated the spontaneous degradation of aspartate and asparagine residues via succinimide intermediates in model peptides in organic co-solvents. We find that the rate of deamidation at asparagine residues is markedly reduced in solvents of low dielectric strength. Theoretical considerations suggest that this decrease in rate is due to the destabilization of the deprotonated peptide bond nitrogen anion that is the postulated attacking species in succinimide formation. This result suggests that asparagine residues in regions with low dielectric constants, such as the interior of a protein or in a membrane bilayer, are protected from this type of degradation reaction. On the other hand, we found little or no effect on the rate of succinimide-mediated isomerization of aspartate residues when subjected to the same changes in dielectric constant. In this case, the destabilization of the attacking peptide bond nitrogen anion may be balanced by increased protonation of the aspartyl side chain carboxyl group, a reaction that results in a superior leaving group. Consequently, any protein structure or conformation that would increase the protonation of an aspartate side chain carboxyl group can be expected to render that residue more labile. These results may help explain why particular aspartate residues have been found to degrade in proteins at rates comparable to those of asparagine residues, even though aspartyl-containing peptides degrade more slowly than corresponding asparaginyl-containing peptides in aqueous solutions.     

Dimer-tetramer transition between solution and crystalline states of streptavidin and avidin mutants

The biotin-binding tetrameric proteins, streptavidin from Streptomyces avidinii and chicken egg white avidin, are excellent models for the study of subunit-subunit interactions of a multimeric protein. Efforts are thus being made to prepare mutated forms of streptavidin and avidin, which would form monomers or dimers, in order to examine their effect on quaternary structure and assembly. In the present communication, we compared the crystal structures of binding site W-->K mutations in streptavidin and avidin. In solution, both mutant proteins are known to form dimers, but upon crystallization, both formed tetramers with the same parameters as the native proteins. All of the intersubunit bonds were conserved, except for the hydrophobic interaction between biotin and the tryptophan that was replaced by lysine. In the crystal structure, the binding site of the mutated apo-avidin contains 3 molecules of structured water instead of the 5 contained in the native protein. The lysine side chain extends in a direction opposite that of the native tryptophan, the void being partially filled by an adjacent lysine residue. Nevertheless, the binding-site conformation observed for the mutant tetramer is an artificial consequence of crystal packing that would not be maintained in the solution-phase dimer. It appears that the dimer-tetramer transition may be concentration dependent, and the interaction among subunits obeys the law of mass action.     

OmpF assembly mutants of Escherichia coli K-12: isolation, characterization, and suppressor analysis.

This paper describes a novel genetic method used to isolate mutations that alter proper assembly of OmpF in the outer membrane. The thermolabile nature of assembly intermediates allowed selection of temperature-sensitive mutations within the ompF gene. A variant allele of ompF (ompF-Dex) was used because it provided a convenient selectable phenotype (Dex+). Assembly mutants were isolated in two steps. First, amber mutations were obtained that mapped in ompF-Dex. This resulted in a Dex- phenotype. Starting with these Dex- strains, Dex+ revertants were isolated. Mutants that displayed a temperature-sensitive Dex+ phenotype were further characterized. Three such mutants possessed a single substitution within ompF that reverted the nonsense codon to a sense codon which replaced W214 with either an E or Q and Y231 with a Q residue in the mature OmpF protein. All three mutant OmpF proteins showed an assembly defect. This defect led to a substantial reduction in the amount of stable OmpF trimers with the concomitant increase of a high-molecular-weight form of OmpF which migrated at the top of the gel. Suppressor mutations were sought that corrected the assembly defect of OmpF. These extragenic suppressor mutations were mapped at 45 min on the Escherichia coli chromosome. The suppressor mutations displayed no allele specificity and were recessive to the wild-type allele. In the presence of a suppressor, mutant stable trimers appeared in an almost normal manner. The appearance of stable trimers concurred with a substantial loss of the high-molecular-weight OmpF species. At this stage, it is not clear whether the high-molecular-weight species of OmpF is a normal assembly intermediate or a dead-end assembly product. The results presented in this study raise the intriguing possibility of a chaperone-like activity for the wild-type suppressor gene product.     

Environmental effects on the protonation states of active site residues in bacteriorhodopsin.

Finite difference solutions of the Poisson-Boltzmann equation are used to calculate the pKa values of the functionally important ionizable groups in bacteriorhodopsin. There are strong charge-charge interactions between the residues in the binding site leading to the possibility of complex titration behavior. Structured water molecules, if they exist in the binding site, can have significant effects on the calculated pKa by strongly stabilizing ionized species. The ionization states of the Schiff base and Asp-85 are found to be strongly coupled. Small environmental changes, which might occur as a consequence of trans-cis isomerization, are capable of causing large shifts in the relative pKa values of these two groups. This provides an explanation for the protonation of Asp-85 and the deprotonation of the Schiff base in the M state of bacteriorhodopsin. The different behavior of Asp-85 and Asp-212 is discussed in this regard.     

Complex interactions between native and invasive fish: the simultaneous effects of multiple negative interactions

We suggest that the ultimate outcome of interactions between native species and invasive species (extinction or coexistence) depends on the number of simultaneous negative interactions (competition and predation), which depends on relative body sizes of the species. Multiple simultaneous interactions may constrain the ability of native species to trade fitness components (i.e., reduced growth for reduced risk of predation) causing a spiral to extinction. We found evidence for five types of interactions between the adults and juveniles of introduced western mosquitofish (Gambusia affinis) and the juveniles of native least chub (Iotichthys phlegethontis). We added ten large (23–28 mm) and seven small (9–13 mm) young-of-the-year (YOY) least chub to replicate enclosures with zero, low, and high densities of mosquitofish in a desert spring ecosystem. Treatments with mosquitofish reduced the average survival of least chub by one-third. No small YOY least chub survived in enclosures with high mosquitofish densities. We also performed two laboratory experiments to determine mortality to predation, aggressiveness, and habitat selection of least chub in the presence of mosquitofish. Mean mortality of least chub due to predation by large mosquitofish was 69.7% over a 3-h trial. Least chub were less aggressive, selected protected habitats (Potamogeton spp.), and were more stationary in the presence of mosquitofish where the dominance hierarchy was large mosquitofish>>large least chubsmall mosquitofish>>small least chub. Least chub juveniles appear to be figuratively caught in a vice. Rapid growth to a size refuge could reduce the risk of predation, but the simultaneous effects of competition decreased least chub growth and prolonged the period when juveniles were vulnerable to mosquitofish predation.     

Ionization states of the complex formed between 2-benzyl-3-phosphonopropionic acid and carboxypeptidase A.

The binding to carboxypeptidase A of two phosphonic acid analogues of 2-benzylsuccinate, 2-DL-2-benzyl-3-phosphonopropionic acid (inhibitor I) and 2-DL-2-benzyl-3-(-O-ethylphosphono)propionic acid (inhibitor II) was studied by observing their 31P resonances when free and bound to the enzyme in the range of pH from 5 to 10. The binding of I by co-ordination to the active-site Zn(II) lowered the highest pKa of I from a value of 7.66(+/- 0.10) to a value of 6.71(+/- 0.17). No titration of any protons on II occurred over the pH range studied. The enzyme-bound inhibitor II also did not titrate over the pH range 6.17-7.60. The pH-dependencies of the apparent inhibition constants for I and II were also investigated by using N-(-2-(furanacryloyl)-L-phenylalanyl-L-phenylalanine as substrate. Two enzymic functional groups with pKa values of 5.90(+/- 0.06) and 9.79(+/- 0.14) must be protonated for binding of inhibitor I, and two groups with pKa values of 6.29(+/- 0.10) and 9.19(+/- 0.15) for binding of inhibitor II. Over the pH range from 6.71 to 7.66, inhibitor I binds to the enzyme in a complex of the enzyme in a more protonated form, and the inhibitor in a less protonated form than the predominant unligated forms at this pH. Mock & Tsay [(1986) Biochemistry 25, 2920-2927] made a similar finding for the binding of L-2-(1-carboxy-2-phenylethyl)-4-phenylazophenol over a pH range of nearly 4 units. The true inhibition constant for the dianionic form of inhibitor I (racemic) was calculated to be 54.0(+/- 5.9) nM and that of the trianionic form to be 5.92(+/- 0.65) nM. The true inhibition constant of the fully ionized II (racemic) was calculated to be 79.8(+/- 6.4) nM.     

Assembly of LamB and OmpF in deep rough lipopolysaccharide mutants of Escherichia coli K-12.

Assembly of the OmpF and LamB proteins was kinetically retarded in deep rough lipopolysaccharide mutants of Escherichia coli K-12. OmpF assembly was affected at the step of conversion of metastable trimers to stable trimers, whereas LamB assembly was influenced both at the monomer-to-metastable trimer and metastable-to-stable trimer steps. These assembly defects were reversed in the presence of the sfaA1 and sfaB3 suppressor alleles, which were isolated by using ompF assembly mutants.     

Resilience in a two-population system: interactions between Allee effects and connectivity

Resilience in ecosystems and resistance to regime shifts has been a major focus in ecological research. How migration and general network dynamics affect the resilience of populations or induce regime shift cascades is a particularly challenging open question in theoretical ecology. We focus on regime shifts in populations with variable-strength Allee effects to demonstrate the effect of migration on resilience in two-population systems with critical transitions. The result is a mathematical model that justifies the assumption that resilience can be averaged across connected populations and suggests several management strategies to either avoid or induce regime shift cascades.     

Ciprofloxacin interactions with bacterial protein OmpF: Modelling of FRET from a multi-tryptophan protein trimer

The outer membrane protein F (OmpF) is known to play an important role in the uptake of fluoroquinolone antibiotics by bacteria. In this study, the degree of binding of the fluoroquinolone antibiotic ciprofloxacin to OmpF in a lipid membrane environment is quantified using a methodology based on Förster resonance energy transfer (FRET). Analysis of the fluorescence quenching of OmpF is complex as each OmpF monomer presents two tryptophans at different positions, thus sensing two different distributions of acceptors in the bilayer plane. Specific FRET formalisms were derived accounting for the different energy transfer contributions to quenching of each type of tryptophan of OmpF, allowing the recovery of upper and lower boundaries for the ciprofloxacin-OmpF binding constant (K_B). log (K_B) was found to lie in the range 3.15-3.62 or 3.58-4.00 depending on the location for the ciprofloxacin binding site assumed in the FRET modelling, closer to the centre or to the periphery of the OmpF trimer, respectively. This methodology is suitable for the analysis of FRET data obtained with similar protein systems and can be readily adapted to different geometries.     

Stabilization of partially folded states of cytochrome c in aqueous surfactant: effects of ionic and hydrophobic interactions

The interaction of submicellar concentrations of sodium dodecyl sulfate (SDS) with horse heart cytochrome c has been found to stabilize two spectroscopically distinct partially folded intermediates at pH 7. The first intermediate is formed by the interaction of SDS with native cytochrome c, and this intermediate retains the majority of the secondary structure while the tertiary structure of the protein is lost. The unfolding of this intermediate with urea leads to the formation of a second intermediate, which is also formed on refolding of the unfolded protein (unfolded by urea) by SDS. The second intermediate retains about 50% of the native secondary structure with no tertiary structure of the protein. The second intermediate was found to be absent at low pH. While induction of helical structure of a protein by SDS in the native condition has been reported earlier, this is possibly the first report of the refolding of a protein in a strongly denaturing condition (in the presence of 10 M urea). The relative contributions of the hydrophobic and the electrostatic interactions of the surfactants with cytochrome c have been determined from the formation of the molten globule species from the acid-induced unfolded protein in the presence of SDS or lauryl maltoside.     

Selection for mutants altered in the expression or export of outer membrane porin OmpF.

Strains in which the lacZ gene (which specifies beta-galactosidase) is fused to a gene encoding an envelope protein often exhibit a phenotype termed overproduction lethality. In such strains, high-level synthesis of the cognate hybrid protein interferes with the process of protein export, and this leads ultimately to cell death. A variation of this phenomenon has been discovered with lacZ fusions to the gene specifying the major outer membrane porin protein OmpF. In this case, we find that lambda transducing phage carrying an ompF-lacZ fusion will not grow on a host strain that constitutively overexpresses ompF. We have exploited this observation to develop a selection for ompF mutants. Using this protocol, we have isolated mutants altered in ompF expression and have identified mutations that block OmpF export. Our results suggest that it should be possible to adapt this selection for use with other genes specifying exported proteins.     

Role of charged residues at the OmpF porin channel constriction probed by mutagenesis and simulation

The channel constriction of OmpF porin, a pore protein in the bacterial outer membrane, is highly charged due to the presence of three arginines (R42, R82, and R132) and two acidic residues (D113 and E117). The influence of these charges on ion conductance, ion selectivity, and voltage gating has been studied with mutants D113N/E117Q, R42A/R82A/R132A/D113N/E117Q, and V18K/G131K, which were designed to remove or add protein charge at the channel constriction. The crystal structures revealed no or only local changes compared to wild-type OmpF, thus allowing a comparative study. The single-channel conductance of the isosteric D113N/E117Q variant was found to be 2-fold reduced, and that of the pentuple mutant was 70% of the wild-type value, despite a considerably larger pore cross section. Ion selectivity was drastically altered by the mutations with cation/anion permeability ratios ranging from 1 to 12. Ion flow through these and eight other mutants, which have been characterized previously, was simulated by Brownian dynamics based on the detailed crystal structures. The calculated ion selectivity and relative channel conductance values agree well with the experimental data. This demonstrates that ion translocation through porin is mainly governed by pore geometry and charge, the two factors that are properly represented in the simulations.