Asp Side Chain Research Articles

A Threonine Stabilizes the NiC and NiR Catalytic Intermediates of [NiFe]-hydrogenase

The heterodimeric [NiFe] hydrogenase from Desulfovibrio fructosovorans catalyzes the reversible oxidation of H2 into protons and electrons. The catalytic intermediates have been attributed to forms of the active site (NiSI, NiR, and NiC) detected using spectroscopic methods under potentiometric but non-catalytic conditions. Here, we produced variants by replacing the conserved Thr-18 residue in the small subunit with Ser, Val, Gln, Gly, or Asp, and we analyzed the effects of these mutations on the kinetic (H2 oxidation, H2 production, and H/D exchange), spectroscopic (IR, EPR), and structural properties of the enzyme. The mutations disrupt the H-bond network in the crystals and have a strong effect on H2 oxidation and H2 production turnover rates. However, the absence of correlation between activity and rate of H/D exchange in the series of variants suggests that the alcoholic group of Thr-18 is not necessarily a proton relay. Instead, the correlation between H2 oxidation and production activity and the detection of the NiC species in reduced samples confirms that NiC is a catalytic intermediate and suggests that Thr-18 is important to stabilize the local protein structure of the active site ensuring fast NiSI-NiC-NiR interconversions during H2 oxidation/production.

Roles of Intramolecular and Intermolecular Hydrogen Bonding in a Three-Water-Assisted Mechanism of Succinimide Formation from Aspartic Acid Residues

Aspartic acid (Asp) residues in peptides and proteins are prone to isomerization to the β-form and racemization via a five-membered succinimide intermediate. These nonenzymatic reactions have relevance to aging and age-related diseases. In this paper, we report a three water molecule-assisted, six-step mechanism for the formation of succinimide from Asp residues found by density functional theory calculations. The first two steps constitute a stepwise iminolization of the C-terminal amide group. This iminolization involves a quintuple proton transfer along intramolecular and intermolecular hydrogen bonds formed by the C-terminal amide group, the side-chain carboxyl group, and the three water molecules. After a conformational change (which breaks the intramolecular hydrogen bond involving the iminol nitrogen) and a reorganization of water molecules, the iminol nitrogen nucleophilically attacks the carboxyl carbon of the Asp side chain to form a five-membered ring. This cyclization is accompanied by a triple proton transfer involving two water molecules, so that a gem-diol tetrahedral intermediate is formed. The last step is dehydration of the gem-diol group catalyzed by one water molecule, and this is the rate-determining step. The calculated overall activation barrier (26.7 kcal mol−1) agrees well with an experimental activation energy.

Structure of a sedoheptulose 7-phosphate cyclase: ValA from Streptomyces hygroscopicus.

Sedoheptulose 7-phosphate cyclases (SH7PCs) encompass three enzymes involved in producing the core cyclitol structures of pseudoglycosides and similar bioactive natural products. One such enzyme is ValA from Streptomyces hygroscopicus subsp. jinggangensis 5008, which makes 2-epi-5-epi-valiolone as part of the biosynthesis of the agricultural antifungal agent validamycin A. We present, as the first SH7PC structure, the 2.1 Å resolution crystal structure of ValA in complex with NAD+ and Zn2+ cofactors. ValA has a fold and active site organization resembling those of the sugar phosphate cyclase dehydroquinate synthase (DHQS) and contains two notable, previously unrecognized interactions between NAD+ and Asp side chains conserved in all sugar phosphate cyclases that may influence catalysis. Because the domains of ValA adopt a nearly closed conformation even though no sugar substrate is present, comparisons with a ligand-bound DHQS provide a model for aspects of substrate binding. One striking active site difference is a loop that adopts a distinct conformation as a result of an Asp → Asn change with respect to DHQS and alters the identity and orientation of a key Arg residue. This and other active site differences in ValA are mostly localized to areas where the ValA substrate differs from that of DHQS. Sequence comparisons with a second SH7PC making a product with distinct stereochemistry lead us to postulate that the product stereochemistry of a given SH7PC is not the result of events taking place during catalysis but is accomplished by selective binding of either the α or β pyranose anomer of the substrate.

Intramolecular cyclization of aspartic acid residues assisted by three water molecules: a density functional theory study

Aspartic acid (Asp) residues in peptides and proteins (l-Asp) are known to undergo spontaneous nonenzymatic reactions to form l-β-Asp, d-Asp, and d-β-Asp residues. The formation of these abnormal Asp residues in proteins may affect their three-dimensional structures and hence their properties and functions. Indeed, the reactions have been thought to contribute to aging and pathologies. Most of the above reactions of the l-Asp residues proceed via a cyclic succinimide intermediate. In this paper, a novel three-water-assisted mechanism is proposed for cyclization of an Asp residue (forming a gem-diol precursor of the succinimide) by the B3LYP/6-31 + G(d,p) density functional theory calculations carried out for an Asp-containing model compound (Ace−Asp−Nme, where Ace = acetyl and Nme = NHCH3). The three water molecules act as catalysts by mediating 'long-range' proton transfers. In the proposed mechanism, the amide group on the C-terminal side of the Asp residue is first converted to the tautomeric iminol form (iminolization). Then, reorientation of a water molecule and a conformational change occur successively, followed by the nucleophilic attack of the iminol nitrogen on the carboxyl carbon of the Asp side chain to form the gem-diol species. A satisfactory agreement was obtained between the calculated and experimental energetics.

Just three water molecules can trigger the undesired nonenzymatic reactions of aspartic acid residues: new insight from a quantum-chemical study

Aspartic acid (Asp) residues in peptides and proteins (L-Asp) can undergo spontaneous, nonenzymatic reactions under physiological conditions by which abnormal L-β-Asp, D-Asp, and/or D-β-Asp residues are formed. These altered Asp residues may affect the three-dimensional structures of the peptides and proteins and hence their properties and functions. In fact, the altered Asp residues are relevant to age-related diseases such as cataract and Alzheimer's disease. Most of the above reactions of the L-Asp residue proceed via a cyclic succinimide intermediate. In this paper, I propose a detailed mechanism of cyclization of an Asp residue (forming a precursor of the succinimide) by the B3LYP/6-31+G(d,p) density functional theory calculations carried out for a small Asp-containing model compound complexed with three water molecules which act as general acid-base catalysts in proton transfers. In the proposed mechanism, the amide group on the C-terminal side of the Asp residue is first converted to the tautomeric iminol form. Then, successive reorientation of a water molecule and conformational change occur followed by the nucleophilic attack of the iminol nitrogen atom on the carboxyl carbon atom of the Asp side chain to form a five-membered ring. A satisfactory agreement was obtained between the calculated and experimental energetics.

Complete integrin headpiece opening in eight steps

Carefully soaking crystals with Arg-Gly-Asp (RGD) peptides, we captured eight distinct RGD-bound conformations of the αIIbβ3 integrin headpiece. Starting from the closed βI domain conformation, we saw six intermediate βI conformations and finally the fully open βI with the hybrid domain swung out in the crystal lattice. The β1-α1 backbone that hydrogen bonds to the Asp side chain of RGD was the first element to move followed by adjacent to metal ion-dependent adhesion site Ca(2+), α1 helix, α1' helix, β6-α7 loop, α7 helix, and hybrid domain. We define in atomic detail how conformational change was transmitted over long distances in integrins, 40 Å from the ligand binding site to the opposite end of the βI domain and 80 Å to the far end of the hybrid domain. During these movements, RGD slid in its binding groove toward αIIb, and its Arg side chain became ordered. RGD concentration requirements in soaking suggested a >200-fold higher affinity after opening. The thermodynamic cycle shows how higher affinity pays the energetic cost of opening.

Adsorption of Arginine–Glycine–Aspartate Tripeptide onto Negatively Charged Rutile (110) Mediated by Cations: The Effect of Surface Hydroxylation

Classical molecular dynamics (MD) simulations were employed to investigate the adsorption behaviors of arginine-glycine-aspartate (RGD) tripeptide onto the negatively charged hydroxylated/nonhydroxylated rutile (110) surfaces, mediated by biologically important cations (Na+ or Ca2+). The simulation results indicate that the inherent nature of the cation determines its binding strength, thereby regulating the adsorption geometry of the peptide. The sparse hydroxyl groups on the nonhydroxylated rutile diminish the probability of H-bond formation between RGD and the surface, resulting in an early desorption of the peptide even with a mediating Na+ ion. In contrast, the negatively charged aspartate (Asp) side chain is bridged to the negatively charged hydroxylated rutile by an inner-sphere Na+ ion, in coordination with the Asp-rutile hydrogen bonds at the anchoring sites. The inner- and outer-sphere Ca2+ ions are demonstrated to be capable of "trapping" RGD on both hydroxylated and nonhydroxylated rutile, in the absence of hydrogen bonds with the surface. The strongly bound inner-sphere mediating Ca2+ ion exerts a "gluing" effect on the Asp side chain, producing a tightly packed RGD-rutile complex, whereas a less localized distribution density of the outer-sphere mediating Ca2+ ion results in a higher mobility of the Asp side chain. The intramolecular interaction is suggested to facilitate the structural stability of RGD adsorbed on the negative rutile in a "horseshoe" configuration.

Structural basis for recognition of 5'-phosphotyrosine adducts by Tdp2.

SummaryThe DNA repair enzyme TDP2 resolves 5′-phosphotyrosyl-DNA adducts, and is responsible for resistance to anti-cancer drugs that target covalent topoisomerase-DNA complexes. TDP2 also participates in key signaling pathways during development and tumorigenesis, and cleaves a protein-RNA linkage during picornavirus replication. The crystal structure of zebrafish TDP2 bound to DNA reveals a deep and narrow basic groove that selectively accommodates the 5′-end of single-stranded DNA in a stretched conformation. The crystal structure of the full-length C. elegans TDP2 shows that this groove can also accommodate an acidic peptide stretch in vitro, with Glu and Asp sidechains occupying the DNA backbone phosphate binding sites. This extensive molecular mimicry suggests a potential mechanism for auto-regulation and how TDP2 may interact with phosphorylated proteins in signaling. Our study provides a framework to interrogate functions of TDP2 and develop inhibitors for chemotherapeutic and antiviral applications.

Comparative analysis of photocaged RGDS peptides for cell patterning

Photocaged RGDS is a cell nonadhesive tetrapeptide that can be activated with light to become cell-adhesive. Such molecules can find useful applications in controlling cell adhesion for biological study, drug development, and in forming dynamic, adhesion-controlled biomaterials. Herein, we prepared RGDS peptide photocaged either on the Arg-Gly backbone amide nitrogen atom (R[-]GDS) or Asp side chain carboxyl (RG[D]S). A critical comparison of the peptides' chemical and physiological properties relevant for biological applications was carried out. It was observed that RG[D]S was synthesized more readily via automated solid-phase synthesis, underwent uncaging with a rate constant 3-fold higher than R[-]GDS, and was more stable in aqueous solution. Automated docking studies were performed to examine the interactions of various caged RGDS peptides with cell surface integrin receptor to identify suitable locations for the photosensitive 2-nitrobenzyl (NB) group for biological applications. A competitive binding ELISA method compared the ability of various peptides to bind to α(V)β(3) cell integrin receptors and the data were found to be consistent with the modeling predictions. Finally, the application of our caged RGDS peptides in controlling cell adhesion to form cell patterns on a hydrogel material was presented.

Effect of Glutamate Side Chain Length on Intrahelical Glutamate–Lysine Ion Pairing Interactions

Ion pairing interactions between oppositely charged amino acids are important for protein structure stability. Despite the apparent electrostatic nature of these interactions, the charged amino acids Lys, Arg, Glu, and Asp have a different number of hydrophobic methylenes linking the charged functionality to the backbone. To investigate the effect of Glu (and Asp) side chain length on ion pairing interactions, a series of 36 monomeric α-helical peptides containing Zbb-Xaa (i, i+3), (i, i+4), and (i, i+5) (Zbb = Aad, Glu, Asp; Xaa = Lys, Orn, Dab, Dap) sequence patterns were studied by circular dichroism (CD) spectroscopy at pH 7 and 2. Peptides with Glu and Aad exhibited similar helicity and pH dependence, whereas peptides with Asp behaved distinctly different. The side chain interaction energetics were derived from the CD data using the nesting block method coupled with modified Lifson-Roig theory. At pH 7, no Zbb-Xaa (i, i+5) interaction was observed, regardless of side chain length (consistent with the helix geometry). Interestingly, only Lys was capable of supporting Zbb-Xaa (i, i+3) interactions, whereas any Xaa side chain length supported Zbb-Xaa (i, i+4) interactions. In particular, the magnitude of both Zbb(-)-Lys (i, i+4) and Zbb(-)-Orn (i, i+4) interaction energies followed the trend Asp > Glu > Aad. Side chain conformational analysis by molecular mechanics calculations showed that the Zbb-Xaa (i, i+3) interactions involved the χ(1) dihedral combination (g+, g+) for the i and i+3 residues, whereas the Zbb-Xaa (i, i+4) interactions were supported by the χ(1) dihedral combination (t, g+) for the i and i+4 residues. These calculated low energy conformers were consistent with conformations of intrahelical Asp-Lys and Glu-Lys salt bridges in a nonredundant protein structure database. These results suggest that Asp and Glu provide natural variation, and lengthening the Glu side chain further to Aad does not furnish additional characteristics that Glu cannot supply.

Structural Elucidation of Dextran Degradation Mechanism by Streptococcus mutans Dextranase Belonging to Glycoside Hydrolase Family 66

Dextranase is an enzyme that hydrolyzes dextran α-1,6 linkages. Streptococcus mutans dextranase belongs to glycoside hydrolase family 66, producing isomaltooligosaccharides of various sizes and consisting of at least five amino acid sequence regions. The crystal structure of the conserved fragment from Gln(100) to Ile(732) of S. mutans dextranase, devoid of its N- and C-terminal variable regions, was determined at 1.6 Å resolution and found to contain three structural domains. Domain N possessed an immunoglobulin-like β-sandwich fold; domain A contained the enzyme's catalytic module, comprising a (β/α)(8)-barrel; and domain C formed a β-sandwich structure containing two Greek key motifs. Two ligand complex structures were also determined, and, in the enzyme-isomaltotriose complex structure, the bound isomaltooligosaccharide with four glucose moieties was observed in the catalytic glycone cleft and considered to be the transglycosylation product of the enzyme, indicating the presence of four subsites, -4 to -1, in the catalytic cleft. The complexed structure with 4',5'-epoxypentyl-α-d-glucopyranoside, a suicide substrate of the enzyme, revealed that the epoxide ring reacted to form a covalent bond with the Asp(385) side chain. These structures collectively indicated that Asp(385) was the catalytic nucleophile and that Glu(453) was the acid/base of the double displacement mechanism, in which the enzyme showed a retaining catalytic character. This is the first structural report for the enzyme belonging to glycoside hydrolase family 66, elucidating the enzyme's catalytic machinery.

NMR Observable-Based Structure Refinement of DAP12-NKG2C Activating Immunoreceptor Complex in Explicit Membranes

NMR observables, such as NOE-based distance measurements, are increasingly being used to characterize membrane protein structures. However, challenges in membrane protein NMR studies often yield a relatively small number of such restraints that can create ambiguities in defining critical side chain-side chain interactions. In the recent solution NMR structure of the DAP12-NKG2C immunoreceptor transmembrane helix complex, five functionally required interfacial residues (two Asps and two Thrs in the DAP12 dimer and one Lys in NKG2C) display a surprising arrangement in which one Asp side chain faces the membrane hydrophobic core. To explore whether these side-chain interactions are energetically optimal, we used the published distance restraints for molecular dynamics simulations in explicit micelles and bilayers. The structures refined by this protocol are globally similar to the published structure, but the side chains of those five residues form a stable network of salt bridges and hydrogen bonds, leaving the Asp side chain shielded from the hydrophobic core, which is also consistent with available experimental observations. Moreover, the simulations show similar short-range interactions between the transmembrane complex and lipid/detergent molecules in micelles and bilayers, respectively. This study illustrates the efficacy of NMR membrane protein structure refinements in explicit membrane systems.

Structure and assembly of Gram-positive bacterial pili: unique covalent polymers

Bacterial pili are long, multi-subunit protein assemblies that extend from bacterial surfaces, mediating adhesion and colonisation. The recently characterised pili expressed by Gram-positive pathogens represent a novel variation; completely covalent polymers in which sortase-mediated isopeptide bonds link successive pilin subunits. Recent structural studies of the component pilins have revealed a common pattern of tandem immunoglobulin (Ig)-like domains, joined end-on-end. This long thin assembly is further stabilised by autocatalytically generated isopeptide bond crosslinks within the domains, joining Lys and Asn(or Asp) side chains. Specialised subunits at the tip and the base complete the assembly, with the tip pilins presenting novel adhesive structures.

Modeling the Interaction between Integrin-Binding Peptide (RGD) and Rutile Surface: The Effect of Cation Mediation on Asp Adsorption

The binding of a negatively charged residue, aspartic acid (Asp) in tripeptide arginine-glycine-aspartic acid, onto a negatively charged hydroxylated rutile (110) surface in aqueous solution, containing divalent (Mg(2+), Ca(2+), or Sr(2+)) or monovalent (Na(+), K(+), or Rb(+)) cations, was studied by molecular dynamics (MD) simulations. The results indicate that ionic radii and charges will significantly affect the hydration, adsorption geometry, and distance of cations from the rutile surface, thereby regulating the Asp/rutile binding mode. The adsorption strength of monovalent cations on the rutile surface in the order Na(+) > K(+) > Rb(+) shows a "reverse" lyotropic trend, while the divalent cations on the same surface exhibit a "regular" lyotropic behavior with decreasing crystallographic radii (the adsorption strength of divalent cations: Sr(2+) > Ca(2+) > Mg(2+)). The Asp side chain in NaCl, KCl, and RbCl solutions remains stably H-bonded to the surface hydroxyls and the inner-sphere adsorbed compensating monovalent cations act as a bridge between the COO(-) group and the rutile, helping to "trap" the negatively charged Asp side chain on the negatively charged surface. In contrast, the mediating divalent cations actively participate in linking the COO(-) group to the rutile surface; thus the Asp side chain can remain stably on the rutile (110) surface, even if it is not involved in any hydrogen bonds with the surface hydroxyls. Inner- and outer-sphere geometries are all possible mediation modes for divalent cations in bridging the peptide to the rutile surface.

Modeling the Interaction between Integrin-Binding Peptide (RGD) and Rutile Surface: The Effect of Na+ on Peptide Adsorption

The dynamics of a single tripeptide Arg-Gly-Asp (RGD) adsorbing onto negatively charged hydroxylated rutile (110) surface in aqueous solution was studied using molecular dynamics (MD) simulations. The results indicate that the adsorbed Na{sup +} ions play an important role in determining the binding geometry of RGD. With an initial 'horseshoe' configuration, the charged side groups (COO{sup -} and NH{sub 2}) of the peptide are able to interact with the surface through direct hydrogen bonds (H bonds) in the very early stage of adsorption. The Na{sup +} ions approach the positively charged Arg side chain, competing with the Arg side chain for adsorption to the negatively charged hydroxyl oxygen. In coordination with the structural adjustment of the peptide, the Arg residue is driven to detach from the rutile surface. In contrast, the Na+ ions in close proximity to the negatively charged Asp side chain contribute to the binding of the COO{sup -} group on the surface, helping the carboxyl oxygen not involved in COO{sup -}-surface H bonds to orientate toward the hydroxyl hydrogens. Once both carboxyl oxygens form enough H bonds with the hydroxyl hydrogens, the redundant ions move toward a more favorable adsorption site.

Low temperature FTIR spectroscopy provides new insights in the pH-dependent proton pathway of proteorhodopsin

In the presented study the low pH photocycle of proteorhodopsin is extensively investigated by means of low temperature FTIR spectroscopy. Besides the already well-known characteristics of the all-trans and 13-cis retinal vibrations the 77K difference spectrum at pH 5.1 shows an additional negative signal at 1744cm−1 which is interpreted as indicator for the L state. The subsequent photocycle steps are investigated at temperatures higher than 200K. The combination of visible and FTIR spectroscopy enabled us to observe that the deprotonation of the Schiff base is linked to the protonation of an Asp or Glu side chain — the new proton acceptor under acidic conditions. The difference spectra of the late intermediates are characterized by large amide I changes and two further bands ((−)1751cm−1/(+)1725cm-1) in the spectral region of the Asp/Glu ν(C=O) vibrations. The band position of the negative signature points to a transient deprotonation of Asp-97. In addition, the pH dependence of the acidic photocycle was investigated. The difference spectra at pH 5.5 show distinct differences connected to changes in the protonation state of key residues. Based on our data we propose a three-state model that explains the complex pH dependence of PR.

Structural and Enzymatic Insights into Caspase-2 Protein Substrate Recognition and Catalysis

Caspase-2, the most evolutionarily conserved member in the human caspase family, may play important roles in stress-induced apoptosis, cell cycle regulation, and tumor suppression. In biochemical assays, caspase-2 uniquely prefers a pentapeptide (such as VDVAD) rather than a tetrapeptide, as required for efficient cleavage by other caspases. We investigated the molecular basis for pentapeptide specificity using peptide analog inhibitors and substrates that vary at the P5 position. We determined the crystal structures of apo caspase-2, caspase-2 in complex with peptide inhibitors VDVAD-CHO, ADVAD-CHO, and DVAD-CHO, and a T380A mutant of caspase-2 in complex with VDVAD-CHO. Two residues, Thr-380 and Tyr-420, are identified to be critical for the P5 residue recognition; mutation of the two residues reduces the catalytic efficiency by about 4- and 40-fold, respectively. The structures also provide a series of snapshots of caspase-2 in different catalytic states, shedding light on the mechanism of capase-2 activation, substrate binding, and catalysis. By comparing the apo and inhibited caspase-2 structures, we propose that the disruption of a non-conserved salt bridge between Glu-217 and the invariant Arg-378 is important for the activation of caspase-2. These findings broaden our understanding of caspase-2 substrate specificity and catalysis.

Structural basis for endosomal recruitment of ESCRT-I by ESCRT-0 in yeast

The ESCRT-0 and ESCRT-I complexes coordinate the clustering of ubiquitinated cargo with intralumenal budding of the endosomal membrane, two essential steps in vacuolar/lysosomal protein sorting from yeast to humans. The 1.85-Å crystal structure of interacting regions of the yeast ESCRT-0 and ESCRT-I complexes reveals that PSDP motifs of the Vps27 ESCRT-0 subunit bind to a novel electropositive N-terminal site on the UEV domain of the ESCRT-I subunit Vps23 centred on Trp16. This novel site is completely different from the C-terminal part of the human UEV domain that binds to P(S/T)AP motifs of human ESCRT-0 and HIV-1 Gag. Disruption of the novel PSDP-binding site eliminates the interaction in vitro and blocks enrichment of Vps23 in endosome-related class E compartments in yeast cells. However, this site is non-essential for sorting of the ESCRT cargo Cps1. Taken together, these results show how a conserved motif/domain pair can evolve to use strikingly different binding modes in different organisms.

Protonation/deprotonation effects on the stability of the Trp-cage miniprotein

The Trp-cage miniprotein is a 20 amino acid peptide that exhibits many of the properties of globular proteins. In this protein, the hydrophobic core is formed by a buried Trp side chain. The folded state is stabilized by an ion pair between aspartic acid and an arginine side chain. The effect of protonating the aspartic acid on the Trp-cage miniprotein folding/unfolding equilibrium is studied by explicit solvent molecular dynamics simulations of the protein in the charged and protonated Asp9 states. Unbiased Replica Exchange Molecular Dynamics (REMD) simulations, spanning a wide temperature range, are carried out to the microsecond time scale, using the AMBER99SB forcefield in explicit TIP3P water. The protein structural ensembles are studied in terms of various order parameters that differentiate the folded and unfolded states. We observe that in the folded state the root mean square distance (rmsd) from the backbone of the NMR structure shows two highly populated basins close to the native state with peaks at 0.06 nm and 0.16 nm, which are consistent with previous simulations using the same forcefield. The fraction of folded replicas shows a drastic decrease because of the absence of the salt bridge. However, significant populations of conformations with the arginine side chain exposed to the solvent, but within the folded basin, are found. This shows the possibility to reach the folded state without formation of the ion pair. We also characterize changes in the unfolded state. The equilibrium populations of the folded and unfolded states are used to characterize the thermodynamics of the system. We find that the change in free energy difference due to the protonation of the Asp amino acid is 3 kJ mol(-1) at 297 K, favoring the charged state, and resulting in ΔpK(1) = 0.5 units for Asp9. We also study the differences in the unfolded state ensembles for the two charge states and find significant changes at low temperature, where the protonated Asp side chain makes multiple hydrogen bonds to the protein backbone.

Modeling the Cu+ Binding in the 1−16 Region of the Amyloid-β Peptide Involved in Alzheimer’s Disease

The coordination of copper to the amyloid-β (1-16) (Aβ) peptide has been investigated because of its relevance for understanding Cu redox activity when the ion is embedded in peptides involved in neurodegenerative diseases. In this work, several reasonable models of Cu(+) coordination were built on the basis of experimental information and investigated by first-principles molecular dynamics simulations in the Car-Parrinello scheme. The propensity of a linear Nδ (His)-Cu-Nδ (His) coordination for Cu(+) is shown by all the models investigated here, with distortions due to weak interactions with the carbonyl O of His 6 and His 13 and with the amide N of His 14. Though the His 6-Cu-His 14 linear coordination is favored in truncated models, the His 13-Cu-His 14 linear coordination is favored by interactions present in the complete solvated and in vacuo models of Cu-Aβ (1-16). These interactions include steric hindrance for the expulsion of His 13, hydrogen bonds between Asp and His side chains and a network of electrostatic interactions stabilizing two separated 1-10 and 11-16 peptide regions. The role of linear His 13-Cu-His 14 coordination in stabilizing Cu(I) and in increasing the Cu(II)/Cu(I) reorganization energy can be therefore modulated by boundary conditions acting on the Aβ (1-16) ligand.