Hydrogen Bond Constraints Research Articles

The Evolutionarily Conserved Family of Cyanovirin-N Homologs: Structures and Carbohydrate Specificity

Solution structures for three members of the recently discovered cyanovirin-N (CV-N) homolog family of lectins have been determined. Cyanovirin-N homologs (CVNHs) from Tuber borchii, Ceratopteris richardii, and Neurospora crassa, representing each of the three phylogenetic groups, were selected. All proteins exhibit the same fold, and the overall structures resemble that of the founding member of the family, CV-N, albeit with noteworthy differences in loop conformation and detailed local structure. Since no data are available regarding the proteins' function or their natural ligands, extensive carbohydrate-binding studies were conducted. We delineated ligand-binding sites on all three proteins by nuclear magnetic resonance and identified which sugars interact by array screening. The number and location of binding sites vary for the three proteins, and different ligand specificities exist. Potential physiological roles for two family members, TbCVNH and NcCVNH, were probed in nutrition deprivation experiments that suggest a possible involvement of these proteins in lifestyle-related responses.

Three-Dimensional Structure of Conotoxin tx3a: An m-1 Branch Peptide of the M-Superfamily

The M-superfamily, one of eight major conotoxin superfamilies found in the venom of the cone snail, contains a Cys framework with disulfide-linked loops labeled 1, 2, and 3 (-CC (1) C (2) C (3) CC-). M-Superfamily conotoxins can be divided into the m-1, -2, -3, and -4 branches, based upon the number of residues located in the third Cys loop between the fourth and fifth Cys residues. Here we provide a three-dimensional solution structure for the m-1 conotoxin tx3a found in the venom of Conus textile. The 15-amino acid peptide, CCSWDVCDHPSCTCC, has disulfide bonds between Cys (1) and Cys (14), Cys (2) and Cys (12), and Cys (7) and Cys (15) typical of the C1-C5, C2-C4, and C3-C6 connectivity pattern seen in m-1 branch peptides. The tertiary structure of tx3a was determined by two-dimensional (1)H NMR in combination with the combined assignment and dynamics algorithm for nuclear magnetic resonance (NMR) applications CYANA program. Input for structure calculations consisted of 62 inter- and intraproton, five phi angle, and four hydrogen bond constraints. The root-mean-square deviation values for the 20 final structures are 0.32 +/- 0.07 and 0.84 +/- 0.11 A for the backbone and heavy atoms, respectively. Surprisingly, the structure of tx3a has a "triple-turn" motif seen in the m-2 branch conotoxin mr3a, which is absent in mr3e, the only other member of the m-1 branch of the M-superfamily whose structure is known. Interestingly, injection of tx3a into mice elicits an excitatory response similar to that of the m-2 branch peptide mr3a, even though the conotoxins have different disulfide connectivity patterns.

Native atomic burials, supplemented by physically motivated hydrogen bond constraints, contain sufficient information to determine the tertiary structure of small globular proteins

We investigate the possibility that atomic burials, as measured by their distances from the structural geometrical center, contain sufficient information to determine the tertiary structure of globular proteins. We report Monte Carlo simulated annealing results of all-atom hard-sphere models in continuous space for four small proteins: the all-beta WW-domain 1E0L, the alpha/beta protein-G 1IGD, the all-alpha engrailed homeo-domain 1ENH, and the alpha + beta engineered monomeric form of the Cro protein 1ORC. We used as energy function the sum over all atoms, labeled by i, of |R(i) - R(i) (*)|, where R(i) is the atomic distance from the center of coordinates, or central distance, and R(i) (*) is the "ideal" central distance obtained from the native structure. Hydrogen bonds were taken into consideration by the assignment of two ideal distances for backbone atoms forming hydrogen bonds in the native structure depending on the formation of a geometrically defined bond, independently of bond partner. Lowest energy final conformations turned out to be very similar to the native structure for the four proteins under investigation and a strong correlation was observed between energy and distance root mean square deviation (DRMS) from the native in the case of all-beta 1E0L and alpha/beta 1IGD. For all alpha 1ENH and alpha + beta 1ORC the overall correlation between energy and DRMS among final conformations was not as high because some trajectories resulted in high DRMS but low energy final conformations in which alpha-helices adopted a non-native mutual orientation. Comparison between central distances and actual accessible surface areas corroborated the implicit assumption of correlation between these two quantities. The Z-score obtained with this native-centric potential in the discrimination of native 1ORC from a set of random compact structures confirmed that it contains a much smaller amount of native information when compared to a traditional contact Go potential but indicated that simple sequence-dependent burial potentials still need some improvement in order to attain a similar discriminability. Taken together, our results suggest that central distances, in conjunction to physically motivated hydrogen bond constraints, contain sufficient information to determine the native conformation of these small proteins and that a solution to the folding problem for globular proteins could arise from sufficiently accurate burial predictions from sequence followed by minimization of a burial-dependent energy function.

The Existence of Multiple Conformers of Interleukin-21 Directs Engineering of a Superpotent Analogue

The high resolution three-dimensional structure of human interleukin (hIL)-21 has been resolved by heteronuclear NMR spectroscopy. Overall, the hIL-21 structure is dominated by a well defined central four-helical bundle, arranged in an up-up-down-down topology, as observed for other cytokines. A segment of the hIL-21 molecule that includes the third helical segment, helix C, is observed to exist in two distinct and interchangeable states. In one conformer, the helix C segment is presented in a regular, alpha-helical conformation, whereas in the other conformer, this segment is largely disordered. A structure-based sequence alignment of hIL-21 with receptor complexes of the related cytokines, interleukin-2 and -4, implied that this particular segment is involved in receptor binding. An hIL-21 analog was designed to stabilize the region around helix C through the introduction of a segment grafted from hIL-4. This novel hIL-21 analog was demonstrated to exhibit a 10-fold increase in potency in a cellular assay.

Physical‐chemical determinants of turn conformations in globular proteins

Globular proteins are assemblies of alpha-helices and beta-strands, interconnected by reverse turns and longer loops. Most short turns can be classified readily into a limited repertoire of discrete backbone conformations, but the physical-chemical determinants of these distinct conformational basins remain an open question. We investigated this question by exhaustive analysis of all backbone conformations accessible to short chain segments bracketed by either an alpha-helix or a beta-strand (i.e., alpha-segment-alpha, beta-segment-beta, alpha-segment-beta, and beta-segment-alpha) in a nine-state model. We find that each of these four secondary structure environments imposes its own unique steric and hydrogen-bonding constraints on the intervening segment, resulting in a limited repertoire of conformations. In greater detail, an exhaustive set of conformations was generated for short backbone segments having reverse-turn chain topology and bracketed between elements of secondary structure. This set was filtered, and only clash-free, hydrogen-bond-satisfied conformers having reverse-turn topology were retained. The filtered set includes authentic turn conformations, observed in proteins of known structure, but little else. In particular, over 99% of the alternative conformations failed to satisfy at least one criterion and were excluded from the filtered set. Furthermore, almost all of the remaining alternative conformations have close tolerances that would be too tight to accommodate side chains longer than a single beta-carbon. These results provide a molecular explanation for the observation that reverse turns between elements of regular secondary can be classified into a small number of discrete conformations.

Insights to Substrate Binding and Processing by West Nile Virus NS3 Protease through Combined Modeling, Protease Mutagenesis, and Kinetic Studies

West Nile Virus is becoming a widespread pathogen, infecting people on at least four continents with no effective treatment for these infections or many of their associated pathologies. A key enzyme that is essential for viral replication is the viral protease NS2B-NS3, which is highly conserved among all flaviviruses. Using a combination of molecular fitting of substrates to the active site of the crystal structure of NS3, site-directed enzyme and cofactor mutagenesis, and kinetic studies on proteolytic processing of panels of short peptide substrates, we have identified important enzyme-substrate interactions that define substrate specificity for NS3 protease. In addition to better understanding the involvement of S2, S3, and S4 enzyme residues in substrate binding, a residue within cofactor NS2B has been found to strongly influence the preference of flavivirus proteases for lysine or arginine at P2 in substrates. Optimization of tetrapeptide substrates for enhanced protease affinity and processing efficiency has also provided important clues for developing inhibitors of West Nile Virus infection.

Solution structure of candoxin, a novel three-finger toxin from the venom of Bungarus candidus

Candoxin, a novel toxin purified from the venom of the Malayan krait (Bungarus candidus), is a 66- residue polypeptide containing five disulfide bridges, and is a reversible antagonist of postjunctional nicotinic acetylcholine receptors of the neuromuscular junction. A family of structures were calculated using a combination of distance geometry and simulated annealing with nOe, hydrogen bond and dihedral angle constraints. After refinement 19 structures, which satisfy the experimental constraints, were obtained. A comparison of each of the final structures with the average structure, gives an RMSD of 1.20 A for backbone atoms. Candoxin has an overall conformation similar to other three-finger toxins with a two-stranded and a three-stranded β-sheets. In spite of the low sequence homology, the tertiary fold of candoxin closely resembles those of erabutoxin b and cobrotoxin, with the exception of the disulfide bridge in the Loop I. Though candoxin lags the most conserved Tyr at the origin of the central loop it posses the most popular 3-dimensional array of residues (W31, R35 and K49), which is reported to be critical for curaremimetic neurotoxicity. The presence of the kink at the tip of the first loop (Loop I) caused by the C6-C11 disulphide bridge, shortens the Loop and hence may isolate this portion of the molecule and prevent its interference with the rest of the molecule, and facilitate specific interaction with receptor/acceptor protein. These observations suggest that candoxin may bind to acetylcholine receptor with lower affinity compared to other neurotoxins and hence could account for its relatively lower toxicity.

Effect of chemical structure on the volume-phase transition in neutral and weakly charged poly(N-alkyl(meth)acrylamide) hydrogels studied by ultrasmall-angle x-ray scattering

Neutral poly(N-isopropylacrylamide) (PIPAAm), poly(N,N-diethylacrylamide) (PDEAAm), and poly(N-isopropylmethacrylamide) (PIPMAm) hydrogels and their weakly charged counterparts prepared by copolymerizing with sodium methacrylate (x(MNa)=0,0.025,0.05) were studied using ultrasmall-angle x-ray scattering. The volume-phase transition in hydrogels was observed as an increase in the inhomogeneity correlation length of the networks. The change in inhomogeneity correlation length was abrupt in neutral PIPAAm and PIPMAm gels with increase in temperature but was continuous in neutral PDEAAm gels. Addition of ionic comonomer to the network backbone suppressed the volume-phase transition in poly(N-alkylacrylamide)s but not in PIPMAm. The observed differences in temperature-induced volume change of these three polymers in water cannot be rationalized based on their relative hydrophobicity and are instead explained by considering the hydrogen-bonding constraints on their thermal fluctuations. Both PIPAAm and PDEAAm undergo volume collapse since their thermal fluctuations are constrained by hydrogen bonding with water to an extent that beyond a critical temperature they seek entropic compensation. Although thermal fluctuations in both PIPAAm and PIPMAm are equally constrained, thermal energy of the latter can be relaxed via the rotation of alpha-methyl groups allowing it greater flexibility. Compared to N-alkylacrylamides, N-alkylmethacrylamide can thus sustain hydrogen bonding to relatively higher temperatures before seeking entropic compensation by undergoing volume collapse.

3D structure of Alzheimer's amyloid-β(1–42) fibrils

Alzheimer's disease is the most fatal neurodegenerative disorder wherein the process of amyloid-beta (Abeta) amyloidogenesis appears causative. Here, we present the 3D structure of the fibrils comprising Abeta(1-42), which was obtained by using hydrogen-bonding constraints from quenched hydrogen/deuterium-exchange NMR, side-chain packing constraints from pairwise mutagenesis studies, and parallel, in-register beta-sheet arrangement from previous solid-state NMR studies. Although residues 1-17 are disordered, residues 18-42 form a beta-strand-turn-beta-strand motif that contains two intermolecular, parallel, in-register beta-sheets that are formed by residues 18-26 (beta1) and 31-42 (beta2). At least two molecules of Abeta(1-42) are required to achieve the repeating structure of a protofilament. Intermolecular side-chain contacts are formed between the odd-numbered residues of strand beta1 of the nth molecule and the even-numbered residues of strand beta2 of the (n - 1)th molecule. This interaction pattern leads to partially unpaired beta-strands at the fibrillar ends, which explains the sequence selectivity, the cooperativity, and the apparent unidirectionality of Abeta fibril growth. It also provides a structural basis for fibrillization inhibitors.

Constrained geometric simulation of diffusive motion in proteins

We describe a new computational method, FRODA (framework rigidity optimized dynamic algorithm), for exploring the internal mobility of proteins. The rigid regions in the protein are first determined, and then replaced by ghost templates which are used to guide the movements of the atoms in the protein. Using random moves, the available conformational phase space of a 100 residue protein can be well explored in approximately 10–100 min of computer time using a single processor. All of the covalent, hydrophobic and hydrogen bond constraints are maintained, and van der Waals overlaps are avoided, throughout the simulation. We illustrate the results of a FRODA simulation on barnase, and show that good agreement is obtained with nuclear magnetic resonance experiments. We additionally show how FRODA can be used to find a pathway from one conformation to another. This directed dynamics is illustrated with the protein dihydrofolate reductase.

The Three-dimensional Solution Structure of Ca 2+-bound S100A1 as Determined by NMR Spectroscopy

S100A1 is an EF-hand-containing Ca(2+)-binding protein that undergoes a conformational change upon binding calcium as is necessary to interact with protein targets and initiate a biological response. To better understand how calcium influences the structure and function of S100A1, the three-dimensional structure of calcium-bound S100A1 was determined by multidimensional NMR spectroscopy and compared to the previously determined structure of apo. In total, 3354 nuclear Overhauser effect-derived distance constraints, 240 dihedral constraints, 160 hydrogen bond constraints, and 362 residual dipolar coupling restraints derived from a series of two-dimensional, three-dimensional, and four-dimensional NMR experiments were used in its structure determination (>21 constraints per residue). As with other dimeric S100 proteins, S100A1 is a symmetric homodimer with helices 1, 1', 4, and 4' associating into an X-type four-helix bundle at the dimer interface. Within each subunit there are four alpha-helices and a short antiparallel beta-sheet typical of two helix-loop-helix EF-hand calcium-binding domains. The addition of calcium did not change the interhelical angle of helices 1 and 2 in the pseudo EF-hand significantly; however, there was a large reorientation of helix 3 in the typical EF-hand. The large conformational change exposes a hydrophobic cleft, defined by residues in the hinge region, the C terminus, and regions of helix 3, which are important for the interaction between S100A1 and a peptide (TRTK-12) derived from the actin-capping protein CapZ.

Conformational analysis of human calcitonin in solution

The solution conformation of human calcitonin in a mixture of 60% water and 40% trifluoroethanol has been determined by the combined use of 1H NMR spectroscopy and distance geometry calculations with a distributed computing technique. 1H NMR spectroscopy provided 195 distance constraints and 13 hydrogen bond constraints. The 20 best converged structures exhibit atomic rmsd of 0.43 A for the backbone atoms from the averaged coordinate position in the region of Asn3-Phe22. The conformation is characterized by a nearly amphiphilic alpha-helix domain that extends from Leu4 in the cyclic region to His20. There are no significant differences observed among the overall structures of a series of calcitonins obtained from ultimobranchial bodies, including those that possess 20- to 50-fold greater activity. Three aromatic amino acid residues, Tyr12, Phe16 and Phe19, form a hydrophobic surface of human calcitonin. Bulky side chains on the surface could interfere with the ligand-receptor interaction thereby causing its low activity, relative to those of other species.

The Nociceptin Pharmacophore Site for Opioid Receptor Binding Derived from the NMR Structure and Bioactivity Relationships

Nociceptin, a 17 amino acid opioid-like peptide that has an inhibitory effect on synaptic transmission in the nervous system, is involved in learning, memory, attention, and emotion and is also implicated in the perception of pain and visual, auditory, and olfactory functions. In this study, we investigated the NMR solution structure of nociceptin in membrane-like environments (trifluoroethanol and SDS micelles) and found it to have a relatively stable helix conformation from residues 4-17 with functionally important N-terminal residues being folded aperidoically on top of the helix. In functional assays for receptor binding and calcium flux, alanine-scanning variants of nociceptin indicated that functionally important residues generally followed helix periodicity, consistent with the NMR structural model. Structure-activity relationships allowed identification of pharmacophore sites that were used in small molecule data base searches, affording hits with demonstrated nociceptin receptor binding affinities.

Solution NMR Structure Investigation for Releasing Mechanism of Neocarzinostatin Chromophore from the Holoprotein

Holo-neocarzinostatin (holo-NCS) is a complex protein carrying the anti-tumor active enediyne ring chromophore by a scaffold consisting of an immunoglobulin-like seven-stranded anti-parallel beta-barrel. Because of the labile chromophore reflecting its extremely strong DNA cleavage activity and complete stabilization in the complex, holo-NCS has attracted much attention in clinical use as well as for drug delivery systems. Despite many structural analyses for holo-NCS, the chromophore-releasing mechanism to trigger prompt attacks on the target DNA is still unclear. We determined the three-dimensional structure of the protein and the internal motion by multinuclear NMR to investigate the releasing mechanism. The internal motion studied by 13C NMR methine relaxation experiments showed that the complex has a rigid structure for its loops as well as the beta-barrel in aqueous solution. This agrees with the refined NMR solution structure, which has good convergence in the loop regions. We also showed that the chromophore displayed a similar internal motion as the protein moiety. The structural comparison between the refined solution structure and x-ray crystal structure indicated characteristic differences. Based on the findings, we proposed the chromophore-releasing mechanism by a three-state equilibrium, which sufficiently describes both the strong binding and the prompt releasing of the chromophore. We demonstrated that we could bridge the dynamic properties and the static structure features with simple kinetic assumptions to solve the biochemical function.

Solution Structure of Prosurvival Mcl-1 and Characterization of Its Binding by Proapoptotic BH3-only Ligands

The B cell lymphoma-2 (Bcl-2) homologs myeloid cell leukemia-1 (Mcl-1) and A1 are prosurvival factors that selectively bind a subset of proapoptotic Bcl homology (BH) 3-only proteins. To investigate the molecular basis of the selectivity, we determined the solution structure of the C-terminal Bcl-2-like domain of Mcl-1. This domain shares features expected of a prosurvival Bcl-2 protein, having a helical fold centered on a core hydrophobic helix and a surface-exposed hydrophobic groove for binding its cognate partners. A number of residues in the binding groove differentiate Mcl-1 from its homologs, and in contrast to other Bcl-2 homologs, Mcl-1 has a binding groove in a conformation intermediate between the open structures characterized by peptide complexes and the closed state observed in unliganded structures. Mutagenesis of potential binding site residues was used to probe the contributions of groove residues to the binding properties of Mcl-1. Although mutations in Mcl-1 had little impact on binding, a single mutation in the BH3-only ligand Bad enabled it to bind both Mcl-1 and A1 while retaining its binding to Bcl-2, Bcl-xL, and Bcl-w. Elucidating the selective action of certain BH3-only ligands is required for delineating their mode of action and will aid the search for effective BH3-mimetic drugs.

Solution Structure of Plant Nonspecific Lipid Transfer Protein-2 from Rice (Oryza sativa)

The three-dimensional structure of rice nonspecific lipid transfer protein (nsLTP2) has been solved for the first time. The structure of nsLTP2 was obtained using 813 distance constraints, 30 hydrogen bond constraints, and 19 dihedral angle constraints. Fifteen of the 50 random simulated annealing structures satisfied all of the constraints and possessed good nonbonded contacts. The novel three-dimensional fold of rice nsLTP2 contains a triangular hydrophobic cavity formed by three prominent helices. The four disulfide bonds required for stabilization of the nsLTP2 structure show a different pattern of cysteine pairing compared with nsLTP1. The C terminus of the protein is very flexible and forms a cap over the hydrophobic cavity. Molecular modeling studies suggested that the hydrophobic cavity could accommodate large molecules with rigid structures, such as sterols. The positively charged residues on the molecular surface of nsLTP2 are structurally similar to other plant defense proteins.

Three-dimensional solution structure of the calcium-signaling protein apo-S100A1 as determined by NMR.

S100A1, a member of the S100 protein family, is an EF-hand containing Ca(2+)-binding protein (93 residues per subunit) with noncovalent interactions at its dimer interface. Each subunit of S100A1 has four alpha-helices and a small antiparallel beta-sheet consistent with two helix-loop-helix calcium-binding domains [Baldiserri et al. (1999) J. Biomol. NMR 14, 87-88]. In this study, the three-dimensional structure of reduced apo-S100A1 was determined by NMR spectroscopy using a total of 2220 NOE distance constraints, 258 dihedral angle constraints, and 168 backbone hydrogen bond constraints derived from a series of 2D, 3D, and 4D NMR experiments. The final structure was found to be globular and compact with the four helices in each subunit aligning to form a unicornate-type four-helix bundle. Intermolecular NOE correlations were observed between residues in helices 1 and 4 from one subunit to residues in helices 1' and 4' of the other subunit, respectively, consistent with the antiparallel alignment of the two subunits to form a symmetric X-type four-helix bundle as found for other members of the S100 protein family. Because of the similarity of the S100A1 dimer interface to that found for S100B, it was possible to calculate a model of the S100A1/B heterodimer. This model is consistent with a number of NMR chemical shift changes observed when S100A1 is titrated into a sample of (15)N-labeled S100B. Helix 3 (and 3') of S100A1 was found to have an interhelical angle of -150 degrees with helix 4 (and 4') in the apo state. This crossing angle is quite different (>50 degrees ) from that typically found in other EF-hand containing proteins such as apocalmodulin and apotroponin C but more similar to apo-S100B, which has an interhelical angle of -166 degrees. As with S100B, it is likely that the second EF-hand of apo-S100A1 reorients dramatically upon the addition of Ca(2+), which can explain the Ca(2+) dependence that S100A1 has for binding several of its biological targets.

Sequestration of the Active Site by Interdomain Shifting: CRYSTALLOGRAPHIC AND SPECTROSCOPIC EVIDENCE FOR DISTINCT CONFORMATIONS OF l-3-HYDROXYACYL-CoA DEHYDROGENASE

Sequestration of the Active Site by Interdomain Shifting: CRYSTALLOGRAPHIC AND SPECTROSCOPIC EVIDENCE FOR DISTINCT CONFORMATIONS OF l-3-HYDROXYACYL-CoA DEHYDROGENASE

Establishing a degree of order: obtaining high-resolution NMR structures from molecular alignment

I would like to thank Eva de Alba for providing structural data as well as Ad Bax and James Baber for their critical comments.

Solution structure of an artificial Fe8S8 ferredoxin: the D13C variant of Bacillus schlegelii Fe7S8 ferredoxin.

The solution structure of the D13C variant of the thermostable Fe7S8 ferredoxin from Bacillus schlegelii has been determined by 1H-NMR spectroscopy in its oxidized form. In a variable-temperature NMR study the D13C variant was as thermostable (up to 90 degrees C) as the wild-type protein (WT). Seventy-five out of 77 amino acid residues and 81% of all theoretically expected proton resonances in the D13C Fe8S8 protein have been assigned. Its structure was determined through torsion angle dynamics calculations with the program DYANA, using 935 meaningful NOEs (from a total of 1251), hydrogen bond constraints, and NMR-derived dihedral angle constraints for the cluster-ligating cysteines. Afterwards, restrained energy minimization and restrained molecular dynamics were applied to each conformer of the family. The final family of 20 structures has RMSD values from the mean structure of 0.055 nm for the backbone atoms and of 0.099 nm for all heavy atoms. The overall folding of the WT is maintained in the mutant, except for the immediate vicinity of the new cysteine, which becomes much more similar to native Fe8S8 proteins. The two residues at positions 11 and 12, which constitute an insertion with respect to all known Fe8S8 proteins, assume a conformation that does not prevent the preceding and following residues from folding like in native Fe8S8 proteins. Clear evidence for the existence of two conformations involving almost half of the amino acid residues was found. The two conformations are structurally indistinguishable. Temperature-dependent NMR experiments show that one of them is thermodynamically more stable than the other.