Mapping Conformational Ensembles of Aβ Oligomers in Molecular Dynamics Simulations

Although the oligomers formed by Aβ peptides appear to be the primary cytotoxic species in Alzheimer's disease, detailed information about their structures appears to be lacking. In this article, we use exhaustive replica exchange molecular dynamics and an implicit solvent united-atom model to study the structural properties of Aβ monomers, dimers, and tetramers. Our analysis suggests that the conformational ensembles of Aβ dimers and tetramers are very similar, but sharply distinct from those sampled by the monomers. The key conformational difference between monomers and oligomers is the formation of β-structure in the oligomers occurring together with the loss of intrapeptide interactions and helix structure. Our simulations indicate that, independent of oligomer order, the Aβ aggregation interface is largely confined to the sequence region 10-23, which forms the bulk of interpeptide interactions. We show that the fractions of β structure computed in our simulations and measured experimentally are in good agreement.     

Association of neighboring β-strands of outer membrane protein A in lipid bilayers revealed by site-directed fluorescence quenching

We present a detailed study on the formation of neighboring β-strands during the folding of a monomeric integral membrane protein of the β-barrel type. β-Strand and β-barrel formations were investigated for the eight-stranded transmembrane domain of outer membrane protein A (OmpA) with single-tryptophan (W), single-cysteine (C) OmpA mutants. Based on the OmpA structure, W and C were introduced in two neighboring β-strands oriented toward the hydrocarbon core of the membrane. Replaced residue pairs were closer to either the periplasmic turns (named cis-side) or the outer loops (named trans-side) of the strand. W_nC_m OmpA mutants containing W at position n and C at position m along the polypeptide chain were labeled at the C by a nitroxyl spin label, which is a short-range fluorescence quencher. To monitor the association of neighboring β-strands, we determined the proximity between fluorescent W and labeled C in OmpA folding experiments by intramolecular fluorescence quenching. Formation of native β-strand contacts in folding experiments required the lipid membrane. Residues in the trans-side of strands β₁, β₂, and β₃, represented by mutants W₁₅C₃₅ (β₁β₂, trans) and W₅₇C₃₅ (β₃β₂, trans), reached close proximity prior to residues in the N(β₁)- and C(β₈)-terminal strands as examined for mutants W₁₅C₁₆₂ (β₁β₈, trans) and W₇C₁₇₀ (β₁β₈, cis). Tryptophan and cysteine converged slightly faster in W₁₅C₁₆₂ (β₁β₈, trans) than in W₇C₁₇₀ (β₁β₈, cis). The last folding step was observed for residues at the cis-ends of strands β₁ and β₂ for the mutant W₇C₄₃ (β₁β₂, cis). The data also demonstrate that the neighboring β-strands associate upon insertion into the hydrophobic core of the lipid bilayer.     

Glimpses of the molecular mechanisms of β2-microglobulin fibril formation in vitro: Aggregation on a complex energy landscape

β₂-microglobulin (β₂m) is a 99-residue protein that aggregates to form amyloid fibrils in dialysis-related amyloidosis. The protein provides a powerful model for exploration of the structural molecular mechanisms of fibril formation from a full-length protein in vitro. Fibrils have been assembled from β₂m under both low pH conditions, where the precursor is disordered, and at neutral pH where the protein is initially natively folded. Here we discuss the roles of sequence and structure in amyloid formation, the current understanding of the structural mechanisms of the early stages of aggregation of β₂m at both low and neutral pH, and the common and distinct features of these assembly pathways.     

Characterization of a novel dextransucrase from Weissella confusa isolated from sourdough

Weissella confusa and Weissella cibaria isolated from wheat sourdoughs produce, from sucrose, linear dextrans due to a single soluble dextransucrase. In this study, the first complete gene sequence encoding dextransucrase from a W. confusa strain (LBAE C39-2) along with the one from a W. cibaria strain (LBAE K39) were reported. Corresponding gene cloning was achieved using specific primers designed on the basis of the draft genome sequence of these species. Deduced amino acid sequence of W. confusa and W. cibaria dextransucrase revealed common structural features of the glycoside hydrolase family 70. Notably, the regions located in the vicinity of the catalytic triad (D, E, D) are highly conserved. However, comparison analysis also revealed that Weissella dextransucrases form a distinct phylogenetic group within glucansucrases of other lactic acid bacteria. We then cloned the W. confusa C39-2 dextransucrase gene and successfully expressed the mature corresponding enzyme in Escherichia coli. The purified recombinant enzyme rDSRC39-2 catalyzed dextran synthesis from sucrose with a K _m of 8.6 mM and a V _max of 20 μmol/mg/min. According to ¹H and ¹³C NMR analysis, the polymer is a linear class 1 dextran with 97.2 % α-(1→6) linkages and 2.8 % α-(1→3) branch linkages, similar to the one produced by W. confusa C39-2 strain. The enzyme exhibited optimum catalytic activity for temperatures ranging from 35 to 40 °C and a pH of 5.4 in 20 mM sodium acetate buffer. This novel dextransucrase is responsible for production of dextran with predominant α-(1→6) linkages that could find applications as food hydrocolloids.     

Construction of an ultra-high density consensus genetic map,and enhancement of the physical map from genome sequencing in <Emphasis Type="Italic">Lupinus angustifolius</Emphasis>

Key message

An ultra-high density genetic map containing 34,574 sequence-defined markers was developed in Lupinus angustifolius. Markers closely linked to nine genes of agronomic traits were identified. A physical map was improved to cover 560.5 Mb genome sequence.

Abstract

Lupin (Lupinus angustifolius L.) is a recently domesticated legume grain crop. In this study, we applied the restriction-site associated DNA sequencing (RADseq) method to genotype an F₉ recombinant inbred line population derived from a wild type × domesticated cultivar (W × D) cross. A high density linkage map was developed based on the W × D population. By integrating sequence-defined DNA markers reported in previous mapping studies, we established an ultra-high density consensus genetic map, which contains 34,574 markers consisting of 3508 loci covering 2399 cM on 20 linkage groups. The largest gap in the entire consensus map was 4.73 cM. The high density W × D map and the consensus map were used to develop an improved physical map, which covered 560.5 Mb of genome sequence data. The ultra-high density consensus linkage map, the improved physical map and the markers linked to genes of breeding interest reported in this study provide a common tool for genome sequence assembly, structural genomics, comparative genomics, functional genomics, QTL mapping, and molecular plant breeding in lupin.

     

Structural Determination of Aβ25-35 Micelles by Molecular Dynamics Simulations

Amyloid-β (Aβ) peptides and other amyloidogenic proteins can form a wide range of soluble oligomers of varied morphologies at the early aggregation stage, and some of these oligomers are biologically relevant to the pathogenesis of Alzheimer's disease. Spherical micelle-like oligomers have been often observed for many different types of amyloids. Here, we report a hybrid computational approach to systematically construct, search, optimize, and rank soluble micelle-like Aβ_25-35 structures with different side-chain packings at the atomic level. Simulations reveal for the first time, to our knowledge, that two Aβ micelles with antiparallel peptide organization and distinct surface hydrophobicity display high structural stability. Stable micelles experience a slow secondary structural transition from turn to α-helix. Energetic analysis coupled with computational mutagenesis reveals that van der Waals and solvation energies play a more pronounced role in stabilizing the micelles, whereas the electrostatic energies present a stable but minor energetic contribution to peptide assemblies. Modeled Aβ micelles with shapes and dimensions similar to those of experimentally derived spherical structures also provide detailed information about the roles of structural dynamics and transition in the formation of amyloid fibrils. The strong binding affinity of our micelles to antibodies implies that micelles may be a biologically relevant species.     

Recombinant Minimalist Spider Wrapping Silk Proteins Capable of Native-Like Fiber Formation

Spider silks are desirable biomaterials characterized by high tensile strength, elasticity, and biocompatibility. Spiders produce different types of silks for different uses, although dragline silks have been the predominant focus of previous studies. Spider wrapping silk, made of the aciniform protein (AcSp1), has high toughness because of its combination of high elasticity and tensile strength. AcSp1 in Argiope trifasciata contains a 200-aa sequence motif that is repeated at least 14 times. Here, we produced in E. coli recombinant proteins consisting of only one to four of the 200-aa AcSp1 repeats, designated W₁ to W₄. We observed that purified W₂, W₃ and W₄ proteins could be induced to form silk-like fibers by shear forces in a physiological buffer. The fibers formed by W₄ were ∼3.4 µm in diameter and up to 10 cm long. They showed an average tensile strength of 115 MPa, elasticity of 37%, and toughness of 34 J cm⁻³. The smaller W₂ protein formed fewer fibers and required a higher protein concentration to form fibers, whereas the smallest W₁ protein did not form silk-like fibers, indicating that a minimum of two of the 200-aa repeats was required for fiber formation. Microscopic examinations revealed structural features indicating an assembly of the proteins into spherical structures, fibrils, and silk-like fibers. CD and Raman spectral analysis of protein secondary structures suggested a transition from predominantly α-helical in solution to increasingly β-sheet in fibers.     

Logical analysis of the mechanism of protein folding III. Prediction of the strong long-range interactions.

It has been found that strong long-range interactions occur in regions having large β-structural potentials. As has been described previously (Nagano, 1974), interactions among regions having both helical and β-structural potentials (αβ-gaβ interactions) are also very important. Accordingly, an idea is presented in this paper that the relative stability of a protein conformation could be estimated by a relatively simple mathematical function of sequence and conformation. The function P(p,q) is called the non-energy part of pseudo-free energy, because minimization of the sum of P(p,q) and energy functions (cf. Levitt, 1974; Warme &; Scheraga, 1974) can be expected to lead to a plausible model of a protein. A merit of the function is that it can help us decide which way to go in manipulating a temporarily built model, e.g. towards a helix-rich protein or towards a β-structure-rich protein. The estimation of P(p,q) as an artificial potential does not use much computer time because only the co-ordinates of the β-carbon atoms (α-carbon atoms if the residue is Gly) are used. It is composed of terms of the long-range interactions P_L and short-range interactions P_S. The term P_L represents the relative strength of helix-helix interactions, helix-β-candidate interactions and β-candidate-β-candidate interactions. It is assumed that both helical and β-structural potentials can be measured as the differences between the predicted function for helix and β-structure, respectively, as defined previously (Nagano, 1973), and the corresponding largest values ever found. A hypothesis that two residues distantly separated in the primary sequence contribute less to the stability of the whole molecule is finally discarded because the true conformation of concanavalin A becomes very unstable compared with its false conformation folded like the main part of subtilisin. The parameters thus determined indicate that the helix-β-candidate interactions are almost as important as the β-candidate-β-candidate interactions. Both helix and loop prediction functions are combined to give the short-range interactions term, P_S, according to whether the region is really helical or not, and to whether it is really looped or not. The function P(p,q) can be used as a criterion for judging whether the predicted conformation is realistic or false, because the parameters can be adjusted to give, within limits, reasonable values of −10 kcal/residue for true conformations and higher than −5 kcal/residue for false conformations.As an application of the present theory of protein folding, the tertiary structure of bacteriophage T4 lysozyme is predicted and presented in Figure 1, prior to the X-ray structure becoming available.     

NMR evidence for forming highly populated helical conformations in the partially folded hNck2 SH3 domain

Recent studies of several proteins implied that the folding of β-proteins may follow a nonhierarchical mechanism in which two major transitions are essential, i.e., the collapse of a random coil to form a nonnative helical intermediate, followed by a transformation into the native β-structure. We report that the first hNck2 SH3 domain, assuming an all-β barrel in the native form, can be reversibly transformed into a stable and nonnative helical state by acid-unfolding. We also conducted extensive NMR and mutagenesis studies that led to two striking findings: 1), NMR analysis reveals that in the helical state formed at pH 2.0, the first and last β-strands in the native form become unstructured, whereas the rest is surprisingly converted into two highly populated helices with a significantly limited backbone motion; and 2), a conserved four-residue sequence is identified on the second β-strand, a mutation of which suddenly renders the SH3 domain into a helical state even at pH 6.5, with NMR conformational and dynamic properties highly similar to those of the wild-type at pH 2.0. This observation implies that the region might contribute key interactions to disrupt the helical state, and to facilitate a further transformation into the native SH3 fold in the second transition.     

Computer analyses of characteristic infrared bands of globular proteins.

Infrared spectra of myoglobin, ribonuclease, lysozyme, α-chymotrypsin, α-lactalbumin, and β-lactoglobulin A were obtained in deuterium oxide solution in units of absorbance versus wavenumber from 1340 to 1750 cm^?1. The spectra were resolved into Gaussian components by means of an iterative computer program. Resolved characteristic absorption peaks for the two infrared active amide I′ components of antiparallel chain-pleated sheets (β-structure) were obtained. The characteristic amide I′ peaks of α-helical regions and apparently unordered regions overlap in D₂O solution. Absorptivity values for the resolved β-structure peak around 1630 cm^?1 were estimated on the basis of the known structure of ribonuclease, lysozyme, and β-chymotrypsin. The β-structure content of β-lactoglobulin was estimated to be ca. 48% of α-lactalbumin ca. 18%, and of α_s-casein close to zero. The results are in general agreement with conclusions drawn from circular dichroism and optical rotatory dispersion studies.     

Tandem genetic duplications in Salmonella typhimurium: amplification of the histidine operon.

Bovine pancreatic phospholipase A₂ (M_r = 14,000) has been crystallized and its three-dimensional structure determined by X-ray diffraction analysis to a resolution of 2.4 Å. Three heavy-atom derivatives were used in the phase calculations with inclusion of the anomalous dispersion differences. The resulting electron density map allowed an easy and unambiguous tracing of the peptide chain. Two of the seven disulfide connections appeared to be different from what was suggested by the earlier chemical and structural work. The bovine phospholipase A₂ structure contains about 50% α-helix and 10% β-structure. The bovine enzyme structure was found to deviate substantially from the previously published porcine prophospholipase structure.     

The 1.8 Å Resolution Structure of Hevamine,a Plant Chitinase/Lysozyme,and Analysis of the Conserved Sequence and Structure Motifs of Glycosyl Hydrolase Family 18

The three-dimensional structure of hevamine, a plant enzyme with chitinase and lysozyme activity, has been refined at 1.8 Å resolution to an R-factor of 14.9% and a freeR-factor of 19.6%. The final model consists of all 273 amino acid residues and 206 ordered water molecules. Two non-prolinecis-peptides were identified, involving Phe32 and Trp255, both of which are implicated in substrate binding.Other glycosyl hydrolase family 18 proteins with known three-dimen sional structure are bacterial chitinase A, endo-β-N-acetylglucosaminidase F₁, endo-β-N-acetylglucosaminidase H, and the two plant proteins concanavalin B and narbonin, which have no known enzymatic activity. All these structures contain a (βα)₈barrel fold, with the two family 18 consensus regions roughly corresponding to the third and fourth barrel strands. This confirms the grouping of these proteins into family 18, which was only based on weak and local sequence similarity. The substrate specificity of the enzymes is determined by the loops following the barrel strands that form the substrate binding site. All enzymes have an aspartic acid and a glutamic acid residue in positions identical with Asp 125 and the catalytic Glu127 of hevamine. The lack of chitinase activity of concanavalin B and narbonin can be explained by the absence of one of these carboxylate groups, and by differences in the loops that form the substrate-binding cleft in hevamine.     

Modeling the Early Stages of Phase Separation in Disordered Elastin-like Proteins

Elastin-like proteins (ELPs) are known to undergo liquid-liquid phase separation reversibly above a concentration-dependent transition temperature. Previous studies suggested that, as temperature increases, ELPs experience an increased propensity for type II β-turns. However, how the ELPs behave below the phase transition temperature itself is still elusive. Here, we investigate the importance of β-turn formation during the early stages of ELP self-association. We examined the behavior of two ELPs, a 150-repeat construct that had been investigated previously (ELP[V₅G₃A₂-150] as well as a new 40-repeat construct (ELP40) suitable for nuclear magnetic resonance measurements. Structural analysis of ELP40 reveals a disordered conformation, and chemical shifts throughout the sequence are insensitive to changes in temperature over 20°C. However, a low population of β-turn conformation cannot be ruled out based on chemical shifts alone. To examine the structural consequences of β-turns in ELPs, a series of structural ensembles of ELP[V₅G₃A₂-150] were generated, incorporating differing amounts of β-turn bias throughout the chain. To mimic the early stages of the phase change, two monomers were paired, assuming preferential interaction at β-turn regions. This approach was justified by the observation that buried hydrophobic turns are commonly observed to interact in the Protein Data Bank. After dimerization, the ensemble-averaged hydrodynamic properties were calculated for each degree of β-turn bias, and the results were compared with analytical ultracentrifugation experiments at various temperatures. We find that the temperature dependence of the sedimentation coefficient $\left(s_{20,\text{w}}^{\text{o}}\right)$ can be reproduced by increasing the β-turn content in the structural ensemble. This analysis allows us to estimate the presence of β-turns and weak associations under experimental conditions. Because disordered proteins frequently exhibit weak biases in secondary structure propensity, these experimentally-driven ensemble calculations may complement existing methods for modeling disordered proteins generally.     

Backbone relaxation coupled to the ionization of internal groups in proteins: a self-guided Langevin dynamics study

Pathways of structural relaxation triggered by ionization of internal groups in staphylococcal nuclease (SNase) were studied through multiple self-guided Langevin dynamics (SGLD) simulations. Circular dichroism, steady-state Trp fluorescence, and nuclear magnetic resonance spectroscopy have shown previously that variants of SNase with internal Glu, Asp, and Lys at positions 66 or 92, and Arg at position 66, exhibit local reorganization or global unfolding when the internal ionizable group is charged. Except for Arg-66, these internal ionizable groups have unusual pK_a values and are neutral at physiological pH. The structural trends observed in the simulations are in general agreement with experimental observations. The I92D variant, which unfolds globally upon ionization of Asp-92, in simulations often exhibits extensive hydration of the protein core, and sometimes also significant perturbations of the β-barrel. In the crystal structure of the V66R variant, the β1 strand from the β-barrel is domain-swapped; in the simulations, the β1 strand is sometimes partially released. The V66K variant, which in solutions shows reorganization of six residues at the C-terminus of helix α1 and perturbations in the β-barrel structure, exhibits fraying of three residues of helix α1 in one simulation, and perturbations and partial unfolding of three β-strands in a few other simulations. In sharp contrast, very small structural changes were observed in simulations of the wild-type protein. The simulations indicate that charging of internal groups frequently triggers penetration of water into the protein interior. The pK_a values of Asp-92 and Arg-66 calculated with continuum methods on SGLD-relaxed structures reached the normal values in most simulations. Detailed analysis of accuracy and performance of SGLD demonstrates that SGLD outperforms LD in sampling of alternative protein conformations without loss of the accuracy and level of detail characteristic of regular LD.     

Design of 11-residue peptides with unusual biophysical properties: induced secondary structure in the absence of water

A series of oligopeptides with β-forming and adhesive properties, were synthesized and analyzed for adhesion shear strength, secondary structure, and association properties. The sequences contained related hydrophobic core segments varying in length from 5 to 12 residues and flanked by di- or tri-lysine segments. Three remarkable peptides consisting of just 11 residues with hydrophobic core sequences of FLIVI, IGSII, and IVIGS flanked by three lysine residues gave the highest dry adhesion shear strength and displayed unusual biophysical properties in the presence and absence of water. KKKFLIVIKKK had its highest adhesion strength at 2% (w/v) at pH 12.0 and showed the highest adhesion strength after exposure to water (water resistance). Both KKKIGSIIKKK and KKKIVIGSKKK, at 4% (w/v) at pH 12.0, displayed nearly identical dry shear strength values to that with the FLIVI core sequence. The peptide with IGSII core, however, displayed a lower water resistance and the latter, IVIGS, showed no water resistance, completely delaminating upon soaking in water. These are the smallest peptides with adhesive properties reported to date and show remarkable adhesion strength even at lower concentrations of 0.2% (w/v), which corresponds to 1.6 mM. The FLIVI containing peptide adopted a β-sheet secondary structure in water while the IGSII- and IVIGS-containing sequences folded similarly only in the absence of water. Analytical ultracentrifugation studies showed that when the FLIVI sequence adopts β-structure in aqueous solution, it associates into a large molecular weight assembly. The random coils of IGSII and IVIGS showed no tendency to associate at any pH.     

四川柏木人工林林下植被生物量与林分结构的关系

森林结构与林下植被生物量的关系是森林持续经营与森林碳计量监测的科学基础,但一直缺乏必要的研究。以四川柏木(Cupressus funebris)人工林为研究对象,揭示林下植被生物量(Wu)、灌木生物量(Ws)和草本生物量(Wh)与林分结构的关系,并试图构建区域性林下植被生物量估测的混合模型。结果表明:(1)乔、灌、草群体共12个结构因子中,灌木群体的平均基径(Ds)、盖度(Cs)、高度(Hs)、体积(Vs)与林下植被生物量关系更紧密,在林下植被生物量模型构建中更有效;(2)多模型拟合与比较表明,柏木林Ws最佳估算模型为Ws=0.0005V1.0411s(R2a=0.762,P0.001,n=40),而Wu的最佳估算模型为ln Wu=0.0158Hs+0.0111Cs-0.5358(R2a=0.695,P0.001,n=40),但对于Wh未能获得较为理想的估算模型(R2a0.410,P0.01,n=40);(3)林分密度(Du)整合进入多元线性模型提高了林下植被生物量的估测精度,ln Wu=a+b Du+c Hs+d Cs(R2a=0.721,P0.001,n=40)。研究为区域性林下生物量估测模型构建提供了新论据。     

Structure based inhibitor design targeting glycogen phosphorylase b. Virtual screening,synthesis, biochemical and biological assessment of novel N-acyl-β-d-glucopyranosylamines

Glycogen phosphorylase (GP) is a validated target for the development of new type 2 diabetes treatments. Exploiting the Zinc docking database, we report the in silico screening of 1888 N-acyl-β-d-glucopyranosylamines putative GP inhibitors differing only in their R groups. CombiGlide and GOLD docking programs with different scoring functions were employed with the best performing methods combined in a ‘consensus scoring’ approach to ranking of ligand binding affinities for the active site. Six selected candidates from the screening were then synthesized and their inhibitory potency was assessed both in vitro and ex vivo. Their inhibition constants’ values, in vitro, ranged from 5 to 377 μM while two of them were effective at causing inactivation of GP in rat hepatocytes at low μM concentrations. The crystal structures of GP in complex with the inhibitors were defined and provided the structural basis for their inhibitory potency and data for further structure based design of more potent inhibitors.     

Structural studies of triazole inhibitors with promising inhibitor effects against antibiotic resistance metallo-β-lactamases

Metallo-β-lactamases (MBLs) are an emerging cause of bacterial antibiotic resistance by hydrolysing all classes of β-lactams except monobactams, and the MBLs are not inhibited by clinically available serine-β-lactamase inhibitors. Two of the most commonly encountered MBLs in clinical isolates worldwide – the New Delhi metallo-β-lactamase (NDM-1) and the Verona integron-encoded metallo-β-lactamase (VIM-2) – are included in this study.A series of several NH-1,2,3-triazoles was prepared by a three-step protocol utilizing Banert cascade reaction as the key step. The inhibitor properties were evaluated in biochemical assays against the MBLs VIM-2, NDM-1 and GIM-1, and VIM-2 showed IC₅₀ values down to nanomolar range. High-resolution crystal structures of four inhibitors in complex with VIM-2 revealed hydrogen bonds from the triazole inhibitors to Arg228 and to the backbone of Ala231 or Asn233, along with hydrophobic interactions to Trp87, Phe61 and Tyr67. The inhibitors show reduced MIC in synergy assays with Pseudomonas aeruginosa and Escherichia coli strains harbouring VIM enzymes. The obtained results will be useful for further structural guided design of MBL inhibitors.     

Body condition (Wr) and reproductive potential of bighead and silver carp hybrids: Postzygotic selection in the Mississippi River Basin

Invasive bighead (Hypophthalmichthys nobilis) and silver carp (H. molitrix) are reproductively isolated in their native range, but form a bimodal, multigenerational hybrid swarm within the Mississippi River Basin (MRB). Despite observed F₁ hybrid superiority in experimental settings, effects of postzygotic selection on bighead and silver carp hybrids have not been tested in a natural system. Individual parent and hybrid genotypes were resolved at 57 species‐specific loci and used to evaluate postzygotic selection for body condition (W_r) and female reproductive potential (presence of spawning stage gonads and gonadosomatic index [GSI]) in the MRB during 2009–2011. Body condition in the Marseilles Reach, Illinois River declined with a decrease in species‐specific allele frequency from 1.0 to 0.4 for each species and early generation hybrids (F₁, F₂, and first‐generation backcross) had lower mean W_r than late generation hybrids (2nd+ generation backcrosses) and parentals. Proportions of stage IV and stage V (spawning stage) female gonads differed between bighead and silver carp, but not among parentals and their early and late generation hybrids within the MRB. Mean GSI values did not differ between parentals and hybrids. Because reproductive potential did not differ between hybrids and parentals, our results suggest that early generation hybrids occur in low frequency either as a factor of poor condition (W_r) and postreproductive survival, infrequent reproductive encounters by parental bighead and silver carp, or selection pressures acting on juvenile or immature life stages. Our results suggest that a combination of genetic and environmental factors may contribute to the postzygotic success of bighead and silver carp hybrids in the Mississippi River Basin.