Prediction of charge-induced molecular alignment: residual dipolar couplings at pH 3 and alignment in surfactant liquid crystalline phases

Investigating protein conformational energy landscapes and atomic resolution dynamics from NMR dipolar couplings: a review

Nuclear magnetic resonance spectroscopy is exquisitely sensitive to protein dynamics. In particular inter-nuclear dipolar couplings, that become measurable in solution when the protein is dissolved in a dilute liquid crystalline solution, report on all conformations sampled up to millisecond timescales. As such they provide the opportunity to describe the Boltzmann distribution present in solution at atomic resolution, and thereby to map the conformational energy landscape in unprecedented detail. The development of analytical methods and approaches based on numerical simulation and their application to numerous biologically important systems is presented.

Comparative studies of disordered proteins with similar sequences: application to Aβ40 and Aβ42

Quantitative comparisons of intrinsically disordered proteins (IDPs) with similar sequences, such as mutant forms of the same protein, may provide insights into IDP aggregation-a process that plays a role in several neurodegenerative disorders. Here we describe an approach for modeling IDPs with similar sequences that simplifies the comparison of the ensembles by utilizing a single library of structures. The relative population weights of the structures are estimated using a Bayesian formalism, which provides measures of uncertainty in the resulting ensembles. We applied this approach to the comparison of ensembles for Aβ40 and Aβ42. Bayesian hypothesis testing finds that although both Aβ species sample β-rich conformations in solution that may represent prefibrillar intermediates, the probability that Aβ42 samples these prefibrillar states is roughly an order of magnitude larger than the frequency in which Aβ40 samples such structures. Moreover, the structure of the soluble prefibrillar state in our ensembles is similar to the experimentally determined structure of Aβ that has been implicated as an intermediate in the aggregation pathway. Overall, our approach for comparative studies of IDPs with similar sequences provides a platform for future studies on the effect of mutations on the structure and function of disordered proteins.

Domain cooperativity in multidomain proteins: what can we learn from molecular alignment in anisotropic media?

Many proteins have modular design with multiple globular domains connected via flexible linkers. As a simple model of such system, we study a tandem construct consisting of two identical SH3 domains and a variable-length Gly/Ser linker. When the linker is short, this construct represents a dumbbell-shaped molecule with limited amount of domain-domain mobility. Due to its elongated shape, this molecule efficiently aligns in steric alignment media. As the length of the linker increases, the two domains become effectively uncoupled and begin to behave as independent entities. Consequently, their degree of alignment drops, approaching that found in the (near-spherical) isolated SH3 domains. To model the dependence of alignment parameters on the length of the interdomain linker, we have generated in silico a series of conformational ensembles representing SH3 tandems with different linker length. These ensembles were subsequently used as input for alignment prediction software PALES. The predicted alignment tensors were compared with the results of experimental measurements using a series of tandem-SH3 samples in PEG/hexanol alignment media. This comparison broadly confirmed the expected trends. At the same time, it has been found that the isolated SH3 domain aligns much stronger than expected. This finding can be attributed to complex morphology of the PEG/hexanol media and/or to weak site-specific interactions between the protein and the media. In the latter case, there are strong indications that electrostatic interactions may play a role. The fact that PEG/hexanol does not behave as a simple steric media should serve as a caution for studies that use PALES as a quantitative prediction tool (especially for disordered proteins). Further progress in this area depends on our ability to accurately model the anisotropic media and its site-specific interactions with protein molecules. Once this ability is improved, it should be possible to use the alignment parameters as a measure of domain-domain cooperativity, thus identifying the situations where two domains transiently interact with each other or become coupled through a partially structured linker.

Ligand-induced changes in the structure and dynamics of Escherichia coli peptide deformylase

Peptide deformylase (PDF) is an enzyme that is responsible for removing the formyl group from nascently synthesized polypeptides in bacteria, attracting much attention as a potential target for novel antibacterial agents. Efforts to develop potent inhibitors of the enzyme have progressed on the basis of classical medicinal chemistry, combinatorial chemistry, and structural approaches, yet the validity of PDF as an antibacterial target hangs, in part, on the ability of inhibitors to selectively target this enzyme in favor of structurally related metallohydrolases. We have used (15)N NMR spectroscopy and isothermal titration calorimetry to investigate the high-affinity interaction of EcPDF with actinonin, a naturally occurring potent EcPDF inhibitor. Backbone amide chemical shifts, residual dipolar couplings, hydrogen-deuterium exchange, and (15)N relaxation reveal structural and dynamic effects of ligand binding in the immediate vicinity of the ligand-binding site as well as at remote sites. A comparison of the crystal structures of free and actinonin-bound EcPDF with the solution data suggests that most of the consequences of the ligand binding to the protein are lost or obscured during crystallization. The results of these studies improve our understanding of the thermodynamic global minimum and have important implications for structure-based drug design.

NMR: prediction of molecular alignment from structure using the PALES software

Orientational restraints such as residual dipolar couplings promise to overcome many of the problems that traditionally limited liquid-state nuclear magnetic resonance spectroscopy. Recently, we developed methods to predict a molecular alignment tensor and thus residual dipolar couplings for a given molecular structure. This provides many new opportunities for the study of the structure and dynamics of proteins, nucleic acids, oligosaccharides and small molecules. This protocol details the use of the software PALES (Prediction of AlignmEnt from Structure) for prediction of an alignment tensor from a known three-dimensional (3D) coordinate file of a solute. The method is applicable to alignment of molecules in many neutral and charged orienting media and takes into account the molecular shape and 3D charge distribution of the molecule.

A unifying probabilistic framework for analyzing residual dipolar couplings

Residual dipolar couplings provide complementary information to the nuclear Overhauser effect measurements that are traditionally used in biomolecular structure determination by NMR. In a de novo structure determination, however, lack of knowledge about the degree and orientation of molecular alignment complicates the analysis of dipolar coupling data. We present a probabilistic framework for analyzing residual dipolar couplings and demonstrate that it is possible to estimate the atomic coordinates, the complete molecular alignment tensor, and the error of the couplings simultaneously. As a by-product, we also obtain estimates of the uncertainty in the coordinates and the alignment tensor. We show that our approach encompasses existing methods for determining the alignment tensor as special cases, including least squares estimation, histogram fitting, and elimination of an explicit alignment tensor in the restraint energy.

Evidence of molecular alignment fluctuations in aqueous dilute liquid crystalline media

Protein dynamics can be studied by NMR measurements of aqueous dilute liquid crystalline samples. However, the measured residual dipolar couplings are sensitive not only to internal fluctuations but to all changes in internuclear vectors relative to the laboratory frame. We show that side-chain fluctuations and bond librations in the ps-ns time scale perturb the molecular shape and charge distribution of a small globular protein sufficiently to cause a noticeable variation in the molecular alignment. The alignment variation disperses the bond vectors of a conformational ensemble even further from the dispersion already caused by internal fluctuations of a protein. Consequently RDC-probed order parameters are lower than those obtained by laboratory frame relaxation measurements.

Theoretical framework for NMR residual dipolar couplings in unfolded proteins

A theoretical framework for the prediction of nuclear magnetic resonance (NMR) residual dipolar couplings (RDCs) in unfolded proteins under weakly aligning conditions is presented. The unfolded polypeptide chain is modeled as a random flight chain while the alignment medium is represented by a set of regularly arranged obstacles. For the case of bicelles oriented perpendicular to the magnetic field, a closed-form analytical result is derived. With the obtained analytical expression the RDCs are readily accessible for any locus along the chain, for chains of differing length, and for varying bicelle concentrations. The two general features predicted by the model are (i) RDCs in the center segments of a polypeptide chain are larger than RDCs in the end segments, resulting in a bell-shaped sequential distribution of RDCs, and (ii) couplings are larger for shorter chains than for longer chains at a given bicelle concentration. Experimental data available from the literature confirm the first prediction of the model, providing a tool for recognizing fully unfolded polypeptide chains. With less certainty experimental data appear to support the second prediction as well. However, more systematic experimental studies are needed in order to validate or disprove the predictions of the model. The presented framework is an important step towards a solid theoretical foundation for the analysis of experimentally measured RDCs in unfolded proteins in the case of alignment media such as polyacrylamide gels and neutral bicelle systems which align biomacromolecules by a steric mechanism. Various improvements and generalizations are possible within the suggested approach.