A glycophorin A-like framework for the dimerization of photosynthetic core complexes

A Rigorous Framework for Calculating Protein-Protein Binding Affinities in Membranes

Calculating the binding free energy of integral transmembrane (TM) proteins is crucial for understanding the mechanisms by which they recognize one another and reversibly associate. The glycophorin A (GpA) homodimer, composed of two α-helical segments, has long served as a model system for studying TM protein reversible association. The present work establishes a methodological framework for calculating the binding affinity of the GpA homodimer in the heterogeneous environment of a membrane. Our investigation carefully considered a variety of protocols, including the appropriate choice of the force field, rigorous standardization reflecting the experimental conditions, sampling algorithm, anisotropic environment, and collective variables, to accurately describe GpA dimerization via molecular dynamics-based approaches. Specifically, two strategies were explored: (i) an unrestrained potential mean force (PMF) calculation, which merely enhances sampling along the separation of the two binding partners without any restraint, and (ii) a so-called "geometrical route", whereby the α-helices are progressively separated with imposed restraints on their orientational, positional, and conformational degrees of freedom to accelerate convergence. Our simulations reveal that the simplified, unrestrained PMF approach is inadequate for the description of GpA dimerization. Instead, the geometrical route, tailored specifically to GpA in a membrane environment, yields excellent agreement with experimental data within a reasonable computational time. A dimerization free energy of -10.7 kcal/mol is obtained, in fairly good agreement with available experimental data. The geometrical route further helps elucidate how environmental forces drive association before helical interactions stabilize it. Our simulations also brought to light a distinct, long-lived spatial arrangement that potentially serves as an intermediate state during dimer formation. The methodological advances in the generalized geometrical route provide a powerful tool for accurate and efficient binding-affinity calculations of intricate TM protein complexes in inhomogeneous environments.

Variation in supramolecular organisation of the photosynthetic membrane of Rhodobacter sphaeroides induced by alteration of PufX

In purple bacteria of the genus Rhodobacter (Rba.), an LH1 antenna complex surrounds the photochemical reaction centre (RC) with a PufX protein preventing the LH1 complex from completely encircling the RC. In membranes of Rba. sphaeroides, RC-LH1 complexes associate as dimers which in turn assemble into longer range ordered arrays. The present work uses linear dichroism (LD) and dark-minus-light difference LD (ΔLD) to probe the organisation of genetically altered RC-LH1 complexes in intact membranes. The data support previous proposals that Rba. capsulatus, and Rba. sphaeroides heterologously expressing the PufX protein from Rba. capsulatus, produce monomeric core complexes in membranes that lack long-range order. Similarly, Rba. sphaeroides with a point mutation in the Gly 51 residue of PufX, which is located on the membrane-periplasm interface, assembles mainly non-ordered RC-LH1 complexes that are most likely monomeric. All the Rba. sphaeroides membranes in their ΔLD spectra exhibited a spectral fingerprint of small degree of organisation implying the possibility of ordering influence of LH1, and leading to an important conclusion that PufX itself has no influence on ordering RC-LH1 complexes, as long-range order appears to be induced only through its role of configuring RC-LH1 complexes into dimers.

Multidimensional umbrella sampling and replica-exchange molecular dynamics simulations for structure prediction of transmembrane helix dimers

Structural information of a transmembrane (TM) helix dimer is useful in understanding molecular mechanisms of important biological phenomena such as signal transduction across the cell membrane. Here, we describe an umbrella sampling (US) scheme for predicting the structure of a TM helix dimer in implicit membrane using the interhelical crossing angle and the TM-TM relative rotation angles as the reaction coordinates. This scheme conducts an efficient conformational search on TM-TM contact interfaces, and its robustness is tested by predicting the structures of glycophorin A (GpA) and receptor tyrosine kinase EphA1 (EphA1) TM dimers. The nuclear magnetic resonance (NMR) structures of both proteins correspond to the global free-energy minimum states in their free-energy landscapes. In addition, using the landscape of GpA as a reference, we also examine the protocols of temperature replica-exchange molecular dynamics (REMD) simulations for structure prediction of TM helix dimers in implicit membrane. A wide temperature range in REMD simulations, for example, 250-1000 K, is required to efficiently obtain a free-energy landscape consistent with the US simulations. The interhelical crossing angle and the TM-TM relative rotation angles can be used as reaction coordinates in multidimensional US and be good measures for conformational sampling of REMD simulations.

Cross-species investigation of the functions of the Rhodobacter PufX polypeptide and the composition of the RC-LH1 core complex

In well-characterised species of the Rhodobacter (Rba.) genus of purple photosynthetic bacteria it is known that the photochemical reaction centre (RC) is intimately-associated with an encircling LH1 antenna pigment protein, and this LH1 antenna is prevented from completely surrounding the RC by a single copy of the PufX protein. In Rba. veldkampii only monomeric RC-LH1 complexes are assembled in the photosynthetic membrane, whereas in Rba. sphaeroides and Rba. blasticus a dimeric form is also assembled in which two RCs are surrounded by an S-shaped LH1 antenna. The present work established that dimeric RC-LH1 complexes can also be isolated from Rba. azotoformans and Rba. changlensis, but not from Rba. capsulatus or Rba. vinaykumarii. The compositions of the monomers and dimers isolated from these four species of Rhodobacter were similar to those of the well-characterised RC-LH1 complexes present in Rba. sphaeroides. Pigment proteins were also isolated from strains of Rba. sphaeroides expressing chimeric RC-LH1 complexes. Replacement of either the Rba. sphaeroides LH1 antenna or PufX with its counterpart from Rba. capsulatus led to a loss of the dimeric form of the RC-LH1 complex, but the monomeric form had a largely unaltered composition, even in strains in which the expression level of LH1 relative to the RC was reduced. The chimeric RC-LH1 complexes were also functional, supporting bacterial growth under photosynthetic conditions. The findings help to tease apart the different functions of PufX in different species of Rhodobacter, and a specific protein structural arrangement that allows PufX to fulfil these three functions is proposed.

Oligomerization state of photosynthetic core complexes is correlated with the dimerization affinity of a transmembrane helix

In the Rhodobacter (Rba.) species of photosynthetic purple bacteria, a single transmembrane α-helix, PufX, is found within the core complex, an essential photosynthetic macromolecular assembly that performs the absorption and the initial processing of light energy. Despite its structural simplicity, many unresolved questions surround PufX, the most important of which is its location within the photosynthetic core complex. One proposed placement of PufX is at the center of a core complex dimer, where two PufX helices associate in the membrane and form a homodimer. Inability for PufX of certain Rba. species to form a homodimer is thought to lead to monomeric core complexes. In the present study, we employ a combination of computational and experimental techniques to test the hypothesized homodimerization of PufX. We carry out a systematic investigation to measure the dimerization affinity of PufX from four Rba. species, Rba. blasticus , Rba. capsulatus , Rba. sphaeroides , and Rba. veldkampii , using a molecular dynamics-based free-energy method, as well as experimental TOXCAT assays. We found that the four PufX helices have substantially different dimerization affinities. Both computational and experimental techniques demonstrate that species with dimeric core complexes have PufX that can potentially form a homodimer, whereas the one species with monomeric core complexes has a PufX with little to no dimerization propensity. Our analysis of the helix-helix interface revealed a number of positions that may be important for PufX dimerization and the formation of a hydrogen-bond network between these GxxxG-containing helices. Our results suggest that the different oligomerization states of core complexes in various Rba. species can be attributed, among other factors, to the different propensity of its PufX helix to homodimerize.

Molecular dynamics simulations of the interactions of kinin peptides with an anionic POPG bilayer

We have performed molecular dynamics simulations of peptide hormone bradykinin (BK) and its fragment des-Arg9-BK in the presence of an anionic lipid bilayer, with an aim toward delineating the mechanism of action related to their bioactivity. Starting from the initial aqueous environment, both of the peptides are quickly adsorbed and stabilized on the cell surface. Whereas BK exhibits a stronger interaction with the membrane and prefers to stay on the interface, des-Arg9-BK, with the loss of C-terminal Arg, penetrates further. The heterogeneous lipid-water interface induces β-turn-like structure in the otherwise inherently flexible peptides. In the membrane-bound state, we observed C-terminal β-turn formation in BK, whereas for des-Arg9-BK, with the deletion of Arg9, turn formation occurred in the middle of the peptide. The basic Arg residues anchor the peptide to the bilayer by strong electrostatic interactions with charged lipid headgroups. Simulations with different starting orientations of the peptides with respect to the bilayer surface lead to the same observations, namely, the relative positioning of the peptides on the membrane surface, deeper penetration of the des-Arg9-BK, and the formation of turn structures. The lipid headgroups adjacent to the bound peptides become substantially tilted, causing bilayer thinning near the peptide contact region and increase the degree of disorder in nearby lipids. Again, because of hydrogen bonding with the peptide, the neighboring lipid's polar heads exhibit considerably reduced flexibility. Corroborating findings from earlier experiments, our results provide important information about how the lipid environment promotes peptide orientation/conformation and how the peptide adapts to the environment.

Formation of salt bridges mediates internal dimerization of myosin VI medial tail domain

The unconventional motor protein, myosin VI, is known to dimerize upon cargo binding to its C-terminal end. It has been shown that one of its tail domains, called the medial tail domain, is a dimerization region. The domain contains an unusual pattern of alternating charged residues and a few hydrophobic residues. To reveal the unknown dimerization mechanism of the medial tail domain, we employed molecular dynamics and single-molecule experimental techniques. Both techniques suggest that the formation of electrostatic-based interhelical salt bridges between oppositely charged residues is a key dimerization factor. For the dimerization to occur, the two identical helices within the dimer do not bind in a symmetric fashion, but rather with an offset of about one helical repeat. Calculations of the dimer-dissociation energy find the contribution of hydrophobic residues to the dimerization process to be minor; they also find that the asymmetric homodimer state is energetically favorable over a state of separate helices.

Experimental evidence that the membrane-spanning helix of PufX adopts a bent conformation that facilitates dimerisation of the Rhodobacter sphaeroides RC-LH1 complex through N-terminal interactions

The PufX polypeptide is an integral component of some photosynthetic bacterial reaction center-light harvesting 1 (RC-LH1) core complexes. Many aspects of the structure of PufX are unresolved, including the conformation of its long membrane-spanning helix and whether C-terminal processing occurs. In the present report, NMR data recorded on the Rhodobacter sphaeroides PufX in a detergent micelle confirmed previous conclusions derived from equivalent data obtained in organic solvent, that the α-helix of PufX adopts a bent conformation that would allow the entire helix to reside in the membrane interior or at its surface. In support of this, it was found through the use of site-directed mutagenesis that increasing the size of a conserved glycine on the inside of the bend in the helix was not tolerated. Possible consequences of this bent helical structure were explored using a series of N-terminal deletions. The N-terminal sequence ADKTIFNDHLN on the cytoplasmic face of the membrane was found to be critical for the formation of dimers of the RC-LH1 complex. It was further shown that the C-terminus of PufX is processed at an early stage in the development of the photosynthetic membrane. A model in which two bent PufX polypeptides stabilise a dimeric RC-LH1 complex is presented, and it is proposed that the N-terminus of PufX from one half of the dimer engages in electrostatic interactions with charged residues on the cytoplasmic surface of the LH1α and β polypeptides on the other half of the dimer.

Applications of the molecular dynamics flexible fitting method

In recent years, cryo-electron microscopy (cryo-EM) has established itself as a key method in structural biology, permitting the structural characterization of large biomolecular complexes in various functional states. The data obtained through single-particle cryo-EM has recently seen a leap in resolution thanks to landmark advances in experimental and computational techniques, resulting in sub-nanometer resolution structures being obtained routinely. The remaining gap between these data and revealing the mechanisms of molecular function can be closed through hybrid modeling tools that incorporate known atomic structures into the cryo-EM data. One such tool, molecular dynamics flexible fitting (MDFF), uses molecular dynamics simulations to combine structures from X-ray crystallography with cryo-EM density maps to derive atomic models of large biomolecular complexes. The structures furnished by MDFF can be used subsequently in computational investigations aimed at revealing the dynamics of the complexes under study. In the present work, recent applications of MDFF are presented, including the interpretation of cryo-EM data of the ribosome at different stages of translation and the structure of a membrane-curvature-inducing photosynthetic complex.

Native architecture of the photosynthetic membrane from Rhodobacter veldkampii

The photosynthetic membrane in purple bacteria contains several pigment-protein complexes that assure light capture and establishment of the chemiosmotic gradient. The bioenergetic tasks of the photosynthetic membrane require the strong interaction between these various complexes. In the present work, we acquired the first images of the native outer membrane architecture and the supramolecular organization of the photosynthetic apparatus in vesicular chromatophores of Rhodobacter (Rb.) veldkampii. Mixed with LH2 (light-harvesting complex 2) rings, the PufX-containing LH1-RC (light-harvesting complex 1--reaction center) core complexes appear as C-shaped monomers, with random orientations in the photosynthetic membrane. Within the LH1 fence surrounding the RC, a remarkable gap that is probably occupied (or partially occupied) by PufX is visualized. Sequence alignment revealed that one specific region in PufX may be essential for PufX-induced core dimerization. In this region of ten amino acids in length all Rhodobacter species had five conserved amino acids, with the exception of Rb. veldkampii. Our findings provide direct evidence that the presence of PufX in Rb. veldkampii does not directly govern the dimerization of LH1-RC core complexes in the native membrane. It is indicated, furthermore, that the high membrane curvature of Rb. veldkampii chromatophores (Rb. veldkampii features equally small vesicular chromatophores alike Rb. sphaeroides) is not due to membrane bending induced by dimeric RC-LH1-PufX cores, as it has been proposed in Rb. sphaeroides.

Energy Transfer Dynamics in an RC-LH1-PufX Tubular Photosynthetic Membrane

Light absorption and the subsequent transfer of excitation energy are the first two steps of the photosynthetic process, carried out by protein-bound pigments, mainly bacteriochlorophylls (BChls), in photosynthetic bacteria. BChls are anchored in light-harvesting (LH) complexes, such as light-harvesting complex I (LH1), which directly associates with the reaction center (RC), forming the RC-LH1 core complex. In Rhodobacter sphaeroides, RC-LH1 core complexes contain an additional protein, PufX, and assemble into dimeric RC-LH1-PufX core complexes. In the absence of light-harvesting complexes II, the former complexes can aggregate into a helically ordered tubular photosynthetic membrane. We examined the excitation transfer dynamics in a single RC-LH1-PufX core complex dimer using the hierarchical equations of motion for dissipative quantum dynamics that accurately, yet computationally costly, treat the coupling between BChls and their protein environment. A widely employed description, generalized Förster theory, was also used to calculate the transfer rates of the same excitonic system in order to verify the accuracy of this computationally cheap method. Additionally, in light of the structural uncertainties in the Rhodobacter sphaeroides RC-LH1-PufX core complex, geometrical alterations were introduced in the BChl organization. It is shown that the energy transfer dynamics is not affected by the considered changes in the BChl organization, and that generalized Förster theory provides accurate transfer rates. An all-atom model for a tubular photosynthetic membrane is then constructed on the basis of electron microscopy data, and the overall energy transfer properties of this membrane are computed.

Dimerisation of the Rhodobacter sphaeroides RC-LH1 photosynthetic complex is not facilitated by a GxxxG motif in the PufX polypeptide

In purple photosynthetic bacteria the initial steps of light energy transduction take place in an RC-LH1 complex formed by the photochemical reaction centre (RC) and the LH1 light harvesting pigment-protein. In Rhodobacter sphaeroides, the RC-LH1 complex assembles in a dimeric form in which two RCs are surrounded by an S-shaped LH1 antenna. There is currently debate over the detailed architecture of this dimeric RC-LH1 complex, with particular emphasis on the location and precise function of a minor polypeptide component termed PufX. It has been hypothesised that the membrane-spanning helical region of PufX contains a GxxxG dimerisation motif that facilitates the formation of a dimer of PufX at the interface of the RC-LH1 dimer, and more specifically that the formation of this PufX dimer seeds assembly of the remaining RC-LH1 dimer (J. Busselez et al., 2007). In the present work this hypothesis was tested by site directed mutagenesis of the glycine residues proposed to form the GxxxG motif. Mutation of these glycines to leucine did not decrease the propensity of the RC-LH1 complex to assemble in a dimeric form, as would be expected from experimental studies of the effect of mutation on GxxxG motifs in other membrane proteins. Indeed increased yields of dimer were seen in two of the glycine-to-leucine mutants constructed. It is concluded that the PufX from Rhodobacter sphaeroides does not contain a genuine GxxxG helix dimerisation motif.

Self-assembly of photosynthetic membranes

Bacterial photosynthetic membranes, also known as chromatophores, are tightly packed with integral membrane proteins that work together to carry out photosynthesis. Chromatophores display a wide range of cellular morphologies; spherical, tubular, and lamellar chromatophores have all been observed in different bacterial species, or with different protein constituents. Through recent computational modeling and simulation, it has been demonstrated that the light-harvesting complexes abundant in chromatophores induce local membrane curvature via multiple mechanisms. These protein complexes assemble to generate a global curvature and sculpt the chromatophores into various cellular-scale architectures.