Computational challenges for macromolecular structure determination by X-ray crystallography and solution NMR-spectroscopy

Quality and bias of protein disorder predictors

Disorder in proteins is vital for biological function, yet it is challenging to characterize. Therefore, methods for predicting protein disorder from sequence are fundamental. Currently, predictors are trained and evaluated using data from X-ray structures or from various biochemical or spectroscopic data. However, the prediction accuracy of disordered predictors is not calibrated, nor is it established whether predictors are intrinsically biased towards one of the extremes of the order-disorder axis. We therefore generated and validated a comprehensive experimental benchmarking set of site-specific and continuous disorder, using deposited NMR chemical shift data. This novel experimental data collection is fully appropriate and represents the full spectrum of disorder. We subsequently analyzed the performance of 26 widely-used disorder prediction methods and found that these vary noticeably. At the same time, a distinct bias for over-predicting order was identified for some algorithms. Our analysis has important implications for the validity and the interpretation of protein disorder, as utilized, for example, in assessing the content of disorder in proteomes.

Common TFIIH recruitment mechanism in global genome and transcription-coupled repair subpathways

Nucleotide excision repair is initiated by two different damage recognition subpathways, global genome repair (GGR) and transcription-coupled repair (TCR). In GGR, XPC detects DNA lesions and recruits TFIIH via interaction with the pleckstrin homology (PH) domain of TFIIH subunit p62. In TCR, an elongating form of RNA Polymerase II detects a lesion on the transcribed strand and recruits TFIIH by an unknown mechanism. Here, we found that the TCR initiation factor UVSSA forms a stable complex with the PH domain of p62 via a short acidic string in the central region of UVSSA, and determined the complex structure by NMR. The acidic string of UVSSA binds strongly to the basic groove of the PH domain by inserting Phe408 and Val411 into two pockets, highly resembling the interaction mechanism of XPC with p62. Mutational binding analysis validated the structure and identified residues crucial for binding. TCR activity was markedly diminished in UVSSA-deficient cells expressing UVSSA mutated at Phe408 or Val411. Thus, a common TFIIH recruitment mechanism is shared by UVSSA in TCR and XPC in GGR.

Bayesian Modeling of Biomolecular Assemblies with Cryo-EM Maps

A growing array of experimental techniques allows us to characterize the three-dimensional structure of large biological assemblies at increasingly higher resolution. In addition to X-ray crystallography and nuclear magnetic resonance in solution, new structure determination methods such cryo-electron microscopy (cryo-EM), crosslinking/mass spectrometry and solid-state NMR have emerged. Often it is not sufficient to use a single experimental method, but complementary data need to be collected by using multiple techniques. The integration of all datasets can only be achieved by computational means. This article describes Inferential structure determination, a Bayesian approach to integrative modeling of biomolecular complexes with hybrid structural data. I will introduce probabilistic models for cryo-EM maps and outline Markov chain Monte Carlo algorithms for sampling model structures from the posterior distribution. I will focus on rigid and flexible modeling with cryo-EM data and discuss some of the computational challenges of Bayesian inference in the context of biomolecular modeling.

PSCD Domains of Pleuralin-1 from the Diatom Cylindrotheca fusiformis: NMR Structures and Interactions with Other Biosilica-Associated Proteins

Diatoms are eukaryotic unicellular algae characterized by silica cell walls and associated with three unique protein families, the pleuralins, frustulins, and silaffins. The NMR structure of the PSCD4 domain of pleuralin-1 from Cylindrotheca fusiformis contains only three short helical elements and is stabilized by five unique disulfide bridges. PSCD4 contains two binding sites for Ca(2+) ions with millimolar affinity. NMR-based interaction studies show an interaction of the domain with native silaffin-1A as well as with α-frustulins. The interaction sites of the two proteins mapped on the PSCD4 structure are contiguous and show only a small overlap. A plausible functional role of pleuralin could be to bind simultaneously silaffin-1A located inside the cell wall and α-frustulin coating the cell wall, thus connecting the interfaces between hypotheca and epitheca at the girdle bands. Restrained molecular dynamics calculations suggest a bead-chain-like structure of the central part of pleuralin-1.

Efficient and Minimal Method to Bias Molecular Simulations with Experimental Data

A primary goal in molecular simulations is to modify the potential energy of a system so that properties of the simulation match experimental data. This is traditionally done through iterative cycles of simulation and reparameterization. An alternative approach is to bias the potential energy so that the system matches experimental data. This can be done while minimally changing the underlying free energy of the molecular simulation. Current minimal biasing methods require replicas, which can lead to unphysical dynamics and introduces new complexity: the choice of replica number and their properties. Here, we describe a new method, called experiment directed simulation that does not require replicas, converges rapidly, can match many data simultaneously, and minimally modifies the potential. The experiment directed simulation method is demonstrated on model systems and a three-component electrolyte simulation. The theory used to derive the method also provides insight into how changing a molecular force-field impacts the expected value of observables in simulation.

Bayesian weighting of statistical potentials in NMR structure calculation

The use of statistical potentials in NMR structure calculation improves the accuracy of the final structure but also raises issues of double counting and possible bias. Because statistical potentials are averaged over a large set of structures, they may not reflect the preferences of a particular structure or data set. We propose a Bayesian method to incorporate a knowledge-based backbone dihedral angle potential into an NMR structure calculation. To avoid bias exerted through the backbone potential, we adjust its weight by inferring it from the experimental data. We demonstrate that an optimally weighted potential leads to an improvement in the accuracy and quality of the final structure, especially with sparse and noisy data. Our findings suggest that no universally optimal weight exists, and that the weight should be determined based on the experimental data. Other knowledge-based potentials can be incorporated using the same approach.

Accelerate Sampling in Atomistic Energy Landscapes Using Topology-Based Coarse-Grained Models

We describe a multiscale enhanced sampling (MSES) method where efficient topology-based coarse-grained models are coupled with all-atom ones to enhance the sampling of atomistic protein energy landscape. The bias from the coupling is removed by Hamiltonian replica exchange, thus allowing one to benefit simultaneously from faster transitions of coarse-grained modeling and accuracy of atomistic force fields. The method is demonstrated by calculating the conformational equilibria of several small but nontrivial β-hairpins with varied stabilities.

Inference of structure ensembles of flexible biomolecules from sparse, averaged data

We present the theoretical foundations of a general principle to infer structure ensembles of flexible biomolecules from spatially and temporally averaged data obtained in biophysical experiments. The central idea is to compute the Kullback-Leibler optimal modification of a given prior distribution with respect to the experimental data and its uncertainty. This principle generalizes the successful inferential structure determination method and recently proposed maximum entropy methods. Tractability of the protocol is demonstrated through the analysis of simulated nuclear magnetic resonance spectroscopy data of a small peptide.

Calibration of Boltzmann distribution priors in Bayesian data analysis

The Boltzmann distribution is commonly used as a prior probability in Bayesian data analysis. Examples include the Ising model in statistical image analysis and the canonical ensemble based on molecular dynamics force fields in protein structure calculation. These models involve a temperature or weighting factor that needs to be inferred from the data. Bayesian inference stipulates to determine the temperature based on the model evidence. This is challenging because the model evidence, a ratio of two high-dimensional normalization integrals, cannot be calculated analytically. We outline a replica-exchange Monte Carlo scheme that allows us to estimate the model evidence by use of multiple histogram reweighting. The method is illustrated for an Ising model and examples in protein structure determination.

Applications and limitations of in silico models in drug discovery

Drug discovery in the late twentieth and early twenty-first century has witnessed a myriad of changes that were adopted to predict whether a compound is likely to be successful, or conversely enable identification of molecules with liabilities as early as possible. These changes include integration of in silico strategies for lead design and optimization that perform complementary roles to that of the traditional in vitro and in vivo approaches. The in silico models are facilitated by the availability of large datasets associated with high-throughput screening, bioinformatics algorithms to mine and annotate the data from a target perspective, and chemoinformatics methods to integrate chemistry methods into lead design process. This chapter highlights the applications of some of these methods and their limitations. We hope this serves as an introduction to in silico drug discovery.

The C-terminal tail of human neuronal calcium sensor 1 regulates the conformational stability of the Ca²⁺₋ activated state

Neuronal calcium sensor 1 (NCS-1) and orthologs are expressed in all organisms from yeast to humans. In the latter, NCS-1 plays an important role in neurotransmitter release and interacts with a plethora of binding partners mostly through a large solvent-exposed hydrophobic crevice. The structural basis behind the multispecific binding profile is not understood. To begin to address this, we applied NMR spectroscopy to determine the solution structure of calcium-bound human NCS-1. The structure in solution demonstrates interdomain flexibility and, in the absence of a binding partner, the C-terminal tail residues occupy the hydrophobic crevice as a ligand mimic. A variant with a C-terminal tail deletion shows lack of a defined structure but maintained cooperative unfolding and dramatically reduced global stability. The results suggest that the C-terminal tail is important for regulating the conformational stability of the Ca(2+)-activated state. Furthermore, a single amino acid mutation that was recently diagnosed in a patient with autistic spectrum disorder was seen to affect the C-terminal tail and binding crevice in NCS-1.

Statistical mechanics analysis of sparse data

Inferential structure determination uses Bayesian theory to combine experimental data with prior structural knowledge into a posterior probability distribution over protein conformational space. The posterior distribution encodes everything one can say objectively about the native structure in the light of the available data and additional prior assumptions and can be searched for structural representatives. Here an analogy is drawn between the posterior distribution and the canonical ensemble of statistical physics. A statistical mechanics analysis assesses the complexity of a structure calculation globally in terms of ensemble properties. Analogs of the free energy and density of states are introduced; partition functions evaluate the consistency of prior assumptions with data. Critical behavior is observed with dwindling restraint density, which impairs structure determination with too sparse data. However, prior distributions with improved realism ameliorate the situation by lowering the critical number of observations. An in-depth analysis of various experimentally accessible structural parameters and force field terms will facilitate a statistical approach to protein structure determination with sparse data that avoids bias as much as possible.

Three-dimensional molecular modeling with single molecule FRET

Single molecule fluorescence energy transfer experiments enable investigations of macromolecular conformation and folding by the introduction of fluorescent dyes at specific sites in the macromolecule. Multiple such experiments can be performed with different labeling site combinations in order to map complex conformational changes or interactions between multiple molecules. Distances that are derived from such experiments can be used for determination of the fluorophore positions by triangulation. When combined with a known structure of the macromolecule(s) to which the fluorophores are attached, a three-dimensional model of the system can be determined. However, care has to be taken to properly derive distance from fluorescence energy transfer efficiency and to recognize the systematic or random errors for this relationship. Here we review the experimental and computational methods used for three-dimensional modeling based on single molecule fluorescence resonance transfer, and describe recent progress in pushing the limits of this approach to macromolecular complexes.

Structure and disorder in an unfolded state under nondenaturing conditions from ensemble models consistent with a large number of experimental restraints

Obtaining detailed structural models of disordered states of proteins under nondenaturing conditions is important for a better understanding of both functional intrinsically disordered proteins and unfolded states of folded proteins. Extensive experimental characterization of the drk N-terminal SH3 domain unfolded state has shown that, although it appears to be highly disordered, it possesses significant nonrandom secondary and tertiary structure. In our previous attempts to generate structural models of the unfolded state using the program ENSEMBLE, we were limited by insufficient experimental restraints and conformational sampling. In this study, we have vastly expanded our experimental restraint set to include (1)H-(15)N residual dipolar couplings, small-angle X-ray scattering measurements, nitroxide paramagnetic relaxation enhancements, O(2)-induced (13)C paramagnetic shifts, hydrogen-exchange protection factors, and (15)N R(2) data, in addition to the previously used nuclear Overhauser effects, amino terminal Cu(2+)-Ni(2+) binding paramagnetic relaxation enhancements, J-couplings, chemical shifts, hydrodynamic radius, and solvent accessibility restraints. We have also implemented a new ensemble calculation methodology that uses iterative conformational sampling and seeks to calculate the simplest possible ensemble models. As a result, we can now generate ensembles that are consistent with much larger experimental data sets than was previously possible. Although highly heterogeneous and having broad molecular size distributions, the calculated drk N-terminal SH3 domain unfolded-state ensembles have very different properties than expected for random or statistical coils and possess significant nonnative alpha-helical structure and both native-like and nonnative tertiary structure.

Improving the species cross-reactivity of an antibody using computational design

The high degree of specificity displayed by antibodies often results in varying potencies against antigen orthologs, which can affect the efficacy of these molecules in different animal models of disease. We have used a computational design strategy to improve the species cross-reactivity of an antibody-based inhibitor of the cancer-associated serine protease MT-SP1. In silico predictions were tested in vitro, and the most effective mutation, T98R, was shown to improve antibody affinity for the mouse ortholog of the enzyme 14-fold, resulting in an inhibitor with a K(I) of 340 pM. This improved affinity will be valuable when exploring the role of MT-SP1 in mouse models of cancer, and the strategy outlined here could be useful in fine-tuning antibody specificity.

Accurate NMR structures through minimization of an extended hybrid energy

The use of generous distance bounds has been the hallmark of NMR structure determination. However, bounds necessitate the estimation of data quality before the calculation, reduce the information content, introduce human bias, and allow for major errors in the structures. Here, we propose a new rapid structure calculation scheme based on Bayesian analysis. The minimization of an extended energy function, including a new type of distance restraint and a term depending on the data quality, results in an estimation of the data quality in addition to coordinates. This allows for the determination of the optimal weight on the experimental information. The resulting structures are of better quality and closer to the X-ray crystal structure of the same molecule. With the new calculation approach, the analysis of discrepancies from the target distances becomes meaningful. The strategy may be useful in other applications-for example, in homology modeling.

Structural biology by NMR: structure, dynamics, and interactions

The function of bio-macromolecules is determined by both their 3D structure and conformational dynamics. These molecules are inherently flexible systems displaying a broad range of dynamics on time-scales from picoseconds to seconds. Nuclear Magnetic Resonance (NMR) spectroscopy has emerged as the method of choice for studying both protein structure and dynamics in solution. Typically, NMR experiments are sensitive both to structural features and to dynamics, and hence the measured data contain information on both. Despite major progress in both experimental approaches and computational methods, obtaining a consistent view of structure and dynamics from experimental NMR data remains a challenge. Molecular dynamics simulations have emerged as an indispensable tool in the analysis of NMR data.

NMR determination of oligonucleotide structure

This unit provides an overview of the use of NMR to determine oligonucleotide structure. It covers basic NMR spectral properties, acquisition of interproton distance restraints and torsion angle restraints, structure refinement, assessment of the quality of the structure obtained. Software programs used in the process are also described.

Error distribution derived NOE distance restraints

Errors and imprecisions in distance restraints derived from NOESY peak volumes are usually accounted for by generous lower and upper bounds on the distances. In this paper, we propose a new form of distance restraints, replacing the subjective bounds by a potential function obtained from the error distribution of the distances. We derived the shape of the potential from molecular dynamics calculations and by comparison of NMR data with X-ray crystal structures. We used complete cross-validation to derive the optimal weight for the data in the calculation. In a model system with synthetic restraints, the accuracy of the structures improved significantly compared to calculations with the usual form of restraints. For experimental data sets, the structures systematically approach the X-ray crystal structures of the same protein. Also standard quality indicators improve compared to standard calculations. The results did not depend critically on the exact shape of the potential. The new approach is less subjective and uses fewer assumptions in the interpretation of NOESY peak volumes as distance restraints than the usual approach. Figures of merit for the structures, such as the RMS difference from the average structure or the RMS difference from the data, are therefore less biased and more meaningful measures of structure quality than with the usual form of restraints.

Weighting of experimental evidence in macromolecular structure determination

The determination of macromolecular structures requires weighting of experimental evidence relative to prior physical information. Although it can critically affect the quality of the calculated structures, experimental data are routinely weighted on an empirical basis. At present, cross-validation is the most rigorous method to determine the best weight. We describe a general method to adaptively weight experimental data in the course of structure calculation. It is further shown that the necessity to define weights for the data can be completely alleviated. We demonstrate the method on a structure calculation from NMR data and find that the resulting structures are optimal in terms of accuracy and structural quality. Our method is devoid of the bias imposed by an empirical choice of the weight and has some advantages over estimating the weight by cross-validation.

Bayesian inference applied to macromolecular structure determination

The determination of macromolecular structures from experimental data is an ill-posed inverse problem. Nevertheless, conventional techniques to structure determination attempt an inversion of the data by minimization of a target function. This approach leads to problems if the data are sparse, noisy, heterogeneous, or difficult to describe theoretically. We propose here to view biomolecular structure determination as an inference rather than an inversion problem. Probability theory then offers a consistent formalism to solve any structure determination problem: We use Bayes' theorem to derive a probability distribution for the atomic coordinates and all additional unknowns. This distribution represents the complete information contained in the data and can be analyzed numerically by Markov chain Monte Carlo sampling techniques. We apply our method to data obtained from a nuclear magnetic resonance experiment and discuss the estimation of theory parameters.

Protein flexibility: multiple molecular dynamics simulations of insulin chain B

Multiple molecular dynamics simulations totaling more than 100 ns were performed on chain B of insulin in explicit solvent at 300 K and 400 K. Despite some individual variations, a comparison of the protein dynamics of each simulation showed similar trends and most structures were consistent with NMR experimental values, even at the elevated temperature. The importance of packing interactions in determining the conformational transitions of the protein was observed, sometimes resulting in conformations induced by localized hydrophobic interactions. The high temperature simulation generated a more diverse range of structures with similar elements of secondary structure and populated conformations to the simulations at room temperature. A broad sampling of the conformational space of insulin chain B illustrated a wide range of conformational states with many transitions at room temperature in addition to the conformational states observed experimentally. The T-state conformation associated with insulin activity was consistently present and a possible mechanism of behavior was suggested.

Inferential structure determination

Macromolecular structures calculated from nuclear magnetic resonance data are not fully determined by experimental data but depend on subjective choices in data treatment and parameter settings. This makes it difficult to objectively judge the precision of the structures. We used Bayesian inference to derive a probability distribution that represents the unknown structure and its precision. This probability distribution also determines additional unknowns, such as theory parameters, that previously had to be chosen empirically. We implemented this approach by using Markov chain Monte Carlo techniques. Our method provides an objective figure of merit and improves structural quality.

Quantitative study of the effects of chemical shift tolerances and rates of SA cooling on structure calculation from automatically assigned NOE data

The calculation of protein structures from nuclear magnetic resonance (NMR) data has been greatly facilitated by improvements in software for the automatic assignment of NOESY spectra. Nevertheless, for larger proteins, resonance overlap may lead to an overwhelming number of assignment options per peak. Although most software for automatic NOESY assignment can deal with a certain level of assignment ambiguity, structure calculations fail when this becomes too high. Reducing the number of assignment options per peak by reducing the chemical shift tolerances can lead to correct assignments being excluded, and thus also to incorrect structures. We have investigated, systematically, for three proteins of different size, the influence of the chemical shift tolerance limits (Delta) and of the number of simulated annealing (SA) cooling steps on the performance of the software ARIA. Large tolerance windows, and the correspondingly high levels of ambiguity, did not cause problems when appropriately slower cooling was used in our SA protocol. In cases where a high percentage of well-converged structures was not achieved, we demonstrate that it is more productive to calculate fewer structures whilst applying slow cooling, than to calculate many structures with fast cooling. In this way, high-quality structures were obtained even for proteins whose NMR spectra showed great degeneracy, and where there was much inconsistency in peak alignment between different samples. The method described herein opens the way to the automated structure determination of larger proteins from NMR data.

Contribution of single tryptophan residues to the fluorescence and stability of ribonuclease Sa

Ribonuclease Sa (RNase Sa) contains no tryptophan (Trp) residues. We have added single Trp residues to RNase Sa at sites where Trp is found in four other microbial ribonucleases, yielding the following variants of RNase Sa: Y52W, Y55W, T76W, and Y81W. We have determined crystal structures of T76W and Y81W at 1.1 and 1.0 A resolution, respectively. We have studied the fluorescence properties and stabilities of the four variants and compared them to wild-type RNase Sa and the other ribonucleases on which they were based. Our results should help others in selecting sites for adding Trp residues to proteins. The most interesting findings are: 1), Y52W is 2.9 kcal/mol less stable than RNase Sa and the fluorescence intensity emission maximum is blue-shifted to 309 nm. Only a Trp in azurin is blue-shifted to a greater extent (308 nm). This blue shift is considerably greater than observed for Trp71 in barnase, the Trp on which Y52W is based. 2), Y55W is 2.1 kcal/mol less stable than RNase Sa and the tryptophan fluorescence is almost completely quenched. In contrast, Trp59 in RNase T1, on which Y55W is based, has a 10-fold greater fluorescence emission intensity. 3), T76W is 0.7 kcal/mol more stable than RNase Sa, indicating that the Trp side chain has more favorable interactions with the protein than the threonine side chain. The fluorescence properties of folded Y76W are similar to those of the unfolded protein, showing that the tryptophan side chain in the folded protein is largely exposed to solvent. This is confirmed by the crystal structure of the T76W which shows that the side chain of the Trp is only approximately 7% buried. 4), Y81W is 0.4 kcal/mol less stable than RNase Sa. Based on the crystal structure of Y81W, the side chain of the Trp is 87% buried. Although all of the Trp side chains in the variants contribute to the unusual positive circular dichroism band observed near 235 nm for RNase Sa, the contribution is greatest for Y81W.

The solution structure of frenatin 3, a neuronal nitric oxide synthase inhibitor from the giant tree frog,Litoria infrafrenata

The peptide frenatin 3 is a major component of the skin secretion of the Australian giant tree frog, Litoria infrafrenata. Frenatin 3 is 22 amino acids in length, and shows neither antimicrobial nor anticancer activity. It inhibits the production of nitric oxide by the enzyme neuronal nitric oxide synthase at a micromolar concentration by binding to its regulatory protein, Ca2+ calmodulin, a protein known to recognize and bind amphipathic alpha-helices. The solution structure of frenatin 3 has been investigated using NMR spectroscopy and restrained molecular dynamics calculations. In trifluoroethanol/water mixtures, the peptide forms an amphipathic alpha-helix over residues 1-14 while the C-terminal eight residues are more flexible and less structured. The flexible region may be responsible for the lack of antimicrobial activity. In water, frenatin 3 exhibits some alpha-helical character in its N-terminal region.

The 0.93Å Crystal Structure of Sphericase: A Calcium-loaded Serine Protease from Bacillus sphaericus

We have previously isolated sphericase (Sph), an extracellular mesophilic serine protease produced by Bacillus sphaericus. The Sph amino acid sequence is highly homologous to two cold-adapted subtilisins from Antarctic bacilli S39 and S41 (76% and 74% identity, respectively). Sph is calcium-dependent, 310 amino acid residues long and has optimal activity at pH 10.0. S41 and S39 have not as yet been structurally analysed. In the present work, we determined the crystal structure of Sph by the Eu/multiwavelength anomalous diffraction method. The structure was extended to 0.93A resolution and refined to a crystallographic R-factor of 9.7%. The final model included all 310 amino acid residues, one disulfide bond, 679 water molecules and five calcium ions. Although Sph is a mesophilic subtilisin, its amino acid sequence is similar to that of the psychrophilic subtilisins, which suggests that the crystal structure of these subtilisins is very similar. The presence of five calcium ions bound to a subtilisin molecule, as found here for Sph, has not been reported for the subtilisin superfamily. None of these calcium-binding sites correlates with the well-known high-affinity calcium-binding site (site I or site A), and only one site has been described previously. This calcium-binding pattern suggests that a reduction in the flexibility of the surface loops of Sph by calcium binding may be responsible for its adaptation to mesophilic organisms.

Relative stability of protein structures determined by X-ray crystallography or NMR spectroscopy: a molecular dynamics simulation study

The relative stability of protein structures determined by either X-ray crystallography or nuclear magnetic resonance (NMR) spectroscopy has been investigated by using molecular dynamics simulation techniques. Published structures of 34 proteins containing between 50 and 100 residues have been evaluated. The proteins selected represent a mixture of secondary structure types including all alpha, all beta, and alpha/beta. The proteins selected do not contain cysteine-cysteine bridges. In addition, any crystallographic waters, metal ions, cofactors, or bound ligands were removed before the systems were simulated. The stability of the structures was evaluated by simulating, under identical conditions, each of the proteins for at least 5 ns in explicit solvent. It is found that not only do NMR-derived structures have, on average, higher internal strain than structures determined by X-ray crystallography but that a significant proportion of the structures are unstable and rapidly diverge in simulations.

The solution structures and activity of caerin 1.1 and caerin 1.4 in aqueous trifluoroethanol and dodecylphosphocholine micelles

The caerin 1 peptides are among the most powerful of the broad-spectrum antibiotic amphibian peptides. Caerin 1.1 has previously been shown to form an amphipathic helix-bend-helix structure in aqueous trifluoroethanol (H. Wong, J. H. Bowie, and J. A. Carver European Journal Biochemistry, 1997, Vol. 247, pp. 545-557) and structure-activity relationship studies indicate that both helices are required for activity, as well as flexibility in the bend region connecting the two. The structure of caerin 1.1 in dodecylphosphocholine micelles was investigated and shown to be very similar to that determined in aqueous trifluoroethanol. Caerin 1.4, which is identical to caerin 1.1, but with serine residues replacing Val5 and Gly7, is less active than caerin 1.1 against most bacterial species but has improved activity against Escherichia coli and Micrococcus luteus. The solution NMR structure of caerin 1.4 was determined in both aqueous trifluoroethanol and dodecylphosphocholine micelles, and was shown to be similar to caerin 1.1. It was concluded that differences in the hydrophobicity and hydrophilic angle of the first helix are probably responsible for the different spectra of antibacterial activity. The similarity of the structures calculated in aqueous trifluoroethanol and dodecylphosphocholine micelles suggests that, for caerin 1.1 and 1.4, these solvent systems are equally as good at representing a membrane environment.

Conformational differences in liganded and unliganded states of Galectin-3

The conformation of the carbohydrate recognition domain of Galectin-3, a lectin known to bind galactose containing oligosaccharides in mammalian systems, has been investigated in the absence of ligand and in the presence of N-acetylactosamine. A new methodology based on the measurement of residual dipolar couplings from NMR spectra has been used to characterize differences in protein structure along the backbone in the presence and absence of ligand, as well as the binding geometry of the ligand itself. The data on the ligand are consistent with the ligand binding geometry found in a crystal structure of the complexed state. However, a significant rearrangement of backbone loops near the binding site appears to occur in the absence of ligand. The implications for ligand specificity and protein functionality are discussed.

nNOS inhibition, antimicrobial and anticancer activity of the amphibian skin peptide, citropin 1.1 and synthetic modifications. The solution structure of a modified citropin 1.1

A large number of bioactive peptides have been isolated from amphibian skin secretions. These peptides have a variety of actions including antibiotic and anticancer activities and the inhibition of neuronal nitric oxide synthase. We have investigated the structure-activity relationship of citropin 1.1, a broad-spectrum antibiotic and anticancer agent that also causes inhibition of neuronal nitric oxide synthase, by making a number of synthetically modified analogues. Citropin 1.1 has been shown previously to form an amphipathic alpha-helix in aqueous trifluoroethanol. The results of the structure-activity studies indicate the terminal residues are important for bacterial activity and increasing the overall positive charge, while maintaining an amphipathic distribution of residues, increases activity against Gram-negative organisms. Anticancer activity generally mirrors antibiotic activity suggesting a common mechanism of action. The N-terminal residues are important for inhibition of neuronal nitric oxide synthase, as is an overall positive charge greater than three. The structure of one of the more active synthetic modifications (A4K14-citropin 1.1) was determined in aqueous trifluoroethanol, showing that this peptide also forms an amphipathic alpha-helix.

Molecular modelling in structural biology

Molecular modelling is a powerful methodology for analysing the three dimensional structure of biological macromolecules. There are many ways in which molecular modelling methods have been used to address problems in structural biology. It is not widely appreciated that modelling methods are often an integral component of structure determination by NMR spectroscopy and X-ray crystallography. In this review we consider some of the numerous ways in which modelling can be used to interpret and rationalise experimental data and in constructing hypotheses that can be tested by experiment. Genome sequencing projects are producing a vast wealth of data describing the protein coding regions of the genome under study. However, only a minority of the protein sequences thus identified will have a clear sequence homology to a known protein. In such cases valuable three-dimensional models of the protein coding sequence can be constructed by homology modelling methods. Threading methods, which used specialised schemes to relate protein sequences to a library of known structures, have been shown to be able to identify the likely protein fold even in cases where there is no clear sequence homology. The number of protein sequences that cannot be assigned to a structural class by homology or threading methods, simply because they belong to a previously unidentified protein folding class, will decrease in the future as collaborative efforts in systematic structure determination begin to develop. For this reason, modelling methods are likely to become increasingly useful in the near future. The role of the blind prediction contests, such as the Critical Assessment of techniques for protein Structure Prediction (CASP), will be briefly discussed. Methods for modelling protein-ligand and protein-protein complexes are also described and examples of their applications given.

Conserved segments 1A and 2B of the intermediate filament dimer: their atomic structures and role in filament assembly

Intermediate filaments (IFs) are key components of the cytoskeleton in higher eukaryotic cells. The elementary IF 'building block' is an elongated coiled-coil dimer consisting of four consecutive alpha-helical segments. The segments 1A and 2B include highly conserved sequences and are critically involved in IF assembly. Based on the crystal structures of three human vimentin fragments at 1.4-2.3 A resolution (PDB entries 1gk4, 1gk6 and 1gk7), we have established the molecular organization of these two segments. The fragment corresponding to segment 1A forms a single, amphipatic alpha-helix, which is compatible with a coiled-coil geometry. While this segment might yield a coiled coil within an isolated dimer, monomeric 1A helices are likely to play a role in specific dimer-dimer interactions during IF assembly. The 2B segment reveals a double-stranded coiled coil, which unwinds near residue Phe351 to accommodate a 'stutter'. A fragment containing the last seven heptads of 2B interferes heavily with IF assembly and also transforms mature vimentin filaments into a new kind of structure. These results provide the first insight into the architecture and functioning of IFs at the atomic level.

Tyrosine hydrogen bonds make a large contribution to protein stability

The aim of this study was to gain a better understanding of the contribution of hydrogen bonds by tyrosine -OH groups to protein stability. The amino acid sequences of RNases Sa and Sa3 are 69 % identical and each contains eight Tyr residues with seven at equivalent structural positions. We have measured the stability of the 16 tyrosine to phenylalanine mutants. For two equivalent mutants, the stability increases by 0.3 kcal/mol (RNase Sa Y30F) and 0.5 kcal/mol (RNase Sa3 Y33F) (1 kcal=4.184 kJ). For all of the other mutants, the stability decreases with the greatest decrease being 3.6 kcal/mol for RNase Sa Y52F. Seven of the 16 tyrosine residues form intramolecular hydrogen bonds and the average decrease in stability for these is 2.0(+/-1.0) kcal/mol. For the nine tyrosine residues that do not form intramolecular hydrogen bonds, the average decrease in stability is 0.4(+/-0.6) kcal/mol. Thus, most tyrosine -OH groups contribute favorably to protein stability even if they do not form intramolecular hydrogen bonds. Generally, the stability changes for equivalent positions in the two proteins are remarkably similar. Crystal structures were determined for two of the tyrosine to phenylalanine mutants of RNase Sa: Y80F (1.2 A), and Y86F (1.7 A). The structures are very similar to that of wild-type RNase Sa, and the hydrogen bonding partners of the tyrosine residues always form intermolecular hydrogen bonds to water in the mutants. These results provide further evidence that the hydrogen bonding and van der Waals interactions of polar groups in the tightly packed interior of folded proteins are more favorable than similar interactions with water in the unfolded protein, and that polar group burial makes a substantial contribution to protein stability.

Solution structure and dynamics of ribonuclease Sa

We have used NMR methods to characterize the structure and dynamics of ribonuclease Sa in solution. The solution structure of RNase Sa was obtained using the distance constraints provided by 2,276 NOEs and the C6-C96 disulfide bond. The 40 resulting structures are well determined; their mean pairwise RMSD is 0.76 A (backbone) and 1.26 A (heavy atoms). The solution structures are similar to previously determined crystal structures, especially in the secondary structure, but exhibit new features: the loop composed of Pro 45 to Ser 48 adopts distinct conformations and the rings of tyrosines 51, 52, and 55 have reduced flipping rates. Amide protons with greatly reduced exchange rates are found predominantly in interior beta-strands and the alpha-helix, but also in the external 3/10 helix and edge beta-strand linked by the disulfide bond. Analysis of (15)N relaxation experiments (R1, R2, and NOE) at 600 MHz revealed five segments, consisting of residues 1-5, 28-31, 46-50, 60-65, 74-77, retaining flexibility in solution. The change in conformation entropy for RNase SA folding is smaller than previously believed, since the native protein is more flexible in solution than in a crystal.

Structural analyses of nucleotide binding to an aminoglycoside phosphotransferase

3',5"-Aminoglycoside phosphotransferase type IIIa [APH(3')-IIIa] is a bacterial enzyme that confers resistance to a range of aminoglycoside antibiotics while exhibiting striking homology to eukaryotic protein kinases (ePK). The structures of APH(3')-IIIa in its apoenzyme form and in complex with the nonhydrolyzable ATP analogue AMPPNP were determined to 3.2 and 2.4 A resolution, respectively. Furthermore, refinement of the previously determined ADP complex was completed. The structure of the apoenzyme revealed alternate positioning of a flexible loop (analogous to the P-loop of ePK's), occupying part of the nucleotide-binding pocket of the enzyme. Despite structural similarity to protein kinases, there was no evidence of domain movement associated with nucleotide binding. This rigidity is due to the presence of more extensive interlobe interactions in the APH(3')-IIIa structure than in the ePK's. Differences between the ADP and AMPPNP complexes are confined to the area of the nucleotide-binding pocket. The position of conserved active site residues and magnesium ions remains unchanged, but there are differences in metal coordination between the two nucleotide complexes. Comparison of the di/triphosphate binding site of APH(3')-IIIa with that of ePK's suggests that the reaction mechanism of APH(3")-IIIa and related aminoglycoside kinases will closely resemble that of eukaryotic protein kinases. However, the orientation of the adenine ring in the binding pocket differs between APH(3')-IIIa and the ePK's by a rotation of approximately 40 degrees. This alternate binding mode is likely a conserved feature among aminoglycoside kinases and could be exploited for the structure-based drug design of compounds to combat antibiotic resistance.

Conformation and dynamics of heparin and heparan sulfate

The glycosaminoglycans heparin and heparan sulfate contain similar structural units in varying proportions providing considerable diversity in sequence and biological function. Both compounds are alternating copolymers of glucosamine with both iduronate- and glucuronate-containing sequences bearing N-sulfate, N-acetyl, and O-sulfate substitution. Protein recognition of these structurally-diverse compounds depends upon substitution pattern, overall molecular shape, and on internal mobility. In this review particular attention is paid to the dynamic aspects of heparin/heparan sulfate conformation. The iduronate residue possesses an unusually flexible pyranose ring conformation. This extra source of internal mobility creates special problems in rationalization of experimental data for these compounds. We present herein the solution-state NMR parameters, fiber diffraction data, crystallographic data, and molecular modeling methods employed in the investigation of heparin and heparan sulfate. Heparin is a useful model compound for the sulfated, protein-binding regions of heparan sulfate. The literature contains a number of solution and solid-state studies of heparin oligo- and polysaccharides for both isolated heparin species and those bound to protein receptors. These studies indicate a diversity of iduronate ring conformations, but a limited range of glycosidic linkage geometries in the repeating disaccharides. In this sense, heparin exhibits a well-defined overall shape within which iduronate ring forms can freely interconvert. Recent work suggests that computational modeling could potentially identify heparin binding sites on protein surfaces.

The antibiotic and anticancer active aurein peptides from the Australian Bell Frogs Litoria aurea and Litoria raniformis the solution structure of aurein 1.2

Seventeen aurein peptides are present in the secretion from the granular dorsal glands of the Green and Golden Bell Frog Litoria aurea, and 16 from the corresponding secretion of the related Southern Bell Frog L. raniformis. Ten of these peptides are common to both species. Thirteen of the aurein peptides show wide-spectrum antibiotic and anticancer activity. These peptides are named in three groups (aureins 1-3) according to their sequences. Amongst the more active peptides are aurein 1.2 (GLFDIIKKIAESF-NH2), aurein 2.2 (GLFDIVKKVVGALGSL-NH2) and aurein 3.1 (GLFDIVKKIAGHIAGSI-NH2). Both L. aurea and L. raniformis have endoproteases that deactivate the major membrane-active aurein peptides by removing residues from both the N- and C-termini of the peptides. The most abundant degradation products have two residues missing from the N-terminal end of the peptide. The solution structure of the basic peptide, aurein 1.2, has been determined by NMR spectroscopy to be an amphipathic alpha-helix with well-defined hydrophilic and hydrophobic regions. Certain of the aurein peptides (e.g. aureins 1.2 and 3.1) show anticancer activity in the NCI test regime, with LC50 values in the 10-5-10-4 M range. The aurein 1 peptides have only 13 amino-acid residues: these are the smallest antibiotic and anticancer active peptides yet reported from an anuran. The longer aurein 4 and 5 peptides, e.g. aurein 4.1 (GLIQTIKEKLKELAGGLVTGIQS-OH) and aurein 5. 1 (GLLDIVTGLLGNLIVDVLKPKTPAS-OH) show neither antibacterial nor anticancer activity.