Automated robust and accurate assignment of protein resonances for solid state NMR

The process of resonance assignment represents a time-consuming and potentially error-prone bottleneck in structural studies of proteins by solid-state NMR (ssNMR). Software for the automation of this process is therefore of high interest. Procedures developed through the last decades for solution-state NMR are not directly applicable for ssNMR due to the inherently lower data quality caused by lower sensitivity and broader lines, leading to overlap between peaks. Recently, the first efforts towards procedures specifically aimed for ssNMR have been realized (Schmidt et al. in J Biomol NMR 56(3):243–254, 2013). Here we present a robust automatic method, which can accurately assign protein resonances using peak lists from a small set of simple 2D and 3D ssNMR experiments, applicable in cases with low sensitivity. The method is demonstrated on three uniformly ¹³C, ¹⁵N labeled biomolecules with different challenges on the assignments. In particular, for the immunoglobulin binding domain B1 of streptococcal protein G automatic assignment shows 100 % accuracy for the backbone resonances and 91.8 % when including all side chain carbons. It is demonstrated, by using a procedure for generating artificial spectra with increasing line widths, that our method, GAMES_ASSIGN can handle a significant amount of overlapping peaks in the assignment. The impact of including different ssNMR experiments is evaluated as well.     

Two-dimensional solid-state NMR applied to a chimeric potassium channel

Solid-state NMR (ssNMR) represents a spectroscopic method to study membrane protein structure and dynamics in lipid bilayers. We present two-dimensional correlation experiments conducted on a fully [13C,15N] labeled version of a chimeric potassium (KcsA-Kv1.3) channel. Data obtained by using two different ion concentrations suggest a structural conservation of the selectivity filter region. SsNMR experiments conducted at two different temperatures point to differential molecular dynamics of the channel.     

ICMRBS founder’s medal 2006: Biological solid-state NMR,methods and applications

Solid-state NMR (ssNMR) provides increasing possibilities to study structure and dynamics of biomolecular systems. Our group has been interested in developing ssNMR-based approaches that are applicable to biomolecules of increasing molecular size and complexity without the need of specific isotope-labelling. Methodological aspects ranging from spectral assignments to the indirect detection of proton–proton contacts in multi-dimensional ssNMR are discussed and applied to (membrane) protein complexes.     

Apolipoprotein E isoform-specific binding to the low-density lipoprotein receptor

Apolipoprotein E (apoE) is a ligand for members of the low-density lipoprotein receptor (LDLR) family and functions in plasma cholesterol homeostasis. A fluorescence-based assay has been employed in molecular studies of receptor-ligand interactions. Competition experiments revealed isoform-specific differences in binding of lipid-associated apoE N terminal (NT) domain to a recombinant soluble LDLR (sLDLR). In a similar manner, lipid--associated-but not lipid-free--full-length apoE3 showed binding activity to sLDLR. The molecular chaperone, receptor-associated protein, inhibited apoE3-NT-phospholipid complex binding to sLDLR. Kinetic studies of apoE3-NT-phospholipid complex interaction with sLDLR revealed time-dependent effects of apoE-NT isoform binding to sLDLR. The results reveal a discerning method for study of the molecular basis of ligand interactions that likely influence receptor function in maintenance of whole body cholesterol homeostasis.     

Two-Dimensional Solid-State NMR Applied to a Chimeric Potassium Channel

Solid-state NMR (ssNMR) represents a spectroscopic method to study membrane protein structure and dynamics in lipid bilayers. We present two-dimensional correlation experiments conducted on a fully [¹³C,¹⁵N] labeled version of a chimeric potassium (KcsA-Kv1.3) channel. Data obtained by using two different ion concentrations suggest a structural conservation of the selectivity filter region. SsNMR experiments conducted at two different temperatures point to differential molecular dynamics of the channel.     

Determination of structural topology of a membrane protein in lipid bilayers using polarization optimized experiments (POE) for static and MAS solid state NMR spectroscopy

The low sensitivity inherent to both the static and magic angle spinning techniques of solid-state NMR (ssNMR) spectroscopy has thus far limited the routine application of multidimensional experiments to determine the structure of membrane proteins in lipid bilayers. Here, we demonstrate the advantage of using a recently developed class of experiments, polarization optimized experiments, for both static and MAS spectroscopy to achieve higher sensitivity and substantial time-savings for 2D and 3D experiments. We used sarcolipin, a single pass membrane protein, reconstituted in oriented bicelles (for oriented ssNMR) and multilamellar vesicles (for MAS ssNMR) as a benchmark. The restraints derived by these experiments are then combined into a hybrid energy function to allow simultaneous determination of structure and topology. The resulting structural ensemble converged to a helical conformation with a backbone RMSD ~0.44 Å, a tilt angle of 24° ± 1°, and an azimuthal angle of 55° ± 6°. This work represents a crucial first step toward obtaining high-resolution structures of large membrane proteins using combined multidimensional oriented solid-state NMR and magic angle spinning solid-state NMR.     

Molecular dynamics simulation of the induced‐fit binding process of DNA aptamer and L‐argininamide

Aptamers are rare functional nucleic acids with binding affinity to and specificity for target ligands. Recent experiments have lead to the proposal of an induced‐fit binding mechanism for L ‐argininamide (Arm) and its binding aptamer. However, at the molecular level, this mechanism between the aptamer and its coupled ligand is still poorly understood. The present study used explicit solvent molecular dynamics (MD) simulations to examine the critical bases involved in aptamer‐Arm binding and the induced‐fit binding process at atomic resolution. The simulation results revealed that the Watson‐Crick pair (G10‐C16), C9, A12, and C17 bases play important roles in aptamer‐Arm binding, and that binding of Arm results in an aptamer conformation optimized through a general induced‐fit process. In an aqueous solution, the mechanism has the following characteristic stages: (a) adsorption stage, the Arm anchors to the binding site of aptamer with strong electrostatic interaction; (b) binding stage, the Arm fits into the binding site of aptamer by hydrogen‐bond formation; and (c) complex stabilization stage, the hydrogen bonding and electrostatic interactions cooperatively stabilize the complex structure. This study provides dynamics information on the aptamer‐ligand induced‐fit binding mechanism. The critical bases in aptamer‐ligand binding may provide a guideline in aptamer design for molecular recognition engineering.     

Mathematical analysis of ligand-induced monomerization and dimerization in the monomer dimer equilibrium of proteins. Application to D-amino acid oxidase.

We have mathematically analyzed ligand-induced monomerization and dimerization in a protein monomer-dimer equilibrium system, in which the monomer has one and the dimer two binding sites. These dimer sites have the same binding constants for the first ligand but may cooperatively interact when one of them is occupied by a ligand molecule. In this system, the apparent dimerization constant and the apparent molecular weight are functions of free ligand concentration, and depend on the intrinsic binding constants of the ligand molecule to the monomer and the dimer. The behavior of these functions is classified into 17 cases according to the values of the three intrinsic binding constants, and some calculated examples are shown graphically for selected parameters. The theory was also applied to D-amino acid oxidase [EC 1.4.3.3], a flavoprotein, and the pH dependence of the apparent dimerization constant and the apparent molecular weight in the presence of ligand, p-aminobenzoate, were studied theoretically using parameters obtained in our previous experiments (5).     

Probabilistic modeling of rosette formation

Rosetting, or forming a cell aggregate between a single target nucleated cell and a number of red blood cells (RBCs), is a simple assay for cell adhesion mediated by specific receptor-ligand interaction. For example, rosette formation between sheep RBC and human lymphocytes has been used to differentiate T cells from B cells. Rosetting assay is commonly used to determine the interaction of Fc gamma-receptors (FcgammaR) expressed on inflammatory cells and IgG coated on RBCs. Despite its wide use in measuring cell adhesion, the biophysical parameters of rosette formation have not been well characterized. Here we developed a probabilistic model to describe the distribution of rosette sizes, which is Poissonian. The average rosette size is predicted to be proportional to the apparent two-dimensional binding affinity of the interacting receptor-ligand pair and their site densities. The model has been supported by experiments of rosettes mediated by four molecular interactions: FcgammaRIII interacting with IgG, T cell receptor and coreceptor CD8 interacting with antigen peptide presented by major histocompatibility molecule, P-selectin interacting with P-selectin glycoprotein ligand 1 (PSGL-1), and L-selectin interacting with PSGL-1. The latter two are structurally similar and are different from the former two. Fitting the model to data enabled us to evaluate the apparent effective two-dimensional binding affinity of the interacting molecular pairs: 7.19x10(-5) microm4 for FcgammaRIII-IgG interaction, 4.66x10(-3) microm4 for P-selectin-PSGL-1 interaction, and 0.94x10(-3) microm4 for L-selectin-PSGL-1 interaction. These results elucidate the biophysical mechanism of rosette formation and enable it to become a semiquantitative assay that relates the rosette size to the effective affinity for receptor-ligand binding.     

Structural characterization of Ca(2+)-ATPase-bound phospholamban in lipid bilayers by solid-state nuclear magnetic resonance (NMR) spectroscopy

Phospholamban (PLN) regulates cardiac contractility by modulation of sarco(endo)plasmic reticulum calcium ATPase (SERCA) activity. While PLN and SERCA1a, an isoform from skeletal muscle, have been structurally characterized in great detail, direct information about the conformation of PLN in complex with SERCA has been limited. We used solid-state NMR (ssNMR) spectroscopy to deduce structural properties of both the A 36F 41A 46 mutant (AFA-PLN) and wild-type PLN (WT-PLN) when bound to SERCA1a after reconstitution in a functional lipid bilayer environment. Chemical-shift assignments in all domains of AFA-PLN provide direct evidence for the presence of two terminal alpha helices connected by a linker region of reduced structural order that differs from previous findings on free PLN. ssNMR experiments on WT-PLN show no significant difference in binding compared to AFA-PLN and do not support the coexistence of a significantly populated dynamic state of PLN after formation of the PLN/SERCA complex. A combination of our spectroscopic data with biophysical and biochemical data using flexible protein-protein docking simulations provides a structural basis for understanding the interaction between PLN and SERCA1a.     

Deuterium off-resonance rotating frame spin-lattice relaxation of macromolecular bound ligands.

Deuterated 3-trimethylsilylpropionic acid binding to bovine serum albumin was used as a model system to examine the feasibility and limitations of using the deuterium off-resonance rotating frame spin-lattice relaxation experiment for the study of equilibrium ligand-binding behavior to proteins. The results of this study demonstrate that the rotational-diffusion behavior of the bound species can be monitored directly, i.e., the observed correlation time of the ligand in the presence of a protein is approximately equal to the correlation time of the ligand in the bound state, provided that the fraction of bound ligand is at least 0.20. The presence of local ligand motion and/or chemical exchange contributions to relaxation in the bound state was inferred from the observation that the correlation time of the bound ligand was somewhat smaller than the correlation time characterizing the overall tumbling of the protein. An approximate value for the fraction of bound ligand was obtained from off-resonance relaxation experiments when supplemental spin-lattice or transverse relaxation times were employed in the analysis. Incorporation of local motion effects for the bound species into the theoretical relaxation formalism enabled the evaluation of an order parameter and an effective correlation time, which in conjunction with a wobbling in a cone model, provided additional information about ligand motion in the bound state.     

Solution- and Adsorbed-State Structural Ensembles Predicted for the Statherin-Hydroxyapatite System

We have developed a multiscale structure prediction technique to study solution- and adsorbed-state ensembles of biomineralization proteins. The algorithm employs a Metropolis Monte Carlo-plus-minimization strategy that varies all torsional and rigid-body protein degrees of freedom. We applied the technique to fold statherin, starting from a fully extended peptide chain in solution, in the presence of hydroxyapatite (HAp) (001), (010), and (100) monoclinic crystals. Blind (unbiased) predictions capture experimentally observed macroscopic and high-resolution structural features and show minimal statherin structural change upon adsorption. The dominant structural difference between solution and adsorbed states is an experimentally observed folding event in statherin's helical binding domain. Whereas predicted statherin conformers vary slightly at three different HAp crystal faces, geometric and chemical similarities of the surfaces allow structurally promiscuous binding. Finally, we compare blind predictions with those obtained from simulation biased to satisfy all previously published solid-state NMR (ssNMR) distance and angle measurements (acquired from HAp-adsorbed statherin). Atomic clashes in these structures suggest a plausible, alternative interpretation of some ssNMR measurements as intermolecular rather than intramolecular. This work demonstrates that a combination of ssNMR and structure prediction could effectively determine high-resolution protein structures at biomineral interfaces.     

Biophysical and molecular docking insight into interaction mechanism and thermal stability of human serum albumin isoforms with a semi-synthetic water-soluble camptothecin analog irinotecan hydrochloride

In the present work, we have examined the binding parameters, thermodynamics, and stability of human serum albumin (HSA) isoforms at pH 7.4 and 9.0, using spectroscopic, calorimetric, and molecular docking methods in the presence of water-soluble camptothecin analog irinotecan hydrochloride (CPT-11). We observed that CPT-11 binds to HSA through a static quenching procedure of ground-state complex formation with N-isoform and B-isoform. Hydrogen bond and hydrophobic interactions are the major governing forces that participating in the formation of protein–drug complex. To determine the binding site of CPT-11 within HSA molecules, we also have performed molecular docking experiments. We explored the CPT-11-mediated stability and modulation of HSA by performing dynamic light scattering (DLS) and differential scanning calorimetry (DSC) experiments. DLS and DSC techniques are used to determine the size and the melting point (T_m) of HSA, which was decreased in the presence of CPT-11. Therefore, CPT-11 plays an important role in HSA stability and protein–ligand interactions. The present study provides valuable information in the field of pharmacokinetics, pharmaco-dynamics, and drug discovery.     

Seeing the light: preassembly and ligand-induced changes of the interferon gamma receptor complex in cells

Our experiments were designed to test the hypothesis that the cell surface interferon gamma receptor chains are preassembled rather than associated by ligand and to assess the molecular changes on ligand binding. To accomplish this, we used fluorescence resonance energy transfer, a powerful spectroscopic technique that has been used to determine molecular interactions and distances between the donor and acceptor. However, current commercial instruments do not provide sufficient sensitivity or the full spectra to provide decisive results of interactions between proteins labeled with blue and green fluorescent proteins in living cells. In our experiments, we used the blue fluorescent protein and green fluorescent protein pair, attached a monochrometer and charge-coupled device camera to a modified confocal microscope, reduced background fluorescence with the use of two-photon excitation, and focused on regions of single cells to provide clear spectra of fluorescence resonance energy transfer. In contrast to the prevailing view, the results demonstrate that the receptor chains are preassociated and that the intracellular domains move apart on binding the ligand interferon gamma. Application of this technology should lead to new rapid methods for high throughput screening and delineation of the interactome of cells.     

Switches of hydrogen bonds during ligand–protein association processes determine binding kinetics

Revealing the processes of ligand–protein associations deepens our understanding of molecular recognition and binding kinetics. Hydrogen bonds (H‐bonds) play a crucial role in optimizing ligand–protein interactions and ligand specificity. In addition to the formation of stable H‐bonds in the final bound state, the formation of transient H‐bonds during binding processes contributes binding kinetics that define a ligand as a fast or slow binder, which also affects drug action. However, the effect of forming the transient H‐bonds on the kinetic properties is little understood. Guided by results from coarse‐grained Brownian dynamics simulations, we used classical molecular dynamics simulations in an implicit solvent model and accelerated molecular dynamics simulations in explicit waters to show that the position and distribution of the H‐bond donor or acceptor of a drug result in switching intermolecular and intramolecular H‐bond pairs during ligand recognition processes. We studied two major types of HIV‐1 protease ligands: a fast binder, xk263, and a slow binder, ritonavir. The slow association rate in ritonavir can be attributed to increased flexibility of ritonavir, which yields multistep transitions and stepwise entering patterns and the formation and breaking of complex H‐bond pairs during the binding process. This model suggests the importance of conversions of spatiotemporal H‐bonds during the association of ligands and proteins, which helps in designing inhibitors with preferred binding kinetics. Copyright © 2014 John Wiley & Sons, Ltd.     

Investigating the binding behaviour of two avidin‐based testosterone binders using molecular recognition force spectroscopy

Molecular recognition force spectroscopy, a biosensing atomic force microscopy technique allows to characterise the dissociation of ligand–receptor complexes at the molecular level. Here, we used molecular recognition force spectroscopy to study the binding capability of recently developed testosterone binders. The two avidin‐based proteins called sbAvd‐1 and sbAvd‐2 are expected to bind both testosterone and biotin but differ in their binding behaviour towards these ligands. To explore the ligand binding and dissociation energy landscape of these proteins, we tethered biotin or testosterone to the atomic force microscopy probe while the testosterone‐binding protein was immobilized on the surface. Repeated formation and rupture of the ligand–receptor complex at different pulling velocities allowed determination of the loading rate dependence of the complex‐rupturing force. In this way, we obtained the molecular dissociation rate (k_off) and energy landscape distances (x_β) of the four possible complexes: sbAvd‐1‐biotin, sbAvd‐1‐testosterone, sbAvd‐2‐biotin and sbAvd‐2‐testosterone. It was found that the kinetic off‐rates for both proteins and both ligands are similar. In contrast, the x_β values, as well as the probability of complex formations, varied considerably. In addition, competitive binding experiments with biotin and testosterone in solution differ significantly for the two testosterone‐binding proteins, implying a decreased cross‐reactivity of sbAvd‐2. Unravelling the binding behaviour of the investigated testosterone‐binding proteins is expected to improve their usability for possible sensing applications. Copyright © 2014 John Wiley & Sons, Ltd.     

Atomistic mechanism of polyphenol amyloid aggregation inhibitors: molecular dynamics study of Curcumin,Exifone, and Myricetin interaction with the segment of tau peptide oligomer

Amyloid fibrils are highly ordered protein aggregates associated with many diseases affecting millions of people worldwide. Polyphenols such as Curcumin, Exifone, and Myricetin exhibit modest inhibition toward fibril formation of tau peptide which is associated with Alzheimer’s disease. However, the molecular mechanisms of this inhibition remain elusive. We investigated the binding of three polyphenol molecules to the protofibrils of an amyloidogenic fragment VQIVYK of tau peptide by molecular dynamics simulations in explicit solvent. We find that polyphenols induce conformational changes in the oligomer aggregate. These changes disrupt the amyloid H bonding, perturbing the aggregate. While the structural evolution of the control oligomer with no ligand is limited to the twisting of the β-sheets without their disassembly, the presence of polyphenol molecule pushes the β-sheets apart, and leads to a loosely packed structure where two of four β-sheets dissociate in each of the three cases considered here. The H-bonding capacity of polyphenols is responsible for the observed behavior. The calculated binding free energies and its individual components enabled better understanding of the binding. Results indicated that the contribution from Van der Waals interactions is more significant than electrostatic contribution to the binding. The findings from this study are expected to assist in the development of aggregation inhibitors. Significant binding between polyphenols and aggregate oligomer identified in our simulations confirms the previous experimental observations in which polyphenols refold the tau peptide without forming covalent bonds.     

Transient ligand docking sites in Cerebratulus lacteus mini-hemoglobin

The monomeric hemoglobin of the nemertean worm Cerebratulus lacteus functions as an oxygen storage protein to maintain neural activity under hypoxic conditions. It shares a large, apolar matrix tunnel with other small hemoglobins, which has been implicated as a potential ligand migration pathway. Here we explore ligand migration and binding within the distal heme pocket, to which the tunnel provides access to ligands from the outside. FTIR/TDS experiments performed at cryogenic temperatures reveal the presence of three transient ligand docking sites within the distal pocket, the primary docking site B on top of pyrrole C and secondary sites C and D. Site C is assigned to a cavity adjacent to the distal portion of the heme pocket, surrounded by the B and E helices. It has an opening to the apolar tunnel and is expected to be on the pathway for ligand entry and exit, whereas site D, circumscribed by TyrB10, GlnE7, and the CD corner, most likely is located on a side pathway of ligand migration. Flash photolysis experiments at ambient temperatures indicate that the rate-limiting step for ligand binding to CerHb is migration through the apolar channel to site C. Movement from C to B and iron-ligand bond formation involve low energy barriers and thus are very rapid processes in the wt protein.     

Ligand Binding Mechanics of Maltose Binding Protein

In the past decade, single-molecule force spectroscopy has provided new insights into the key interactions stabilizing folded proteins. A few recent studies probing the effects of ligand binding on mechanical protein stability have come to quite different conclusions. While some proteins seem to be stabilized considerably by a bound ligand, others appear to be unaffected. Since force acts as a vector in space, it is conceivable that mechanical stabilization by ligand binding is dependent on the direction of force application. In this study, we vary the direction of the force to investigate the effect of ligand binding on the stability of maltose binding protein (MBP). MBP consists of two lobes connected by a hinge region that move from an open to a closed conformation when the ligand maltose binds. Previous mechanical experiments, where load was applied to the N and C termini, have demonstrated that MBP is built up of four building blocks (unfoldons) that sequentially detach from the folded structure. In this study, we design the pulling direction so that force application moves the two MBP lobes apart along the hinge axis. Mechanical unfolding in this geometry proceeds via an intermediate state whose boundaries coincide with previously reported MBP unfoldons. We find that in contrast to N-C-terminal pulling experiments, the mechanical stability of MBP is increased by ligand binding when load is applied to the two lobes and force breaks the protein-ligand interactions directly. Contour length measurements indicate that MBP is forced into an open conformation before unfolding even if ligand is bound. Using mutagenesis experiments, we demonstrate that the mechanical stabilization effect is due to only a few key interactions of the protein with its ligand. This work illustrates how varying the direction of the applied force allows revealing important details about the ligand binding mechanics of a large protein.     

An application of CIFAP for predicting the binding affinity of Chk1 inhibitors derived from 2‐aminothiazole‐4‐carboxamide

Investigation of protein‐ligand interactions obtained from experiments has a crucial part in the design of newly discovered and effective drugs. Analyzing the data extracted from known interactions could help scientists to predict the binding affinities of promising ligands before conducting experiments. The objective of this study is to advance the CIFAP (compressed images for affinity prediction) method, which is relevant to a protein‐ligand model, identifying 2D electrostatic potential images by separating the binding site of protein‐ligand complexes and using the images for predicting the computational affinity information represented by pIC₅₀ values. The CIFAP method has 2 phases, namely, data modeling and prediction. In data modeling phase, the separated 3D structure of the binding pocket with the ligand inside is fitted into an electrostatic potential grid box, which is then compressed through 3 orthogonal directions into three 2D images for each protein‐ligand complex. Sequential floating forward selection technique is performed for acquiring prediction patterns from the images. In the prediction phase, support vector regression (SVR) and partial least squares regression are used for testing the quality of the CIFAP method for predicting the binding affinity of 45 CHK1 inhibitors derived from 2‐aminothiazole‐4‐carboxamide. The results show that the CIFAP method using both support vector regression and partial least squares regression is very effective for predicting the binding affinities of CHK1‐ligand complexes with low‐error values and high correlation. As a future work, the results could be improved by working on the pose of the ligands inside the grid.