Electrostatic steering at acetylcholine binding sites

Rational pore engineering reveals the relative contribution of enzymatic sites and self-assembly towards rapid ferroxidase activity and mineralization: impact of electrostatic guiding and cage-confinement in bacterioferritin

The self-assembled ferritin protein nanocage plays a pivotal role during oxidative stress, iron metabolism, and host-pathogen interaction by executing rapid iron uptake, oxidation and its safe-storage. Self-assembly creates a nanocompartment and various pores/channels for the uptake of charged substrates (Fe2+) and develops a concentration gradient across the protein shell. This phenomenon fuels rapid ferroxidase activity by an upsurge in the substrate concentration at the catalytic sites. However, it is difficult to segregate the relative contributions of the catalytic sites and self-assembly towards rapid ferroxidase/mineralization activity owing to the inherent self-assembly propensity of ferritins. In the current work, 3-fold pore electrostatics of bacterioferritin from Mycobacterium tuberculosis were rationally altered by site-directed mutagenesis to generate self-assembled (E121A and E121Q) and assembly-defective (E121K and E121F) variants. In comparison to the autoxidation of Fe2+ in buffer, the assembly-defective variants exhibited significantly faster ferroxidase/mineralization activity and O2 consumption kinetics due to their functional catalytic sites, but failed to level-up with the self-assembled variants even at 100-fold higher Fe2+ concentration. Only the self-assembled variants exhibited cooperativity in iron oxidation, maintained biomineral solubility, and protected DNA against the Fenton reaction. This report highlights the concerted effect of self-assembly and ferroxidase sites that propels the rapid Fe2+ uptake, its oxidation and biomineralization in bacterioferritin. The findings also establish the importance of electrostatic guiding and nanoconfinement offered by ferritin self-assembly towards its enzymatic activity and antioxidative properties. Moreover, this work identifies the key electrostatic interactions ("hot-spots") at the subunit contact points that control the cage/pore formation, impart cage stability and influence ferritin's natural functions. Manipulation of hot-spot residues can be further extended towards the encapsulation of cargo, for various bio-medical applications, by strategically inducing its disassembly and subsequent reassembly through adjustments in ionic strength. This would bypass the need for extreme/harsh reaction conditions and minimize the loss of cargo/protein.

Glycosylation: mechanisms, biological functions and clinical implications

Protein post-translational modification (PTM) is a covalent process that occurs in proteins during or after translation through the addition or removal of one or more functional groups, and has a profound effect on protein function. Glycosylation is one of the most common PTMs, in which polysaccharides are transferred to specific amino acid residues in proteins by glycosyltransferases. A growing body of evidence suggests that glycosylation is essential for the unfolding of various functional activities in organisms, such as playing a key role in the regulation of protein function, cell adhesion and immune escape. Aberrant glycosylation is also closely associated with the development of various diseases. Abnormal glycosylation patterns are closely linked to the emergence of various health conditions, including cancer, inflammation, autoimmune disorders, and several other diseases. However, the underlying composition and structure of the glycosylated residues have not been determined. It is imperative to fully understand the internal structure and differential expression of glycosylation, and to incorporate advanced detection technologies to keep the knowledge advancing. Investigations on the clinical applications of glycosylation focused on sensitive and promising biomarkers, development of more effective small molecule targeted drugs and emerging vaccines. These studies provide a new area for novel therapeutic strategies based on glycosylation.

On the possibility of the existence of orienting hydrodynamic steering effects in the kinetics of receptor-ligand association

In the vast majority of biologically relevant cases of receptor-ligand complex formation, the binding site of the receptor is a small part of its surface, and moreover, formation of a biologically active complex often requires a specific orientation of the ligand relative to the binding site. Before the formation of the initial form of the complex, only long-range, electrostatic and hydrodynamic interactions can act between the ligand approaching the binding site and the receptor. In this context, the question arises whether as a result of these interactions, there is a pre-orientation of the ligand towards the binding site, which to some extent would accelerate the formation of the complex. The role of electrostatic interactions in the orientation of the ligand relative to the binding site of the receptor is well documented. The analogous role of hydrodynamic interactions, although assessed as very significant by Brune and Kim (PNAS 91, 2930-2934, (1994)), is still debatable. In this article, I present the current state of knowledge on this subject and consider the possibilities of demonstrating the orienting effect of hydrodynamic interactions in the processes of receptor-ligand association, in an experimental way supported by computer simulations.

Co-translational formation of disulfides guides folding of the SARS-CoV-2 receptor binding domain

Many secreted proteins, including viral proteins, contain multiple disulfide bonds. How disulfide formation is coupled to protein folding in the cell remains poorly understood at the molecular level. Here, we combine experiment and simulation to address this question as it pertains to the SARS-CoV-2 receptor binding domain (RBD). We show that the RBD can only refold reversibly if its native disulfides are present before folding. But in their absence, the RBD spontaneously misfolds into a nonnative, molten-globule-like state that is structurally incompatible with complete disulfide formation and that is highly prone to aggregation. Thus, the RBD native structure represents a metastable state on the protein's energy landscape with reduced disulfides, indicating that nonequilibrium mechanisms are needed to ensure native disulfides form before folding. Our atomistic simulations suggest that this may be achieved via co-translational folding during RBD secretion into the endoplasmic reticulum. Namely, at intermediate translation lengths, native disulfide pairs are predicted to come together with high probability, and thus, under suitable kinetic conditions, this process may lock the protein into its native state and circumvent highly aggregation-prone nonnative intermediates. This detailed molecular picture of the RBD folding landscape may shed light on SARS-CoV-2 pathology and molecular constraints governing SARS-CoV-2 evolution.

Differences in IgG autoantibody Fab glycosylation across autoimmune diseases

BACKGROUND

Increased prevalence of autoantibody Fab glycosylation has been demonstrated for several autoimmune diseases.

OBJECTIVES

To study whether elevated Fab glycosylation is a common feature of autoimmunity, this study investigated Fab glycosylation levels on serum IgG and its subclasses for autoantibodies associated with a range of different B cell-mediated autoimmune diseases, including rheumatoid arthritis, myasthenia gravis subtypes, pemphigus vulgaris, antineutrophil cytoplasmic antibody-associated vasculitis, systemic lupus erythematosus, anti-glomerular basement membrane glomerulonephritis, thrombotic thrombocytopenic purpura, and Guillain-Barré syndrome.

METHODS

The level of Fab glycosylated IgG antibodies was assessed by lectin affinity chromatography and autoantigen-specific immunoassays.

RESULTS

In 6 of 10 autoantibody responses, in 5 of 8 diseases, the investigators found increased levels of Fab glycosylation on IgG autoantibodies that varied from 86% in rheumatoid arthritis to 26% in systemic lupus erythematosus. Elevated autoantibody Fab glycosylation was not restricted to IgG₄, which is known to be prone to Fab glycosylation, but was also present in IgG₁. When autoimmune diseases with a chronic disease course were compared with more acute autoimmune illnesses, increased Fab glycosylation was restricted to the chronic diseases. As a proxy for chronic autoantigen exposure, the investigators determined Fab glycosylation levels on antibodies to common latent herpes viruses, as well as to glycoprotein 120 in individuals who are chronically HIV-1-infected. Immunity to these viral antigens was not associated with increased Fab glycosylation levels, indicating that chronic antigen-stimulation as such does not lead to increased Fab glycosylation levels.

CONCLUSIONS

These data indicate that in chronic but not acute B cell-mediated autoimmune diseases, disease-specific autoantibodies are enriched for Fab glycans.

A conserved arginine with non-conserved function is a key determinant of agonist selectivity in α7 nicotinic ACh receptors

BACKGROUND AND PURPOSE

The α7 and α4β2* ("*" denotes possibly assembly with another subunit) nicotinic acetylcholine receptors (nAChRs) are the most abundant nAChRs in the mammalian brain. These receptors are the most targeted nAChRs in drug discovery programmes for brain disorders. However, the development of subtype-specific agonists remains challenging due to the high degree of sequence homology and conservation of function in nAChRs. We have developed C(10) variants of cytisine, a partial agonist of α4β2 nAChR that has been used for smoking cessation. The C(10) methyl analogue used in this study displays negligible affinity for α7 nAChR, while retaining high affinity for α4β2 nAChR.

EXPERIMENTAL APPROACH

The structural underpinning of the selectivity of 10-methylcytisine for α7 and α4β2 nAChRs was investigated using molecular dynamic simulations, mutagenesis and whole-cell and single-channel current recordings.

KEY RESULTS

We identified a conserved arginine in the β3 strand that exhibits a non-conserved function in nAChRs. In α4β2 nAChR, the arginine forms a salt bridge with an aspartate residue in loop B that is necessary for receptor expression, whereas in α7 nAChR, this residue is not stabilised by electrostatic interactions, making its side chain highly mobile. This lack of constrain produces steric clashes with agonists and affects the dynamics of residues involved in agonist binding and the coupling network.

CONCLUSION AND IMPLICATIONS

We conclude that the high mobility of the β3-strand arginine in the α7 nAChR influences agonist binding and possibly gating network and desensitisation. The findings have implications for rational design of subtype-selective nAChR agents.

Applications of in Silico Methods for Design and Development of Drugs Targeting Protein-Protein Interactions

BACKGROUND

Identification of Protein-Protein Interactions (PPIs) is a major challenge in modern molecular biology and biochemistry research, due to the unquestionable role of proteins in cells, biological process and pathological states. Over the past decade, the PPIs have evolved from being considered a highly challenging field of research to being investigated and examined as targets for pharmacological intervention.

OBJECTIVE

Comprehension of protein interactions is crucial to known how proteins come together to build signalling pathways, to carry out their functions, or to cause diseases, when deregulated. Multiplicity and great amount of PPIs structures offer a huge number of new and potential targets for the treatment of different diseases.

METHODS

Computational techniques are becoming predominant in PPIs studies for their effectiveness, flexibility, accuracy and cost. As a matter of fact, there are effective in silico approaches which are able to identify PPIs and PPI site. Such methods for computational target prediction have been developed through molecular descriptors and data-mining procedures.

RESULTS

In this review, we present different types of interactions between protein-protein and the application of in silico methods for design and development of drugs targeting PPIs. We described computational approaches for the identification of possible targets on protein surface and to detect of stimulator/ inhibitor molecules.

CONCLUSION

A deeper study of the most recent bioinformatics methodologies for PPIs studies is vital for a better understanding of protein complexes and for discover new potential PPI modulators in therapeutic intervention.

Lack of Dependence of the Sizes of the Mesoscopic Protein Clusters on Electrostatics

Protein-rich clusters of steady submicron size and narrow size distribution exist in protein solutions in apparent violation of the classical laws of phase equilibrium. Even though they contain a minor fraction of the total protein, evidence suggests that they may serve as essential precursors for the nucleation of ordered solids such as crystals, sickle-cell hemoglobin polymers, and amyloid fibrils. The cluster formation mechanism remains elusive. We use the highly basic protein lysozyme at nearly neutral and lower pH as a model and explore the response of the cluster population to the electrostatic forces, which govern numerous biophysical phenomena, including crystallization and fibrillization. We tune the strength of intermolecular electrostatic forces by varying the solution ionic strength I and pH and find that despite the weaker repulsion at higher I and pH, the cluster size remains constant. Cluster responses to the presence of urea and ethanol demonstrate that cluster formation is controlled by hydrophobic interactions between the peptide backbones, exposed to the solvent after partial protein unfolding that may lead to transient protein oligomers. These findings reveal that the mechanism of the mesoscopic clusters is fundamentally different from those underlying the two main classes of ordered protein solid phases, crystals and amyloid fibrils, and partial unfolding of the protein chain may play a significant role.

Recent advances in the understanding of two-step nucleation of protein crystals

The two-step mechanism of nucleation of crystals in solutions posits that the formation of crystal nuclei occurs within structures of extended lifetimes, in which the nucleating solute is at high concentration. The validity of this mechanism has been demonstrated for proteins, small-molecule organic and inorganic materials, colloids, and polymers. Due to large molecule sizes, proteins are an ideal system to study the details of this nucleation pathway, in particular the formation mechanisms of the nucleation precursors and the associated physico-chemical rules. The precursors of protein crystal nuclei are protein-rich clusters of sizes ∼100 nm that contain 10 000–100 000 molecules and occupy less than 10⁻³of the total solution volume. Here we demonstrate, using oblique illumination microscopy, the liquid nature of the clusters of the protein lysozyme and reveal their inhomogeneous structure. We test a hypothesis put forth by theory that clusters primarily consist of transient protein oligomers. For this, we explore how varying the strength of the Coulomb interaction affects the cluster characteristics. We find that the cluster’s size is insensitive to variations of pH and ionic strength. In contrast, the addition of urea, a chaotropic agent that leads to protein unfolding, strongly decreases the cluster size. Shear stress, a known protein denaturant, induced by bubbling of the solutions with an inert gas, elicits a similar response. These observations support partial protein unfolding, followed by dimerization, as the mechanism of cluster formation. The amide hydrogen–deuterium exchange, monitored by nuclear magnetic resonance, highlights that lysozyme conformational flexibility is a condition for the formation of the protein-rich clusters and facilitates the nucleation of protein crystals.

Acetylcholine promotes binding of α-conotoxin MII at α3 β2 nicotinic acetylcholine receptors

α-Conotoxin MII (α-CTxMII) is a 16-residue peptide with the sequence GCCSNPVCHLEHSNLC, containing Cys2-Cys8 and Cys3-Cys16 disulfide bonds. This peptide, isolated from the venom of the marine cone snail Conus magus, is a potent and selective antagonist of neuronal nicotinic acetylcholine receptors (nAChRs). To evaluate the impact of channel-ligand interactions on ligand-binding affinity, homology models of the heteropentameric α3β2-nAChR were constructed. The models were created in MODELLER with the aid of experimentally characterized structures of the Torpedo marmorata-nAChR (Tm-nAChR, PDB ID: 2BG9) and the Aplysia californica-acetylcholine binding protein (Ac-AChBP, PDB ID: 2BR8) as templates for the α3- and β2-subunit isoforms derived from rat neuronal nAChR primary amino acid sequences. Molecular docking calculations were performed with AutoDock to evaluate interactions of the heteropentameric nAChR homology models with the ligands acetylcholine (ACh) and α-CTxMII. The nAChR homology models described here bind ACh with binding energies commensurate with those of previously reported systems, and identify critical interactions that facilitate both ACh and α-CTxMII ligand binding. The docking calculations revealed an increased binding affinity of the α3β2-nAChR for α-CTxMII with ACh bound to the receptor, and this was confirmed through two-electrode voltage clamp experiments on oocytes from Xenopus laevis. These findings provide insights into the inhibition and mechanism of electrostatically driven antagonist properties of the α-CTxMIIs on nAChRs.

Attractive hydration forces in DNA-dendrimer interactions on the nanometer scale

The energetic contribution of attractive hydration forces arising from water ordering is an interesting but often neglected aspect of macromolecular interactions. Ordering effects of water can bring about cooperativity in many intermolecular transactions, in both the short and long range. Given its high charge density, this is of particular importance for DNA. For instance, in nanotechnology, highly charged dendrimers are used for DNA compaction and transfection. Hypothesizing that water ordering and hydration forces should be maximal for DNA complexes that show charge complementarity (positive-negative), we present here an analysis of water ordering from molecular dynamics simulations and free energy calculations of the interaction between DNA and a nanoparticle with a high positive charge density. Our results indicate not only that complexation of the dendrimer with DNA affects the local water structure but also that ordered water molecules facilitate long-range interactions between the molecules. This contributes significantly to the free energy of binding of dendrimers to DNA and extends the interaction well beyond the electrostatic range of the DNA. Such water effects are of potentially substantial importance in cases when molecules appear to recognize each other across sizable distances, or for which kinetic rates are too fast to be due to pure diffusion. Our results are in good agreement with experiments on the role of solvent in DNA condensation by multivalent cations and exemplify a microscopic realization of mean-field phenomenological theories for hydration forces between mesoscopic surfaces.

On the role of electrostatics in protein-protein interactions

The role of electrostatics in protein-protein interactions and binding is reviewed in this paper. A brief outline of the computational modeling, in the framework of continuum electrostatics, is presented and the basic electrostatic effects occurring upon the formation of the complex are discussed. The effect of the salt concentration and pH of the water phase on protein-protein binding free energy is demonstrated which indicates that the increase of the salt concentration tends to weaken the binding, an observation that is attributed to the optimization of the charge-charge interactions across the interface. It is pointed out that the pH-optimum (pH of optimal binding affinity) varies among the protein-protein complexes, and perhaps is a result of their adaptation to particular subcellular compartments. The similarities and differences between hetero- and homo-complexes are outlined and discussed with respect to the binding mode and charge complementarity.

A distal point mutation in the streptavidin-biotin complex preserves structure but diminishes binding affinity: experimental evidence of electronic polarization effects?

We have identified a distal point mutation in streptavidin that causes a 1000-fold reduction in biotin binding affinity without disrupting the equilibrium complex structure. The F130L mutation creates a small cavity occupied by a water molecule; however, all neighboring side chain positions are preserved, and protein-biotin hydrogen bonds are unperturbed. Molecular dynamics simulations reveal a reduced mobility of biotin binding residues but no observable destabilization of protein-ligand interactions. Our combined structural and computational studies suggest that the additional water molecule may affect binding affinity through an electronic polarization effect that impacts the highly cooperative hydrogen bonding network in the biotin binding pocket.

Structural answers and persistent questions about how nicotinic receptors work

The electron diffraction structure of nicotinic acetylcholine receptor (nAChR) from Torpedo marmorata and the X-ray crystallographic structure of acetylcholine binding protein (AChBP) are providing new answers to persistent questions about how nAChRs function as biophysical machines and as participants in cellular and systems physiology. New high-resolution information about nAChR structures might come from advances in crystallography and NMR, from extracellular domain nAChRs as high fidelity models, and from prokaryotic nicotinoid proteins. At the level of biophysics, structures of different nAChRs with different pharmacological profiles and kinetics will help describe how agonists and antagonists bind to orthosteric binding sites, how allosteric modulators affect function by binding outside these sites, how nAChRs control ion flow, and how large cytoplasmic domains affect function. At the level of cellular and systems physiology, structures of nAChRs will help characterize interactions with other cellular components, including lipids and trafficking and signaling proteins, and contribute to understanding the roles of nAChRs in addiction, neurodegeneration, and mental illness. Understanding nAChRs at an atomic level will be important for designing interventions for these pathologies.

A paramagnetic molecular voltmeter

We have developed a general electron paramagnetic resonance (EPR) method to measure electrostatic potential at spin labels on proteins to millivolt accuracy. Electrostatic potential is fundamental to energy-transducing proteins like myosin, because molecular energy storage and retrieval is primarily electrostatic. Quantitative analysis of protein electrostatics demands a site-specific spectroscopic method sensitive to millivolt changes. Previous electrostatic potential studies on macromolecules fell short in sensitivity, accuracy and/or specificity. Our approach uses fast-relaxing charged and neutral paramagnetic relaxation agents (PRAs) to increase nitroxide spin label relaxation rate solely through collisional spin exchange. These PRAs were calibrated in experiments on small nitroxides of known structure and charge to account for differences in their relaxation efficiency. Nitroxide longitudinal (R(1)) and transverse (R(2)) relaxation rates were separated by applying lineshape analysis to progressive saturation spectra. The ratio of measured R(1) increases for each pair of charged and neutral PRAs measures the shift in local PRA concentration due to electrostatic potential. Voltage at the spin label is then calculated using the Boltzmann equation. Measured voltages for two small charged nitroxides agree with Debye-Hückel calculations. Voltage for spin-labeled myosin fragment S1 also agrees with calculation based on the pK shift of the reacted cysteine.

Nicotinic acetylcholine receptor channel electrostatics determined by diffusion-enhanced luminescence energy transfer

The electrostatic potentials within the pore of the nicotinic acetylcholine receptor (nAChR) were determined using lanthanide-based diffusion-enhanced fluorescence energy transfer experiments. Freely diffusing Tb3+ -chelates of varying charge constituted a set of energy transfer donors to the acceptor, crystal violet, a noncompetitive antagonist of the nAChR. Energy transfer from a neutral Tb3+ -chelate to nAChR-bound crystal violet was reduced 95% relative to the energy transfer to free crystal violet. This result indicated that crystal violet was strongly shielded from solvent when bound to the nAChR. Comparison of energy transfer from positively and negatively charged chelates indicate negative electrostatic potentials of -25 mV in the channel, measured in low ionic strength, and -10 mV measured in physiological ionic strength. Debye-Hückel analyses of potentials determined at various ionic strengths were consistent with 1-2 negative charges within 8 A of the crystal violet binding site. To complement the energy transfer experiments, the influence of pH and ionic strength on the binding of [3H]phencyclidine were determined. The ionic strength dependence of binding affinity was consistent with -3.3 charges within 8 A of the binding site, according to Debye-Hückel analysis. The pH dependence of binding had an apparent pKa of 7.2, a value indicative of a potential near -170 mV if the titratable residues are constituted of aspartates and glutamates. It is concluded that long-range potentials are small and likely contribute little to selectivity or conductance whereas close interactions are more likely to contribute to electrostatic stabilization of ions and binding of noncompetitive antagonists within the channel.

Physical activity in relation to cancer of the colon and rectum in a cohort of male smokers

We examined the association between occupational and leisure physical activity and colorectal cancer in a cohort of male smokers. Among the 29,133 men aged 50-69 years in the Alpha-Tocopherol, Beta-Carotene Cancer Prevention study,152 colon and 104 rectal cancers were documented during up to 12 years of follow-up. For colon cancer, compared with sedentary workers, men in light occupational activity had a relative risk (RR) of 0.60 [95% confidence interval (CI), 0.34-1.04], whereas those in moderate/heavy activity had an RR of 0.45 (CI, 0.26-0.78; P for trend, 0.003). Subsite analysis revealed a significant association for moderate/heavy occupational activity in the distal colon (RR, 0.21; CI, 0.09-0.51) but not in the proximal colon (RR, 0.87; CI, 0.40-1.92). There was no significant association between leisure activity and colon cancer (active versus sedentary; RR, 0.82; CI, 0.59-1.13); however, the strongest inverse association was found among those most active in both work and leisure (RR, 0.33; CI, 0.16-0.71). For rectal cancer, there were risk reductions for those in light (RR, 0.71; CI, 0.36-1.37) and moderate/heavy occupational activity (RR, 0.50; CI, 0.26-0.97; P for trend, 0.04), and no association for leisure activity. These data provide evidence for a protective role of physical activity in colon and rectal cancer.