Observation of the closing of individual hydrogen bonds during TFE-induced helix formation in a peptide

Helix formation of an S-peptide analog, comprising the first 20 residues of Ribonuclease A and two additional N-terminal residues, was studied by measuring hydrogen bond (H-bond) (h3)J(NC') scalar couplings as a function of 2,2,2-trifluoroethanol (TFE) concentration. The (h3)J(NC') couplings give direct evidence for the closing of individual backbone N-H***O = C H-bonds during the TFE-induced formation of secondary structure. Whereas no (h3)J(NC') correlations could be detected without TFE, alpha-helical (i,i +4) H-bond correlations were observed for the amides of residues A5 to M15 in the presence of TFE. The analysis of individual coupling constants indicates that alpha-helix formation starts at the center of the S-peptide around residue E11 and proceeds gradually from there to both peptide ends as the TFE concentration is increased. At 60% to 90% TFE, well-formed alpha-helical H-bonds were observed for the amides hydrogens of residues K9 to Q13, whereas H-bonds of residues T5 to A8, H14, and M15 are affected by fraying. No intramolecular backbone H-bonds are present at and beyond the putative helix stop signal D16. As the (h3)J(NC') constants represent ensemble averages and the dependence of (h3)J(NC') on H-bond lengths is very steep, the size of the individual (h3)J(NC') coupling constants can be used as a measure for the population of a closed H-bond. These individual populations are in agreement with results derived from the Lifson-Roig theory for coil-to-helix transitions. The present work shows that the closing of individual H-bonds during TFE-induced helix formation can be monitored by changes in the size of H-bond scalar couplings.     

Quantification of H/D Isotope Effects on Protein Hydrogen-bonds by h3JNC′ and 1JNC′ Couplings and Peptide Group 15N and 13C′ Chemical Shifts

The effect of hydrogen/deuterium exchange on protein hydrogen bond coupling constants (h3)J(NC') has been investigated in the small globular protein ubiquitin. The couplings across deuterated or protonated hydrogen bonds were measured by a long-range quantitative HA(CACO)NCO experiment. The analysis is combined with a determination of the H(N)/D(N) isotope effect on the amide group (1)J(NC') couplings and the (15)N and (13)C' chemical shifts. On average, H-bond deuteration exchange weakens (h3)J(NC') and strengthens (1)J(NC') couplings. A correlation is found between the size of the (15)N isotope shift, the (15)N chemical shift, and the (h3)J(NC') coupling constants. The data are consistent with a reduction of donor-acceptor overlap as expected from the classical Ubbelohde effect and the common understanding that H(N)/D(N) exchange leads to a shortening of the N-hydron bond length.     

Peptaibol zervamicin IIb structure and dynamics refinement from transhydrogen bond J couplings

Zervamicin IIB (Zrv-IIB) is a channel-forming peptaibol antibiotic of fungal origin. The measured transhydrogen bond (3h)J(NC') couplings in methanol solution heaving average value of -0.41 Hz indicate that the stability of the Zrv-IIB helix in this milieu is comparable to the stability of helices in globular proteins. The N-terminus of the peptide forms an alpha-helix, whereas 3(10)-helical hydrogen bonds stabilize the C-terminus. However, two weak transhydrogen bond peaks are observed in a long-range HNCO spectrum for HN Aib(12). Energy calculations using the Empirical Conformation Energy Program for Peptides (ECEPP)/2 force field and the implicit solvent model show that the middle of the peptide helix accommodates a bifurcated hydrogen bond that is simultaneously formed between HN Aib(12) and CO Leu(8) and CO Aib(9). Several lowered (3h)J(NC') on a polar face of the helix correlate with the conformational exchange process observed earlier and imply dynamic distortions of a hydrogen bond pattern with the predominant population of a properly folded helical structure. The refined structure of Zrv-IIB on the basis of the observed hydrogen bond pattern has a small ( approximately 20 degrees ) angle of helix bending that is virtually identical to the angle of bending in dodecylphosphocholine (DPC) micelles, indicating the stability of a hinge region in different environments. NMR parameters ((1)HN chemical shifts and transpeptide bond (1)J(NC') couplings) sensitive to hydrogen bonding along with the solvent accessible surface area of carbonyl oxygens indicate a large polar patch on the convex side of the helix formed by three exposed backbone carbonyls of Aib(7), Aib(9), and Hyp(10) and polar side chains of Hyp(10), Gln(11), and Hyp(13). The unique structural features, high helix stability and the enhanced polar patch, set apart Zrv-IIB from other peptaibols (for example, alamethicin) and possibly underlie its biological and physiological properties.     

Molecular details of Itk activation by prolyl isomerization and phospholigand binding: the NMR structure of the Itk SH2 domain bound to a phosphopeptide

The Src homology 2 (SH2) domain of interleukin-2 tyrosine kinase (Itk) is a critical component of the regulatory apparatus controlling the activity of this immunologically important enzyme. To gain insight into the structural features associated with the activated form of Itk, we have solved the NMR structure of the SH2 domain bound to a phosphotyrosine-containing peptide (pY) and analyzed changes in trans-hydrogen bond scalar couplings ((3h)J(NC')) that result from pY binding. Isomerization of a single prolyl imide bond in this domain is responsible for simultaneous existence of two distinct SH2 conformers. Prolyl isomerization directs ligand recognition: the trans conformer preferentially binds pY. The structure of the SH2/pY complex provides insight into the ligand specificity; the BG loop in the ligand-free trans SH2 conformer is pre-arranged for optimal contacts with the pY+3 residue of the ligand. Analysis of (3h)J(NC') couplings arising from hydrogen bonds has revealed propagation of structural changes from the pY binding pocket to the CD loop containing conformationally heterogeneous proline as well as to the alphaB helix, on the opposite site of the domain. These findings offer a structural framework for understanding the roles of prolyl isomerization and pY binding in Itk regulation.     

Ligand-induced strain in hydrogen bonds of the c-Src SH3 domain detected by NMR

Changes in the molecular conformation of proteins can result from a variety of perturbations, and can play crucial roles in the regulation of biological activity. A new solution NMR method has been applied to monitor ligand-induced changes in hydrogen bond geometry in the chicken c-Src SH3 domain. The structural response of this domain to ligand binding has been investigated by measuring trans-hydrogen bond (15)N-(13)C' scalar couplings in the free state and when bound to the high affinity class I ligand RLP2, containing residues RALPPLPRY. A comparison between hydrogen bonds in high resolution X-ray structures of this domain and those observed via (h3)J(NC') couplings in solution shows remarkable agreement. Two backbone-to-side-chain hydrogen bonds are observed in solution, and each appears to play a role in stabilization of loop structure. Reproducible ligand-induced changes in trans-hydrogen bond scalar couplings are observed across the domain that translate into changes in hydrogen bond length ranging between 0.02 to 0.12 A. The observed changes can be rationalized by an induced fit mechanism in which hydrogen bonds across the protein participate in a compensatory response to forces imparted at the protein-ligand interface. Upon ligand binding, mutual intercalation of the two Leu-Pro segments of the ligand between three aromatic side-chains protruding from the SH3 surface wedges apart secondary structural elements within the SH3 domain. This disruption is transmitted in a domino-like effect across the domain through networks of hydrogen bonded peptide planes. The unprecedented resolution obtained demonstrates the ability to characterize subtle structural rearrangements within a protein upon perturbation, and represents a new step in the endeavor to understand how hydrogen bonds contribute to the stabilization and function of biological macromolecules.     

Solvent-induced differentiation of protein backbone hydrogen bonds in calmodulin

In apo and holoCaM, almost half of the hydrogen bonds (H-bonds) at the protein backbone expected from the corresponding NMR or X-ray structures were not detected by h3JNC' couplings. The paucity of the h3JNC' couplings was considered in terms of dynamic features of these structures. We examined a set of seven proteins and found that protein-backbone H-bonds form two groups according to the h3JNC' couplings measured in solution. H-bonds that have h3JNC' couplings above the threshold of 0.2 Hz show distance/angle correlation among the H-bond geometrical parameters, and appear to be supported by the backbone dynamics in solution. The other H-bonds have no such correlation; they populate the water-exposed and flexible regions of proteins, including many of the CaM helices. The observed differentiation in a dynamical behavior of backbone H-bonds in apo and holoCaM appears to be related to protein functions.     

Peptide bond distortions from planarity: new insights from quantum mechanical calculations and peptide/protein crystal structures

By combining quantum-mechanical analysis and statistical survey of peptide/protein structure databases we here report a thorough investigation of the conformational dependence of the geometry of peptide bond, the basic element of protein structures. Different peptide model systems have been studied by an integrated quantum mechanical approach, employing DFT, MP2 and CCSD(T) calculations, both in aqueous solution and in the gas phase. Also in absence of inter-residue interactions, small distortions from the planarity are more a rule than an exception, and they are mainly determined by the backbone ψ dihedral angle. These indications are fully corroborated by a statistical survey of accurate protein/peptide structures. Orbital analysis shows that orbital interactions between the σ system of C(α) substituents and the π system of the amide bond are crucial for the modulation of peptide bond distortions. Our study thus indicates that, although long-range inter-molecular interactions can obviously affect the peptide planarity, their influence is statistically averaged. Therefore, the variability of peptide bond geometry in proteins is remarkably reproduced by extremely simplified systems since local factors are the main driving force of these observed trends. The implications of the present findings for protein structure determination, validation and prediction are also discussed.     

Rapid structural fluctuations of the free HIV protease flaps in solution: Relationship to crystal structures and comparison with predictions of dynamics calculations

Crystal structures have shown that the HIV-1 protease flaps, domains that control access to the active site, are closed when the active site is occupied by a ligand. Although flap structures ranging from closed to semi-open are observed in the free protease, crystal structures reveal that even the semi-open flaps block access to the active site, indicating that the flaps are mobile in solution. The goals of this paper are to characterize the secondary structure and fast (sub-ns) dynamics of the flaps of the free protease in solution, to relate these results to X-ray structures and to compare them with predictions of dynamics calculations. To this end we have obtained nearly complete backbone and many sidechain signal assignments of a fully active free-protease construct that is stabilized against autoproteolysis by three point mutations. The secondary structure of this protein was characterized using the chemical shift index, measurements of (3h)J(NC') couplings across hydrogen bonds, and NOESY connectivities. Analysis of these measurements indicates that the protease secondary structure becomes irregular near the flap tips, residues 49-53. Model-free analysis of (15)N relaxation parameters, T(1), T(2) (T(1rho)) and (15)N-[(1)H] NOE, shows that residues in the flap tips are flexible on the sub-ns time scale, in contrast with previous observations on the inhibitor-bound protease. These results are compared with theoretical predictions of flap dynamics and the possible biological significance of the sub-ns time scale dynamics of the flap tips is discussed.     

A simple and sensitive experiment for measurement of JCC couplings between backbone carbonyl and methyl carbons in isotopically enriched proteins

Summary A simple 2D difference experiment is described that allows quantitative measurement of ¹³C–¹³C J couplings between backbone carbonyl and side-chain carbons. Precise ³J_CC values were measured from data recorded in just 2 h for a 1-mM solution of the 20-kD complex between the protein calmodulin and a 26-residue synthetic peptide. The J couplings aid in determining the ¹ angles of valine, isoleucine and threonine residues, and in making stereospecific assignments of the Val C methyl groups. Error analysis indicates that the uncertainty in the derived J couplings is generally less than ca. 0.3 Hz.     

Solution structure of Ca(2+)-calmodulin reveals flexible hand-like properties of its domains.

The solution structure of Ca(2+)-ligated calmodulin is determined from residual dipolar couplings measured in a liquid crystalline medium and from a large number of heteronuclear J couplings for defining side chains. Although the C-terminal domain solution structure is similar to the X-ray crystal structure, the EF hands of the N-terminal domain are considerably less open. The substantial differences in interhelical angles correspond to negligible changes in short interproton distances and, therefore, cannot be identified by comparison of NOEs and X-ray data. NOE analysis, however, excludes a two-state equilibrium in which the closed apo conformation is partially populated in the Ca(2+)-ligated state. The difference between the crystal and solution structures of Ca(2+)-calmodulin indicates considerable backbone plasticity within the domains of calmodulin, which is key to their ability to bind a wide range of targets. In contrast, the vast majority of side chains making up the target binding surface are locked into the same chi(1) rotameric states as in complexes with target peptide.     

Impact of N-terminal acetylation of α-synuclein on its random coil and lipid binding properties

N-Terminal acetylation of α-synuclein (aS), a protein implicated in the etiology of Parkinson's disease, is common in mammals. The impact of this modification on the protein's structure and dynamics in free solution and on its membrane binding properties has been evaluated by high-resolution nuclear magnetic resonance and circular dichroism (CD) spectroscopy. While no tetrameric form of acetylated aS could be isolated, N-terminal acetylation resulted in chemical shift perturbations of the first 12 residues of the protein that progressively decreased with the distance from the N-terminus. The directions of the chemical shift changes and small changes in backbone (3)J(HH) couplings are consistent with an increase in the α-helicity of the first six residues of aS, although a high degree of dynamic conformational disorder remains and the helical structure is sampled <20% of the time. Chemical shift and (3)J(HH) data for the intact protein are virtually indistinguishable from those recorded for the corresponding N-terminally acetylated and nonacetylated 15-residue synthetic peptides. An increase in α-helicity at the N-terminus of aS is supported by CD data on the acetylated peptide and by weak medium-range nuclear Overhauser effect contacts indicative of α-helical character. The remainder of the protein has chemical shift values that are very close to random coil values and indistinguishable between the two forms of the protein. No significant differences in the fibrillation kinetics were observed between acetylated and nonacetylated aS. However, the lipid binding properties of aS are strongly impacted by acetylation and exhibit distinct behavior for the first 12 residues, indicative of an initiation role for the N-terminal residues in an "initiation-elongation" process of binding to the membrane.     

NMR of hydrogen bonding in cold-shock protein A and an analysis of the influence of crystallographic resolution on comparisons of hydrogen bond lengths

Hydrogen bonding in cold-shock protein A of Escherichia coli has been investigated using long-range HNCO spectroscopy. Nearly half of the amide protons involved in hydrogen bonds in solution show no measurable protection from exchange in water, cautioning against a direct correspondence between hydrogen bonding and hydrogen exchange protection. The N to O atom distance across a hydrogen bond, R(NO), is related to the size of the (3h)J(NC') trans hydrogen bond coupling constant and the amide proton chemical shift. Both NMR parameters show poorer agreement with the 2.0-A resolution X-ray structure of the cold-shock protein studied by NMR than with a 1.2-A resolution X-ray structure of a homologous cold-shock protein from the thermophile B. caldolyticus. The influence of crystallographic resolution on comparisons of hydrogen bond lengths was further investigated using a database of 33 X-ray structures of ribonuclease A. For highly similar structures, both hydrogen bond R(NO) distance and Calpha coordinate root mean square deviations (RMSD) show systematic increases as the resolution of the X-ray structure used for comparison decreases. As structures diverge, the effects of coordinate errors on R(NO) distance and Calpha coordinate root mean square deviations become progressively smaller. The results of this study are discussed with regard to the influence of data precision on establishing structure similarity relationships between proteins.     

Exact Solutions for Internuclear Vectors and Backbone Dihedral Angles from NH Residual Dipolar Couplings in Two Media, and their Application in a Systematic Search Algorithm for Determining Protein Backbone Structure

We have derived a quartic equation for computing the direction of an internuclear vector from residual dipolar couplings (RDCs) measured in two aligning media, and two simple trigonometric equations for computing the backbone (phi,psi) angles from two backbone vectors in consecutive peptide planes. These equations make it possible to compute, exactly and in constant time, the backbone (phi,psi) angles for a residue from RDCs in two media on any single backbone vector type. Building upon these exact solutions we have designed a novel algorithm for determining a protein backbone substructure consisting of alpha-helices and beta-sheets. Our algorithm employs a systematic search technique to refine the conformation of both alpha-helices and beta-sheets and to determine their orientations using exclusively the angular restraints from RDCs. The algorithm computes the backbone substructure employing very sparse distance restraints between pairs of alpha-helices and beta-sheets refined by the systematic search. The algorithm has been demonstrated on the protein human ubiquitin using only backbone NH RDCs, plus twelve hydrogen bonds and four NOE distance restraints. Further, our results show that both the global orientations and the conformations of alpha-helices and beta-strands can be determined with high accuracy using only two RDCs per residue. The algorithm requires, as its input, backbone resonance assignments, the identification of alpha-helices and beta-sheets as well as sparse NOE distance and hydrogen bond restraints.     

Stress and strain in staphylococcal nuclease.

Protein molecules generally adopt a tertiary structure in which all backbone and side chain conformations are arranged in local energy minima; however, in several well-refined protein structures examples of locally strained geometries, such as cis peptide bonds, have been observed. Staphylococcal nuclease A contains a single cis peptide bond between residues Lys 116 and Pro 117 within a type VIa beta-turn. Alternative native folded forms of nuclease A have been detected by NMR spectroscopy and attributed to a mixture of cis and trans isomers at the Lys 116-Pro 117 peptide bond. Analyses of nuclease variants K116G and K116A by NMR spectroscopy and X-ray crystallography are reported herein. The structure of K116A is indistinguishable from that of nuclease A, including a cis 116-117 peptide bond (92% populated in solution). The overall fold of K116G is also indistinguishable from nuclease A except in the region of the substitution (residues 112-117), which contains a predominantly trans Gly 116-Pro 117 peptide bond (80% populated in solution). Both Lys and Ala would be prohibited from adopting the backbone conformation of Gly 116 due to steric clashes between the beta-carbon and the surrounding residues. One explanation for these results is that the position of the ends of the residue 112-117 loop only allow trans conformations where the local backbone interactions associated with the phi and psi torsion angles are strained. When the 116-117 peptide bond is cis, less strained backbone conformations are available. Thus the relaxation of the backbone strain intrinsic to the trans conformation compensates for the energetically unfavorable cis X-Pro peptide bond. With the removal of the side chain from residue 116 (K116G), the backbone strain of the trans conformation is reduced to the point that the conformation associated with the cis peptide bond is no longer favorable.     

Identification of hydrogen bonds to the Rieske cluster through the weakly coupled nitrogens detected by electron spin echo envelope modulation spectroscopy

The interaction of the reduced[2Fe-2S] cluster of isolated Rieske fragment from the bc1 complex of Rhodobacter sphaeroides with nitrogens (14N and 15N) from the local protein environment has been studied by X- and S-band pulsed EPR spectroscopy. The two-dimensional electron spin echo envelope modulation spectra of uniformly 15N-labeled protein show two well resolved cross-peaks with weak couplings of approximately 0.3-0.4 and 1.1 MHz in addition to couplings in the range of 6-8 MHz from two coordinating Ndelta of histidine ligands. The quadrupole coupling constants for weakly coupled nitrogens determined from S-band electron spin echo envelope modulation spectra identify them as Nepsilon of histidine ligands and peptide nitrogen (Np), respectively. Analysis of the line intensities in orientation-selected S-band spectra indicated that Np is the backbone N-atom of Leu-132 residue. The hyperfine couplings from Nepsilon and Np demonstrate the predominantly isotropic character resulting from the transfer of unpaired spin density onto the 2s orbitals of the nitrogens. Spectra also show that other peptide nitrogens in the protein environment must carry a 5-10 times smaller amount of spin density than the Np of Leu-132 residue. The appearance of the excess unpaired spin density on the Np of Leu-132 residue indicates its involvement in hydrogen bond formation with the bridging sulfur of the Rieske cluster. The configuration of the hydrogen bond therefore provides a preferred path for spin density transfer. Observation of similar splittings in the 15N spectra of other Rieske-type proteins and [2Fe-2S] ferredoxins suggests that a hydrogen bond between the bridging sulfur and peptide nitrogen is a common structural feature of [2Fe-2S] clusters.     

Structure and conformation of linear peptides. V. Structure of L-prolyl-glycyl-glycine

The tripeptide, L-prolyl-glycyl-glycine, crystallizes in the trigonal space group P3(2), with a = b = 8.682(2) A, c = 12.008(2) and Z = 3. The structure was solved by direct methods and refined to an R-value of 0.07 for 727 reflections (I greater than 1.0 sigma). The molecule exists as a zwitterion in the crystal. The peptide units are trans and show significant deviations from planarity (omega 1 = 169.7 degrees, omega 2 = -170.1 degrees). The peptide backbone adopts a left-handed helical conformation similar to that of polyglycine II and polyproline II.     

Strategy for supplementing structure calculations using limited data with hydrophobic distance restraints

Introductory biochemistry texts often note that the fold of a protein is completely defined when the dihedral angles phi and psi are known for each amino acid. This assertion was examined with torsion angle dynamics and simulated annealing (TAD/SA) calculations of protein G using only dihedral angle restraints. When all dihedral angles were restrained to within 1 degrees of the values of the X-ray structure, the TAD/SA structures gave a backbone root mean square deviation to the target of 4 A. Factors that contributed to divergence from the correct solution include deviations of peptide bonds from planarity, internal conflicts resulting from the nonuniform energies of different phi, psi combinations, and relaxation to extended conformations in the absence of long-range constraints. Simulations including hydrogen-bond restraints showed that even a few long-range contacts constrain the fold better than a complete set of accurate dihedral restraints. A procedure is described for TAD/SA calculations using hydrogen-bond restraints, idealized dihedral restraints for residues in regular secondary structures, and "hydrophobic distance restraints" derived from the positions of hydrophobic residues in the amino acid sequence. The hydrogen-bond restraints are treated as inviolable, whereas violated hydrophobic restraints are removed following reduction of restraint upper bounds from 2 to 1 times the predicted radius of gyration. The strategy was tested with simulated restraints from X-ray structures of proteins from different fold classes and NMR data for cold shock protein A that included only backbone chemical shifts and hydrogen bonds obtained from a long-range HNCO experiment.     

Assignment of the natural abundance 13C spectrum of proteins using 13C 1H-detected heteronuclear multiple-bond correlation NMR spectroscopy: structural information and stereospecific assignments from two- and three-bond carbon-hydrogen coupling constants.

Proton-detected heteronuclear multiple-bond 1H-13C correlations (HMBC) previously have been used for assignment purposes in a variety of isotopically enriched proteins. In the present study it is demonstrated that the technique yields an almost complete assignment of the natural abundance 13C spectrum of the protein basic pancreatic trypsin inhibitor (BPTI). In addition, the technique permits additional 1H assignments to be made for this well-studied protein. The intensities of observed correlations permit rough estimates to be made of 2J(C,H) and 3J(C,H) coupling constants. These couplings can be used for conformational studies of both the side chains and the backbone. Intra- and interresidue coupling between C alpha H and the carbonyl carbon provides information about the backbone angles psi and phi. Side-chain conformations can be determined from both two- and three-bond carbon-hydrogen coupling constants. The present study of BPTI together with its known high-precision solution structure yields an experimental correlation between resonance intensities and secondary structure. The spectra show the potential of the method in analyzing 13C NMR spectra of nonenriched proteins. The method yields 13C NMR chemical shifts, which are versatile parameters to be used to monitor structural changes, titrations, etc.     

Automated extraction of backbone deuteration levels from amide H/2H mass spectrometry experiments

A Fourier deconvolution method has been developed to explicitly determine the amount of backbone amide deuterium incorporated into protein regions or segments by hydrogen/deuterium (H/D) exchange with high-resolution mass spectrometry. Determination and analysis of the level and number of backbone amide exchanging in solution provide more information about the solvent accessibility of the protein than do previous centroid methods, which only calculate the average deuterons exchanged. After exchange, a protein is digested into peptides as a way of determining the exchange within a local area of the protein. The mass of a peptide upon deuteration is a sum of the natural isotope abundance, fast exchanging side-chain hydrogens (present in MALDI-TOF H/2H data) and backbone amide exchange. Removal of the components of the isotopic distribution due to the natural isotope abundances and the fast exchanging side-chains allows for a precise quantification of the levels of backbone amide exchange, as is shown by an example from protein kinase A. The deconvoluted results are affected by overlapping peptides or inconsistent mass envelopes, and evaluation procedures for these cases are discussed. Finally, a method for determining the back exchange corrected populations is presented, and its effect on the data is discussed under various circumstances.     

An improved algorithm for MFR fragment assembly

A method for generating protein backbone models from backbone only NMR data is presented, which is based on molecular fragment replacement (MFR). In a first step, the PDB database is mined for homologous peptide fragments using experimental backbone-only data i.e. backbone chemical shifts (CS) and residual dipolar couplings (RDC). Second, this fragment library is refined against the experimental restraints. Finally, the fragments are assembled into a protein backbone fold using a rigid body docking algorithm using the RDCs as restraints. For improved performance, backbone nuclear Overhauser effects (NOEs) may be included at that stage. Compared to previous implementations of MFR-derived structure determination protocols this model-building algorithm offers improved stability and reliability. Furthermore, relative to CS-ROSETTA based methods, it provides faster performance and straightforward implementation with the option to easily include further types of restraints and additional energy terms.