The role of dynamics in allosteric regulation

Physical mechanisms of the Sec machinery operation

The Sec complex, composed of a motor protein SecA and a channel SecYEG, is an ATP-driven molecular machine for the transport of proteins across the plasma membrane in bacteria. Today, there is a consensus about a general "rough" model of the complex activation and operation, which, however, lacks understanding of the physical mechanisms behind it. Molecular dynamics simulations were employed to address a way of allosteric activation, conformational transition of SecYEG from the closed to the open state, and driving forces of protein transport. We found that binding of SecA (in the ATP-bound state) and the protein signal sequence leads to a transmembrane helix rearrangment that weakens contacts inside the hydrophobic core of SecYEG and provides a driving force for plug opening. The conformational transitions are enabled by a delicate interplay between hydrophobic forces on one side and PEES (proton motive force, external - due to binding with the translocation partners - entropic, and solvent-induced) on the other side. In the open state, SecYEG still provides a barrier for bulky residues that contributes to the driving forces of transport. Other important contributions come from SecA and the membrane potential acting in different stages of protein transport to guarantee a nearly constant driving force. Given that the different forces act on different types of residues, the suggested mechanisms taken together provide a directional motion for any substrate, thereby maximizing the efficiency of the Sec machinery.

Structure and Dynamics of the CCCH-Type Tandem Zinc Finger Domain of POS-1 and Implications for RNA Binding Specificity

CCCH-type tandem zinc finger (TZF) motifs are found in many RNA-binding proteins involved in regulating mRNA stability, translation, and splicing. In Caenorhabditis elegans, several RNA-binding proteins that regulate embryonic development and cell fate determination contain CCCH TZF domains, including POS-1. Previous biochemical studies have shown that despite high levels of sequence conservation, POS-1 recognizes a broader set of RNA sequences compared to the human homologue tristetraprolin. However, the molecular basis of these differences remains unknown. In this study, we refined the consensus RNA sequence and determined the differing binding specificities of the two zinc fingers of POS-1. We also determined the solution structure and characterized the internal dynamics of the TZF domain of POS-1. From the structure, we identified unique features that define the RNA binding specificity of POS-1. We also observed that the TZF domain of POS-1 is in equilibrium between interconverting conformations. Transitions between these conformations require internal motions involving many residues with correlated dynamics in each ZF. We propose that the correlated dynamics are necessary to allow allosteric communication between the nucleotide-binding pockets observed in the N-terminal ZF. Our study shows that both the structure and conformational plasticity of POS-1 are important in ensuring recognition of its RNA binding targets.

Utilizing Molecular Dynamics Simulations, Machine Learning, Cryo-EM, and NMR Spectroscopy to Predict and Validate Protein Dynamics

Protein dynamics play a crucial role in biological function, encompassing motions ranging from atomic vibrations to large-scale conformational changes. Recent advancements in experimental techniques, computational methods, and artificial intelligence have revolutionized our understanding of protein dynamics. Nuclear magnetic resonance spectroscopy provides atomic-resolution insights, while molecular dynamics simulations offer detailed trajectories of protein motions. Computational methods applied to X-ray crystallography and cryo-electron microscopy (cryo-EM) have enabled the exploration of protein dynamics, capturing conformational ensembles that were previously unattainable. The integration of machine learning, exemplified by AlphaFold2, has accelerated structure prediction and dynamics analysis. These approaches have revealed the importance of protein dynamics in allosteric regulation, enzyme catalysis, and intrinsically disordered proteins. The shift towards ensemble representations of protein structures and the application of single-molecule techniques have further enhanced our ability to capture the dynamic nature of proteins. Understanding protein dynamics is essential for elucidating biological mechanisms, designing drugs, and developing novel biocatalysts, marking a significant paradigm shift in structural biology and drug discovery.

De novo design of allosterically switchable protein assemblies

Allosteric modulation of protein function, wherein the binding of an effector to a protein triggers conformational changes at distant functional sites, plays a central part in the control of metabolism and cell signalling1-3. There has been considerable interest in designing allosteric systems, both to gain insight into the mechanisms underlying such 'action at a distance' modulation and to create synthetic proteins whose functions can be regulated by effectors4-7. However, emulating the subtle conformational changes distributed across many residues, characteristic of natural allosteric proteins, is a significant challenge8,9. Here, inspired by the classic Monod-Wyman-Changeux model of cooperativity10, we investigate the de novo design of allostery through rigid-body coupling of peptide-switchable hinge modules11 to protein interfaces12 that direct the formation of alternative oligomeric states. We find that this approach can be used to generate a wide variety of allosterically switchable systems, including cyclic rings that incorporate or eject subunits in response to peptide binding and dihedral cages that undergo effector-induced disassembly. Size-exclusion chromatography, mass photometry13 and electron microscopy reveal that these designed allosteric protein assemblies closely resemble the design models in both the presence and absence of peptide effectors and can have ligand-binding cooperativity comparable to classic natural systems such as haemoglobin14. Our results indicate that allostery can arise from global coupling of the energetics of protein substructures without optimized side-chain-side-chain allosteric communication pathways and provide a roadmap for generating allosterically triggerable delivery systems, protein nanomachines and cellular feedback control circuitry.

Dynamic negative allosteric effect: regulation of catalysis via multicomponent rotor speed

Nanorotor R1 (420 kHz) was assembled from five components utilizing three orthogonal interactions. Post-modification at the distal position generated the advanced six component rotor R2 (45 kHz). The decrease in R2 speed leads to the inhibition of a three-component reaction by reducing catalyst release.

Membrane mimetic-dependence of GPCR energy landscapes

We leveraged variable-temperature 19F-NMR spectroscopy to compare the conformational equilibria of the human A2A adenosine receptor (A2AAR), a class A G protein-coupled receptor (GPCR), across a range of temperatures ranging from lower temperatures typically employed in 19F-NMR experiments to physiological temperature. A2AAR complexes with partial agonists and full agonists showed large increases in the population of a fully active conformation with increasing temperature. NMR data measured at physiological temperature were more in line with functional data. This was pronounced for complexes with partial agonists, where the population of active A2AAR was nearly undetectable at lower temperature but became evident at physiological temperature. Temperature-dependent behavior of complexes with either full or partial agonists exhibited a pronounced sensitivity to the specific membrane mimetic employed. Cellular signaling experiments correlated with the temperature-dependent conformational equilibria of A2AAR in lipid nanodiscs but not in some detergents, underscoring the importance of the membrane environment in studies of GPCR function.

Structural basis for the allosteric behaviour and substrate specificity of Lactococcus lactis Prolidase

Prolidase (EC 3.4.13.9) is an enzyme that specifically hydrolyzes Xaa-Pro dipeptides into free amino acids. We previously studied kinetic behaviours and solved the crystal structure of wild-type (WT) Lactococcus lactis prolidase (Llprol), showing that this homodimeric enzyme has unique characteristics: allosteric behaviour and substrate inhibition. In this study, we focused on solving the crystal structures of three Llprol mutants (D36S, H38S, and R293S) which behave differently in v-S plots. The D36S and R293S Llprol mutants do not show allosteric behaviour, and the Llprol mutant H38S has allosteric behaviour comparable to the WT enzyme (Hill constant 1.52 and 1.58, respectively). The crystal structures of Llprol variants suggest that the active site of Llprol formed with amino acid residues from both monomers, i.e., located in an interfacial area of dimer. The comparison between the structure models of Llprol indicated that the two monomers in the dimers of Llprol variants have different relative positions among Llprol variants. They showed different interatomic distances between the amino acid residues bridging the two monomers and varied sizes of the solvent-accessible interface areas in each Llprol variant. These observations indicated that Llprol could adapt to different conformational states with distinctive substrate affinities. It is strongly speculated that the domain movements required for productive substrate binding are restrained in allosteric Llprol (WT and H38S). At low substrate concentrations, only one out of the two active sites at the dimer interface could accept substrate; as a result, the asymmetrical activated dimer leads to allosteric behaviour.

Perspectives on Computational Enzyme Modeling: From Mechanisms to Design and Drug Development

Understanding enzyme mechanisms is essential for unraveling the complex molecular machinery of life. In this review, we survey the field of computational enzymology, highlighting key principles governing enzyme mechanisms and discussing ongoing challenges and promising advances. Over the years, computer simulations have become indispensable in the study of enzyme mechanisms, with the integration of experimental and computational exploration now established as a holistic approach to gain deep insights into enzymatic catalysis. Numerous studies have demonstrated the power of computer simulations in characterizing reaction pathways, transition states, substrate selectivity, product distribution, and dynamic conformational changes for various enzymes. Nevertheless, significant challenges remain in investigating the mechanisms of complex multistep reactions, large-scale conformational changes, and allosteric regulation. Beyond mechanistic studies, computational enzyme modeling has emerged as an essential tool for computer-aided enzyme design and the rational discovery of covalent drugs for targeted therapies. Overall, enzyme design/engineering and covalent drug development can greatly benefit from our understanding of the detailed mechanisms of enzymes, such as protein dynamics, entropy contributions, and allostery, as revealed by computational studies. Such a convergence of different research approaches is expected to continue, creating synergies in enzyme research. This review, by outlining the ever-expanding field of enzyme research, aims to provide guidance for future research directions and facilitate new developments in this important and evolving field.

T-cell virtuosity in ''knowing thyself"

Major Histocompatibility Complex (MHC) I and II and the αβ T-cell antigen receptor (TCRαβ) govern fundamental traits of adaptive immunity. They form a membrane-borne ligand-receptor system weighing host proteome integrity to detect contamination by nonself proteins. MHC-I and -II exhibit the "MHC-fold", which is able to bind a large assortment of short peptides as proxies for self and nonself proteins. The ensuing varying surfaces are mandatory ligands for Ig-like TCRαβ highly mutable binding sites. Conserved molecular signatures guide TCRαβ ligand binding sites to focus on the MHC-fold (MHC-restriction) while leaving many opportunities for its most hypervariable determinants to contact the peptide. This riveting molecular strategy affords many options for binding energy compatible with specific recognition and signalling aimed to eradicated microbial pathogens and cancer cells. While the molecular foundations of αβ T-cell adaptive immunity are largely understood, uncertainty persists on how peptide-MHC binding induces the TCRαβ signals that instruct cell-fate decisions. Solving this mystery is another milestone for understanding αβ T-cells' self/nonself discrimination. Recent developments revealing the innermost links between TCRαβ structural dynamics and signalling modality should help dissipate this long-sought-after enigma.

Dynamic correlations: exact and approximate methods for mutual information

MOTIVATION

Proteins are dynamic entities that undergo conformational changes critical for their functions. Understanding the communication pathways and information transfer within proteins is crucial for elucidating allosteric interactions in their mechanisms. This study utilizes mutual information (MI) analysis to probe dynamic allostery. Using two cases, Ubiquitin and PLpro, we have evaluated the accuracy and limitations of different approximations including the exact anisotropic and isotropic models, multivariate Gaussian model, isotropic Gaussian model, and the Gaussian Network Model (GNM) in revealing allosteric interactions.

RESULTS

Our findings emphasize the required trajectory length for capturing accurate mutual information profiles. Long molecular dynamics trajectories, 1 ms for Ubiquitin and 100 µs for PLpro are used as benchmarks, assuming they represent the ground truth. Trajectory lengths of approximately 5 µs for Ubiquitin and 1 µs for PLpro marked the onset of convergence, while the multivariate Gaussian model accurately captured mutual information with trajectories of 5 ns for Ubiquitin and 350 ns for PLpro. However, the isotropic Gaussian model is less successful in representing the anisotropic nature of protein dynamics, particularly in the case of PLpro, highlighting its limitations. The GNM, however, provides reasonable approximations of long-range information exchange as a minimalist network model based on a single crystal structure. Overall, the optimum trajectory lengths for effective Gaussian approximations of long-time dynamic behavior depend on the inherent dynamics within the protein's topology. The GNM, by showcasing dynamics across relatively diverse time scales, can be used either as a standalone method or to gauge the adequacy of MD simulation lengths.

AVAILABILITY AND IMPLEMENTATION

Mutual information codes are available at https://github.com/kemaldemirtas/prc-MI.git.

QMrebind: incorporating quantum mechanical force field reparameterization at the ligand binding site for improved drug-target kinetics through milestoning simulations

Understanding the interaction of ligands with biomolecules is an integral component of drug discovery and development. Challenges for computing thermodynamic and kinetic quantities for pharmaceutically relevant receptor-ligand complexes include the size and flexibility of the ligands, large-scale conformational rearrangements of the receptor, accurate force field parameters, simulation efficiency, and sufficient sampling associated with rare events. Our recently developed multiscale milestoning simulation approach, SEEKR2 (Simulation Enabled Estimation of Kinetic Rates v.2), has demonstrated success in predicting unbinding (koff) kinetics by employing molecular dynamics (MD) simulations in regions closer to the binding site. The MD region is further subdivided into smaller Voronoi tessellations to improve the simulation efficiency and parallelization. To date, all MD simulations are run using general molecular mechanics (MM) force fields. The accuracy of calculations can be further improved by incorporating quantum mechanical (QM) methods into generating system-specific force fields through reparameterizing ligand partial charges in the bound state. The force field reparameterization process modifies the potential energy landscape of the bimolecular complex, enabling a more accurate representation of the intermolecular interactions and polarization effects at the bound state. We present QMrebind (Quantum Mechanical force field reparameterization at the receptor-ligand binding site), an ORCA-based software that facilitates reparameterizing the potential energy function within the phase space representing the bound state in a receptor-ligand complex. With SEEKR2 koff estimates and experimentally determined kinetic rates, we compare and interpret the receptor-ligand unbinding kinetics obtained using the newly reparameterized force fields for model host-guest systems and HSP90-inhibitor complexes. This method provides an opportunity to achieve higher accuracy in predicting receptor-ligand koff rate constants.

Membrane Mimetic-Dependence of GPCR Energy Landscapes

Protein function strongly depends on temperature, which is related to temperature-dependent changes in the equilibria of protein conformational states. We leveraged variable-temperature 19F-NMR spectroscopy to interrogate the temperature dependence of the conformational landscape of the human A2A adenosine receptor (A2AAR), a class A GPCR. Temperature-induced changes in the conformational equilibria of A2AAR in lipid nanodiscs were markedly dependent on the efficacy of bound drugs. While antagonist complexes displayed only modest changes as the temperature rose, both full and partial agonist complexes exhibited substantial increases in the active state population. Importantly, the temperature-dependent response of complexes with both full and partial agonists exhibited a pronounced sensitivity to the specific membrane mimetic employed. In striking contrast to observations within lipid nanodiscs, in detergent micelles the active state population exhibited different behavior for A2AAR complexes with both full and partial agonists. This underscores the importance of the protein environment in understanding the thermodynamics of GPCR activation.

Hard-Soft Cooperativity of an Allosteric Tweezer

Drawing inspiration from allosteric proteins, a zigzag-shaped π-conjugation was structurally engineered into a tweezer-like ionophore having multiple disparate binding sites. When a soft metal ion binds to the central tridentate ligand motif, the rigid backbone folds, bringing two macrocyclic arms into close proximity. Stabilized by a coordinating anion, the tweezer-like conformation of the resulting metalloligand recruits a hard cation to form a sandwich-like complex with a remarkably enhanced binding affinity and selectivity.

IMPDH1-associated autosomal dominant retinitis pigmentosa: natural history of novel variant Lys314Gln and a comprehensive literature search

BACKGROUND

Inosine monophosphate dehydrogenase (IMPDH) is a key regulatory enzyme in the de novo synthesis of the purine base guanine. Mutations in the inosine monophosphate dehydrogenase 1 gene (IMPDH1) are causative for RP10 autosomal dominant retinitis pigmentosa (adRP). This study reports a novel variant in a family with IMPDH1-associated retinopathy. We also performed a comprehensive review of all reported IMPDH1 disease causing variants with their associated phenotype.

MATERIALS AND METHODS

Multimodal imaging and functional studies documented the phenotype including best-corrected visual acuity (BCVA), fundus photograph, fundus autofluorescence (FAF), full field electroretinogram (ffERG), optical coherence tomography (OCT) and visual field (VF) data were collected. A literature search was performed in the PubMed and LOVD repositories.

RESULTS

We report 3 cases from a 2-generation family with a novel heterozygous likely pathogenic variant p. (Lys314Gln) (exon 10). The ophthalmic phenotype showed diffuse outer retinal atrophy with mild pigmentary changes with sparse pigmentary changes. FAF showed early macular involvement with macular hyperautofluorescence (hyperAF) surrounded by hypoAF. Foveal ellipsoid zone island can be found in the youngest patient but not in the older ones. The literature review identified a further 56 heterozygous, 1 compound heterozygous, and 2 homozygous variant. The heterozygous group included 43 missense, 3 in-frame, 1 nonsense, 2 frameshift, 1 synonymous, and 6 intronic variants. Exon 10 was noted as a hotspot harboring 18 variants.

CONCLUSIONS

We report a novel IMPDH1 variant. IMPDH1-associated retinopathy presents most frequently in the first decade of life with early macular involvement.

Electron Tunneling in Biology: When Does it Matter?

Electrons can tunnel between cofactor molecules positioned along biological electron transport chains up to a distance of ≃ 20 Å on the millisecond time scale of enzymatic turnover. This tunneling range determines the design of biological energy chains facilitating the cross-membrane transport of electrons. Tunneling distance and cofactors' redox potentials become the main physical parameters affecting the rate of electron transport. In addition, universal charge-transport properties are assigned to all proteins, making protein identity, flexibility, and dynamics insignificant. This paradigm is challenged by dynamical models of electron transfer, showing that the electron hopping rate is constant within the crossover distance R* ≃ 12 Å, followed with an exponential falloff at longer distances. If this hypothesis is fully confirmed, natural and man-made energy chains for electron transport should be best designed by placing redox cofactors near the crossover distance R*. Protein flexibility and dynamics affect the magnitude of the maximum hopping rate within the crossover distance. Changes in protein flexibility between forward and backward transitions contribute to vectorial charge transport. For biological energy chains, charge transport through proteins is not defined by universal parameters, and protein identity matters.

An NMR portrait of functional and dysfunctional allosteric cooperativity in cAMP-dependent protein kinase A

The cAMP-dependent protein kinase A (PKA) is the archetypical eukaryotic kinase. The catalytic subunit (PKA-C) structure is highly conserved among the AGC-kinase family. PKA-C is a bilobal enzyme with a dynamic N-lobe, harbouring the Adenosine-5'-triphosphate (ATP) binding site and a more rigid helical C-lobe. The substrate-binding groove resides at the interface of the two lobes. A distinct feature of PKA-C is the positive binding cooperativity between nucleotide and substrate. Several PKA-C mutations lead to the development of adenocarcinomas, myxomas, and other rare forms of liver tumours. Nuclear magnetic resonance (NMR) spectroscopy shows that these mutations disrupt the allosteric communication between the two lobes, causing a drastic decrease in binding cooperativity. The loss of cooperativity correlates with changes in substrate fidelity and reduced kinase affinity for the endogenous protein kinase inhibitor (PKI). The similarity between PKI and the inhibitory sequence of the kinase regulatory subunits suggests that the overall mechanism of regulation of the kinase may be disrupted. We surmise that a reduced or obliterated cooperativity may constitute a common trait for both orthosteric and allosteric mutations of PKA-C that may lead to dysregulation and disease.

The "violin model": Looking at community networks for dynamic allostery

Allosteric regulation of proteins continues to be an engaging research topic for the scientific community. Models describing allosteric communication have evolved from focusing on conformation-based descriptors of protein structural changes to appreciating the role of internal protein dynamics as a mediator of allostery. Here, we explain a "violin model" for allostery as a contemporary method for approaching the Cooper-Dryden model based on redistribution of protein thermal fluctuations. Based on graph theory, the violin model makes use of community network analysis to functionally cluster correlated protein motions obtained from molecular dynamics simulations. This Review provides the theory and workflow of the methodology and explains the application of violin model to unravel the workings of protein kinase A.

Ligands selectively tune the local and global motions of neurotensin receptor 1 (NTS1)

Nuclear magnetic resonance (NMR) studies have revealed that fast methyl sidechain dynamics can report on entropically-driven allostery. Yet, NMR applications have been largely limited to the super-microsecond motional regimes of G protein-coupled receptors (GPCRs). We use ¹³C^ε-methionine chemical shift-based global order parameters to test if ligands affect the fast dynamics of a thermostabilized GPCR, neurotensin receptor 1 (NTS₁). We establish that the NTS₁ solution ensemble includes substates with lifetimes on several, discrete timescales. The longest-lived states reflect those captured in agonist- and inverse agonist-bound crystal structures, separated by large energy barriers. We observe that the rapid fluctuations of individual methionine residues, superimposed on these long-lived states, respond collectively with the degree of fast, global dynamics correlating with ligand pharmacology. This approach lends confidence to interpreting spectra in terms of local structure and methyl dihedral angle geometry. The results suggest a role for sub-microsecond dynamics and conformational entropy in GPCR ligand discrimination.

Analysis of coordinated NMR chemical shifts to map allosteric regulatory networks in proteins

The exquisite sensitivity of the NMR chemical shift to local environment makes it an ideal probe to assess atomic level perturbations in proteins of all sizes and structural compositions. Recent advances in solution and solid-state NMR spectroscopy of biomolecules have leveraged the chemical shift to report on short- and long-range couplings between individual amino acids to establish "networks" of residues that form the basis of allosteric pathways that transmit chemical signals through the protein matrix to induce functional responses. The simple premise that thermodynamically and functionally coupled regions of a protein (i.e. active and allosteric sites) should be reciprocally sensitive to structural or dynamic perturbations has enabled NMR spectroscopy, the premier method for molecular resolution of protein structural fluctuations, to occupy a place at the forefront of investigations into protein allostery. Here, we detail several key methods of NMR chemical shift analysis to extract mechanistic information about long-range chemical signaling in a protein, focusing on practical methodological aspects and the circumstances under which a given approach would be relevant. We also detail some of the experimental considerations that should be made when applying these methods to specific protein systems.

Allosteric rescue of catalytically impaired ATP phosphoribosyltransferase variants links protein dynamics to active-site electrostatic preorganisation

ATP phosphoribosyltransferase catalyses the first step of histidine biosynthesis and is controlled via a complex allosteric mechanism where the regulatory protein HisZ enhances catalysis by the catalytic protein HisG_S while mediating allosteric inhibition by histidine. Activation by HisZ was proposed to position HisG_S Arg56 to stabilise departure of the pyrophosphate leaving group. Here we report active-site mutants of HisG_S with impaired reaction chemistry which can be allosterically restored by HisZ despite the HisZ:HisG_S interface lying ~20 Å away from the active site. MD simulations indicate HisZ binding constrains the dynamics of HisG_S to favour a preorganised active site where both Arg56 and Arg32 are poised to stabilise leaving-group departure in WT-HisG_S. In the Arg56Ala-HisG_S mutant, HisZ modulates Arg32 dynamics so that it can partially compensate for the absence of Arg56. These results illustrate how remote protein-protein interactions translate into catalytic resilience by restoring damaged electrostatic preorganisation at the active site.

Substrate Induced Movement of the Metal Cofactor between Active and Resting State

Regulation of enzyme activity is vital for living organisms. In metalloenzymes, far-reaching rearrangements of the protein scaffold are generally required to tune the metal cofactor's properties by allosteric regulation. Here structural analysis of hydroxyketoacid aldolase from Sphingomonas wittichii RW1 (SwHKA) revealed a dynamic movement of the metal cofactor between two coordination spheres without protein scaffold rearrangements. In its resting state configuration (M²⁺ _R ), the metal constitutes an integral part of the dimer interface within the overall hexameric assembly, but sterical constraints do not allow for substrate binding. Conversely, a second coordination sphere constitutes the catalytically active state (M²⁺ _A ) at 2.4 Å distance. Bidentate coordination of a ketoacid substrate to M²⁺ _A affords the overall lowest energy complex, which drives the transition from M²⁺ _R to M²⁺ _A . While not described earlier, this type of regulation may be widespread and largely overlooked due to low occupancy of some of its states in protein crystal structures.

Viruses, immunity and evolution

In this article, the evolution of viruses is analyzed in terms of their complexity. It is shown that the evolution of viruses is a partially directed process. The participation of viruses and mobile genetic elements in the evolution of other organisms by integration into the genome is also an a priori directed process. The high variability of genomes (including the genes of antibodies), which differs by orders of magnitude for various viruses and their hosts, is not a random process but is the result of the action of a molecular genetic control system. Herein, a model of partially directed evolution of viruses is proposed. Throughout the life cycle of viruses, there is an interaction of complex biologically important molecules that cannot be explained on the basis of classic laws. The interaction of a virus with a cell is essentially a quantum event, including selective long-range action. Such an interaction can be interpreted as the "remote key-lock" principle. In this article, a model of the interaction of biologically important viral molecules with cellular molecules based on nontrivial quantum interactions is proposed. Experiments to test the model are also proposed.

Emerging Methods and Applications to Decrypt Allostery in Proteins and Nucleic Acids

Many large protein-nucleic acid complexes exhibit allosteric regulation. In these systems, the propagation of the allosteric signaling is strongly coupled to conformational dynamics and catalytic function, challenging state-of-the-art analytical methods. Here, we review established and innovative approaches used to elucidate allosteric mechanisms in these complexes. Specifically, we report network models derived from graph theory and centrality analyses in combination with molecular dynamics (MD) simulations, introducing novel schemes that implement the synergistic use of graph theory with enhanced simulations methods and ab-initio MD. Accelerated MD simulations are used to construct "enhanced network models", describing the allosteric response over long timescales and capturing the relation between allostery and conformational changes. "Ab-initio network models" combine graph theory with ab-initio MD and quantum mechanics/molecular mechanics (QM/MM) simulations to describe the allosteric regulation of catalysis by following the step-by-step dynamics of biochemical reactions. This approach characterizes how the allosteric regulation changes from reactants to products and how it affects the transition state, revealing a tense-to-relaxed allosteric regulation along the chemical step. Allosteric models and applications are showcased for three paradigmatic examples of allostery in protein-nucleic acid complexes: (i) the nucleosome core particle, (ii) the CRISPR-Cas9 genome editing system and (iii) the spliceosome. These methods and applications create innovative protocols to determine allosteric mechanisms in protein-nucleic acid complexes that show tremendous promise for medicine and bioengineering.

Mutagenic Activation of Glutathione Peroxidase-4: Approaches toward Rational Design of Allosteric Drugs

Glutathione peroxidase 4 (GPX4) reduces lipid hydroperoxides in lipid membranes, effectively inhibiting iron-dependent cell death or ferroptosis. The upregulation of the enzyme by the mutations at residues D21 and D23 has been suggested to be associated with higher protein activity, which confers more protection against neurodegenerative diseases such as Alzheimer's, Parkinson's, and Huntington's diseases. Therefore, it has become an attractive target for treating and preventing neurodegenerative diseases. However, identifying means of mimicking the beneficial effects of these mutations distant from the active site constitutes a formidable challenge in moving toward therapeutics. In this study, we explore using molecular dynamics simulations to computationally map the conformational and energetic landscape of the wild-type GPX4 protein and three mutant variants to identify the allosteric networks of the enzyme. We present the conformational dynamic profile providing the desired signature behavior of the enzyme. We also discuss the implications of these findings for drug design efforts.

Conformational dynamics modulating electron transfer

Diffusional dynamics of the donor-acceptor distance are responsible for the appearance of a new time scale of diffusion over the distance of electronic tunneling in electron-transfer reactions. The distance dynamics compete with the medium polarization dynamics in the dynamics-controlled electron-transfer kinetics. The pre-exponential factor of the electron-transfer rate constant switches, at the crossover distance, between a distance-independent, dynamics-controlled plateau and exponential distance decay. The crossover between two regimes is controlled by an effective relaxation time slowed down by a factor exponentially depending on the variance of the donor-acceptor displacement. Flexible donor-acceptor complexes must show a greater tendency for dynamics-controlled electron transfer. Energy chains based on electron transport are best designed by placing the redox cofactors near the crossover distance.

In Search of a Dynamical Vocabulary: A Pipeline to Construct a Basis of Shared Traits in Large-Scale Motions of Proteins

The paradigmatic sequence–structure–dynamics–function relation in proteins is currently well established in the scientific community; in particular, a large effort has been made to probe the first connection, indeed providing convincing evidence of its strength and rationalizing it in a quantitative and general framework. In contrast, however, the role of dynamics as a link between structure and function has eluded a similarly clear-cut verification and description. In this work, we propose a pipeline aimed at building a basis for the quantitative characterization of the large-scale dynamics of a set of proteins, starting from the sole knowledge of their native structures. The method hinges on a dynamics-based clusterization, which allows a straightforward comparison with structural and functional protein classifications. The resulting basis set, obtained through the application to a group of related proteins, is shown to reproduce the salient large-scale dynamical features of the dataset. Most interestingly, the basis set is shown to encode the fluctuation patterns of homologous proteins not belonging to the initial dataset, thus highlighting the general applicability of the pipeline used to build it.

Prediction of Protein Allosteric Signalling Pathways and Functional Residues Through Paths of Optimised Propensity

Allostery commonly refers to the mechanism that regulates protein activity through the binding of a molecule at a different, usually distal, site from the orthosteric site. The omnipresence of allosteric regulation in nature and its potential for drug design and screening render the study of allostery invaluable. Nevertheless, challenges remain as few computational methods are available to effectively predict allosteric sites, identify signalling pathways involved in allostery, or to aid with the design of suitable molecules targeting such sites. Recently, bond-to-bond propensity analysis has been shown successful at identifying allosteric sites for a large and diverse group of proteins from knowledge of the orthosteric sites and its ligands alone by using network analysis applied to energy-weighted atomistic protein graphs. To address the identification of signalling pathways, we propose here a method to compute and score paths of optimised propensity that link the orthosteric site with the identified allosteric sites, and identifies crucial residues that contribute to those paths. We showcase the approach with three well-studied allosteric proteins: h-Ras, caspase-1, and 3-phosphoinositide-dependent kinase-1 (PDK1). Key residues in both orthosteric and allosteric sites were identified and showed agreement with experimental results, and pivotal signalling residues along the pathway were also revealed, thus providing alternative targets for drug design. By using the computed path scores, we were also able to differentiate the activity of different allosteric modulators.

Prediction and confirmation of a switch-like region within the N-terminal domain of hSIRT1

Many proteins display conformational changes resulting from allosteric regulation. Often only a few residues are crucial in conveying these structural and functional allosteric changes. These regions that undergo a significant change in structure upon receiving an input signal, such as molecular recognition, are defined as switch-like regions. Identifying these key residues within switch-like regions can help elucidate the mechanism of allosteric regulation and provide guidance for synthetic regulation. In this study, we combine a novel computational workflow with biochemical methods to identify a switch-like region in the N-terminal domain of human SIRT1 (hSIRT1), a lysine deacetylase that plays important roles in regulating cellular pathways. Based on primary sequence, computational methods predicted a region between residues 186-193 in hSIRT1 to exhibit switch-like behavior. Mutations were then introduced in this region and the resulting mutants were tested for allosteric reactions to resveratrol, a known hSIRT1 allosteric regulator. After fine-tuning the mutations based on comparison of known secondary structures, we were able to pinpoint M193 as the residue essential for allosteric regulation, likely by communicating the allosteric signal. Mutation of this residue maintained enzyme activity but abolished allosteric regulation by resveratrol. Our findings suggest a method to predict switch-like regions in allosterically regulated enzymes based on the primary sequence. If further validated, this could be an efficient way to identify key residues in enzymes for therapeutic drug targeting and other applications.

Prediction of Allosteric Communication Pathways in Proteins

MOTIVATION

Allostery in proteins is an essential phenomenon in biological processes. In this paper, we present a computational model to predict paths of maximum information transfer between active and allosteric sites. In this information theoretic study, we use mutual information as the measure of information transfer, where transition probability of information from one residue to its contacting neighbors is proportional to the magnitude of mutual information between the two residues. Starting from a given residue and using a Hidden Markov Model, we successively determine the neighboring residues that eventually lead to a path of optimum information transfer. The Gaussian approximation of mutual information between residue pairs is adopted. The limits of validity of this approximation are discussed in terms of a nonlinear theory of mutual information and its reduction to the Gaussian form.

RESULTS

Predictions of the model are tested on six widely studied cases, CheY Bacterial Chemotaxis, B-cell Lymphoma extra-large Bcl-xL, Human proline isomerase cyclophilin A (CypA), Dihydrofolate reductase DHFR, HRas GTPase, and Caspase-1. The communication transmission rendering the propagation of local fluctuations from the active sites throughout the structure in multiple paths correlate well with the known experimental data. Distinct paths originating from the active site may likely represent a multi functionality such as involving more than one allosteric site and/or preexistence of some other functional states. Our model is computationally fast and simple, and can give allosteric communication pathways, which are crucial for the understanding and control of protein functionality.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Identification of Core Allosteric Sites through Temperature- and Nucleus-Invariant Chemical Shift Covariance

Allosteric regulation is essential to control biological function. In addition, allosteric sites offer a promising venue for selective drug targeting. However, accurate mapping of allosteric sites remains challenging since allostery relies on often subtle, yet functionally relevant, structural and dynamical changes. A viable approach proposed to overcome such challenge is the chemical shift covariance analysis (CHESCA). Although CHESCA offers an exhaustive map of allosteric networks, it is critical to define the core allosteric sites to be prioritized in subsequent functional studies or the design of allosteric drugs. Here we propose two new CHESCA-based methodologies, called temperature CHESCA (T-CHESCA) and CLASS-CHESCA, aimed at narrowing down allosteric maps to the core allosteric residues. Both T- and CLASS-CHESCAs rely on the invariance of core inter-residue correlations to changes in the chemical shifts of the active and inactive conformations interconverting in fast exchange. In the T-CHESCA the chemical shifts of such states are modulated through temperature changes, while in the CLASS-CHESCA through variations in the spin-active nuclei involved in pairwise correlations. The T- and CLASS-CHESCAs as well as complete-linkage CHESCA were applied to the cAMP-binding domain of the exchange protein directly activated by cAMP (EPAC), which serves as a prototypical allosteric switch. Residues consistently identified by the three CHESCA methods were found in previously identified EPAC allosteric core sites. Hence, the T-, CLASS- and CL-CHESCA provide a toolset to establish allosteric site hierarchy and triage allosteric sites to be further analyzed by mutations and functional assays. Furthermore, the core allosteric networks selectively revealed through T- and CLASS-CHESCA are expected to facilitate the mechanistic understanding of disease-related mutations and the design of selective allosteric modulators.

Quantitative Analysis of Dynamic Allostery

Dynamic allostery refers to one important class of allosteric regulation that does not involve noticeable conformational changes upon effector binding. In recent years, many "quasi"-dynamic allosteric proteins have been found to only experience subtle conformational changes during allosteric regulation. However, as enthalpic and entropic contributions are coupled to each other and even tiny conformational changes could bring in noticeable free energy changes, a quantitative description is essential to understand the contribution of pure dynamic allostery. Here, by developing a unified anisotropic elastic network model (uANM) considering both side-chain information and ligand heavy atoms, we quantitatively estimated the contribution of pure dynamic allostery in a dataset of known allosteric proteins by excluding the conformational changes upon ligand binding. We found that the contribution of pure dynamic allostery is generally small (much weaker than previously expected) and robustly exhibits an allosteric activation effect, which exponentially decays with the distance between the substrate and the allosteric ligand. We further constructed toy models to study the determinant factors of dynamic allostery in monomeric and oligomeric proteins using the uANM. Analysis of the toy models revealed that a short distance, a small angle between the two ligands, strong protein-ligand interactions, and weak protein internal interactions lead to strong dynamic allostery. Our study provides a quantitative estimation of pure dynamic allostery and facilitates the understanding of dynamic-allostery-controlled biological processes and the design of allosteric drugs and proteins.

Identification of allosteric fingerprints of alpha-synuclein aggregates in matrix metalloprotease-1 and substrate-specific virtual screening with single molecule insights

Alpha-synuclein (aSyn) has implications in pathological protein aggregations in neurodegeneration. Matrix metalloproteases (MMPs) are broad-spectrum proteases and cleave aSyn, leading to aggregation. Previous reports showed that allosteric communications between the two domains of MMP1 on collagen fibril and fibrin depend on substrates, activity, and ligands. This paper reports quantification of allostery using single molecule measurements of MMP1 dynamics on aSyn-induced aggregates by calculating Forster Resonance Energy Transfer (FRET) between two dyes attached to the catalytic and hemopexin domains of MMP1. The two domains of MMP1 prefer open conformations that are inhibited by a single point mutation E219Q of MMP1 and tetracycline, an MMP inhibitor. A two-state Poisson process describes the interdomain dynamics, where the two states and kinetic rates of interconversion between them are obtained from histograms and autocorrelations of FRET values. Since a crystal structure of aSyn-bound MMP1 is unavailable, binding poses were predicted by molecular docking of MMP1 with aSyn using ClusPro. MMP1 dynamics were simulated using predicted binding poses and compared with the experimental interdomain dynamics to identify an appropriate pose. The selected aSyn-MMP1 binding pose near aSyn residue K45 was simulated and analyzed to define conformational changes at the catalytic site. Allosteric residues in aSyn-bound MMP1 exhibiting strong correlations with the catalytic motif residues were compared with allosteric residues in free MMP1, and aSyn-specific residues were identified. The allosteric residues in aSyn-bound MMP1 are K281, T283, G292, G327, L328, E329, R337, F343, G345, N346, Y348, G353, Q354, D363, Y365, S366, S367, F368, P371, R372, V374, K375, A379, F391, A394, R399, M414, F419, V426, and C466. Shannon entropy was defined to quantify MMP1 dynamics. Virtual screening was performed against a site on selected aSyn-MMP1 binding poses, which showed that lead molecules differ between free MMP1 and substrate-bound MMP1. Also, identifying aSyn-specific allosteric residues in MMP1 enabled further selection of lead molecules. In other words, virtual screening needs to take substrates into account for potential substrate-specific control of MMP1 activity in the future. Molecular understanding of interactions between MMP1 and aSyn-induced aggregates may open up the possibility of degrading aggregates by targeting MMPs.

Structural and dynamic insights into the HNH nuclease of divergent Cas9 species

CRISPR-Cas9 is a widely used biochemical tool with applications in molecular biology and precision medicine. The RNA-guided Cas9 protein uses its HNH endonuclease domain to cleave the DNA strand complementary to its endogenous guide RNA. In this study, novel constructs of HNH from two divergent organisms, G. stearothermophilus (GeoHNH) and S. pyogenes (SpHNH) were engineered from their respective full-length Cas9 proteins. Despite low sequence similarity, the X-ray crystal structures of these constructs reveal that the core of HNH surrounding the active site is conserved. Structure prediction of the full-length GeoCas9 protein using Phyre2 and AlphaFold2 also showed that the crystallographic construct of GeoHNH represents the structure of the domain within the full-length GeoCas9 protein. However, significant differences are observed in the solution dynamics of structurally conserved regions of GeoHNH and SpHNH, the latter of which was shown to use such molecular motions to propagate the DNA cleavage signal. Indeed, molecular simulations show that the intradomain signaling pathways, which drive SpHNH function, are non-specific and poorly formed in GeoHNH. Taken together, these outcomes suggest mechanistic differences between mesophilic and thermophilic Cas9 species.

Intrinsically disordered proteins: Ensembles at the limits of Anfinsen's dogma

Intrinsically disordered proteins (IDPs) are proteins that lack rigid 3D structure. Hence, they are often misconceived to present a challenge to Anfinsen's dogma. However, IDPs exist as ensembles that sample a quasi-continuum of rapidly interconverting conformations and, as such, may represent proteins at the extreme limit of the Anfinsen postulate. IDPs play important biological roles and are key components of the cellular protein interaction network (PIN). Many IDPs can interconvert between disordered and ordered states as they bind to appropriate partners. Conformational dynamics of IDPs contribute to conformational noise in the cell. Thus, the dysregulation of IDPs contributes to increased noise and "promiscuous" interactions. This leads to PIN rewiring to output an appropriate response underscoring the critical role of IDPs in cellular decision making. Nonetheless, IDPs are not easily tractable experimentally. Furthermore, in the absence of a reference conformation, discerning the energy landscape representation of the weakly funneled IDPs in terms of reaction coordinates is challenging. To understand conformational dynamics in real time and decipher how IDPs recognize multiple binding partners with high specificity, several sophisticated knowledge-based and physics-based in silico sampling techniques have been developed. Here, using specific examples, we highlight recent advances in energy landscape visualization and molecular dynamics simulations to discern conformational dynamics and discuss how the conformational preferences of IDPs modulate their function, especially in phenotypic switching. Finally, we discuss recent progress in identifying small molecules targeting IDPs underscoring the potential therapeutic value of IDPs. Understanding structure and function of IDPs can not only provide new insight on cellular decision making but may also help to refine and extend Anfinsen's structure/function paradigm.

Complete Model of Vinculin Suggests the Mechanism of Activation by Helical Super-Bundle Unfurling

To shed light onto the activation mechanism of vinculin, we carried out a detailed refinement of chicken vinculin and compared it to the human protein which is greater than 95% identical. Refinement resulted in a complete and significantly improved model. This model includes important elements such as a pro-rich strap region (PRR) and C-terminus. The conformation of the PRR stabilized by its inter- and intra-molecular contacts shows a dynamic, but relatively stable motif that constitutes a docking platform for multiple molecules. The contact of the C-terminus with the PRR suggests that phosphorylation of Tyr1065 might control activation and membrane binding. Improved electron densities showed the presence of large solvent molecules such as phosphates/sulfates and a head-group of PIP2. The improved model allowed for a computational stability analysis to be performed by the program Corex/Best which located numerous hot-spots of increased and decreased stability. Proximity of the identified binding sites for regulatory partners involved in inducing or suppressing the activation of vinculin to the unstable elements sheds new light onto the activation pathway and differential activation. This stability analysis suggests that the activation pathway proceeds by unfurling of the super-bundle built from four bundles of helices without separation of the Vt region (840-1066) from the head. According to our mechanism, when activating proteins bind at the strap region a separation of N and C terminal bundles occurs, followed by unfurling of the super-bundle and flattening of the general shape of the molecule, which exposes the interaction sites for binding of auxiliary proteins.

Recent advances in externally controlled ring-opening polymerisations

Switchable catalysis is a powerful tool in the polymer chemist's toolbox as it allows on demand access to a variety of polymer architectures. Switchable catalysts operate by the generation of a species which is chemically distinct in behaviour and structure to the precursor. This difference in catalytic activity has been exploited to allow spatiotemporal control over polymerisations in the synthesis of (co)polymers. Although switchable methodologies have been applied to other polymerisation mechanisms for quite some time, for ring opening polymerisation (ROP) reactions it is a relatively young area of research. Despite its infancy, the field is accelerating rapidly. Here, we review recent developments for selected external stimuli for ROP, including redox chemistry, light, allosteric and mechanical control. Furthermore, a brief review on switch catalysis involving exogeneous gases will also be provided, although this area differs from traditional switchable catalysis techniques. An outlook on the future of switchable catalysis is also provided.

A Free-Energy Landscape Analysis of Calmodulin Obtained from an NMR Data-Utilized Multi-Scale Divide-and-Conquer Molecular Dynamics Simulation

Calmodulin (CaM) is a multifunctional calcium-binding protein, which regulates a variety of biochemical processes. CaM acts through its conformational changes and complex formation with its target enzymes. CaM consists of two globular domains (N-lobe and C-lobe) linked by an extended linker region. Upon calcium binding, the N-lobe and C-lobe undergo local conformational changes, followed by a major conformational change of the entire CaM to wrap the target enzyme. However, the regulation mechanisms, such as allosteric interactions, which regulate the large structural changes, are still unclear. In order to investigate the series of structural changes, the free-energy landscape of CaM was obtained by multi-scale divide-and-conquer molecular dynamics (MSDC-MD). The resultant free-energy landscape (FEL) shows that the Ca²⁺ bound CaM (holo-CaM) would take an experimentally famous elongated structure, which can be formed in the early stage of structural change, by breaking the inter-domain interactions. The FEL also shows that important interactions complete the structural change from the elongated structure to the ring-like structure. In addition, the FEL might give a guiding principle to predict mutational sites in CaM. In this study, it was demonstrated that the movement process of macroscopic variables on the FEL may be diffusive to some extent, and then, the MSDC-MD is suitable to the parallel computation.

SBMOpenMM: A Builder of Structure-Based Models for OpenMM

Molecular dynamics (MD) simulations have become a standard tool to correlate the structure and function of biomolecules and significant advances have been made in the study of proteins and their complexes. A major drawback of conventional MD simulations is the difficulty and cost of obtaining converged results, especially when exploring potential energy surfaces containing considerable energy barriers. This limits the wide use of MD calculations to determine the thermodynamic properties of biomolecular processes. Indeed, this is true when considering the conformational entropy of such processes, which is ultimately critical in assessing the simulations' convergence. Alternatively, a wide range of structure-based models (SBMs) has been used in the literature to unravel the basic mechanisms of biomolecular dynamics. These models introduce simplifications that focus on the relevant aspects of the physical process under study. Because of this, SBMs incorporate the need to modify the force field definition and parameters to target specific biophysical simulations. Here we introduce SBMOpenMM, a Python library to build force fields for SBMs, that uses the OpenMM framework to create and run SBM simulations. The code is flexible and user-friendly and profits from the high customizability and performance provided by the OpenMM platform.

Dimeric allostery mechanism of the plant circadian clock photoreceptor ZEITLUPE

In Arabidopsis thaliana, the Light-Oxygen-Voltage (LOV) domain containing protein ZEITLUPE (ZTL) integrates light quality, intensity, and duration into regulation of the circadian clock. Recent structural and biochemical studies of ZTL indicate that the protein diverges from other members of the LOV superfamily in its allosteric mechanism, and that the divergent allosteric mechanism hinges upon conservation of two signaling residues G46 and V48 that alter dynamic motions of a Gln residue implicated in signal transduction in all LOV proteins. Here, we delineate the allosteric mechanism of ZTL via an integrated computational approach that employs atomistic simulations of wild type and allosteric variants of ZTL in the functional dark and light states, together with Markov state and supervised machine learning classification models. This approach has unveiled key factors of the ZTL allosteric mechanisms, and identified specific interactions and residues implicated in functional allosteric changes. The final results reveal atomic level insights into allosteric mechanisms of ZTL function that operate via a non-trivial combination of population-shift and dynamics-driven allosteric pathways.

Allosteric activation of T cell antigen receptor signaling by quaternary structure relaxation

The mechanism of T cell antigen receptor (TCR-CD3) signaling remains elusive. Here, we identify mutations in the transmembrane region of TCRβ or CD3ζ that augment peptide T cell antigen receptor (pMHC)-induced signaling not explicable by enhanced ligand binding, lateral diffusion, clustering, or co-receptor function. Using a biochemical assay and molecular dynamics simulation, we demonstrate that the gain-of-function mutations loosen the interaction between TCRαβ and CD3ζ. Similar to the activating mutations, pMHC binding reduces TCRαβ cohesion with CD3ζ. This event occurs prior to CD3ζ phosphorylation and at 0°C. Moreover, we demonstrate that soluble monovalent pMHC alone induces signaling and reduces TCRαβ cohesion with CD3ζ in membrane-bound or solubilised TCR-CD3. Our data provide compelling evidence that pMHC binding suffices to activate allosteric changes propagating from TCRαβ to the CD3 subunits, reconfiguring interchain transmembrane region interactions. These dynamic modifications could change the arrangement of TCR-CD3 boundary lipids to license CD3ζ phosphorylation and initiate signal propagation.

•

Mutations in TCRβ and CD3ζ TMRs that reduce their interaction augment signaling

•

pMHC and anti-CD3 binding to TCR-CD3 induce similar quaternary structure relaxation

•

Soluble monovalent pMHC alone signals and reduces TCRαβ cohesion with CD3ζ

•

Allosteric changes in TCR-CD3 dynamics instigate T cell activation

Coarse-Grained Modeling of Multiple Pathways in Conformational Transitions of Multi-Domain Proteins

Large-scale conformational transitions in multi-domain proteins are often essential for their functions. To investigate the transitions, it is necessary to explore multiple potential pathways, which involve different intermediate structures. Here, we present a multi-basin (MB) coarse-grained (CG) structure-based Go̅ model for describing transitions in proteins with more than two moving domains. This model is an extension of our dual-basin Go̅ model in which system-dependent parameters are determined systematically using the multistate Bennett acceptance ratio method. In the MB Go̅ model for multi-domain proteins, we assume that intermediate structures may have partial inter-domain native contacts. This approach allows us to search multiple transition pathways that involve distinct intermediate structures using the CG molecular dynamics (MD) simulations. We apply this scheme to an enzyme, adenylate kinase (AdK), which has three major domains and can move along two different pathways. Using the optimized mixing parameters for each pathway, AdK shows frequent transitions between the Open, Closed, and the intermediate basins and samples a wide variety of conformations within each basin. The explored multiple transition pathways could be compared with experimental data and examined in more detail by atomistic MD simulations.

Structural resolution of switchable states of a de novo peptide assembly

De novo protein design is advancing rapidly. However, most designs are for single states. Here we report a de novo designed peptide that forms multiple α-helical-bundle states that are accessible and interconvertible under the same conditions. Usually in such designs amphipathic α helices associate to form compact structures with consolidated hydrophobic cores. However, recent rational and computational designs have delivered open α-helical barrels with functionalisable cavities. By placing glycine judiciously in the helical interfaces of an α-helical barrel, we obtain both open and compact states in a single protein crystal. Molecular dynamics simulations indicate a free-energy landscape with multiple and interconverting states. Together, these findings suggest a frustrated system in which steric interactions that maintain the open barrel and the hydrophobic effect that drives complete collapse are traded-off. Indeed, addition of a hydrophobic co-solvent that can bind within the barrel affects the switch between the states both in silico and experimentally.

So far most of the de novo designed proteins are for single states only. Here, the authors present the de novo design and crystal structure determination of a coiled-coil peptide that assembles into multiple, distinct conformational states under the same conditions and further characterise its properties with biophysical experiments, NMR and MD simulations.

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.

Synchronous and Asynchronous Response in Dynamically Perturbed Proteins

We present a dynamic perturbation-response model of proteins based on the Gaussian Network Model, where a residue is perturbed periodically, and the dynamic response of other residues is determined. The model shows that periodic perturbation causes a synchronous response in phase with the perturbation and an asynchronous response that is out of phase. The asynchronous component results from the viscous effects of the solvent and other dispersive factors in the system. The model is based on the solution of the Langevin equation in the presence of solvent, noise, and perturbation. We introduce several novel ideas: The concept of storage and loss compliance of the protein and their dependence on structure and frequency; the amount of work lost and the residues that contribute significantly to the lost work; new dynamic correlations that result from perturbation; causality, that is, the response of j when i is perturbed is not equal to the response of i when j is perturbed. As examples, we study two systems, namely, bovine rhodopsin and the class of nanobodies. The general results obtained are (i) synchronous and asynchronous correlations depend strongly on the frequency of perturbation, their magnitude decreases with increasing frequency, (ii) time-delayed mean-squared fluctuations of residues have only synchronous components. Asynchronicity is present only in cross correlations, that is, correlations between different residues, (iii) perturbation of loop residues leads to a large dissipation of work, (iv) correlations satisfy the hypothesis of pre-existing pathways according to which information transfer by perturbation rides on already existing equilibrium correlations in the system, (v) dynamic perturbation can introduce a selective response in the system, where the perturbation of each residue excites different sets of responding residues, and (vi) it is possible to identify nondissipative residues whose perturbation does not lead to dissipation in the protein. Despite its simplicity, the model explains several features of allosteric manipulation.

Exploring Large Domain Motions in Proteins Using Atomistic Molecular Dynamics with Enhanced Conformational Sampling

Conformational transitions in multidomain proteins are essential for biological functions. The Apo conformations are typically open and flexible, while the Holo states form more compact conformations stabilized by protein-ligand interactions. Unfortunately, the atomically detailed mechanisms for such open-closed conformational changes are difficult to be accessed experimentally as well as computationally. To simulate the transitions using atomistic molecular dynamics (MD) simulations, efficient conformational sampling algorithms are required. In this work, we propose a new approach based on generalized replica-exchange with solute tempering (gREST) for exploring the open-closed conformational changes in multidomain proteins. Wherein, selected surface charged residues in a target protein are defined as the solute region in gREST simulation and the solute temperatures are different in replicas and exchanged between them to enhance the domain motions. This approach is called gREST selected surface charged residues (gREST_SSCR) and is applied to the Apo and Holo states of ribose binding protein (RBP) in solution. The conformational spaces sampled with gREST_SSCR are much wider than those with the conventional MD, sampling open-closed conformational changes while maintaining RBP domains’ stability. The free-energy landscapes of RBP in the Apo and Holo states are drawn along with twist and hinge angles of the two moving domains. The inter-domain salt-bridges that are not observed in the experimental structures are also important in the intermediate states during the conformational changes.

Transient association between proteins elicits alteration of dynamics at sites far away from interfaces

Proteins are known to undergo structural changes upon binding to partner proteins. However, the prevalence, extent, location, and function of change in protein dynamics due to transient protein-protein interactions is not well documented. Here, we have analyzed a dataset of 58 protein-protein complexes of known three-dimensional structure and structures of their corresponding unbound forms to evaluate dynamics changes induced by binding. Fifty-five percent of cases showed significant dynamics change away from the interfaces. This change is not always accompanied by an observed structural change. Binding of protein partner is found to alter inter-residue communication within the tertiary structure in about 90% of cases. Also, residue motions accessible to proteins in unbound form were not always maintained in the bound form. Further analyses revealed functional roles for the distant site where dynamics change was observed. Overall, the results presented here strongly suggest that alteration of protein dynamics due to binding of a partner protein commonly occurs.

Visualizing the Dynamics of a Protein Folding Machinery: The Mechanism of Asymmetric ATP Processing in Hsp90 and its Implications for Client Remodelling

The Hsp90 chaperone system interacts with a wide spectrum of client proteins, forming variable and dynamic multiprotein complexes that involve the intervention of cochaperone partners. Recent results suggest that the role of Hsp90 complexes is to establish interactions that suppress unwanted client activities, allow clients to be protected from degradation and respond to biochemical signals. Cryo-electron microscopy (cryoEM) provided the first key molecular picture of Hsp90 in complex with a kinase, Cdk4, and a cochaperone, Cdc37. Here, we use a combination of molecular dynamics (MD) simulations and advanced comparative analysis methods to elucidate key aspects of the functional dynamics of the complex, with different nucleotides bound at the N-terminal Domain of Hsp90. The results reveal that nucleotide-dependent structural modulations reverberate in a striking asymmetry of the dynamics of Hsp90 and identify specific patterns of long-range coordination between the nucleotide binding site, the client binding pocket, the cochaperone and the client. Our model establishes a direct atomic-resolution cross-talk between the ATP-binding site, the client region that is to be remodeled and the surfaces of the Cdc37-cochaperone.

Mechanisms of ligand binding

Many processes in chemistry and biology involve interactions of a ligand with its molecular target. Interest in the mechanism governing such interactions has dominated theoretical and experimental analysis for over a century. The interpretation of molecular recognition has evolved from a simple rigid body association of the ligand with its target to appreciation of the key role played by conformational transitions. Two conceptually distinct descriptions have had a profound impact on our understanding of mechanisms of ligand binding. The first description, referred to as induced fit, assumes that conformational changes follow the initial binding step to optimize the complex between the ligand and its target. The second description, referred to as conformational selection, assumes that the free target exists in multiple conformations in equilibrium and that the ligand selects the optimal one for binding. Both descriptions can be merged into more complex reaction schemes that better describe the functional repertoire of macromolecular systems. This review deals with basic mechanisms of ligand binding, with special emphasis on induced fit, conformational selection, and their mathematical foundations to provide rigorous context for the analysis and interpretation of experimental data. We show that conformational selection is a surprisingly versatile mechanism that includes induced fit as a mathematical special case and even captures kinetic properties of more complex reaction schemes. These features make conformational selection a dominant mechanism of molecular recognition in biology, consistent with the rich conformational landscape accessible to biological macromolecules being unraveled by structural biology.

A Time-Dependent Quantum Approach to Allostery and a Comparison With Light-Harvesting in Photosynthetic Phenomenon

The allosteric effect is one of the most important processes in regulating the function of proteins, and the elucidation of this phenomenon plays a significant role in understanding emergent behaviors in biological regulation. In this process, a perturbation, generated by a ligand in a part of the macromolecule (the allosteric site), moves along this system and reaches a specific (active) site, dozens of Ångströms away, with a great efficiency. The dynamics of this perturbation in the macromolecule can model precisely the allosteric process. In this article, we will be studying the general characteristics of allostery, using a time-dependent quantum approach to obtain rules that apply to this kind of process. Considering the perturbation as a wave that moves within the molecular system, we will characterize the allosteric process with three of the properties of this wave in the active site: (1) t_a, the characteristic time for reaching that site, (2) A_a, the amplitude of the wave in this site, and (3) B_a, its corresponding spectral broadening. These three parameters, together with the process mechanism and the perturbation efficiency in the process, can describe the phenomenon. One of the main purposes of this paper is to link the parameters t_a, A_a, and B_a and the perturbation efficiency to the characteristics of the system. There is another fundamental process for life that has some characteristics similar to allostery: the light-harvesting (LH) process in photosynthesis. Here, as in allostery, two distant macromolecular sites are involved—two sites dozens of Ångströms away. In both processes, it is particularly important that the perturbation is distributed efficiently without dissipating in the infinite degrees of freedom within the macromolecule. The importance of considering quantum effects in the LH process is well documented in literature, and the quantum coherences are experimentally proven by time-dependent spectroscopic techniques. Given the existing similarities between these two processes in macromolecules, in this work, we suggest using Quantum Mechanics (QM) to study allostery.