Antigen binding forces of individually addressed single-chain Fv antibody molecules

Antibody single-chain Fv fragment (scFv) molecules that are specific for fluorescein have been engineered with a C-terminal cysteine for a directed immobilization on a flat gold surface. Individual scFv molecules can be identified by atomic force microscopy. For selected molecules the antigen binding forces are then determined by using a tip modified with covalently immobilized antigen. An scFv mutant of 12% lower free energy for ligand binding exhibits a statistically significant 20% lower binding force. This strategy of covalent immobilization and measuring well separated single molecules allows the characterization of ligand binding forces in molecular repertoires at the single molecule level and will provide a deeper insight into biorecognition processes.

Equilibrium constants and free energies in unfolding of proteins in urea solutions

A novel thermodynamic approach to the reversible unfolding of proteins in aqueous urea solutions has been developed based on the premise that urea ligands are bound cooperatively to the macromolecule. When successive stoichiometric binding constants have values larger than expected from statistical effects, an equation for moles of bound urea can be derived that contains imaginary terms. For a very steep unfolding curve, one can then show that the fraction of protein unfolded, B̄, depends on the square of the urea concentration, U, and is given by \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} \begin{equation*}\bar {B}=\frac{{\mathit{A}}^{{\mathit{2}}}_{{\mathit{1}}}{\mathit{e}}^{{\mathrm{{\lambda}}}n\bar {B}}{\mathit{U}}^{{\mathit{2}}}}{{\mathrm{1\hspace{.167em}+\hspace{.167em}}}{\mathit{A}}^{{\mathrm{2}}}_{{\mathrm{1}}}{\mathit{e}}^{{\mathrm{{\lambda}}}\bar {B}}{\mathit{U}}^{{\mathrm{2}}}}{\mathrm{.}}\end{equation*}\end{document} A12 is the binding constant as B̄→ 0, and λ is a parameter that reflects the augmentation in affinities of protein for urea as the moles bound increases to the saturation number, n. This equation provides an analytic expression that reproduces the unfolding curve with good precision, suggests a simple linear graphical procedure for evaluating A12 and λ, and leads to the appropriate standard free energy changes. The calculated ΔG° values reflect the coupling of urea binding with unfolding of the protein. Some possible implications of this analysis to protein folding in vivo are described.

Differences in water release for the binding of EcoRI to specific and nonspecific DNA sequences.

The free energy difference between complexes of the restriction nuclease EcoRI with nonspecific DNA and with the enzyme's recognition sequence is linearly dependent on the water chemical potential of the solution, set using several very different solutes, ranging from glycine and glycerol to triethylene glycol and sucrose. This osmotic dependence indicates that the nonspecific complex sequesters some 110 waters more than the specific complex with the recognition sequence. The insensitivity of the difference in number of waters released to the solute identity further indicates that this water is sequestered in a space that is sterically inaccessible to solutes, most likely at the protein-DNA interface of the nonspecific complex. Calculations based on the structure of the specific complex suggest that the apposing DNA and protein surfaces in the nonspecific complex retain approximately a full hydration layer of water.

DNA-strand exchange promoted by RecA protein in the absence of ATP: implications for the mechanism of energy transduction in protein-promoted nucleic acid transactions.

DNA-strand exchange promoted by Escherichia coli RecA protein normally requires the presence of ATP and is accompanied by ATP hydrolysis, thereby implying a need for ATP hydrolysis. Previously, ATP hydrolysis was shown not to be required; here we demonstrate furthermore that a nucleoside triphosphate cofactor is not required for DNA-strand exchange. A gratuitous allosteric effector consisting of the noncovalent complex of ADP and aluminum fluoride, ADP.AIF4-, can both induce the high-affinity DNA-binding state of RecA protein and support the homologous pairing and exchange of up to 800-900 bp of DNA. These results demonstrate that induction of the functionally active, high-affinity DNA-binding state of RecA protein is needed for RecA protein-promoted DNA-strand exchange and that there is no requirement for a high-energy nucleotide cofactor for the exchange of DNA strands. Consequently, the free energy needed to activate the DNA substrates for DNA-strand exchange is not derived from ATP hydrolysis. Instead, the needed free energy is derived from ligand binding and is transduced to the DNA via the associated ligand-induced structural transitions of the RecA protein-DNA complex; ATP hydrolysis simply destroys the effector ligand. This concept has general applicability to the mechanism of energy transduction by proteins.

Interface Energy Measurement of MgO and ZnO: Understanding the Thermodynamic Stability of Nanoparticles

Nanomaterials have triggered excitement in both fundamental science and technological applications in several fields However, the same characteristic high interface area that is responsible for their unique properties causes unconventional instability, often leading to local collapsing during application Thermodynamically, this can be attributed to an increased contribution of the interface to the free energy, activating phenomena such as sintering and grain growth The lack of reliable interface energy data has restricted the development of conceptual models to allow the control of nanoparticle stability on a thermodynamic basis. Here we introduce a novel and accessible methodology to measure interface energy of nanoparticles exploiting the heat released during sintering to establish a quantitative relation between the solid solid and solid vapor interface energies. We exploited this method in MgO and ZnO nanoparticles and determined that the ratio between the solid solid and solid vapor interface energy is 11 for MgO and 0.7 for ZnO. We then discuss that this ratio is responsible for a thermodynamic metastable state that may prevent collapsing of nanoparticles and, therefore, may be used as a tool to design long-term stable nanoparticles.

A computational analysis of substrate binding strength by phosphorylase kinase and protein kinase A

Protein kinases exhibit various degrees of substrate specificity. The large number of different protein kinases in the eukaryotic proteomes makes it impractical to determine the specificity of each enzyme experimentally. To test if it were possible to discriminate potential substrates from non-substrates by simple computational techniques, we analysed the binding enthalpies of modelled enzyme-substrate complexes and attempted to correlate it with experimental enzyme kinetics measurements. The crystal structures of phosphorylase kinase and cAMP-dependent protein kinase were used to generate models of the enzyme with a series of known peptide substrates and non-substrates, and the approximate enthalpy of binding assessed following energy minimization. We show that the computed enthalpies do not correlate closely with kinetic measurements, but the method can distinguish good substrates from weak substrates and non-substrates. Copyright (C) 2002 John Wiley Sons, Ltd.

In silico analysis of microRNAS targeting the HLA-G 3 ` untranslated region alleles and haplotypes

It has been reported that microRNAs (miRNA) may have allele-specific targeting for the 3` untranslated region (3` UTR) of the HLA-G locus. In a previous study, we reported 11 3`UTR haplotypes encompassing the 14-bp insertion/deletion polymorphism and seven SNPs (+3003 T/C, +3010 C/G, +3027 C/A, +3035 C/T, +3142 C/G, +3187A/G,and +3196 C/G), of which only the +3142 C/G SNP has been reported to influence the binding of miRNAs. Using bioinformatics analyses, we identified putative miRNA-binding sites considering the haplotypes encompassing these eight polymorphic sites, and we ranked the lowest free energies that could potentially lead to an mRNA degradation or translational repression. When a specific haplotype or a particular SNP was associated with a miRNA-binding site, we defined a free energy difference of 4 kcal/mol between alleles to classify them energetically distant. The best results were obtained for the miR-513a-5p, miR-518c*, miR-1262 and miR-92a-1*, miR-92a-2*, miR-661, miR-1224-5p, and miR-433 miRNAs, all influencing one or more of the +3003, +3010, +3027, and +3035 SNPs. The miR-2110, miR-93, miR-508-5p, miR-331-5p, miR-616, miR-513b, and miR-589* miRNAs targeted the 14-bp fragment region, and miR-148a, miR-19a*, miR-152, mir-148b,and miR-218-2 also influenced the +3142C/G polymorphism. These results suggest that these miRNAs might play a relevant role on the HLA-G expression pattern. (C) 2009 Published by Elsevier Inc. on behalf of American Society for Histocompatibility and Immunogenetics.

Crystal structures of free, IMP-, and GMP-bound Escherichia coli hypoxanthine phosphoribosyltransferase

Crystal structures have been determined for free Escherichia coli hypoxanthine phosphoribosyltransferase (HPRT) (2.9 Angstrom resolution) and for the enzyme in complex with the reaction products, inosine 5'-monophosphate (IMP) and guanosine 5-monophosphate (GMP) (2.8 Angstrom resolution). Of the known 6-oxopurine phosphoribosyltransferase (PRTase) structures, E. coli HPRT is most similar in structure to that of Tritrichomonas foetus HGXPRT, with a rmsd for 150 Calpha atoms of 1.0 Angstrom. Comparison of the free and product bound structures shows that the side chain of Phe156 and the polypeptide backbone in this vicinity move to bind IMP or GMP. A nonproline cis peptide bond, also found in some other 6-oxopurine PRTases, is observed between Leu46 and Arg47 in both the free and complexed structures. For catalysis to occur, the 6-oxopurine PRTases have a requirement for divalent metal ion, Usually Mg2+ in vivo. In the free structure, a Mg2+, is coordinated to the side chains of Glu103 and Asp104. This interaction may be important for stabilization of the enzyme before catalysis. E. coli HPRT is unique among the known 6-oxopurine PRTases in that it exhibits a marked preference for hypoxanthine as substrate over both xanthine and guanine. The structures suggest that its substrate specificity is due to the modes of binding of the bases. In E. coli HPRT, the carbonyl oxygen of Asp 163 would likely form a hydrogen bond with the 2-exocyclic nitrogen of guanine (in the HPRT-guanine-PRib-PP-Mg2+ complex). However, hypoxanthine does not have a 2-exocyclic atom and the HPRT-IMP structure suggests that hypoxanthine is likely to occupy a different position in the purine-binding pocket.

Effects of molecular crowding by saccharides on alpha-chymotrypsin dimerization

Given the importance of protein complexes as therapeutic targets, it is necessary to understand the physical chemistry of these interactions under the crowded conditions that exist in cells. We have used sedimentation equilibrium to quantify the enhancement of the reversible homodimerization of alpha-chymotrypsin by high concentrations of the osmolytes glucose, sucrose, and raffinose. In an attempt to rationalize the ostuolyte-mediated stabilization of the a-chymotrypsin homodimer, we have used models based on binding interactions (transfer-free energy analysis) and steric interactions (excluded volume theory) to predict the stabilization. Although transfer-free energy analysis predicts reasonably well the relatively small stabilization observed for complex formation between cytochrome c and cytochrome c peroxidase, as well as that between bobtail quail lysozyme and a monoclonal Fab fragment, it underestimates the sugar-mediated stabilization of the alpha-chymotrypsin dimer. Although predictions based on excluded volume theory overestimate the stabilization, it would seem that a major determinant in the observed stabilization of the a-chymotrypsin homodimer is the thermodynamic nonideality arising from molecular crowding by the three small sugars.

Metal ions modulate the folding and stability of the tumor suppressor protein S100A2

Febs Journal (2009)276:1776-1786

Combined simulation and mutagenesis analyses reveal the involvement of key residues for peroxisome proliferator-activated receptor alpha helix 12 dynamic behavior.

The dynamic properties of helix 12 in the ligand binding domain of nuclear receptors are a major determinant of AF-2 domain activity. We investigated the molecular and structural basis of helix 12 mobility, as well as the involvement of individual residues with regard to peroxisome proliferator-activated receptor alpha (PPARalpha) constitutive and ligand-dependent transcriptional activity. Functional assays of the activity of PPARalpha helix 12 mutants were combined with free energy molecular dynamics simulations. The agreement between the results from these approaches allows us to make robust claims concerning the mechanisms that govern helix 12 functions. Our data support a model in which PPARalpha helix 12 transiently adopts a relatively stable active conformation even in the absence of a ligand. This conformation provides the interface for the recruitment of a coactivator and results in constitutive activity. The receptor agonists stabilize this conformation and increase PPARalpha transcription activation potential. Finally, we disclose important functions of residues in PPARalpha AF-2, which determine the positioning of helix 12 in the active conformation in the absence of a ligand. Substitution of these residues suppresses PPARalpha constitutive activity, without changing PPARalpha ligand-dependent activation potential.

Distributed all-atom simulations of intrinsically unstructured proteins

The Computational Biophysics Group at the Universitat Pompeu Fabra (GRIB-UPF) hosts two unique computational resources dedicated to the execution of large scale molecular dynamics (MD) simulations: (a) the ACMD molecular-dynamics software, used on standard personal computers with graphical processing units (GPUs); and (b) the GPUGRID. net computing network, supported by users distributed worldwide that volunteer GPUs for biomedical research. We leveraged these resources and developed studies, protocols and open-source software to elucidate energetics and pathways of a number of biomolecular systems, with a special focus on flexible proteins with many degrees of freedom. First, we characterized ion permeation through the bactericidal model protein Gramicidin A conducting one of the largest studies to date with the steered MD biasing methodology. Next, we addressed an open problem in structural biology, the determination of drug-protein association kinetics; we reconstructed the binding free energy, association, and dissaciociation rates of a drug like model system through a spatial decomposition and a Makov-chain analysis. The work was published in the Proceedings of the National Academy of Sciences and become one of the few landmark papers elucidating a ligand-binding pathway. Furthermore, we investigated the unstructured Kinase Inducible Domain (KID), a 28-peptide central to signalling and transcriptional response; the kinetics of this challenging system was modelled with a Markovian approach in collaboration with Frank Noe’s group at the Freie University of Berlin. The impact of the funding includes three peer-reviewed publication on high-impact journals; three more papers under review; four MD analysis components, released as open-source software; MD protocols; didactic material, and code for the hosting group.

Purification and ligand binding of a soluble class I major histocompatibility complex molecule consisting of the first three domains of H-2Kd fused to beta 2-microglobulin expressed in the baculovirus-insect cell system.

A recombinant baculovirus encoding a single-chain murine major histocompatibility complex class I molecule in which the first three domains of H-2Kd are fused to beta 2-microglobulin (beta 2-m) via a 15-amino acid linker has been isolated and used to infect lepidopteran cells. A soluble, 391-amino acid single-chain H-2Kd (SC-Kd) molecule of 48 kDa was synthesized and glycosylated in insect cells and could be purified in the absence of detergents by affinity chromatography using the anti-H-2Kd monoclonal antibody SF1.1.1.1. We tested the ability of SC-Kd to bind antigenic peptides using a direct binding assay based on photoaffinity labeling. The photoreactive derivative was prepared from the H-2Kd-restricted Plasmodium berghei circumsporozoite protein (P.b. CS) peptide 253-260 (YIPSAEKI), a probe that we had previously shown to be unable to bind to the H-2Kd heavy chain in infected cells in the absence of co-expressed beta 2-microglobulin. SC-Kd expressed in insect cells did not require additional mouse beta 2-m to bind the photoprobe, indicating that the covalently attached beta 2-m could substitute for the free molecule. Similarly, binding of the P.b. CS photoaffinity probe to the purified SC-Kd molecule was unaffected by the addition of exogenous beta 2-m. This is in contrast to H-2KdQ10, a soluble H-2Kd molecule in which beta 2-m is noncovalently bound to the soluble heavy chain, whose ability to bind the photoaffinity probe is greatly enhanced in the presence of an excess of exogenous beta 2-m. The binding of the probe to SC-Kd was allele-specific, since labeling was selectively inhibited only by antigenic peptides known to be presented by the H-2Kd molecule.

Quantifying intracellular protein binding thermodynamics during mechanotransduction based on FRET spectroscopy.

Mechanical force modulates myriad cellular functions including migration, alignment, proliferation, and gene transcription. Mechanotransduction, the transmission of mechanical forces and its translation into biochemical signals, may be mediated by force induced protein conformation changes, subsequently modulating protein signaling. For the paxillin and focal adhesion kinase interaction, we demonstrate that force-induced changes in protein complex conformation, dissociation constant, and binding Gibbs free energy can be quantified by lifetime-resolved fluorescence energy transfer microscopy combined with intensity imaging calibrated by fluorescence correlation spectroscopy. Comparison with in vitro data shows that this interaction is allosteric in vivo. Further, spatially resolved imaging and inhibitor assays show that this protein interaction and its mechano-sensitivity are equal in the cytosol and in the focal adhesions complexes indicating that the mechano-sensitivity of this interaction must be mediated by soluble factors but not based on protein tyrosine phosphorylation.

The influence of the symmetry energy on the giant monopole resonance of neutron-rich nuclei analyzed in Thomas-Fermi theory

We analyze the influence of the density dependence of the symmetry energy on the average excitation energy of the isoscalar giant monopole resonance (GMR) in stable and exotic neutron-rich nuclei by applying the relativistic extended Thomas-Fermi method in scaling and constrained calculations. For the effective nuclear interaction, we employ the relativistic mean field model supplemented by an isoscalar-isovector meson coupling that allows one to modify the density dependence of the symmetry energy without compromising the success of the model for binding energies and charge radii. The semiclassical estimates of the average energy of the GMR are known to be in good agreement with the results obtained in full RPA calculations. The present analysis is performed along the Pb and Zr isotopic chains. In the scaling calculations, the excitation energy is larger when the symmetry energy is softer. The same happens in the constrained calculations for nuclei with small and moderate neutron excess. However, for nuclei of large isospin the constrained excitation energy becomes smaller in models having a soft symmetry energy. This effect is mainly due to the presence of loosely-bound outer neutrons in these isotopes. A sharp increase of the estimated width of the resonance is found in largely neutron-rich isotopes, even for heavy nuclei, which is enhanced when the symmetry energy of the model is soft. The results indicate that at large neutron numbers the structure of the low-energy region of the GMR strength distribution changes considerably with the density dependence of the nuclear symmetry energy, which may be worthy of further characterization in RPA calculations of the response function.