1. Sample Preparation Apo-sample: Dissolve approximately 1-2 mg of the peptide in 450-500 µl of deuterated NMR grade solvent (the preferable solutions for biological samples are 10% D2O in water or d6-DMSO11-12 if the sample does not dissolve in water or reacts with it). Copper-reacted sample: Dissolve approximately 1-2 mg of the peptide with an equimolar amount of metal salt in 450-500 µl of the NMR solvent. Filter the solution using a Sinter glass, filtration paper or any other technique that suits the compounds under investigation and does not adsorb them. This is essential to remove any metallic particles, which will affect the homogeneity. Transfer the solution to an NMR tube and close it. Make sure that the tube quality fits the magnet strength used (for more details see equipment table). If the sample is sensitive to oxygen, these steps must be done under an inert environment and the tube sealed to prevent oxidation. 2. NMR Data Collection and Processing13 Record 1D 1H spectra of the apo- and copper-reacted samples and compare. The copper-containing peptide spectrum should show a significant change of chemical shift in the amide region and the peaks became resolved since the apo-peptide is flexible and shows an average of conformations, but upon reacting with copper, the bound-peptide amides have a single structure (Figure 2). If the spectrum is unchanged, the reaction is assumed to have failed. If the sample turns green, the copper will have oxidized in the atmosphere, which may change its binding properties. In both these cases, the sample will need to be remade. Find optimal temperature conditions for the NMR measurements that give minimal overlap of the amide residues with narrow line-widths. Set up COSY14, TOCSY15 and NOESY16, or ROESY17 NMR experiments under identical conditions (temperature, pH, etc.) and run in sequence. COSY (Figure 3) and TOCSY (Figure 4) experiments will detect through-bond interactions of adjacent and longer range neighboring proton-proton interactions, respectively. COSY spectrum gives information on HN-Hα correlations while TOCSY spectrum gives, in addition information on HN and other H aliphatic protons of correlated residues. NOESY and ROESY (Figure 5) experiments detect through-space proximity with no dependence on bonds to give structural information. The choice between the two is according to the effective size of the compound and the field strength. Some such combinations give theoretical zero signals (Figure 6) and a ROESY experiment will have to be run to detect NOE interactions10,18. If the sample is in water, the experiments need to have a component of water suppression. This is most efficiently accomplished using gradients19. In this case, it is advantageous to rerun the sample, once assigned, in pure D2O to obtain interactions to and among the Hα hydrogens. Optimize the duration of the TOCSY and NOESY or ROESY mixing times by running a number of short experiments to find maximal signal build-up with minimal loss to spin diffusion. Make sure that this is the linear region of the buildup curve. Values common for peptides at 600 MHz fields include mixing times of around 100 msec for TOCSY and 150 msec for NOESY or ROESY experiments. Run 1D experiments in between each experiment to make sure that the sample composition remains constant throughout data acquisition (Figure 7). Process the spectra using appropriate apodization functions to obtain maximal resolution with minimal loss of signal strength, adding zero filling in the t1 dimension. Further correct the baseline in the F2 then F1 dimension with a quadratic polynomial function. Calibrate the chemical shift of the spectra carefully according to the residual water or DMSO peaks in the respective solvents according to their shift from known standards (TMSP or TMS, respectively). 3. Peak Assignment and Integration Using SPARKY20 Prepare a set of COSY and TOCSY spectra overlaid on the NOESY or ROESY spectrum (Figure 8). Assign all NOE peaks in the spectrum according to the sequential assignment methodology developed by Wüthrich21. Start with assigning peaks that overlap TOCSY signals in the fingerprint region, as this will facilitate subsequent peak assignment entries with the SPARKY program (Figure 9). Register the assignment chemical shift (Figure 10), and 3JHNHα values (Figure 11). Translate peaks into distance constraints and 3JHNHα values into dihedral angles. This can be done by integrating the peaks from within the program and translating them into distance constraints using an interaction of known distance (such as the distance between two protons of methylene group or aromatic hydrogen atoms of an aromatic amino acid). If the peaks overlap, and integration methods cannot be used, they can be labeled as strong-medium-weak-very weak by visual estimation, and these designations can be translated into distances up to 2.5, 3.5, 4.5, and 5.5 Å, respectively. This to be most effective in work with peptides. 3JHNHα values are indicative of dihedral angles according to the Karplus equation where most given coupling values can result from more than one dihedral angle. These can be indicative of secondary structure (helices or sheets), random coil or conformational averages thereof. Therefore, it is important to use the constraints carefully, or to use them as confirmation of conformational results. 4. Molecular Mechanics Calculations to Generate Structure Ensemble using XPLOR22 Import the distance constraints and dihedral angles with the correct format for XPLOR. (Figure 12, inter-residual constraints only). XPLOR will search conformational space to find structures that adhere to canonical geometric chemical constraints, such as bond lengths, angles, chirality, etc. in addition to the experimentally found distance constraints, to generate an ensemble in which none of the above parameters are violated. This will constitute the starting ensemble. The first structure determination run will be without using any constraints to the metal. This will show which residues participate in metal binding without any bias. Skip to step 4.4. In addition to the NMR-derived distance constraints, add metal binding constraints to the residues determined. Add appropriate parameters to describe the metal ion and its topology. The appropriate physical information (mass, bond lengths with other atoms, angles, and nonbonding repulsion parameters) must be put into the parameter file. Add the binding information to the topology file: Which bonds are formed and broken and which formal charges are changed as a result of binding. Create a small ensemble of usually 10 members. Introduce the constraints gradually, starting with backbone-to-backbone interactions, then backbone-to-sidechain, then sidechain-to-sidechain interactions, and going from intra-residual, sequential to increasingly longer range interactions. This facilitates identifying mistakes in assignment. Introduce the NOE energy as a square-well potential with a constant force constant of 50 kcal/mol•Å2. Simulated annealing should be done with 1,500 3-fsec steps at 1,500 K and 3,000 1-fsec steps during cooling to 500 K. Finally, minimize the structures using conjugate gradient energy minimization for 4,000 iterations; these are variables that can be changed within XPLOR. Create a final ensemble of usually 50 members. Perform steps from step 4.4 on the entire ensemble. Report the number of each type of NOE interaction that were found, as shown in Figure 13. Create an ensemble of structures that adhere to canonical chemical geometry and the empiric NMR-derived constraints. Report the total number of conformations, the number of these that have violations of the NMR-derived constraints, and the RMSD of the entire ensemble, both backbone, and all heavy atom RMSD values. 5. Structure Analysis The solved ensemble represents the conformational space adopted by the peptide during the NMR measurement. The local flexibility may change in different regions of the molecule and this can be ascribed to structural or functional reasons (the purpose of structural analysis is to determine the structural basis for the stability and mechanism of action). Import all conformations of the structure to the MolMol program23 to create a starting ensemble (Figure 14). Examine the ensemble to determine the local stability of the molecule. Determine the backbone and sidechain RMSD values by selecting subsequent four-residue regions along the sequence and having the program calculate the RMSD to the lowest energy structure or the mean. Determine which regions of the molecule show local stability (Figure 15) by plotting the local RMSD as a function of sequence. Overlay the ensemble along this region of the molecule and use this ensemble for further analysis (Figure 16). Choose low-energy conformations that adhere to the NMR-derived constraints (nonviolated). These will form the "low energy ensemble". Report the number of conformations in the ensemble, the criteria for choosing them and the RMSD values as defined in step 4.6.1. If the metal binding mode was already determined, continue to the next step. Otherwise: Analyze the low energy ensemble and determine which residue sidechains are in correct proximity of each other to be able to bind the metal ion. Once these have been determined, go back to step 3.3 and redo the analysis including the copper-binding data (Figure 17 for the ensemble and Figure 18 for the low energy conformer). Examine the ensemble and determine local secondary structure within the molecule using the default search parameters of the MolMol program. Secondary structures are held by hydrogen-bonding and indicate stable regions of the molecule. Determine the secondary structure of each conformer (Figure 19). Create a table of the percentage of conformers of the ensemble that show each secondary structure. Determine whether the regions of secondary structure overlap the regions that were determined to be the stable regions by local-RMSD. This is often the case. Determine distances and hydrogen bonding within the molecule. Import the ensemble into Chimera24. Use Chimera to determine intramolecular distances between atoms suspected for metal binding. Calculate the average distances in the ensemble. Use Chimera to determine intramolecular hydrogen bonds in each conformer of the ensemble (Figure 20), using the default relaxed hydrogen-bond parameters (by 20%) of the program. Create a table of the percentage of conformers of the ensemble that show each hydrogen-bond. Hydrogen bonds that recur in a majority of the conformations indicate a stable bond that adds to the stability of the structure. Determine electrostatic interactions within the molecule. Identify interactions between charged sidechains. Create a table of the percentage of conformers of the ensemble that show each electrostatic interaction. Calculate the electrostatic potential distribution using the Amber25 force field within the Delphi26 program. Choose a single representative conformation from the ensemble on which to perform the electrostatic potential calculations. Derive electrostatic potential using the Poisson-Boltzman equation and a full Coulombic calculation. Required parameters include the ionic strength of the solution, the dielectric values of the solution (80 for water; 40 for DMSO); the internal peptide (usually around 5.0); the size of the grid (usually 65´65´65 points); and the minimal distance between the peptide and the edge of the grid (10 Å). The results are usually quite robust with respect to these values as the main result is governed by the peptide itself (Figure 21). Examine the electrostatic potential mapped onto the Van der Waals surface of the peptide (Figure 22). Examine the electrostatic potential iso-surfaces, which describe positions having the same potential. Find iso-surfaces that suggest biological relevance, such as extended positive or negative surfaces that might attract or repel other proteins or ligands, or distributions that indicate significance, such as amphipathic distribution or patches of differing properties (Figure 23). Sum up all the structural findings to see how they reinforce each other. Repeat steps 4.3-5.7.4 where relevant.