In Silico Screening of Mutational Effects on Transmembrane Helix Dimerization: Insights from Rigid-Body Docking and Molecular Dynamics Simulations

Integrated structural modeling and super-resolution imaging resolve GPCR oligomers

Formation of G protein-coupled receptors (GPCRs) dimers and higher order oligomers represents a key mechanism in pleiotropic signaling, yet how individual protomers function within oligomers remains poorly understood. For the Class A/rhodopsin subfamily of glycoprotein hormone receptors (GpHRs), di/oligomerization has been demonstrated to play a significant role in regulating its signaling activity at a cellular and physiological level and even pathophysiologically. Here we will describe and discuss the developments in our understanding of GPCR oligomerization, in both health and disease, from the study of this unique and complex subfamily of GPCRs with light on the luteinizing hormone receptor (LHR). Focus will be put on the results of an approach relying on the combination of atomistic modeling by protein-protein docking with super-resolution imaging. The latter could resolve single LHR molecules to ~8nm resolution in functional asymmetric dimers and oligomers, using dual-color photoactivatable dyes and localization microscopy (PD-PALM). Structural modeling of functionally asymmetric LHR trimers and tetramers strongly aligned with PD-PALM-imaged spatial arrangements, identifying multiple possible helix interfaces mediating inter-protomer associations. Diverse spatial and structural assemblies mediating GPCR oligomerization may acutely fine-tune the cellular signaling profile.

Correlation between the molecular effects of mutations at the dimer interface of alanine-glyoxylate aminotransferase leading to primary hyperoxaluria type I and the cellular response to vitamin B6

Primary hyperoxaluria type I (PH1) is a rare disease caused by the deficit of liver alanine-glyoxylate aminotransferase (AGT). AGT prevents oxalate formation by converting peroxisomal glyoxylate to glycine. When the enzyme is deficient, progressive calcium oxalate stones deposit first in the urinary tract and then at the systemic level. Pyridoxal 5'-phosphate (PLP), the AGT coenzyme, exerts a chaperone role by promoting dimerization, as demonstrated by studies at protein and cellular level. Thus, variants showing a destabilized dimeric structure should, in principle, be responsive to vitamin B₆, a precursor of PLP. However, models to predict the extent of responsiveness of each variant are missing. We examined the effects of pathogenic interfacial mutations by combining bioinformatic predictions with molecular and cellular studies on selected variants (R36H, G42E, I56N, G63R, and G216R), in both their holo- (i.e., with bound PLP) and apo- (i.e., without bound PLP) form. We found that all variants displayed structural alterations mainly related to the apoform and consisting of an altered tertiary and quaternary structure. G216R also shows a strongly reduced catalytic efficiency. Moreover, all but G216R respond to vitamin B₆, as shown by their increased specific activity and expression level in a cellular disease model. A global analysis of data unraveled a possible inverse correlation between the degree of destabilization/misfolding induced by a mutation and the extent of B6 responsiveness. These results provide a first explanation of factors influencing B6 response in PH1, a model possibly valuable for other rare diseases caused by protein deficits.

Extension of a protein docking algorithm to membranes and applications to amyloid precursor protein dimerization

Novel adjustments are introduced to the docking algorithm, DOCK/PIERR, for the purpose of predicting structures of transmembrane protein complexes. Incorporating knowledge about the membrane environment is shown to significantly improve docking accuracy. The extended version of DOCK/PIERR is shown to perform comparably to other leading docking packages. This membrane version of DOCK/PIERR is applied to the prediction of coiled-coil homodimer structures of the transmembrane region of the C-terminal peptide of amyloid precursor protein (C99). Results from MD simulation of the C99 homodimer in POPC bilayer and docking are compared. Docking results are found to capture key aspects of the homodimer ensemble, including the existence of three topologically distinct conformers. Furthermore, the extended version of DOCK/PIERR is successful in capturing the effects of solvation in membrane and micelle. Specifically, DOCK/PIERR reproduces essential differences in the homodimer ensembles simulated in POPC bilayer and DPC micelle, where configurational entropy and surface curvature effects bias the handedness and topology of the homodimer ensemble.

Two Dimensional Window Exchange Umbrella Sampling for Transmembrane Helix Assembly

The method of window exchange umbrella sampling molecular dynamics (WEUSMD) with a pre-optimized parameter set was recently used to obtain the most probable conformations and the energetics of transmembrane (TM) helix assembly of a generic TM sequence. When applied to glycophorin A TM domain (GpA-TM) using the restraint potentials along the helix-helix distance, however, tight interfacial packing of GpA-TM resulted in insufficient conformational sampling at short helix-helix separation. This sampling issue is addressed by extending the WEUSMD into two dimensions with the restraint potentials along the helix-helix distance and crossing angle. The two-dimensional WEUSMD results demonstrate that the incomplete sampling in the one-dimensional WEUSMD arises from high barriers along the crossing angle between the GpA-TM helices. Together with the faster convergence in both the assembled conformations and the potential of mean force, the 2D-WEUSMD can be a general and efficient approach in computational studies of TM helix assembly.

Prediction, refinement, and persistency of transmembrane helix dimers in lipid bilayers using implicit and explicit solvent/lipid representations: microsecond molecular dynamics simulations of ErbB1/B2 and EphA1

All-atom simulations are carried out on ErbB1/B2 and EphA1 transmembrane helix dimers in lipid bilayers starting from their solution/DMPC bicelle NMR structures. Over the course of microsecond trajectories, the structures remain in close proximity to the initial configuration and satisfy the majority of experimental tertiary contact restraints. These results further validate CHARMM protein/lipid force fields and simulation protocols on Anton. Separately, dimer conformations are generated using replica exchange in conjunction with an implicit solvent and lipid representation. The implicit model requires further improvement, and this study investigates whether lengthy all-atom molecular dynamics simulations can alleviate the shortcomings of the initial conditions. The simulations correct many of the deficiencies. For example, excessive helix twisting is eliminated over a period of hundreds of nanoseconds. The helix tilt, crossing angles, and dimer contacts approximate those of the NMR-derived structure, although the detailed contact surface remains off-set for one of two helices in both systems. Hence, even microsecond simulations are not long enough for extensive helix rotations. The alternate structures can be rationalized with reference to interaction motifs and may represent still sought after receptor states that are important in ErbB1/B2 and EphA1 signaling.

Transmembrane helix assembly by window exchange umbrella sampling

A method of window exchange umbrella sampling molecular dynamics simulation is employed for transmembrane helix assembly. An analytical expression for the average acceptance probability between neighboring windows is derived and combined with the first passage time optimization method to predetermine a parameter set in an optimal range. With the parameter set, the method provides a substantially more efficient sampling of helix-helix interfaces together with the potential of mean force along the helix-helix distance of a transmembrane helix-dimer model, compared to the umbrella sampling method.

The dimerization interface of the glycoprotein Ibβ transmembrane domain corresponds to polar residues within a leucine zipper motif

Experiments with the transmembrane (TM) domains of the glycoprotein (GP) Ib-IX complex have indicated that the associations between the TM domains of these subunits play an important role in the proper assembly of the complex. As a first step toward understanding these associations, we previously found that the Ibβ TM domain dimerized strongly in Escherichia coli cell membranes and led to Ibβ TM-CYTO (cytoplasmic domain) dimerization in the SDS-PAGE assay, while neither Ibα nor IX TM-CYTO was able to dimerize. In this study, we used the TOXCAT assay to probe the Ibβ TM domain dimerization interface by Ala- and Leu-scanning mutagenesis. Our results show that this interface is based on a leucine zipper-like heptad repeat pattern of amino acids. Mutating either one of polar residues Gln129 or His139 to Leu or Ala disrupted Ibβ TM dimerization dramatically, indicating that polar residues might form part of the leucine zipper-based dimerization interface. Furthermore, these specific mutational effects in the TOXCAT assay were confirmed in the thiol-disulfide exchange and SDS-PAGE assays. The computational modeling studies further revealed that the most likely leucine zipper interface involves hydrogen bonding of Gln129 and electrostatic interaction of the His139 side chain. Correlation of computer modeling results with experimental mutagenesis studies on the Ibβ TM domain may provide insights for understanding the role of the association of TM domains on the assembly of GP Ib-IX complex.

Involvement of the recoverin C-terminal segment in recognition of the target enzyme rhodopsin kinase

NCS (neuronal Ca2+ sensor) proteins belong to a family of calmodulin-related EF-hand Ca2+-binding proteins which, in spite of a high degree of structural similarity, are able to selectively recognize and regulate individual effector enzymes in a Ca2+-dependent manner. NCS proteins vary at their C-termini, which could therefore serve as structural control elements providing specific functions such as target recognition or Ca2+ sensitivity. Recoverin, an NCS protein operating in vision, regulates the activity of rhodopsin kinase, GRK1, in a Ca2+-dependent manner. In the present study, we investigated a series of recoverin forms that were mutated at the C-terminus. Using pull-down assays, surface plasmon resonance spectroscopy and rhodopsin phosphorylation assays, we demonstrated that truncation of recoverin at the C-terminus significantly reduced the affinity of recoverin for rhodopsin kinase. Site-directed mutagenesis of single amino acids in combination with structural analysis and computational modelling of the recoverin–kinase complex provided insight into the protein–protein interface between the kinase and the C-terminus of recoverin. Based on these results we suggest that Phe3 from the N-terminal helix of rhodopsin kinase and Lys192 from the C-terminal segment of recoverin form a cation–π interaction pair which is essential for target recognition by recoverin. Taken together, the results of the present study reveal a novel rhodopsin-kinase-binding site within the C-terminal region of recoverin, and highlights its significance for target recognition and regulation.

Fast predictions of thermodynamics and kinetics of protein-protein recognition from structures: from molecular design to systems biology

The increasing call for an overall picture of the interactions between the components of a biological system that give rise to the observed function is often summarized by the expression systems biology. Both the interpretative and predictive capabilities of holistic models of biochemical systems, however, depend to a large extent on the level of physico-chemical knowledge of the individual molecular interactions making up the network. This review is focused on the structure-based quantitative characterization of protein-protein interactions, ubiquitous in any biochemical pathway. Recently developed, fast and effective computational methods are reviewed, which allow the assessment of kinetic and thermodynamic features of the association-dissociation processes of protein complexes, both in water soluble and membrane environments. The performance and the accuracy of fast and semi-empirical structure-based methods have reached comparable levels with respect to the classical and more elegant molecular simulations. Nevertheless, the broad accessibility and lower computational cost provide the former methods with the advantageous possibility to perform systems-level analyses including extensive in silico mutagenesis screenings and large-scale structural predictions of multiprotein complexes.

In silico screening of mutational effects on enzyme-proteic inhibitor affinity: a docking-based approach

Background

Molecular recognition between enzymes and proteic inhibitors is crucial for normal functioning of many biological pathways. Mutations in either the enzyme or the inhibitor protein often lead to a modulation of the binding affinity with no major alterations in the 3D structure of the complex.

Results

In this study, a rigid body docking-based approach has been successfully probed in its ability to predict the effects of single and multiple point mutations on the binding energetics in three enzyme-proteic inhibitor systems. The only requirement of the approach is an accurate structural model of the complex between the wild type forms of the interacting proteins, with the assumption that the architecture of the mutated complexes is almost the same as that of the wild type and no major conformational changes occur upon binding. The method was applied to 23 variants of the ribonuclease inhibitor-angiogenin complex, to 15 variants of the barnase-barstar complex, and to 8 variants of the bovine pancreatic trypsin inhibitor-β Trypsin system, leading to thermodynamic and kinetic estimates consistent with in vitro data. Furthermore, simulations with and without explicit water molecules at the protein-protein interface suggested that they should be included in the simulations only when their positions are well defined both in the wild type and in the mutants and they result to be relevant for the modulation of mutational effects on the association process.

Conclusion

The correlative models built in this study allow for predictions of mutational effects on the thermodynamics and kinetics of association of three substantially different systems, and represent important extensions of our computational approach to cases in which it is not possible to estimate the absolute free energies. Moreover, this study is the first example in the literature of an extensive evaluation of the correlative weights of the single components of the ZDOCK score on the thermodynamics and kinetics of binding of protein mutants compared to the native state.

Finally, the results of this study corroborate and extend a previously developed quantitative model for in silico predictions of absolute protein-protein binding affinities spanning a wide range of values, i.e. from -10 up to -21 kcal/mol.

The computational approach is simple and fast and can be used for structure-based design of protein-protein complexes and for in silico screening of mutational effects on protein-protein recognition.