Factor H related proteins modulate complement activation on kidney cells

Complement activation at a particular location is determined by the balance of activating and inhibitory proteins. Factor H is a key regulator of the alternative pathway of complement, and genetic or acquired impairments in Factor H are associated with glomerular injury. The human Factor H-related proteins (FHRs) comprise a family of five proteins that are structurally related to Factor H. Variations in the genes or expression levels of the FHRs are also associated with glomerular disease, although the mechanisms of glomerular protection/injury are incompletely understood. To explore the role of the FHRs on complement regulation/dysregulation in the kidney, we expressed and purified recombinant murine FHRs (FHRs A, B, C and E). These four distinct FHRs contain binding regions with high amino acid sequence homology to binding regions within Factor H, but we observed different interactions of the FHRs with Factor H binding ligands, including heparin and C3d. There was differential binding of the FHRs to the resident kidney cell types (mesangial, glomerular endothelial, podocytes, and tubular epithelial). All four FHRs caused complement dysregulation on kidney cell surfaces in vitro, although the magnitude of the effect differed among the FHRs and also varied among the different kidney cells. However, only FHR E caused glomerular complement dysregulation when injected in vivo but did not exacerbate injury when injected into mice with ischemic acute kidney injury, an alternative pathway-mediated model. Thus, our experiments demonstrate that the FHRs have unique, and likely context-dependent, effects on the different cell types within the kidney.

HIV-1 Neutralizing Antibody Signatures and Application to Epitope-Targeted Vaccine Design

Eliciting HIV-1-specific broadly neutralizing antibodies (bNAbs) remains a challenge for vaccine development, and the potential of passively delivered bNAbs for prophylaxis and therapeutics is being explored. We used neutralization data from four large virus panels to comprehensively map viral signatures associated with bNAb sensitivity, including amino acids, hypervariable region characteristics, and clade effects across four different classes of bNAbs. The bNAb signatures defined for the variable loop 2 (V2) epitope region of HIV-1 Env were then employed to inform immunogen design in a proof-of-concept exploration of signature-based epitope targeted (SET) vaccines. V2 bNAb signature-guided mutations were introduced into Env 459C to create a trivalent vaccine, and immunization of guinea pigs with V2-SET vaccines resulted in increased breadth of NAb responses compared with Env 459C alone. These data demonstrate that bNAb signatures can be utilized to engineer HIV-1 Env vaccine immunogens capable of eliciting antibody responses with greater neutralization breadth.

•

HIV-1 bNAb sensitivity signatures from 4 large virus panels mapped across 4 Ab classes

•

Non-contact hypervariable region characteristics are critical for bNAb sensitivity

•

HIV-1 Env 459C used alone as a vaccine can elicit modest tier 2 NAbs in guinea pigs

•

V2 bNAb signature-guided modifications in 459C enhanced neutralization breadth

An assessment of commonly employed satellite-based remote sensors for mapping mangrove species in Mexico using an NDVI-based classification scheme

Optimizing the classification accuracy of a mangrove forest is of utmost importance for conservation practitioners. Mangrove forest mapping using satellite-based remote sensing techniques is by far the most common method of classification currently used given the logistical difficulties of field endeavors in these forested wetlands. However, there is now an abundance of options from which to choose in regards to satellite sensors, which has led to substantially different estimations of mangrove forest location and extent with particular concern for degraded systems. The objective of this study was to assess the accuracy of mangrove forest classification using different remotely sensed data sources (i.e., Landsat-8, SPOT-5, Sentinel-2, and WorldView-2) for a system located along the Pacific coast of Mexico. Specifically, we examined a stressed semiarid mangrove forest which offers a variety of conditions such as dead areas, degraded stands, healthy mangroves, and very dense mangrove island formations. The results indicated that Landsat-8 (30 m per pixel) had  the lowest overall accuracy at 64% and that WorldView-2 (1.6 m per pixel) had the highest at 93%. Moreover, the SPOT-5 and the Sentinel-2 classifications (10 m per pixel) were very similar having accuracies of 75 and 78%, respectively. In comparison to WorldView-2, the other sensors overestimated the extent of Laguncularia racemosa and underestimated the extent of Rhizophora mangle. When considering such type of sensors, the higher spatial resolution can be particularly important in mapping small mangrove islands that often occur in degraded mangrove systems.

HIV-1 ENVELOPE. Effect of the cytoplasmic domain on antigenic characteristics of HIV-1 envelope glycoprotein

A major goal for HIV-1 vaccine development is the production of an immunogen to mimic native, functional HIV-1 envelope trimeric spikes (Env) on the virion surface. We lack a reliable description of a native, functional trimer, however, because of inherent instability and heterogeneity in most preparations. We describe here two conformationally homogeneous Envs derived from difficult-to-neutralize primary isolates. All their non-neutralizing epitopes are fully concealed and independent of their proteolytic processing. Most broadly neutralizing antibodies (bnAbs) recognize these native trimers. Truncation of their cytoplasmic tail has little effect on membrane fusion, but it diminishes binding to trimer-specific bnAbs while exposing non-neutralizing epitopes. These results yield a more accurate antigenic picture than hitherto possible of a genuinely untriggered and functional HIV-1 Env; they can guide effective vaccine development.

Design of lipid nanocapsule delivery vehicles for multivalent display of recombinant Env trimers in HIV vaccination

Immunization strategies that elicit antibodies capable of neutralizing diverse virus strains will likely be an important part of a successful vaccine against HIV. However, strategies to promote robust humoral responses against the native intact HIV envelope trimer structure are lacking. We recently developed chemically cross-linked lipid nanocapsules as carriers of molecular adjuvants and encapsulated or surface-displayed antigens, which promoted follicular helper T-cell responses and elicited high-avidity, durable antibody responses to a candidate malaria antigen. To apply this system to the delivery of HIV antigens, Env gp140 trimers with terminal his-tags (gp140T-his) were anchored to the surface of lipid nanocapsules via Ni-NTA-functionalized lipids. Initial experiments revealed that the large (409 kDa), heavily glycosylated trimers were capable of extracting fluid phase lipids from the membranes of nanocapsules. Thus, liquid-ordered and/or gel-phase lipid compositions were required to stably anchor trimers to the particle membranes. Trimer-loaded nanocapsules combined with the clinically relevant adjuvant monophosphoryl lipid A primed high-titer antibody responses in mice at antigen doses ranging from 5 μg to as low as 100 ng, whereas titers dropped more than 50-fold over the same dose range when soluble trimer was mixed with a strong oil-in-water adjuvant comparator. Nanocapsule immunization also broadened the number of distinct epitopes on the HIV trimer recognized by the antibody response. These results suggest that nanocapsules displaying HIV trimers in an oriented, multivalent presentation can promote key aspects of the humoral response against Env immunogens.

Detection of complement activation using monoclonal antibodies against C3d

During complement activation the C3 protein is cleaved, and C3 activation fragments are covalently fixed to tissues. Tissue-bound C3 fragments are a durable biomarker of tissue inflammation, and these fragments have been exploited as addressable binding ligands for targeted therapeutics and diagnostic agents. We have generated cross-reactive murine monoclonal antibodies against human and mouse C3d, the final C3 degradation fragment generated during complement activation. We developed 3 monoclonal antibodies (3d8b, 3d9a, and 3d29) that preferentially bind to the iC3b, C3dg, and C3d fragments in solution, but do not bind to intact C3 or C3b. The same 3 clones also bind to tissue-bound C3 activation fragments when injected systemically. Using mouse models of renal and ocular disease, we confirmed that, following systemic injection, the antibodies accumulated at sites of C3 fragment deposition within the glomerulus, the renal tubulointerstitium, and the posterior pole of the eye. To detect antibodies bound within the eye, we used optical imaging and observed accumulation of the antibodies within retinal lesions in a model of choroidal neovascularization (CNV). Our results demonstrate that imaging methods that use these antibodies may provide a sensitive means of detecting and monitoring complement activation-associated tissue inflammation.

Delineation of the complement receptor type 2-C3d complex by site-directed mutagenesis and molecular docking

The interactions between the complement receptor type 2 (CR2) and the C3 complement fragments C3d, C3dg, and iC3b are essential for the initiation of a normal immune response. A crystal-derived structure of the two N-terminal short consensus repeat (SCR1-2) domains of CR2 in complex with C3d has previously been elucidated. However, a number of biochemical and biophysical studies targeting both CR2 and C3d appear to be in conflict with these structural data. Previous mutagenesis and heteronuclear NMR spectroscopy studies directed toward the C3d-binding site on CR2 have indicated that the CR2-C3d cocrystal structure may represent an encounter/intermediate or nonphysiological complex. With regard to the CR2-binding site on C3d, mutagenesis studies by Isenman and coworkers [Isenman, D. E., Leung, E., Mackay, J. D., Bagby, S. & van den Elsen, J. M. H. (2010). Mutational analyses reveal that the staphylococcal immune evasion molecule Sbi and complement receptor 2 (CR2) share overlapping contact residues on C3d: Implications for the controversy regarding the CR2/C3d cocrystal structure. J. Immunol. 184, 1946-1955] have implicated an electronegative "concave" surface on C3d in the binding process. This surface is discrete from the CR2-C3d interface identified in the crystal structure. We generated a total of 18 mutations targeting the two (X-ray crystallographic- and mutagenesis-based) proposed CR2 SCR1-2 binding sites on C3d. Using ELISA analyses, we were able to assess binding of mutant forms of C3d to CR2. Mutations directed toward the concave surface of C3d result in substantially compromised CR2 binding. By contrast, targeting the CR2-C3d interface identified in the cocrystal structure and the surrounding area results in significantly lower levels of disruption in binding. Molecular modeling approaches used to investigate disparities between the biochemical data and the X-ray structure of the CR2-C3d cocrystal result in highest-scoring solutions in which CR2 SCR1-2 is docked within the concave surface of C3d.

Biophysical investigations of complement receptor 2 (CD21 and CR2)-ligand interactions reveal amino acid contacts unique to each receptor-ligand pair

Human complement receptor type 2 (CR2 and CD21) is a cell membrane receptor, with 15 or 16 extracellular short consensus repeats (SCRs), that promotes B lymphocyte responses and bridges innate and acquired immunity. The most distally located SCRs, SCR1-2, mediate the interaction of CR2 with its four known ligands (C3d, EBV gp350, IFNalpha, and CD23). To ascertain specific interacting residues on CR2, we utilized NMR studies wherein gp350 and IFNalpha were titrated into (15)N-labeled SCR1-2, and chemical shift changes indicative of specific inter-molecular interactions were identified. With backbone assignments made, the chemical shift changes were mapped onto the crystal structure of SCR1-2. With regard to gp350, the binding region of CR2 is primarily focused on SCR1 and the inter-SCR linker, specifically residues Asn(11), Arg(13), Ala(22), Arg(28), Ser(32), Arg(36), Lys(41), Lys(57), Tyr(64), Lys(67), Tyr(68), Arg(83), Gly(84), and Arg(89). With regard to IFNalpha, the binding is similar to the CR2-C3d interaction with specific residues being Arg(13), Tyr(16), Arg(28), Ser(42), Lys(48), Lys(50), Tyr(68), Arg(83), Gly(84), and Arg(89). We also report thermodynamic properties of each ligand-receptor pair determined using isothermal titration calorimetry. The CR2-C3d interaction was characterized as a two-mode binding interaction with K(d) values of 0.13 and 160 microm, whereas the CR2-gp350 and CR2-IFNalpha interactions were characterized as single site binding events with affinities of 0.014 and 0.035 microm, respectively. The compilation of chemical binding maps suggests specific residues on CR2 that are uniquely important in each of these three binding interactions.

Evaluating the condition of a mangrove forest of the Mexican Pacific based on an estimated leaf area index mapping approach

Given the alarming global rates of mangrove forest loss it is important that resource managers have access to updated information regarding both the extent and condition of their mangrove forests. Mexican mangroves in particular have been identified as experiencing an exceptional high annual rate of loss. However, conflicting studies, using remote sensing techniques, of the current state of many of these forests may be hindering all efforts to conserve and manage what remains. Focusing on one such system, the Teacapán-Agua Brava-Las Haciendas estuarine-mangrove complex of the Mexican Pacific, an attempt was made to develop a rapid method of mapping the current condition of the mangroves based on estimated LAI. Specifically, using an AccuPAR LP-80 Ceptometer, 300 indirect in situ LAI measurements were taken at various sites within the black mangrove (Avicennia germinans) dominated forests of the northern section of this system. From this sample, 225 measurements were then used to develop linear regression models based on their relationship with corresponding values derived from QuickBird very high resolution optical satellite data. Specifically, regression analyses of the in situ LAI with both the normalized difference vegetation index (NDVI) and the simple ration (SR) vegetation index revealed significant positive relationships [LAI versus NDVI (R (2) = 0.63); LAI versus SR (R (2) = 0.68)]. Moreover, using the remaining sample, further examination of standard errors and of an F test of the residual variances indicated little difference between the two models. Based on the NDVI model, a map of estimated mangrove LAI was then created. Excluding the dead mangrove areas (i.e. LAI = 0), which represented 40% of the total 30.4 km(2) of mangrove area identified in the scene, a mean estimated LAI value of 2.71 was recorded. By grouping the healthy fringe mangrove with the healthy riverine mangrove and by grouping the dwarf mangrove together with the poor condition mangrove, mean estimated LAI values of 4.66 and 2.39 were calculated, respectively. Given that the former healthy group only represents 8% of the total mangrove area examined, it is concluded that this mangrove system, considered one of the most important of the Pacific coast of the Americas, is currently experiencing a considerable state of degradation. Furthermore, based on the results of this investigation it is suggested that this approach could provide resource managers and scientists alike with a very rapid and effective method for monitoring the state of remaining mangrove forests of the Mexican Pacific and, possibly, other areas of the tropics.

Mapping of the C3d ligand binding site on complement receptor 2 (CR2/CD21) using nuclear magnetic resonance and chemical shift analysis

Complement receptor 2 (CR2, CD21) is a cell membrane protein, with 15 or 16 extracellular short consensus repeats (SCRs), that promotes B lymphocyte responses and bridges innate and acquired immunity. The most distally located SCRs (SCR1-2) mediate the interaction of CR2 with its four known ligands (C3d, Epstein-Barr virus gp350, interferon-alpha, and CD23). Inhibitory monoclonal antibodies against SCR1-2 block binding of all ligands. To develop ligand-specific inhibitors that would also assist in identifying residues unique to each receptor-ligand interaction, phage were selected from randomly generated libraries by panning with recombinant SCR1-2, followed by specific ligand-driven elution. Derived peptides were tested by competition ELISA. One peptide, C3dp1 (APQHLSSQYSRT) exhibited ligand-specific inhibition at midmicromolar IC(50). C3d was titrated into (15)N-labeled SCR1-2, which revealed chemical shift changes indicative of specific intermolecular interactions. With backbone assignments made, the chemical shift changes were mapped onto the crystal structure of SCR1-2. With regard to C3d, the binding surface includes regions of SCR1, SCR2, and the inter-SCR linker, specifically residues Arg(13), Tyr(16), Arg(28), Tyr(29), Ser(32), Thr(34), Lys(48), Asp(56), Lys(57), Tyr(68), Arg(83), Gly(84), Asn(101), Asn(105), and Ser(109). SCR1 and SCR2 demonstrated distinct binding modes. The CR2 binding surface incorporating SCR1 is inconsistent with a previous x-ray CR2-C3d co-crystal analysis but consistent with mutagenesis, x-ray neutron scattering, and inhibitory monoclonal antibody epitope mapping. Titration with C3dp1 yielded chemical shift changes (Arg(13), Tyr(16), Thr(34), Lys(48), Asp(56), Lys(57), Tyr(68), Arg(83), Gly(84), Asn(105), and Ser(109)) overlapping with C3d, indicating that C3dp1 interacts at the same CR2 site as C3d.

Intrinsic amino acid side-chain hydrophilicity/hydrophobicity coefficients determined by reversed-phase high-performance liquid chromatography of model peptides: comparison with other hydrophilicity/hydrophobicity scales

An accurate determination of the intrinsic hydrophilicity/hydrophobicity of amino acid side-chains in peptides and proteins is fundamental in understanding many area of research, including protein folding and stability, peptide and protein function, protein-protein interactions and peptide/protein oligomerization, as well as the design of protocols for purification and characterization of peptides and proteins. Our definition of intrinsic hydrophilicity/hydrophobicity of side-chains is the maximum possible hydrophilicity/hydrophobicity of side-chains in the absence of any nearest-neighbor effects and/or any conformational effects of the polypeptide chain that prevent full expression of side-chain hydrophilicity/hydrophobicity. In this review, we have compared an experimentally derived intrinsic side-chain hydrophilicity/hydrophobicity scale generated from RP-HPLC retention behavior of de novo designed synthetic model peptides at pH 2 and pH 7 with other RP-HPLC-derived scales, as well as scales generated from classic experimental and calculation-based methods of octanol/water partitioning of Nalpha-acetyl-amino-acid amides or free energy of transfer of free amino acids. Generally poor correlation was found with previous RP-HPLC-derived scales, likely due to the random nature of the peptide mixtures in terms of varying peptide size, conformation and frequency of particular amino acids. In addition, generally poor correlation with the classical approaches served to underline the importance of the presence of a polypeptide backbone when generating intrinsic values. We have shown that the intrinsic scale determined here is in full agreement with the structural characteristics of amino acid side-chains.

HPLC analysis and purification of peptides

High-performance liquid chromatography (HPLC) has proved extremely versatile over the past 25 yr for the isolation and purification of peptides varying widely in their sources, quantity and complexity. This article covers the major modes of HPLC utilized for peptides (size-exclusion, ion-exchange, and reversed-phase), as well as demonstrating the potential of a novel mixed-mode hydrophilic interaction/cation-exchange approach developed in this laboratory. In addition to the value of these HPLC modes for peptide separations, the value of various HPLC techniques for structural characterization of peptides and proteins will be addressed, e.g., assessment of oligomerization state of peptides/proteins by size-exclusion chromatography and monitoring the hydrophilicity/hydrophobicity of amphipathic alpha-helical peptides, a vital precursor for the development of novel antimicrobial peptides. The value of capillary electrophoresis for peptide separations is also demonstrated. Preparative reversed-phase chromatography purification protocols for sample loads of up to 200 mg on analytical columns and instrumentation are introduced for both peptides and recombinant proteins.

Requirements for prediction of peptide retention time in reversed-phase high-performance liquid chromatography: hydrophilicity/hydrophobicity of side-chains at the N- and C-termini of peptides are dramatically affected by the end-groups and location

The value of reversed-phase high-performance liquid chromatography (RP-HPLC) and the field of proteomics would be greatly enhanced by accurate prediction of retention times of peptides of known composition. The present study investigates the hydrophilicity/hydrophobicity of amino acid side-chains at the N- and C-termini of peptides while varying the functional end-groups at the termini. We substituted all 20 naturally occurring amino acids at the N- and C-termini of a model peptide sequence, where the functional end-groups were N(alpha)-acetyl-X- and N(alpha)-amino-X- at the N-terminus and -X-C(alpha)-carboxyl and -X-C(alpha)-amide at the C-terminus. Amino acid coefficients were subsequently derived from the RP-HPLC retention behaviour of these peptides and compared to each other as well as to coefficients determined in the centre of the peptide chain (internal coefficients). Coefficients generated from residues substituted at the C-terminus differed most (between the -X-C(alpha)-carboxyl and -X-C(alpha)-amide peptide series) for hydrophobic side-chains. A similar result was seen for the N(alpha)-acetyl-X- and N(alpha)-amino-X- peptide series, where the largest differences in coefficient values were observed for hydrophobic side-chains. Coefficients derived from substitutions at the C-terminus for hydrophobic amino acids were dramatically different compared to internal coefficients for hydrophobic side-chains, ranging from 17.1 min for Trp to 4.8 min for Cys. In contrast, coefficients derived from substitutions at the N-terminus showed relatively small differences from the internal coefficients. Subsequent prediction of peptide retention time, within an error of just 0.4 min, was achieved by a predictive algorithm using a combination of internal coefficients and coefficients for the C-terminal residues. For prediction of peptide retention time, the sum of the coefficients must include internal and terminal coefficients.

Quantitation of the nearest-neighbour effects of amino acid side-chains that restrict conformational freedom of the polypeptide chain using reversed-phase liquid chromatography of synthetic model peptides with L- and D-amino acid substitutions

Side-chain backbone interactions (or "effects") between nearest neighbours may severely restrict the conformations accessible to a polypeptide chain and thus represent the first step in protein folding. We have quantified nearest-neighbour effects (i to i+1) in peptides through reversed-phase liquid chromatography (RP-HPLC) of model synthetic peptides, where L- and D-amino acids were substituted at the N-terminal end of the peptide sequence, adjacent to a L-Leu residue. These nearest-neighbour effects (expressed as the difference in retention times of L- and D-peptide diastereomers at pHs 2 and 7) were frequently dramatic, depending on the type of side-chain adjacent to the L-Leu residue, albeit such effects were independent of mobile phase conditions. No nearest-neighbour effects were observed when residue i is adjacent to a Gly residue. Calculation of minimum energy conformations of selected peptides supported the view that, whether a L- or D-amino acid is substituted adjacent to L-Leu, its orientation relative to this bulky Leu side-chain represents the most energetically favourable configuration. We believe that such energetically favourable, and different, configurations of L- and D-peptide diastereomers affect their respective interactions with a hydrophobic stationary phase, which are thus quantified by different RP-HPLC retention times. Side-chain hydrophilicity/hydrophobicity coefficients were generated in the presence of these nearest-neighbour effects and, despite the relative difference in such coefficients generated from peptides substituted with L- or D-amino acids, the relative difference in hydrophilicity/hydrophobicity between different amino acids in the L- or D-series is maintained. Overall, our results demonstrate that such nearest-neighbour effects can clearly restrict conformational space of an amino acid side-chain in a polypeptide chain.

Determination of intrinsic hydrophilicity/hydrophobicity of amino acid side chains in peptides in the absence of nearest-neighbor or conformational effects

Understanding the hydrophilicity/hydrophobicity of amino acid side chains in peptides/proteins is one the most important aspects of biology. Though many hydrophilicity/hydrophobicity scales have been generated, an "intrinsic" scale has yet to be achieved. "Intrinsic" implies the maximum possible hydrophilicity/hydrophobicity of side chains in the absence of nearest-neighbor or conformational effects that would decrease the full expression of the side-chain hydrophilicity/hydrophobicity when the side chain is in a polypeptide chain. Such a scale is the fundamental starting point for determining the parameters that affect side-chain hydrophobicity and for quantifying such effects in peptides and proteins. A 10-residue peptide sequence, Ac-X-G-A-K-G-A-G-V-G-L-amide, was designed to enable the determination of the intrinsic values, where position X was substituted by all 20 naturally occurring amino acids and norvaline, norleucine, and ornithine. The coefficients were determined by reversed-phase high-performance liquid chromatography using six different mobile phase conditions involving different pH values (2, 5, and 7), ion-pairing reagents, and the presence and absence of different salts. The results show that the intrinsic hydrophilicity/hydrophobicity of amino acid side chains in peptides (proteins) is independent of pH, buffer conditions, or whether C(8) or C(18) reversed-phase columns were used for 17 side chains (Gly, Ala, Cys, Pro, Val, nVal, Leu, nLeu, Ile, Met, Tyr, Phe, Trp, Ser, Thr, Asn, and Gln) and dependent on pH and buffer conditions, including the type of salt or ion-pairing reagent for potentially charged side chains (Orn, Lys, His, Arg, Asp, and Glu).

Mapping disturbances in a mangrove forest using multi-date landsat TM imagery

To evaluate the accounts of local fishermen, Landsat TM images (1986, 1993, 1999) were examined to assess potential losses in the mangrove forests of the Teacapán-Agua Brava lagoon system, Mexico. A binary change mask derived from image differencing of a band 4/3 ratio was employed to calculate any changes within this forested wetland. The results indicate that by 1986 approximately 18% (or 86 km2) of the mangrove area under study was either dead or in poor condition. The majority of this damage had occurred in the eastern section of the Agua Brava basin, which coincides, with the reports of the elderly fishermen. Examination of aerial photographs from 1970 revealed no adverse impacts in this area and would suggest, as postulated by the fishermen and other scientists, that modifications in environmental conditions following the opening of a canal, Cuautlá canal, in 1972 may have initiated the large-scale mortality. Although these areas of impact are still developing, the results from the satellite data indicate that the majority of the more recent changes are occurring elsewhere in the system. Obvious in the 1999 satellite data, but not so in the 1993, are large areas of mangrove degradation in the northern section of the Teacapán region. In the Agua Brava basin, the more recent transformations are appearing on the western side of the basin. Since long-term records of environmental conditions are absent, it is difficult to determine why these latest changes are occurring or even if the earlier losses were the result of the canal. Potential agents of change that have recently been observed include a hurricane, a second canal, and the uncontrolled expansion of the Cuautlá canal since 1994.