Allosteric Signal Research Articles

Integrated Biophysical Modeling of the SARS-CoV-2 Spike Protein Binding and Allosteric Interactions with Antibodies.

Structural and biochemical studies of the severe acute respiratory syndrome (SARS)-CoV-2 spike glycoproteins and complexes with highly potent antibodies have revealed multiple conformation-dependent epitopes highlighting conformational plasticity of spike proteins and capacity for eliciting specific binding and broad neutralization responses. In this study, we used coevolutionary analysis, molecular simulations, and perturbation-based hierarchical network modeling of the SARS-CoV-2 spike protein complexes with a panel of antibodies targeting distinct epitopes to explore molecular mechanisms underlying binding-induced modulation of dynamics and allosteric signaling in the spike proteins. Through coevolutionary analysis of the SARS-CoV-2 spike proteins, we identified highly coevolving hotspots and functional clusters that enable a functional cross-talk between distant allosteric regions in the SARS-CoV-2 spike complexes with antibodies. Coarse-grained and all-atom molecular dynamics simulations combined with mutational sensitivity mapping and perturbation-based profiling of the SARS-CoV-2 receptor-binding domain (RBD) complexes with CR3022 and CB6 antibodies enabled a detailed validation of the proposed approach and an extensive quantitative comparison with the experimental structural and deep mutagenesis scanning data. By combining in silico mutational scanning, perturbation-based modeling, and network analysis of the SARS-CoV-2 spike trimer complexes with H014, S309, S2M11, and S2E12 antibodies, we demonstrated that antibodies can incur specific and functionally relevant changes by modulating allosteric propensities and collective dynamics of the SARS-CoV-2 spike proteins. The results provide a novel insight into regulatory mechanisms of SARS-CoV-2 S proteins showing that antibody-escaping mutations can preferentially target structurally adaptable energy hotspots and allosteric effector centers that control functional movements and allosteric communication in the complexes.

The Speed of Allosteric Signaling Within a Single-Domain Protein.

While much is known about different allosteric regulation mechanisms, the nature of the allosteric signal and the time scale on which it propagates remains elusive. The PDZ3 domain from postsynaptic density-95 protein is a small protein domain with a terminal third α-helix, i.e., the α3-helix, which is known to be allosterically active. By cross-linking the allosteric helix with an azobenzene moiety, we obtained a photocontrollable PDZ3 variant. Photoswitching triggers its allosteric transition, resulting in a change in binding affinity of a peptide to the remote binding pocket. Using time-resolved infrared and UV/vis spectroscopy, we follow the allosteric signal transduction and reconstruct the timeline in which the allosteric signal propagates through the protein within 200 ns.

Exploring the Allosteric Territory of Protein Function.

While the pervasiveness of allostery in proteins is commonly accepted, we further show the generic nature of allosteric mechanisms by analyzing here transmembrane ion-channel viroporin 3a and RNA-dependent RNA polymerase (RdRp) from SARS-CoV-2 along with metabolic enzymes isocitrate dehydrogenase 1 (IDH1) and fumarate hydratase (FH) implicated in cancers. Using the previously developed structure-based statistical mechanical model of allostery (SBSMMA), we share our experience in analyzing the allosteric signaling, predicting latent allosteric sites, inducing and tuning targeted allosteric response, and exploring the allosteric effects of mutations. This, yet incomplete list of phenomenology, forms a complex and unique allosteric territory of protein function, which should be thoroughly explored. We propose a generic computational framework, which not only allows one to obtain a comprehensive allosteric control over proteins but also provides an opportunity to approach the fragment-based design of allosteric effectors and drug candidates. The advantages of allosteric drugs over traditional orthosteric compounds, complemented by the emerging role of the allosteric effects of mutations in the expansion of the cancer mutational landscape and in the increased mutability of viral proteins, leave no choice besides further extensive studies of allosteric mechanisms and their biomedical implications.

Light-Induced Change of Arginine Conformation Modulates the Rate of Adenosine Triphosphate to Cyclic Adenosine Monophosphate Conversion in the Optogenetic System Containing Photoactivated Adenylyl Cyclase.

We report the first computational characterization of an optogenetic system composed of two photosensing BLUF (blue light sensor using flavin adenine dinucleotide) domains and two catalytic adenylyl cyclase (AC) domains. Conversion of adenosine triphosphate (ATP) to the reaction products, cyclic adenosine monophosphate (cAMP) and pyrophosphate (PPi), catalyzed by ACs initiated by excitation in photosensing domains has emerged in the focus of modern optogenetic applications because of the request in photoregulated enzymes that modulate cellular concentrations of signaling messengers. The photoactivated AC from the soil bacterium Beggiatoa sp. (bPAC) is an important model showing a considerable increase in the ATP to cAMP conversion rate in the catalytic domain after the illumination of the BLUF domain. The 1 μs classical molecular dynamics simulations reveal that the activation of the BLUF domain leading to tautomerization of Gln49 in the chromophore-binding pocket results in switching of the position of the side chain of Arg278 in the active site of AC. Allosteric signal transmission pathways between Gln49 from BLUF and Arg278 from AC were revealed by the dynamical network analysis. The Gibbs energy profiles of the ATP → cAMP + PPi reaction computed using QM(DFT(ωB97X-D3/6-31G**))/MM(CHARMM) molecular dynamics simulations for both Arg278 conformations in AC clarify the reaction mechanism. In the light-activated system, the corresponding arginine conformation stabilizes the pentacoordinated phosphorus of the α-phosphate group in the transition state, thus lowering the activation energy. Simulations of the bPAC system with the Tyr7Phe replacement in the BLUF demonstrate occurrence of both arginine conformations in an equal ratio, explaining the experimentally observed intermediate catalytic activity of the bPAC-Y7F variant as compared with the dark and light states of the wild-type bPAC.

New Mechanistic Insights on Carbon Nanotubes' Nanotoxicity Using Isolated Submitochondrial Particles, Molecular Docking, and Nano-QSTR Approaches.

Simple SummaryCarbon nanotubes are revolutionary materials with applications in a lot of different areas. However, there is a rising concern regarding unlikely toxicity effects these materials may trigger. Due to this, the main aim of this paper is to develop a comprehensive approach to study toxicity effect of carbon nanotubes on the mitochondria F0F1-ATPase. We have employed a combination of experimental and computational study. In so doing, we have combined in vitro inhibition responses in submitochondrial particles with docking elastic network models, fractal surface analysis, and Nano-quantitative structure toxicity relationship models (Nano-QSTR models). Results show that this method may be used for the fast prediction of the nanotoxicity induced by single walled carbon nanotubes (SWCNT), avoiding time- and money-consuming techniques, and may open new avenues toward to the better understanding and prediction of new nanotoxicity mechanisms.Single-walled carbon nanotubes can induce mitochondrial F0F1-ATPase nanotoxicity through inhibition. To completely characterize the mechanistic effect triggering the toxicity, we have developed a new approach based on the combination of experimental and computational study, since the use of only one or few techniques may not fully describe the phenomena. To this end, the in vitro inhibition responses in submitochondrial particles (SMP) was combined with docking, elastic network models, fractal surface analysis, and Nano-QSTR models. In vitro studies suggest that inhibition responses in SMP of F0F1-ATPase enzyme were strongly dependent on the concentration assay (from 3 to 5 µg/mL) for both pristine and COOH single-walled carbon nanotubes types (SWCNT). Besides, both SWCNTs show an interaction inhibition pattern mimicking the oligomycin A (the specific mitochondria F0F1-ATPase inhibitor blocking the c-ring F0 subunit). Performed docking studies denote the best crystallography binding pose obtained for the docking complexes based on the free energy of binding (FEB) fit well with the in vitro evidence from the thermodynamics point of view, following an affinity order such as: FEB (oligomycin A/F0-ATPase complex) = −9.8 kcal/mol > FEB (SWCNT-COOH/F0-ATPase complex) = −6.8 kcal/mol ~ FEB (SWCNT-pristine complex) = −5.9 kcal/mol, with predominance of van der Waals hydrophobic nano-interactions with key F0-ATPase binding site residues (Phe 55 and Phe 64). Elastic network models and fractal surface analysis were performed to study conformational perturbations induced by SWCNT. Our results suggest that interaction may be triggering abnormal allosteric responses and signals propagation in the inter-residue network, which could affect the substrate recognition ligand geometrical specificity of the F0F1-ATPase enzyme in order (SWCNT-pristine > SWCNT-COOH). In addition, Nano-QSTR models have been developed to predict toxicity induced by both SWCNTs, using results of in vitro and docking studies. Results show that this method may be used for the fast prediction of the nanotoxicity induced by SWCNT, avoiding time- and money-consuming techniques. Overall, the obtained results may open new avenues toward to the better understanding and prediction of new nanotoxicity mechanisms, rational drug design-based nanotechnology, and potential biomedical application in precision nanomedicine.

Structural Communication between the E. coli Chaperones DnaK and Hsp90.

The 70 kDa and 90 kDa heat shock proteins Hsp70 and Hsp90 are two abundant and highly conserved ATP-dependent molecular chaperones that participate in the maintenance of cellular homeostasis. In Escherichia coli, Hsp90 (Hsp90) and Hsp70 (DnaK) directly interact and collaborate in protein remodeling. Previous work has produced a model of the direct interaction of both chaperones. The locations of the residues involved have been confirmed and the model has been validated. In this study, we investigate the allosteric communication between Hsp90 and DnaK and how the chaperones couple their conformational cycles. Using elastic network models (ENM), normal mode analysis (NMA), and a structural perturbation method (SPM) of asymmetric and symmetric DnaK-Hsp90, we extract biologically relevant vibrations and identify residues involved in allosteric signaling. When one DnaK is bound, the dominant normal modes favor biological motions that orient a substrate protein bound to DnaK within the substrate/client binding site of Hsp90 and release the substrate from the DnaK substrate binding domain. The presence of one DnaK molecule stabilizes the entire Hsp90 protomer to which it is bound. Conversely, the symmetric model of DnaK binding results in steric clashes of DnaK molecules and suggests that the Hsp90 and DnaK chaperone cycles operate independently. Together, this data supports an asymmetric binding of DnaK to Hsp90.

Study on functional sites in human multidrug resistance protein 1 (hMRP1).

Human multidrug resistance protein 1 (hMRP1) is an important member of the ATP-binding cassette (ABC) transporter superfamily. It can extrude a variety of anticancer drugs and physiological organic anions across the plasma membrane, which is activated by substrate binding, and is accompanied by large-scale cooperative movements between different domains. Currently, it remains unclear completely about how the specific interactions between hMRP1 and its substrate are and which critical residues are responsible for allosteric signal transduction. To the end, we first construct an inward-facing state of hMRP1 using homology modeling method, and then dock substrate proinflammatory agent leukotriene C4 (LTC4) to hMRP1 pocket. The result manifests LTC4 interacts with two parts of hMRP1 pocket, namely the positively charged pocket (P pocket) and hydrophobic pocket (H pocket), similar to its binding mode with bMRP1 (bovine MRP1). Additionally, we use the Gaussian network model (GNM)-based thermodynamic method proposed by us to identify the key residues whose perturbations markedly alter their binding free energy. Here the conventional GNM is improved with covalent/non-covalent interactions and secondary structure information considered (denoted as sscGNM). In the result, sscGNM improves the flexibility prediction, especially for the nucleotide binding domains with rich kinds of secondary structures. The 46 key residue clusters located in different subdomains are identified which are highly consistent with experimental observations. Furtherly, we explore the long-range cooperation within the transporter. This study is helpful for strengthening the understanding of the work mechanism in ABC exporters and can provide important information to scientists in drug design studies.

Information Propagation in Time through Allosteric Signaling

Naively, one expects the information communicated by an enzyme downstream within a signaling network, in which the enzyme is embedded, to grow monotonically with the enzyme's rate of product formation. However, here we observe that this does not necessarily hold true for allosterically regulated enzymes, often observed in signaling networks. In particular, we find that the mutual information between the catalytic sites of an allosterically regulated enzyme and a receiver protein downstream in the signaling pathway depends on the transition kinetics between the different allosterically regulated states of the enzyme and their respective rates of product formation.

Unravelling Allosteric Cross-Talk between Co-Activator Peptide and Ligand Binding Site in Glucocorticoid Receptor

Glucocorticoid receptor (GR) is a nuclear receptor that controls critical biological processes by regulating the transcription of specific genes. There is a known interconnection between the ligand and co-regulator binding sites of the GR ligand binding domain that is crucial for the functional differentiation of the downstream gene transcription. However, the molecular mechanisms underlying such an allosteric control remain elusive. The integration of molecular dynamics simulations, bioinformatic analysis and biophysical measurements led us to detect a network of dynamically interconnected evolutionary conserved residues that mediates the allosteric signal propagation between the co-regulator docking site and the ligand binding pocket. Co-regulator peptides alter the functionality of this network, explaining why the allosteric response depends on the peptide sequence. These findings provide useful insights for the comprehension of GR allosteric regulation that should help the design of novel drugs.

Allosteric signaling in C-linker and cyclic nucleotide-binding domain of HCN2 channels

Opening of hyperpolarization-activated cyclic nucleotide-modulated (HCN) channels is controlled by membrane hyperpolarization and binding of cyclic nucleotides to the tetrameric cyclic nucleotide-binding domain (CNBD), attached to the C-linker (CL) disk. Confocal patch-clamp fluorometry revealed pronounced cooperativity of ligand binding among protomers. However, by which pathways allosteric signal transmission occurs remained elusive. Here, we investigate how changes in the structural dynamics of the CL-CNBD of mouse HCN2 upon cAMP binding relate to inter- and intrasubunit signal transmission. Applying a rigidity-theory-based approach, we identify two intersubunit and one intrasubunit pathways that differ in allosteric coupling strength between cAMP-binding sites or toward the CL. These predictions agree with results from electrophysiological and patch-clamp fluorometry experiments. Our results map out distinct routes within the CL-CNBD that modulate different cAMP-binding responses in HCN2 channels. They signify that functionally relevant submodules may exist within and across structurally discernable subunits in HCN channels.

Dynamics of the HD regulatory subdomain of PARP-1; substrate access and allostery in PARP activation and inhibition.

PARP-1 is a key early responder to DNA damage in eukaryotic cells. An allosteric mechanism links initial sensing of DNA single-strand breaks by PARP-1’s F1 and F2 domains via a process of further domain assembly to activation of the catalytic domain (CAT); synthesis and attachment of poly(ADP-ribose) (PAR) chains to protein sidechains then signals for assembly of DNA repair components. A key component in transmission of the allosteric signal is the HD subdomain of CAT, which alone bridges between the assembled DNA-binding domains and the active site in the ART subdomain of CAT. Here we present a study of isolated CAT domain from human PARP-1, using NMR-based dynamics experiments to analyse WT apo-protein as well as a set of inhibitor complexes (with veliparib, olaparib, talazoparib and EB-47) and point mutants (L713F, L765A and L765F), together with new crystal structures of the free CAT domain and inhibitor complexes. Variations in both dynamics and structures amongst these species point to a model for full-length PARP-1 activation where first DNA binding and then substrate interaction successively destabilise the folded structure of the HD subdomain to the point where its steric blockade of the active site is released and PAR synthesis can proceed.

Chemical systems biology reveals mechanisms of glucocorticoid receptor signaling.

Glucocorticoids display remarkable anti-inflammatory activity, but their use is limited by on-target adverse effects including insulin resistance and skeletal muscle atrophy. We used a chemical systems biology approach, ligand class analysis, to examine ligands designed to modulate glucocorticoid receptor activity through distinct structural mechanisms. These ligands displayed diverse activity profiles, providing the variance required to identify target genes and coregulator interactions that were highly predictive of their effects on myocyte glucose disposal and protein balance. Their anti-inflammatory effects were linked to glucose disposal but not muscle atrophy. This approach also predicted selective modulation in vivo, identifying compounds that were muscle-sparing or anabolic for protein balance and mitochondrial potential. Ligand class analysis defined the mechanistic links between the ligand-receptor interface and ligand-driven physiological outcomes, a general approach that can be applied to any ligand-regulated allosteric signaling system.

Dynamic Network Modeling of Allosteric Interactions and Communication Pathways in the SARS-CoV-2 Spike Trimer Mutants: Differential Modulation of Conformational Landscapes and Signal Transmission via Cascades of Regulatory Switches.

The rapidly growing body of structural and biochemical studies of the SARS-CoV-2 spike glycoprotein has revealed a variety of distinct functional states with radically different arrangements of the receptor-binding domain, highlighting a remarkable function-driven conformational plasticity and adaptability of the spike proteins. In this study, we examined molecular mechanisms underlying conformational and dynamic changes in the SARS-CoV-2 spike mutant trimers through the lens of dynamic analysis of allosteric interaction networks and atomistic modeling of signal transmission. Using an integrated approach that combined coarse-grained molecular simulations, protein stability analysis, and perturbation-based modeling of residue interaction networks, we examined how mutations in the regulatory regions of the SARS-CoV-2 spike protein can differentially affect dynamics and allosteric signaling in distinct functional states. The results of this study revealed key functional regions and regulatory centers that govern collective dynamics, allosteric interactions, and control signal transmission in the SARS-CoV-2 spike proteins. We found that the experimentally confirmed regulatory hotspots that dictate dynamic switching between conformational states of the SARS-CoV-2 spike protein correspond to the key hinge sites and global mediating centers of the allosteric interaction networks. The results of this study provide a novel insight into allosteric regulatory mechanisms of SARS-CoV-2 spike proteins showing that mutations at the key regulatory positions can differentially modulate distribution of states and determine topography of signal communication pathways operating through state-specific cascades of control switch points. This analysis provides a plausible strategy for allosteric probing of the conformational equilibrium and therapeutic intervention by targeting specific hotspots of allosteric interactions and communications in the SARS-CoV-2 spike proteins.

Atomic-scale insights into allosteric inhibition and evolutional rescue mechanism of Streptococcus thermophilus Cas9 by the anti-CRISPR protein AcrIIA6

CRISPR-Cas systems are prokaryotic adaptive immunity against invading phages and plasmids. Phages have evolved diverse protein inhibitors of CRISPR-Cas systems, called anti-CRISPR (Acr) proteins, to neutralize this CRISPR machinery. In response, bacteria have co-evolved Cas variants to escape phage’s anti-CRISPR strategies, called anti-anti-CRISPR systems. Here we explore the anti-CRISPR allosteric inhibition and anti-anti-CRISPR rescue mechanisms between Streptococcus thermophilus Cas9 (St1Cas9) and the anti-CRISPR protein AcrIIA6 at the atomic level, by generating mutants of key residues in St1Cas9. Extensive unbiased molecular dynamics simulations show that the functional motions of St1Cas9 in the presence of AcrIIA6 differ substantially from those of St1Cas9 alone. AcrIIA6 binding triggers a shift of St1Cas9 conformational ensemble towards a less catalytically competent state; this state significantly compromises protospacer adjacent motif (PAM) recognition and nuclease activity by altering interdependently conformational dynamics and allosteric signals among nuclease domains, PAM-interacting (PI) regions, and AcrIIA6 binding motifs. Via in vitro DNA cleavage assays, we further elucidate the rescue mechanism of efficiently escaping AcrIIA6 inhibition harboring St1Cas9 triple mutations (G993K/K1008M/K1010E) in the PI domain and identify the evolutionary landscape of such mutational escape within species. Our results provide mechanistic insights into Acr proteins as natural brakes for the CRISPR-Cas systems and a promising potential for the design of allosteric Acr peptidomimetics.

In silico study of intrinsic dynamics of full-length apo-ACE2 and RBD-ACE2 complex

The key step for SARS-CoV-2 to infect human cells is the membrane fusion triggered by the binding of the viral extracellular Spike protein to the human extracellular receptor, the angiotensin-converting enzyme 2 (ACE2). Although the Cryo-electron microscopy (Cryo-EM) uncovered the static atomic details of ACE2 homodimers, there is still a lack of research on the kinetic and thermodynamic properties of these full-length structures. This information is helpful to understand and interpret the role of ACE2 in the cell entry of SARS-CoV-2. In order to obtain this information, we performed microsecond-scale conventional and accelerated molecular dynamics (MD) simulations of full-length all-atomic systems of the RBD-ACE2 complex, the normal and torsional conformations of the apo-ACE2 homodimer. The comparative analysis of these systems showed that there were differences in their allosteric signal pathways and motion trends. These results may be helpful to further explore the cell entry mechanism of SARS-CoV-2. Moreover, the binding free energy and hydrogen bond distribution analysis of RBD-ACE2 binding interface provided the binding motifs that may be critical to allosteric signal transmission and RBD binding. These multi-conformational binding motifs can be used as targets or templates for the inhibitor design of the cell entry of SARS-CoV-2.

Revealing the Dynamic Allosteric Changes Required for Formation of the Cysteine Synthase Complex by Hydrogen-Deuterium Exchange MS

CysE and CysK, the last two enzymes of the cysteine biosynthetic pathway, engage in a bienzyme complex, cysteine synthase, with yet incompletely characterized three-dimensional structure and regulatory function. Being absent in mammals, the two enzymes and their complex are attractive targets for antibacterial drugs. We have used hydrogen/deuterium exchange MS to unveil how complex formation affects the conformational dynamics of CysK and CysE. Our results support a model where CysE is present in solution as a dimer of trimers, and each trimer can bind one CysK homodimer. When CysK binds to one CysE monomer, intratrimer allosteric communication ensures conformational and dynamic symmetry within the trimer. Furthermore, a long-range allosteric signal propagates through CysE to induce stabilization of the interface between the two CysE trimers, preparing the second trimer for binding the second CysK with a nonrandom orientation. These results provide new molecular insights into the allosteric formation of the cysteine synthase complex and could help guide antibacterial drug design.

Structural analyses of the PKA RIIβ holoenzyme containing the oncogenic DnaJB1-PKAc fusion protein reveal protomer asymmetry and fusion-induced allosteric perturbations in fibrolamellar hepatocellular carcinoma.

When the J-domain of the heat shock protein DnaJB1 is fused to the catalytic (C) subunit of cAMP-dependent protein kinase (PKA), replacing exon 1, this fusion protein, J-C subunit (J-C), becomes the driver of fibrolamellar hepatocellular carcinoma (FL-HCC). Here, we use cryo-electron microscopy (cryo-EM) to characterize J-C bound to RIIβ, the major PKA regulatory (R) subunit in liver, thus reporting the first cryo-EM structure of any PKA holoenzyme. We report several differences in both structure and dynamics that could not be captured by the conventional crystallography approaches used to obtain prior structures. Most striking is the asymmetry caused by the absence of the second cyclic nucleotide binding (CNB) domain and the J-domain in one of the RIIβ:J-C protomers. Using molecular dynamics (MD) simulations, we discovered that this asymmetry is already present in the wild-type (WT) RIIβ2C2 but had been masked in the previous crystal structure. This asymmetry may link to the intrinsic allosteric regulation of all PKA holoenzymes and could also explain why most disease mutations in PKA regulatory subunits are dominant negative. The cryo-EM structure, combined with small-angle X-ray scattering (SAXS), also allowed us to predict the general position of the Dimerization/Docking (D/D) domain, which is essential for localization and interacting with membrane-anchored A-Kinase-Anchoring Proteins (AKAPs). This position provides a multivalent mechanism for interaction of the RIIβ holoenzyme with membranes and would be perturbed in the oncogenic fusion protein. The J-domain also alters several biochemical properties of the RIIβ holoenzyme: It is easier to activate with cAMP, and the cooperativity is reduced. These results provide new insights into how the finely tuned allosteric PKA signaling network is disrupted by the oncogenic J-C subunit, ultimately leading to the development of FL-HCC.

Allosteric Pathways Originating at Cysteine Residues in Regulators of G-Protein Signaling Proteins

Regulators of G-protein signaling (RGS) proteins play a central role in modulating signaling via G-protein coupled receptors (GPCRs). Specifically, RGS proteins bind to activated Gα subunits in G-proteins, accelerate the GTP hydrolysis, and thereby rapidly dampen GPCR signaling. Therefore, covalent molecules targeting conserved cysteine residues among RGS proteins have emerged as potential candidates to inhibit the RGS/Gα protein-protein interaction and enhance GPCR signaling. Although these inhibitors bind to conserved cysteine residues among RGS proteins, we have previously suggested [J. Am. Chem. Soc. 2018;140:3454–3460] that their potencies and specificities are related to differential protein dynamics among RGS proteins. Using data from all-atom molecular dynamics simulations, we reveal these differences in dynamics of RGS proteins by partitioning the protein structural space into a network of communities that allow allosteric signals to propagate along unique pathways originating at inhibitor binding sites and terminating at the RGS/Gα protein-protein interface.

Long-Range Charge Reorganization as an Allosteric Control Signal in Proteins.

A new mechanism of allostery in proteins, based on charge rather than structure, is reported. We demonstrate that dynamic redistribution of charge within a protein can control its function and affect its interaction with a binding partner. In particular, the association of an antibody with its target protein antigen is studied. Dynamic charge shifting within the antibody during its interaction with the antigen is enabled by its binding to a metallic surface that serves as a source for electrons. The kinetics of antibody–antigen association are enhanced when charge redistribution is allowed, even though charge injection happens at a position far from the antigen binding site. This observation points to charge-reorganization allostery, which should be operative in addition or parallel to other mechanisms of allostery, and may explain some current observations on protein interactions.

Exploring the conformational binding mechanism of fibrinogen induced by interactions with penicillin β-lactam antibiotic drugs

Herein, we present an integrated computational and experimental study to tackle the interactions between recognized β-lactam antibiotics (cloxacillin and dicloxacillin) with the fibrinogen blood plasma protein. For this purpose, molecular docking simulation with elastic network based on collective low-frequency normal modes and perturbation response scanning maps were proposed to evaluate the conformational binding mechanism of fibrinogen under the unbound and bound states with the cited β-lactam antibiotics. Aiming to theoretically explore the hidden biochemical mechanisms and structural attributes leading failures in therapy success with β-lactam antibiotics. The computational results pointing that despite these conformational differences, both antibiotics exhibit very similar affinity-based free energies of binding as FEB (cloxacillin/E-region) = − 8.7 kcal/mol and FEB (dicloxacillin/E-region) = −7.7 kcal/mol. We theoretically suggest that the semi-synthetic incorporation of an additional halogen CL-atom in the dicloxacillin, respect to cloxacillin molecule, and its relative docking-pose orientation in the fibrinogen E-region could significantly reduce the appearance of potential fibrinolytic off-target effects usually associated to parenterally administered β-lactam antibiotics. Besides, the performed interactions diagrams revealed that the most relevant antibiotic binding interactions with the fibrinogen E-region (pocket 1) are mainly based on hydrophobic (C···C)-backbone-side-chain non-covalent interactions, acceptor/donor interactions with critical regulatory E-region residues SER50:Q > SER50:N associated to allosteric modulation based long-distance-based perturbations (dicloxacillin > > cloxacillin) in the E-region (Q-chain > N-chain) with remarkable conformational rigidification by decreasing the intrinsic collectivity, and leading different pattern of perturbations as allosteric signal propagation in the intrinsic conformational dynamics under bound state from both β-lactam antibiotics. An experimental validation was carried out by using calorimetric (ITC and DSC) and spectroscopic (Raman and fluorescence) methods. These methods corroborated the computational results, adding quantitative information to explain the binding process. Finally, the obtained results open new perspectives for the “de novo rational drug-design” of new derivatives of β-lactam antibiotics with high pharmacodynamic selectivity/specificity to avoid side-effects toward to achieve optimal benefit/risk rates beyond the antibiotic drug resistance phenomena, favoring the implementation of rigorous criteria for a more personalized antibiotic therapy.