Van der Waals interactions involving proteins.

Van der Waals (dispersion) forces contribute to interactions of proteins with other molecules or with surfaces, but because of the structural complexity of protein molecules, the magnitude of these effects is usually estimated based on idealized models of the molecular geometry, e.g., spheres or spheroids. The calculations reported here seek to account for both the geometric irregularity of protein molecules and the material properties of the interacting media. Whereas the latter are found to fall in the generally accepted range, the molecular shape is shown to cause the magnitudes of the interactions to differ significantly from those calculated using idealized models, with important consequences. First, the roughness of the molecular surface leads to much lower average interaction energies for both protein-protein and protein-surface cases relative to calculations in which the protein molecule is approximated as a sphere. These results indicate that a form of steric stabilization may be an important effect in protein solutions. Underlying this behavior is appreciable orientational dependence, one reflection of which is that molecules of complementary shape are found to exhibit very strong attractive dispersion interactions. Although this has been widely discussed previously in the context of molecular recognition processes, the broader implications of these phenomena may also be important at larger molecular separations, e.g., in the dynamics of aggregation, precipitation, and crystal growth.     

Molecular origins of osmotic second virial coefficients of proteins.

The thermodynamic properties of protein solutions are determined by the molecular interactions involving both solvent and solute molecules. A quantitative understanding of the relationship would facilitate more systematic procedures for manipulating the properties in a process environment. In this work the molecular basis for the osmotic second virial coefficient, B22, is studied; osmotic effects are critical in membrane transport, and the value of B22 has also been shown to correlate with protein crystallization behavior. The calculations here account for steric, electrostatic, and short-range interactions, with the structural and functional anisotropy of the protein molecules explicitly accounted for. The orientational dependence of the protein interactions is seen to have a pronounced effect on the calculations; in particular, the relatively few protein-protein configurations in which the apposing surfaces display geometric complementarity contribute disproportionately strongly to B22. The importance of electrostatic interactions is also amplified in these high-complementarity configurations. The significance of molecular recognition in determining B22 can explain the correlation with crystallization behavior, and it suggests that alteration of local molecular geometry can help in manipulating protein solution behavior. The results also have implications for the role of protein interactions in biological self-organization.     

An In Silico Analysis of Troponin I Mutations in Hypertrophic Cardiomyopathy of Indian Origin

Hypertrophic Cardiomyopathy (HCM) is an autosomal dominant disorder of the myocardium which is hypertrophied resulting in arrhythmias and heart failure leading to sudden cardiac death (SCD). Several sarcomeric proteins and modifier genes have been implicated in this disease. Troponin I, being a part of the Troponin complex (troponin I, troponin C, troponin T), is an important gene for sarcomeric function. Four mutations (1 novel) were identified in Indian HCM cases, namely, Pro82Ser, Arg98Gln, Arg141Gln and Arg162Gln in Troponin I protein, which are in functionally significant domains. In order to analyse the effect of the mutations on protein stability and protein-protein interactions within the Troponin complex, an in silico study was carried out. The freely available X-ray crystal structure (PDB ID: 1JIE) was used as the template to model the protein followed by loop generation and development of troponin complex for both the troponin I wild type and four mutants (NCBI ID: PRJNA194382). The structural study was carried out to determine the effect of mutation on the structural stability and protein-protein interactions between three subunits in the complex. These mutations, especially the arginine to glutamine substitutions were found to result in local perturbations within the troponin complex by creating/removing inter/intra molecular hydrogen bonds with troponin T and troponin C. This has led to a decrease in the protein stability and loss of important interactions between the three subunits. It could have a significant impact on the disease progression when coupled with allelic heterogeneity which was observed in the cases carrying these mutations. However, this can be further confirmed by functional studies on protein levels in the identified cases.     

Modeling implicit reorganization in continuum descriptions of protein-protein interactions

The determination of free energies that govern protein-protein recognition is essential for a detailed molecular understanding of biological specificity. Continuum models of macromolecular interactions, in which the solvent is treated by an implicit representation and the proteins are treated semi-microscopically, are computationally tractable for estimating free energies, yet many questions remain concerning their accuracy. This article reports a continuum analysis of the free-energy changes underlying the binding of 31 interfacial alanine substitutions of two complexes of the antihen egg white lysozyme (HEL) antibody D1.3 bound with HEL or the antibody E5.2. Two implicit schemes for modeling the effects of protein and solvent relaxation were examined, in which the protein environment was treated as either homogeneous with a "protein dielectric constant" of epsilon(p) = 4 or inhomogeneous, with epsilon(p) = 4 for neutral residues and epsilon(p) = 25 for ionized residues. The results showed that the nonuniform dielectric model reproduced the experimental differences better, with an average absolute error of +/-1.1 kcal/mol, compared with +/-1.4 kcal/mol for the uniform model. More importantly, the error for charged residues in the nonuniform model is +/-0.8 kcal/mol and is nearly half of that corresponding to the uniform model. Several substitutions were clearly problematic in determining qualitative trends and probably required explicit structural reorganization at the protein-protein interface.     

Comprehensive and quantitative mapping of energy landscapes for protein-protein interactions by rapid combinatorial scanning

A novel, quantitative saturation (QS) scanning strategy was developed to obtain a comprehensive data base of the structural and functional effects of all possible mutations across a large protein-protein interface. The QS scan approach was applied to the high affinity site of human growth hormone (hGH) for binding to its receptor (hGHR). Although the published structure-function data base describing this system is probably the most extensive for any large protein-protein interface, it is nonetheless too sparse to accurately describe the nature of the energetics governing the interaction. Our comprehensive data base affords a complete view of the binding site and provides important new insights into the general principles underlying protein-protein interactions. The hGH binding interface is highly adaptable to mutations, but the nature of the tolerated mutations challenges generally accepted views about the evolutionary and biophysical pressures governing protein-protein interactions. Many substitutions that would be considered chemically conservative are not tolerated, while conversely, many non-conservative substitutions can be accommodated. Furthermore, conservation across species is a poor predictor of the chemical character of tolerated substitutions across the interface. Numerous deviations from generally accepted expectations indicate that mutational tolerance is highly context dependent and, furthermore, cannot be predicted by our current knowledge base. The type of data produced by the comprehensive QS scan can fill the gaps in the structure-function matrix. The compilation of analogous data bases from studies of other protein-protein interactions should greatly aid the development of computational methods for explaining and designing molecular recognition.     

Do electrostatic interactions destabilize protein–nucleic acid binding?

The negatively charged phosphates of nucleic acids are often paired with positively charged residues upon binding proteins. It was thus counter-intuitive when previous Poisson-Boltzmann (PB) calculations gave positive energies from electrostatic interactions, meaning that they destabilize protein-nucleic acid binding. Our own PB calculations on protein-protein binding have shown that the sign and the magnitude of the electrostatic component are sensitive to the specification of the dielectric boundary in PB calculations. A popular choice for the boundary between the solute low dielectric and the solvent high dielectric is the molecular surface; an alternative is the van der Waals (vdW) surface. In line with results for protein-protein binding, in this article, we found that PB calculations with the molecular surface gave positive electrostatic interaction energies for two protein-RNA complexes, but the signs are reversed when the vdW surface was used. Therefore, whether destabilizing or stabilizing effects are predicted depends on the choice of the dielectric boundary. The two calculation protocols, however, yielded similar salt effects on the binding affinity. Effects of charge mutations differentiated the two calculation protocols; PB calculations with the vdW surface had smaller deviations overall from experimental data.     

Molecular Interactions and Residues Involved in Force Generation in the T4 Viral DNA Packaging Motor

Many viruses utilize molecular motors to package their genomes into preformed capsids. A striking feature of these motors is their ability to generate large forces to drive DNA translocation against entropic, electrostatic, and bending forces resisting DNA confinement. A model based on recently resolved structures of the bacteriophage T4 motor protein gp17 suggests that this motor generates large forces by undergoing a conformational change from an extended to a compact state. This transition is proposed to be driven by electrostatic interactions between complementarily charged residues across the interface between the N- and C-terminal domains of gp17. Here we use atomistic molecular dynamics simulations to investigate in detail the molecular interactions and residues involved in such a compaction transition of gp17. We find that although electrostatic interactions between charged residues contribute significantly to the overall free energy change of compaction, interactions mediated by the uncharged residues are equally if not more important. We identify five charged residues and six uncharged residues at the interface that play a dominant role in the compaction transition and also reveal salt bridging, van der Waals, and solvent hydrogen-bonding interactions mediated by these residues in stabilizing the compact form of gp17. The formation of a salt bridge between Glu309 and Arg494 is found to be particularly crucial, consistent with experiments showing complete abrogation in packaging upon Glu309Lys mutation. The computed contributions of several other residues are also found to correlate well with single-molecule measurements of impairments in DNA translocation activity caused by site-directed mutations.     

Solvent accessibility in native and isolated domain environments: general features and implications to interface predictability

A non-redundant database of 4536 structural domains, comprising more than 790,000 residues, has been used for the calculation of their solvent accessibility in the native protein environment and then in the isolated domain environment. Nearly 140,000 (18%) residues showed a change in accessible surface area in the above two conditions. General features of this change under these two circumstances have been pointed out. Propensities of these interfacing amino acid residues have been calculated and their variation for different secondary structure types has been analyzed. Actual amount of surface area lost by different secondary structures is higher in the case of helix and strands compared to coil and other conformations. Overall change in surface area in hydrophobic and uncharged residues is higher than that in charged residues. An attempt has been made to know the predictability of interface residues from sequence environments. This analysis and prediction results have significant implications towards determining interacting residues in proteins and for the prediction of protein-protein, protein-ligand, protein-DNA and similar interactions.     

Role of charges and solvent on the conformational properties of poly(galacturonic acid) chains: a molecular dynamics study

Pectin shares with many other polysaccharides an intrinsic chemical and physical complexity. The widespread industrial applications have made it one of the most studied polysaccharides. This work presents a theoretical model of poly(galacturonic acid), the major constituent of pectin, suitable to study its structural and dynamical properties. In particular, the effects of solvent and charge status are studied. The dynamics is shown to be severely affected by the presence of charged groups on each residue, making the charged chain much more rigid than the uncharged one. A key structural property for a semirigid polymer, the asymptotic persistence length, is calculated for relatively short charged and uncharged chains in molecular water solvent using a new method. The influence of charge on structural properties of poly(galacturonic acid) is shown to be strong and solvent-dependent. In fact, a large difference is found between continuum solvent adiabatic map calculations and molecular dynamics with explicit solvent, with the latter showing a much larger persistence length.     

Using affinity chromatography to engineer and characterize pH‐dependent protein switches

Conformational changes play important roles in the regulation of many enzymatic reactions. Specific motions of side chains, secondary structures, or entire protein domains facilitate the precise control of substrate selection, binding, and catalysis. Likewise, the engineering of allostery into proteins is envisioned to enable unprecedented control of chemical reactions and molecular assembly processes. We here study the structural effects of engineered ionizable residues in the core of the glutathione‐S‐transferase to convert this protein into a pH‐dependent allosteric protein. The underlying rational of these substitutions is that in the neutral state, an uncharged residue is compatible with the hydrophobic environment. In the charged state, however, the residue will invoke unfavorable interactions, which are likely to induce conformational changes that will affect the function of the enzyme. To test this hypothesis, we have engineered a single aspartate, cysteine, or histidine residue at a distance from the active site into the protein. All of the mutations exhibit a dramatic effect on the protein's affinity to bind glutathione. Whereas the aspartate or histidine mutations result in permanently nonbinding or binding versions of the protein, respectively, mutant GST50C exhibits distinct pH‐dependent GSH‐binding affinity. The crystal structures of the mutant protein GST50C under ionizing and nonionizing conditions reveal the recruitment of water molecules into the hydrophobic core to produce conformational changes that influence the protein's active site. The methodology described here to create and characterize engineered allosteric proteins through affinity chromatography may lead to a general approach to engineer effector‐specific allostery into a protein structure.     

Effect of mutations involving charged residues on the stability of staphylococcal nuclease: a continuum electrostatics study

A continuum electrostatics model is used to calculate the relative stabilities of 117 mutants of staphylococcal nuclease (SNase) involving the mutation of a charged residue to an uncharged residue. The calculations are based on the crystallographic structure of the wild-type protein and attempt to take implicitly into account the effect of the mutations in the denatured state by assuming a linear relationship between the free energy changes caused by the mutation in the native and denatured states. A good correlation (linear correlation coefficient of approximately 0.8) is found with published experimental relative stabilities of these mutants. The results suggest that in the case of SNase (i) charged residues contribute to the stability of the native state mainly through electrostatic interactions, and (ii) native-like electrostatic interactions may persist in the denatured state. The continuum electrostatics method is only moderately sensitive to model parameters and leads to quasi-predictive results for the relative mutant stabilities (error of 2-3 kJ mol(-1) or of the order of k(B)T), except for mutants in which a charged residue is mutated to glycine.     

Structural and Functional Restraints on the Occurrence of Single Amino Acid Variations in Human Proteins

Human genetic variation is the incarnation of diverse evolutionary history, which reflects both selectively advantageous and selectively neutral change. In this study, we catalogue structural and functional features of proteins that restrain genetic variation leading to single amino acid substitutions. Our variation dataset is divided into three categories: i) Mendelian disease-related variants, ii) neutral polymorphisms and iii) cancer somatic mutations. We characterize structural environments of the amino acid variants by the following properties: i) side-chain solvent accessibility, ii) main-chain secondary structure, and iii) hydrogen bonds from a side chain to a main chain or other side chains. To address functional restraints, amino acid substitutions in proteins are examined to see whether they are located at functionally important sites involved in protein-protein interactions, protein-ligand interactions or catalytic activity of enzymes. We also measure the likelihood of amino acid substitutions and the degree of residue conservation where variants occur. We show that various types of variants are under different degrees of structural and functional restraints, which affect their occurrence in human proteome.     

Simultaneous optimization of enzyme activity and quaternary structure by directed evolution

Natural evolution has produced efficient enzymes of enormous structural diversity. We imitated this natural process in the laboratory to augment the efficiency of an engineered chorismate mutase with low activity and an unusual hexameric topology. By applying two rounds of DNA shuffling and genetic selection, we obtained a 400-fold more efficient enzyme, containing three non-active-site mutations. Detailed biophysical characterization of the evolved variant suggests that it exists predominantly as a trimer in solution, but is otherwise similarly stable as the parent hexamer. The dramatic structural and functional effects achieved by a small number of seemingly innocuous substitutions highlights the utility of directed evolution for modifying protein-protein interactions to produce novel quaternary states with optimized activities.     

Mutagenesis of the cysteine-rich clip domain in the Drosophila patterning protease, Snake

A common motif found in invertebrate serine proteases involved in immunity and development is the clip domain, proposed to regulate catalytic activity or protein-protein interactions within proteolytic cascades. Snake functions in a cascade that patterns the Drosophila embryo, and provides an accessible model for exploring the structural requirements for clip domain function. We tested Snake zymogens bearing charged-to-alanine mutations in the clip domain for their ability to rescue embryos lacking endogenous Snake and for their interactions by S2 cell co-transfection with upstream Gastrulation Defective and downstream Easter in the protease cascade. Of 13 single and multiple substitutions, one double mutant in a predicted protruding region exhibited a severe defect in embryonic rescue but showed only minimal defects in the co-transfection assay. We discuss implications of these and other results for potential biological roles of the Snake clip domain and for use of the in vitro assay in predicting protease behavior.     

Photoaffinity labeling combined with mass spectrometric approaches as a tool for structural proteomics

Protein chemistry, such as crosslinking and photoaffinity labeling, in combination with modern mass spectrometric techniques, can provide information regarding protein-protein interactions beyond that normally obtained from protein identification and characterization studies. While protein crosslinking can make tertiary and quaternary protein structure information available, photoaffinity labeling can be used to obtain structural data about ligand-protein interaction sites, such as oligonucleotide-protein, drug-protein and protein-protein interaction. In this article, we describe mass spectrometry-based photoaffinity labeling methodologies currently used and discuss their current limitations. We also discuss their potential as a common approach to structural proteomics for providing 3D information regarding the binding region, which ultimately will be used for molecular modeling and structure-based drug design.     

Dissecting the key recognition features of the MS2 bacteriophage translational repression complex.

The MS2 RNA operator capsid offers an unparalleled opportunity to study sequence-specific protein-protein and RNA-protein interactions in molecular detail. RNA molecules encompassing the minimal translational operator recognition elements can be soaked into crystals of RNA-free coat protein shells, allowing the RNA to access the interior of the capsids and make contact with the operator binding sites. Correct interpretation of these structural studies depends critically on functional analysis in solution to confirm that the interactions seen in the crystal are not an artefact of the unusual approach used to generate the RNA-protein complexes. Here we present a series of in vivo and in vitro functional assays, using coat proteins carrying single amino acid substitutions at residues which either interact with the operator RNA or are involved in stabilizing the conformation of the FG loop, the site of the major quasi-equivalent conformational change. Variant operator RNAs have been assayed for coat protein affinity in vitro. The results reveal the robustness of the operator-coat protein interaction and the requirement for both halves of a protein dimer to contact RNA in order to achieve tight binding. They also suggest that there may be a direct link between the conformation of the FG loop and RNA binding.     

How does protein phosphorylation regulate photosynthesis?

Phosphorylation of light-harvesting antenna proteins redirects absorbed light energy between reaction centres of photosynthetic membranes. A generally accepted explanation for this is that electrostatic forces drive the more negatively charged, phosphorylated antenna proteins between membrane domains that differ in surface charge. However, structural studies on soluble phosphoproteins indicate that phosphorylated amino acid side chains have specific effects on molecular recognition, by ligand blocking or by intramolecular interactions which alter protein structure. These studies suggest alternative mechanisms for phosphorylation in control of pairwise protein-protein interactions in biological membranes. Thus, in photosynthesis, the surface charge model is only one possible interpretation.     

Characterization of Disease-Associated Mutations in Human Transmembrane Proteins

Transmembrane protein coding genes are commonly associated with human diseases. We characterized disease causing mutations and natural polymorphisms in transmembrane proteins by mapping missense genetic variations from the UniProt database on the transmembrane protein topology listed in the Human Transmembrane Proteome database. We found characteristic differences in the spectrum of amino acid changes within transmembrane regions: in the case of disease associated mutations the non-polar to non-polar and non-polar to charged amino acid changes are equally frequent. In contrast, in the case of natural polymorphisms non-polar to charged amino acid changes are rare while non-polar to non-polar changes are common. The majority of disease associated mutations result in glycine to arginine and leucine to proline substitutions. Mutations to positively charged amino acids are more common in the center of the lipid bilayer, where they cause more severe structural and functional anomalies. Our analysis contributes to the better understanding of the effect of disease associated mutations in transmembrane proteins, which can help prioritize genetic variations in personal genomic investigations.     

Variable sensitivity to substitutions in the N-terminal heptad repeat of Mason-Pfizer monkey virus transmembrane protein

The transmembrane protein of Mason-Pfizer monkey virus contains two heptad repeats that are predicted to form amphipathic alpha-helices that mediate the conformational change necessary for membrane fusion. To analyze the relative sensitivity of the predicted hydrophobic face of the N-terminal heptad repeat to the insertion of uncharged, polar, and charged substitutions, mutations that introduced alanine, serine, or glutamic acid into positions 436, 443, 450, and 457 of the envelope protein were examined. Novel systems using Tat protein and the GHOST cell line were developed to test and quantitate the effects of the mutations on Env-mediated fusion and infectivity of the virus. While no single amino acid change at any of the positions interfered significantly with the synthesis, processing, or transport to the plasma membrane of glycoprotein complexes, 9 of the 12 nonconservative mutations in these residues completely abolished fusion activity and virus infectivity. Mutations in the central positions (443 and 450) of the heptad repeat region were the most detrimental to Env function, and even single alanine substitutions in these positions dramatically altered the fusogenicity of the protein. These results demonstrate that this N-terminal heptad repeat plays a critical role in Env-mediated membrane fusion and highlight the key function of central hydrophobic residues in this process and the sensitivity of all positions to charge substitutions.