Abstract

DX600, a small peptide with 26 residues, is a potent, highly selective inhibitor of angiotensin converting enzyme 2 (ACE2). A range of NMR methods including TOCSY and ROESY yield an assignment of its proton spectrum in water and constraints on its conformation. Constrained molecular dynamics simulations of solvated DX600 show that the peptide's most abundant conformer adopts a predominantly random coil conformation. Constrained by the disulfide bond, its backbone defines an overhand knot with frayed ends.

1. Introduction

Angiotensin converting enzyme (ACE), a dipeptidase, is a key component of the renin-angiotensin system that regulates blood pressure [1, 2]. It catalyzes the conversion of inactive angiotensin I to angiotensin II, a potent vasoconstrictor, and the hydrolysis of bradykinin, a vasodilator. Consequently the treatment of hypertension makes extensive use of potent, selective inhibitors of ACE such as captopril and trandolapril [3]. Given the medical importance of the inhibition of ACE, the drug-design community has produced several X-ray structures of ACE-inhibitor complexes [4] and NMR solution-phase structures of the free inhibitors [5, 6].

ACE2 is a closely related enzyme which exhibits considerable sequence and structural homology with ACE [7, 8]. Both are peptidases with zinc at the active center. However, there are differences. Whereas ACE is widely distributed in the body, ACE2 is found principally in the heart, kidney, and testis. It cleaves the terminal leucine residue from the decapeptide angiotensin I and does not play a major role in the control of blood pressure. Knockout studies with mice indicate that ACE2 affects cardiac function [9].

To date, no inhibitors of ACE2 have entered clinical practice [3]. However, inhibitors of ACE2 have been developed as tools for exploring and modulating its biological function [10, 11]. Huang et al. screened peptide libraries displayed on phage and identified six highly potent inhibitors of ACE2 which do not inhibit ACE [12]. DX600, the most potent in the set with (the reciprocal of the equilibrium constant for the binding of the inhibitor, e.g., DX600, to the enzyme) equal to 2.8 nM, has the sequence Gly1-Asp2-Tyr3-Ser4-His5-Cys6-Ser7-Pro8-Leu9-Arg10-Tyr11-Tyr12-Pro13-Trp14-Trp15-Lys16-Cys17-Thr18-Tyr19-Pro20-Asp21-Pro22-Glu23-Gly24-Gly25-Gly26. It is acetylated at the amino terminus and amidated at the carboxyl terminus and a disulfide bond links the two cysteine residues, Cys6 and Cys17. DX600’s properties as an inhibitor mark it as a candidate for an NMR study and its unusual composition makes it particularly attractive. The small peptide contains one disulfide bond, six aromatic residues (Tyr3, Tyr11, Tyr12, Trp14, Trp15, Tyr19), four prolines, and six residues with charged side chains (Asp2, His5, Arg10, Lys16, Asp21, Glu23). This is the first study of the NMR spectroscopy and molecular modeling of DX600. No X-ray structure has been published.

The 1D proton spectra of DX600 in D2O and H2O (Figure 1) show that it is present as several conformers with one conformer predominating. We have employed a broad range of 1D and 2D methods [13] and a combination of proton and carbon-13 spectra to assign the proton spectrum of the most abundant conformer of DX600 in water and to generate NMR constraints on its conformation. An analysis of the NMR results and constrained molecular dynamics (MD) simulations of solvated DX600 show the random-coil character of its most abundant conformer.

2. Materials and Methods

2.1. Sample Preparation

DX600 was purchased from Bachem Americas, Inc. and used without further purification. The supplier confirmed the sample’s purity as >96% via HPLC. 0.65 and 1.3 mM solutions were prepared by dissolving samples of the peptide in D2O and 90% H2O/10% D2O (volume %). The 99.99 atom-% D2O used to prepare the solutions was obtained from Aldrich. The pH of the resulting solutions at 25°C was 3.6.

2.2. NMR Experiment

All NMR measurements were made on a Bruker DPX400 spectrometer equipped with an Avance console, a variable-temperature, proton/carbon, inverse detection probe, and hardware for pulse shaping and gradients. At its magnet field strength of 9.4 Tesla, the proton and carbon-13 Larmor frequencies are 400.13 and 100.62 MHz, respectively. Data acquisition and processing were performed with Version 3.0 of XWINNMR.

The water peak was used as an internal standard for the proton chemical shifts. The position of the water peak with respect to TSP was determined in a separate experiment with an aqueous solution of TSP. This approach of external referencing was validated at the end of the study. TSP was added to the sample and the chemical shifts measured directly with respect to TSP equaled the values obtained by the external referencing. Following the IUPAC recommendation, carbon-13 referencing was based on the carbon-13 and proton values and the position of the TSP singlet in the proton spectrum [14].

Solvent suppression was applied in the acquisition of all proton spectra. In the case of the 100% D2O solutions, an inversion-recovery technique was effective in suppressing the HDO signal without saturating nearby peaks. The preparation period consists of a composite 180° pulse, that is, 90°y-180°x-90°y, followed by a delay, and finally a detection 90° pulse. The delay is set at the null of the HDO peak following the inverting 180° pulse. The delay, 2.0 s at 25°C, is long compared with of the DX600 so its magnetization fully recovers when HDO is at its null. Because of radiation damping, this approach does not work for the 90% H2O solutions. In this case, a standard Watergate pulse sequence is used as the preparation pulse [15].

Standard methods were used for the acquisition of the 1D spectra and the 2D COSY, TOCSY, 2D J resolved, HMQC, and HMBC spectra [15]. Temperature coefficients of chemical shifts were obtained from TOCSY spectra measured at 25°C and 45°C and 1D proton measurements at 2, 25, and 40°C. The d2 delay for the HMQC and HMBC methods was 3.57 ms, corresponding to of 140 Hz. The d6 HMBC delay was 60 ms, corresponding to of 8.3 Hz. The acquisition and processing parameters are tabulated in Table 1. The pulse sequence of Hwang and Shaka was used to acquire the ROESY spectra [16]. The ROESY mixing times were 250 and 500 ms. In order to increase the signal-to-noise ratio of weak nOe’s involving well-resolved signals, we used the NOE difference method of Gordon and Wüthrich to perform a series of 1D experiments [17]. The standard method was modified to include Watergate solvent suppression. The experiment consists of a selective 180° pulse (a Gaussian shaped pulse of length 38 ms), a delay of 250 ms for the buildup of transfer of magnetization, and a Watergate detection pulse. The experiment is repeated with the application of the 180° pulse set well off resonance, for example, −8 ppm. The nOe’s appear in the difference of the two spectra. The method of Stejskal and Tanner modified by Wu et al. involving bipolar gradient pulses and the LED pulse sequence (BP-LDE) was used to determine the translational diffusion constant [18]. Longsworth accurately determined of HDO in D2O [19]. The gradient coil was calibrated by running the method and observing the HDO signal in 99.99 atom% D2O. The following parameters were used in the diffusion measurements: TD, 8 k; d1, 2 s; , 0.4 s; (gradient period), 2 ms; , 5 ms; , 0.1 ms; gradient current, between 1% and 99%.

Upper bounds on distances between protons were extracted from the 500 ms ROESY spectrum using the methods used in our work on erythromycin [20, 21]. Values of  Hz were converted into lower and upper bounds on , the distance between the amide and alpha protons [21, 22]. The difference between the lower and upper bounds, 2.951 Å and 3.002 Å, is proportionate to the uncertainty in the coupling constants and corresponds to an uncertainty of 30° in the torsional angle.

2.3. Molecular Modeling

All modeling calculations were run on a Dell Precision 370 PC running under Red Hat Enterprise Linux Version 5 and using Versions 7.3 and 8.1 of SYBYL (Tripos, Inc., 1669 S. Hanley Road, St. Louis, MO). Molecular surfaces were calculated using MOLCAD. Random search parameters were set to the default values with the following exceptions: Bump Factors, 0.02; Ring Bond Closure, 10 Å; energy cutoff, 10 kcal/mole. The force field for all calculations was MMFF94s, Halgren’s modification of the Merck Molecular Force Field for peptides [23]. The integration time for all MD simulations was 1 fs. With the exception of the temperature and the length of the simulation, SYBYL default parameters were used for all MD simulations and energy minimizations. SYBYL provides a range of algorithms for energy minimization. Our energy minimization calculations employed the Powell algorithm with a Simplex initial optimization. Other available algorithms such as conjugate gradient and steepest descent were examined, but they were not used in the final calculations as they were less efficient. The MD simulations employed a canonical ensemble, for example, at constant volume.

The structure of DX600 was drawn with the Biopolymers module of SYBYL, thus guaranteeing the correct atom types and charges for the side chains: negative for Asp and Glu, positive for His, Lys, and Arg. These charge assignments did not change during the MD simulations. After formation of the disulfide bond, the torsional and distance constraints were added as quadratic penalty functions. The constraints were applied in all the simulations. In the MD simulations, a molecule of DX600 was embedded in a rectangular box containing 1651 water molecules and the energy of the system of water and solute molecules was minimized. In the simulations the interactions between the solute and the water molecules were defined by the MMFF94s force field. Periodic boundary conditions were employed in all MD and energy minimization calculations of hydrated DX600. The global-minimum structure generated by the in vacuo conformational search was a seed in 30 cycles consisting of the following steps: heating at 800 K for 10000 fs, cooling at 298 K for 1500 fs, and 200 cycles of energy minimization. Energy minimization of the structures generated by the final 11 cycles was then performed to convergence, that is, until the average gradient reached 0.05 kcal/mol-Å. This final stage of optimization involved 6000 to 24000 additional cycles. The structures examined in Section 3.3 are based on the results of the final 11 cycles. Several pieces of evidence indicate equilibration of the structure. SYBYL’s energy minimization module colors the atoms to indicate where forces act to change the structure. During the final stages of energy minimization, the color of the water molecules indicated the absence of gradients on the solvent molecules. Two distances, the separation of the alpha carbons on Ser4 and Glu23 and the distance between the center of the aromatic ring on Tyr3 and the amide nitrogen on Asp2, did not vary by more than 0.1 Å during the final energy minimization. SYBYL’s Monitor Distance module ruled out bad contacts between the peptide and the water. It displays violations where the distance between a pair of nonbonded atoms is less than the sum of the van der Waals radii; no violations were found in the final structures.

The in vacuo conformational search involved thousands of cycles of the SYBYL Random Search algorithm. In each cycle, randomly chosen torsional angles are changed and the resulting structure was minimized. This strategy which generated the global minimum and a small family of structures close in energy to the global minimum was very successful in conformational searches of derivatives of erythromycin. A detailed discussion of the methodology can be found in our papers on the NMR structures of these drugs [21, 22].

3. Results and Discussion

3.1. Analysis of the Survey Spectra

The proton survey spectra of DX600 dissolved in D2O and H2O (Figure 1) indicate a predominantly random-coil conformation for the most abundant conformer. The alpha protons of most residues are found within 0.3 ppm of their random coil value. The small peptide has seven aromatic residues and ring current effects are expected to make extensive contributions to the chemical shifts. For example, the chemical shifts of three of the four downfield shifted alpha protons and the upfield shifted methyl protons of Leu9 are affected by nearby aromatic residues. Exchange experiments show the absence of extensive secondary structure. Serial measurements of the proton spectrum made immediately after dissolving the compound in D2O yielded no amide protons with one sample. Measurements with a second sample showed peaks for four amide protons that disappeared with a half life of ca. 15 minutes. They were later assigned to Trp14, Lys16, Tyr19, and Asp21. The temperature coefficients for the chemical shifts of these four amide protons, , are −3.0, −8.2, −0.6, and −5.5 ppb/°C, respectively. The combination of measurable exchange kinetics and a small , which is diagnostic for hydrogen bonding [24], indicates that the amide protons of Trp14 and Tyr19 form weak hydrogen bonds. However, one can conclude from the rapid exchange with D2O that most amide protons are exposed to the solvent and/or unprotected by strong hydrogen bonds. DX600 differs from mast cell degranulating peptide (MCD), another small peptide. Eight amide protons in MCD, a peptide with 22 residues including one proline and two disulfide bridges, exchange slowly enough that they can be observed 30 min. after dissolving the peptide [25]. In this case, the amide protons are stabilized by incorporation in an alpha-helix. The rapid exchange in DX600 indicates the absence of elements of secondary structure such as alpha-helix and beta-sheet.

The survey spectra clearly show the presence of conformers of the most abundant species. Pairs of peaks in the ratio of 2–3 : 1 are seen at 9.9, 8.6, and 1.2 ppm. These are assigned, respectively, to the N1H proton of Trp15, the 2H proton of His5, and the methyl protons of Leu9. The TOCSY spectrum also provides additional evidence for conformers. Furthermore minor cross peaks assigned to Glu23/Gly24 connectivities appear in the ROESY spectrum. The 2D cross peaks assigned to conformers other than the most abundant species correlate positions in the 1D spectrum with peaks integrating to a mole fraction less than 0.3. 1D spectra at temperatures up to 60°C and the ROESY spectrum at 45°C show no evidence of exchange broadening and cross peaks due to conformational exchange. The presence of conformers and their slow exchange can be attributed to the prolines and the disulfide bridge in DX600 and very likely involves cis- and trans-proline peptide bonds. The presence of cis- and trans-proline is a prominent feature of the NMR spectra of MCD [25].

The peaks in the proton spectrum of DX600 are much broader than one would expect for a peptide of its size. Qualitatively the spectrum of MCD is much better resolved. We attribute the result to the random coil conformation. DX600’s overall shape, which is extended and asymmetric rather than spherical, leads to a decrease in the rotational diffusion constant and a corresponding decrease in . Application of the BP-LDE method at 25°C yielded a translational diffusion constant of  m2/s in D2O. A value for the spurious near spherical structure generated by the initial in vacuo modeling can be estimated using the well known relation from hydrodynamic theory between and the molecular radius : where is the viscosity of solvent [26]. Using the molecular volume of the spherical form of DX600, one calculates  m2/s, a value greater than the observed value by a factor of roughly 2.

The nonspherical, random coil conformation of DX600 and the resulting decreased values of and have several consequences for the measurement of NMR spectra. The ROESY spectrum is sparse and the cross-peaks are weak (Figure 2). The ROESY spectrum yields no long-range constraints and only four medium range constraints. Furthermore, the HMBC spectrum, in comparison with the HMQC spectrum, is sparse. A crucial component of the HMBC method is a filter tuned to geminal and vicinal proton-carbon coupling. Small coupling constants require a long filter time, for example, 60 ms, during which the magnetization decays to nearly zero. As a result the HMBC spectrum is dominated by spins attached to relatively mobile residues or side chains; these spins have a longer effective and their magnetization survives the filter. A short is a problem as it frustrates the use of HMBC in the assignment of the proton spectrum. However, the HMBC spectrum identifies the mobile regions of the molecule.

3.2. Assignment of the Proton Spectrum of the Most Abundant Conformer

The sequence-specific assignment of the aqueous solution of DX600 rests on a collection of techniques conducted at 25°C. The assignment was validated by additional measurements at 45°C. Anchor points for the assignment are the readily assigned spin systems: the acetyl group at the amino terminus and the single copies of threonine, leucine, arginine, glutamic acid, and lysine. These amino acids have unique spin systems that can be completely sorted out in the TOCSY spectrum. An additional tag for the lysine spin system is the positively charged terminal ammonium group which yields a prominent, positive cross peak in the ROESY spectrum due to chemical exchange. The spin systems of the 26 residues, identified by the TOCSY spectrum in H2O (Figure 3) and the COSY spectrum in D2O, are linked together by the ROESY spectrum (Figure 2). Weak nOe’s were obtained by the NOE difference experiment. The assignment was further supported by rOe’s between the beta and ring protons of the aromatic residues and HMBC cross-peaks for the mobile side chains of Leu9, Thr18, Glu23, and the aromatic residues. The evidence supporting the assignment is summarized in Table 2. Redundancy of the evidence is important given the sparseness of the ROESY spectrum, the weakness of many cross-peaks, the absence of several cross-peaks, and the existence of conformers. The proton assignments for the most abundant conformer are given in Table 3. The assignment of the residues Tyr19 through Glu23 required special care because of the absence of ROESY cross-peaks for and connectivities. The key to their assignment, summarized in Table 4, is a cluster of nOe’s in which the protons of Asp21 play a central role. Satisfying the data in Table 4 leads to cis-peptide bonds for Pro20, Pro22, and Glu23. The assignment of the alpha protons of Asp21 and Pro22 was validated by the HMQC spectrum.

3.3. The Conformation of DX60

Forty-six distance constraints extracted from the ROESY spectrum and the NOE difference experiments were used in the modeling calculations (Table 5). This set consists of 12 intraresidue constraints, 30 sequential constraints, but only 4 medium-range constraints and no long-range constraints. Half of these constraints involve the more flexible aromatic rings. The few violations of the constraints in the modeling calculations did not exceed 0.05 Å. Constraints on were obtained in the three cases—Tyr11, Cys17, and Thr18—where the value of was large enough (≥9 Hz) to enable an unambiguous interpretation of a nearly s-trans-dihedral angle.

We began our modeling of DX600’s conformation with an in vacuo approach that has worked well for drugs such as erythromycin: an exhaustive search of conformational space and a constrained minimization [2022]. However, this method yielded a result that disagreed with the NMR evidence. The set of low energy conformers generated by in vacuo modeling yielded a well-defined, nearly spherical structure. The magnitude of the monopole-monopole interactions between the six charged side chains argues that the role of the solvent and its large dielectric constant must be explicitly included in the modeling calculations. Therefore, in a second series of calculations, one molecule of DX600 was imbedded in a rectangular box of water molecules. A series of MD simulations for aquated DX600 yielded a predominantly random coil conformation for DX600 in water. Figure 4 shows the overlay of the 11 energy-minimized structures obtained in the final MD cycles. The average RMSD using all backbone atoms is 3.5 Å and using the backbone atoms on residues Ser4 through Glu23 is 2.6 Å. Overall, the backbone of DX600 forms an overhand knot with frayed ends. The poor overlay of the structures at the ends shows that the peptide does not have a well-defined conformation at the amino and carboxy termini. Except for the turn required by the disulfide bond between the two cysteines, the peptide adopts an extended, nonspherical conformation which explains the small translational diffusion constant and the short value of . Given the size of the peptide and its extended conformation, most residues are exposed to the solvent. This exposure explains the relative flexibility of the leucine, threonine, glutamate, and aromatic side chains. The six charged side chains extend into the solvent so that peptide-water Coulombic interactions make a significant contribution to DX600’s conformation.

4. Conclusions

Although DX600 has a predominantly random-coil structure, the constraint of the disulfide bond yields the formation of a loop. Several lessons were learned from this study. A complete assignment of the proton spectrum of the most abundant conformer was only achieved by combining many sources of information. The modeling of a small peptide with a large fraction of charged residues must be handled with care. The development of a structural model from the NMR data required explicit solvation. In vacuo modeling calculations led to an incorrect result.

Acknowledgment

Financial support for this project was provided by the Frederic J. Robbins Fund for Chemistry at Pomona College.