Poisson-Boltzmann (PB) models are a fast and common tool for studying electrostatic processes in proteins, particularly their ionization equilibrium (protonation and/or reduction), often yielding quite good results when compared with more detailed models. Yet, they are conceptually very simple and necessarily approximate, their empirical character being most evident when it comes to the choice of the dielectric constant assigned to the protein region. The present study analyzes several factors affecting the ability of PB-based methods to model protein ionization equilibrium. We give particular attention to a suggestion made by Warshel and co-workers (e.g., Sham et al. J. Phys. Chem. B 1997, 101, 4458) of using different protein dielectric constants for computing the individual (site) and the pairwise (site-site) terms of the ionization free energies. Our prediction of pK(a) values for several proteins indicates that no advantage is obtained by such a procedure, even for sites that are buried and/or display large pK(a) shifts relative to the solution values. In particular, the present methodology gives the best predictions using a dielectric constant around 20, for shifted/buried and nonshifted/exposed sites alike. The similarities and differences between the PB model and Warshel's PDLD/S model are discussed, as well as the reasons behind their apparently discrepant results. The present PB model is shown to predict also good reduction potentials in redox proteins.

The genome of Desulfovibrio vulgaris Hildenborough (DvH) encodes for six hydrogenases (Hases), making it an interesting organism to study the role of these proteins in sulphate respiration. In this work we address the role of the [NiFeSe] Hase, found to be the major Hase associated with the cytoplasmic membrane. The purified enzyme displays interesting catalytic properties, such as a very high H(2) production activity, which is dependent on the presence of phospholipids or detergent, and resistance to oxygen inactivation since it is isolated aerobically in a Ni(II) oxidation state. Evidence was obtained that the [NiFeSe] Hase is post-translationally modified to include a hydrophobic group bound to the N-terminal, which is responsible for its membrane association. Cleavage of this group originates a soluble, less active form of the enzyme. Sequence analysis shows that [NiFeSe] Hases from Desulfovibrionacae form a separate family from the [NiFe] enzymes of these organisms, and are more closely related to [NiFe] Hases from more distant bacterial species that have a medial [4Fe4S](2+/1+) cluster, but not a selenocysteine. The interaction of the [NiFeSe] Hase with periplasmic cytochromes was investigated and is similar to the [NiFe](1) Hase, with the Type I cytochrome c (3) as the preferred electron acceptor. A model of the DvH [NiFeSe] Hase was generated based on the structure of the Desulfomicrobium baculatum enzyme. The structures of the two [NiFeSe] Hases are compared with the structures of [NiFe] Hases, to evaluate the consensual structural differences between the two families. Several conserved residues close to the redox centres were identified, which may be relevant to the higher activity displayed by [NiFeSe] Hases.

ATP-Binding Cassette (ABC) transporters are ubiquitous membrane proteins that use energy from ATP binding or/and hydrolysis to actively transport allocrites across membranes. In this study, we identify ATP-hydrolysis induced conformational changes in a complete ABC exporter (Sav1866) from Staphylococcus aureaus, using molecular dynamics (MD) simulations. By performing MD simulations on the ATP and ADP+IP bound states, we identify the conformational consequences of hydrolysis, showing that the major rearrangements are not restricted to the NBDs, but extend to the transmembrane domains (TMDs) external regions. For the first time, to our knowledge, we see, within the context of a complete transporter, NBD dimer opening in the ADP+IP state in contrast with all ATP-bound states. This opening results from the dissociation of the ABC signature motif from the nucleotide. In addition, in both states, we observe the opening of a gate entrance in the intracellular loop region leading to the exposure of the TMDs internal cavity to the cytoplasm. To see if this opening was large enough to allow allocrite transport, the adiabatic energy profile for doxorubicin passage was determined. For both states, this profile, although an approximation, is overall downhill from the cytoplasmatic to the extracellular side, and the local energy barriers along the TMDs are relatively small, evidencing the exporter nature of Sav1866. The major difference between states is an energy barrier located in the cytoplasmic gate region, which becomes reduced upon hydrolysis, suggesting that allocrite passage is facilitated, and evidencing a possible molecular mechanism for the active transport in these proteins.

Changeux et al. recently suggested that the SARS-CoV-2 spike protein may interact with nicotinic acetylcholine receptors (nAChRs), and that such interactions may be involved in pathology and infectivity. This hypothesis is based on the fact that the SARS-CoV-2 spike protein contains a sequence motif similar to known nAChR antagonists. Here, we use molecular simulations of validated atomically detailed structures of nAChRs, and of the spike, to investigate the possible binding of the Y674-R685 region of the spike to nAChRs. We examine the binding of the Y674-R685 loop to three nAChRs, namely the human α4β2 and α7 subtypes and the muscle-like αβγδ receptor from Tetronarce californica . Our results predict that Y674-R685 has affinity for nAChRs. The region of the spike responsible for binding contains a PRRA motif, a four-residue insertion not found in other SARS-like coronaviruses. The conformational behaviour of the bound Y674-R685 is highly dependent on the receptor subtype: it adopts extended conformations in the α4β2 and α7 complexes, but is more compact when bound to the muscle-like receptor. In the α4β2 and αβγδ complexes, the interaction of Y674-R685 with the receptors forces the loop C region to adopt an open conformation, similar to other known nAChR antagonists. In contrast, in the α7 complex, Y674-R685 penetrates deeply into the binding pocket where it forms interactions with the residues lining the aromatic box, namely with TrpB, TyrC1 and TyrC2. Estimates of binding energy suggest that Y674-R685 forms stable complexes with all three nAChR subtypes. Analyses of simulations of the glycosylated spike show that the Y674-R685 region is accessible for binding. We suggest a potential binding orientation of the spike protein with nAChRs, in which they are in a non-parallel arrangement to one another.

ABC transporters are a large and important family of membrane proteins involved in substrate transport across the membrane. The transported substrates are quite diverse, ranging from monatomic ions to large biomolecules. Consequently, some ABC transporters are involved in biomedically relevant situations, from genetic diseases to multidrug resistance. The most conserved domains in ABC transporters are the nucleotide binding domains (NBDs), which form a dimer responsible for the binding and hydrolysis of ATP, concomitantly with substrate translocation. To elucidate how ATP hydrolysis structurally affects the NBD dimer, and consequently the transporter, we performed a molecular dynamics study on the NBD dimer of the HlyB ABC exporter. We have observed a change in the contact surface between the monomers after hydrolysis, even though we have not seen dimer opening in any of the five 100 ns simulations. We have also identified specific regions that respond to ATP hydrolysis, in particular the X-loop motif of ABC exporters, which has been shown to be in contact with the coupling helices of the transmembrane domains (TMDs). We propose that this motif is an important part of the NBD-TMD communication in ABC exporters. Through nonequilibrium analysis, we have also identified gradual conformational changes within a short time scale after ATP hydrolysis.

Despite the rapid advances in the study of ABC transporters, many fundamental questions linked to ATP binding/hydrolysis and its relation to the transport cycle remain unanswered. In particular, it is still neither clear nor consensual how the ATP energy is used by the nucleotide binding domains (NBDs) to produce mechanical work and drive the substrate translocation. The major conformational changes in the NBDs following ATP hydrolysis during the transport cycle and the role played by the conserved family motifs in harnessing the energy associated with nucleotide hydrolysis are yet unknown. Additionally, the way energy is transmitted from the catalytic to the membrane domains, in order to drive substrate translocation, is also a fundamental question that remains unanswered. Due to the high structure similarities of the NBD architecture throughout the whole ABC family, it is likely that the mechanism of ATP binding, hydrolysis, and communication with the transmembrane domains is similar in all family members, independently of the nature of the transported substrate. In this work, we focused our attention on the consequences of ATP hydrolysis in the NBDs, especially on the structural changes that occur during this process. For that, we use molecular dynamics simulation techniques taking as a starting point the X-ray structure of the MJ0796 dimer from Methanococcus jannaschii. Several potential intermediate states of the ATP hydrolytic cycle are investigated, each consisting of different combinations of nucleotide-bound forms. The results obtained allowed us to identify the conformational rearrangements induced by hydrolysis on the catalytic subunits, as well as the residues involved in this reorganization. The major changes are localized at specific regions of the protein, namely, involving segments 11-19 and 93-124. Additionally, our results together with the knowledge of complete ABC transporter X-ray structures suggest a possible NBD:TMD signal transmission interface.

Highlights d Simulations identify conformational changes in the a4b2 nAChR caused by nicotine d Nicotine induces structural changes in the extracellular and transmembrane domains d Nonequilibrium simulations reveal a binding pocket to ion channel communication d The results suggest a general signal transduction mechanism for the Cys-loop family Authors

Efficient access to C(10) of (À)-cytisine via C-H activation provides access to enantiomerically pure nicotinic acetylcholine receptor ligands that target the highaffinity nicotine a4b2 subtype with enhanced selectivity. These C(10) cytisine variants retain a partial agonist profile at the a4b2 subtype but, critically, display negligible activity at the a7 receptor subtype. Using computational methods, Gallagher and colleagues link receptor selectivity to key protein residues associated with, as well as beyond, the immediate ligand binding site.