Dynamics, flexibility and ligand-induced conformational changes in biological macromolecules: a computational approach

An open source in silico workflow to assist in the design of fusion proteins

Fusion proteins have the potential to become the new norm for targeted therapeutic treatments. Highly specific payload delivery can be achieved by combining custom targeting moieties, such as VHH domains, with active parts of proteins that have a particular activity not naturally targeted to the intended cells. Conversely, novel drug products may make use of the highly specific targeting properties of naturally occurring proteins and combine them with custom payloads. When designing such a product, there is rarely a known structure for the final construct which makes it difficult to assess molecular behaviour that may ultimately impact therapeutic outcome. Considering the time and cost of expressing a construct, optimising the purification procedure, obtaining sufficient quantities for biophysical characterisation, and performing structural studies in vitro, there is an enormous benefit to conduct in silico studies ahead of wet lab work. By following a repeatable, streamlined, and fast workflow of molecular dynamics assessment, it is possible to eliminate low-performing candidates from costly experimental work. There are, however, many aspects to consider when designing a novel fusion protein and it is crucial not to overlook some elements. In this work, we suggest a set of user-friendly, open-source methods which can be used to screen fusion protein candidates from the sequence alone. We used the light chain and translocation domain of botulinum toxin A (BoNT/A) fused with a selected VHH domain, termed here LC-HN-VHH, as a case study for a general approach to designing, modelling, and simulating fusion proteins. Its behaviour in silico correlated well with initial in vitro work, with SEC HPLC showing multiple protein states in solution and a dynamic protein shifting between these states over time without loss of material.

Investigating the effects of cysteine-118 oxidation on G12D KRas structure and dynamics: insights from MD simulations

Mutations of Ras proteins are believed to be among the most prominent causes of cancer. There is increasing evidence that the activity of Ras may be controlled by the redox state of cysteine residues located within the NKCD motif. This redox signaling is critical to both physiological and pathological processes and occurs when C118 is oxidized in a reversible manner. In this study, we used atomistic molecular dynamics simulations and Markov state models to investigate the structural and conformational effects of C118 oxidation on the oncogenic mutant KRas(G12D). While both mutants share common features and exhibit some distinct conformational states and fluctuations, we have found that the oxidized variant KRas(G12D/C118SOH) is more dynamic than the unoxidized counterpart, particularly in the switch II region. Additionally, C118 oxidation is found to alter the structure of the nucleotide-binding site and the switch regions as well as perturb the conformational equilibrium between Ras active and inactive states. These conformational preferences may alter the affinity to different effectors, resulting in selective downstream activation. Our results are anticipated to help future drug development efforts aimed at KRAS-related anticancer treatment.Communicated by Ramaswamy H. Sarma.

The complexities of ligand/receptor interactions: Exploring the role of molecular vibrations and quantum tunnelling

Molecular vibrations and quantum tunneling may link ligand binding to the function of pharmacological receptors. The well-established lock-and-key model explains a ligand's binding and recognition by a receptor; however, a general mechanism by which receptors translate binding into activation, inactivation, or modulation remains elusive. The Vibration Theory of Olfaction was proposed in the 1930s to explain this subset of receptor-mediated phenomena by correlating odorant molecular vibrations to smell, but a mechanism was lacking. In the 1990s, inelastic electron tunneling was proposed as a plausible mechanism for translating molecular vibration to odorant physiology. More recently, studies of ligands' vibrational spectra and the use of deuterated ligand analogs have provided helpful information to study this admittedly controversial hypothesis in metabotropic receptors other than olfactory receptors. In the present work, based in part on published experiments from our laboratory using planarians as an experimental organism, I will present a rationale and possible experimental approach for extending this idea to ligand-gated ion channels.

Molecular Docking and Molecular Dynamics Simulations in Related to Leishmania donovani: An Update and Literature Review

Leishmaniasis, a disease caused by Leishmania parasites and transmitted via sandflies, presents in two main forms: cutaneous and visceral, the latter being more severe. With 0.7 to 1 million new cases each year, primarily in Brazil, diagnosing remains challenging due to diverse disease manifestations. Traditionally, the identification of Leishmania species is inferred from clinical and epidemiological data. Advances in disease management depend on technological progress and the improvement of parasite identification programs. Current treatments, despite the high incidence, show limited efficacy due to factors like cost, toxicity, and lengthy regimens causing poor adherence and resistance development. Diagnostic techniques have improved but a significant gap remains between scientific progress and application in endemic areas. Complete genomic sequence knowledge of Leishmania allows for the identification of therapeutic targets. With the aid of computational tools, testing, searching, and detecting affinity in molecular docking are optimized, and strategies that assess advantages among different options are developed. The review focuses on the use of molecular docking and molecular dynamics (MD) simulation for drug development. It also discusses the limitations and advancements of current treatments, emphasizing the importance of new techniques in improving disease management.

Elucidation of critical pH-dependent structural changes in Botulinum Neurotoxin E

Botulinum Neurotoxins (BoNT) are the most potent toxins currently known. However, they also have therapeutic applications for an increasing number of motor related conditions due to their specificity, and low diffusion into the system. Although the start- and end- points for the BoNT mechanism of action are well-studied, a critical step remains poorly understood. It is theorised that BoNTs undergo a pH-triggered conformational shift, activating the neurotoxin by priming it to form a transmembrane (TM) channel. To test this hypothesis, we combined molecular dynamics (MD) simulations and small-angle x-ray scattering (SAXS), revealing a new conformation of serotype E (BoNT/E). This conformation was exclusively observed in simulations below pH 5.5, as determined by principal component analysis (PCA), and its theoretical SAXS profile matched an experimental SAXS profile obtained at pH 4. Additionally, a localised secondary structural change was observed in MD simulations below pH 5.5, in a region previously identified as instrumental for membrane insertion for serotype A (BoNT/A). These changes were found at a critical pH value for BoNTs in vivo, and may be relevant for their therapeutic use.

Structure-based identification of potential SARS-CoV-2 main protease inhibitors

To address coronavirus disease (COVID-19), currently, no effective drug or vaccine is available. In this regard, molecular modeling approaches are highly useful to discover potential inhibitors of the main protease (Mpro) enzyme of SARS-CoV-2. Since, the Mpro enzyme plays key roles in mediating viral replication and transcription; therefore, it is considered as an attractive drug target to control SARS-CoV-2 infection. By using structure-based drug design, pharmacophore modeling, and virtual high throughput drug screening combined with docking and all-atom molecular dynamics simulation approach, we have identified five potential inhibitors of SARS-CoV-2 Mpro. MD simulation studies revealed that compound 54035018 binds to the Mpro with high affinity (ΔGbind -37.40 kcal/mol), and the complex is more stable in comparison with other protein-ligand complexes. We have identified promising leads to fight COVID-19 infection as these compounds fulfill all drug-likeness properties. However, experimental and clinical validations are required for COVID-19 therapy.Communicated by Ramaswamy H. Sarma.

Probing Structural Perturbation of Biomolecules by Extracting Cryo-EM Data Heterogeneity

Single-particle cryogenic electron microscopy (cryo-EM) has become an indispensable tool to probe high-resolution structural detail of biomolecules. It enables direct visualization of the biomolecules and opens a possibility for averaging molecular images to reconstruct a three-dimensional Coulomb potential density map. Newly developed algorithms for data analysis allow for the extraction of structural heterogeneity from a massive and low signal-to-noise-ratio (SNR) cryo-EM dataset, expanding our understanding of multiple conformational states, or further implications in dynamics, of the target biomolecule. This review provides an overview that briefly describes the workflow of single-particle cryo-EM, including imaging and data processing, and new methods developed for analyzing the data heterogeneity to understand the structural variability of biomolecules.

The role of C-terminal helix in the conformational transition of an arginine binding protein

•

Probe the role of C-ter. helix (CTH) in conformational transition of TmArgBP.

•

Presence of CTH almost doubles the barrier to access the closed-state.

•

In the absence of CTH, the protein can fluctuate between the two conformations.

•

CTH not only constraints the open-state conformation but also guides in accessing it.

The thermotoga maritima arginine binding protein (TmArgBP) is a periplasmic binding protein that has a short helix at the C-terminal end (CTH), which is swapped between the two chains. We apply a coarse-grained structure-based model (SBM) and all-atom MD simulation on this protein to understand the mechanism and the role of CTH in the conformational transition. When the results of SBM simulations of TmArgBP in the presence and absence of CTH are compared, we find that CTH is strategically located at the back of the binding pocket restraining the open-state conformation thereby disengaging access to the closed-state. We also ran all-atom MD simulations of open-state TmArgBP with and without CTH and discovered that in the absence of CTH the protein could reach the closed-state within 250 ns, while in its presence, the protein remained predominantly in its open-state conformation. In the simulation started from unliganded closed-state conformation without CTH, the protein exhibited multiple transitions between the two states, suggesting CTH as an essential structural element to stabilize the open-state conformation. In another simulation that began with an unliganded closed-state conformation with CTH, the protein was able to access the open-state. In this simulation the CTH was observed to reorient itself to interact with the protein emphasizing its role in assisting the conformational change. Based on our findings, we believe that CTH not only acts as a structural element that constraints the protein in its open-state but it may also guide the protein back to its open-state conformation upon ligand unbinding.

Pseudokinases repurpose flexibility signatures associated with the protein kinase fold for noncatalytic roles

The bilobal protein kinase-like fold in pseudokinases lack one or more catalytic residues, conserved in canonical protein kinases, and are considered enzymatically deficient. Tertiary structures of pseudokinases reveal that their loops topologically equivalent to activation segments of kinases adopt contracted configurations, which is typically extended in active conformation of kinases. Herein, anisotropic network model based normal mode analysis (NMA) was conducted on 51 active conformation structures of protein kinases and 26 crystal structures of pseudokinases. Our observations indicate that although backbone fluctuation profiles are similar for individual kinase-pseudokinase families, low intensity mean square fluctuations in pseudo-activation segment and other sub-structures impart rigidity to pseudokinases. Analyses of collective motions from functional modes reveal that pseudokinases, compared to active kinases, undergo distinct conformational transitions using the same structural fold. All-atom NMA of protein kinase-pseudokinase pairs from each family, sharing high amino acid sequence identities, yielded distinct community clusters, partitioned by residues exhibiting highly correlated fluctuations. It appears that atomic fluctuations from equivalent activation segments guide community membership and network topologies for respective kinase and pseudokinase. Our findings indicate that such adaptations in backbone and side-chain fluctuations render pseudokinases competent for catalysis-independent roles.

Ligand Binding Introduces Significant Allosteric Shifts in the Locations of Protein Fluctuations

Allostery is usually considered to be a mechanism for transmission of signals associated with physical or dynamic changes in some part of a protein. Here, we investigate the changes in fluctuations across the protein upon ligand binding based on the fluctuations computed with elastic network models. These results suggest that binding reduces the fluctuations at the binding site but increases fluctuations at remote sites, but not to fully compensating extents. If there were complete conservation of entropy, then only the enthalpies of binding would matter and not the entropies; however this does not appear to be the case. Experimental evidence also suggests that energies and entropies of binding can compensate but that the extent of compensation varies widely from case to case. Our results do however always show transmission of an allosteric signal to distant locations where the fluctuations are increased. These fluctuations could be used to compute entropies to improve evaluations of the thermodynamics of binding. We also show the allosteric relationship between peptide binding in the GroEL trans-ring that leads directly to the release of GroES from the GroEL-GroES cis-ring. This finding provides an example of how calculating these changes to protein dynamics induced by the binding of an allosteric ligand can regulate protein function and mechanism.

Molecular flexibility in computational protein design: an algorithmic perspective

Computational protein design (CPD) is a powerful technique for engineering new proteins, with both great fundamental implications and diverse practical interests. However, the approximations usually made for computational efficiency, using a single fixed backbone and a discrete set of side chain rotamers, tend to produce rigid and hyper-stable folds that may lack functionality. These approximations contrast with the demonstrated importance of molecular flexibility and motions in a wide range of protein functions. The integration of backbone flexibility and multiple conformational states in CPD, in order to relieve the inaccuracies resulting from these simplifications and to improve design reliability, are attracting increased attention. However, the greatly increased search space that needs to be explored in these extensions defines extremely challenging computational problems. In this review, we outline the principles of CPD and discuss recent effort in algorithmic developments for incorporating molecular flexibility in the design process.

Acute Intermittent Porphyria: An Overview of Therapy Developments and Future Perspectives Focusing on Stabilisation of HMBS and Proteostasis Regulators

Acute intermittent porphyria (AIP) is an autosomal dominant inherited disease with low clinical penetrance, caused by mutations in the hydroxymethylbilane synthase (HMBS) gene, which encodes the third enzyme in the haem biosynthesis pathway. In susceptible HMBS mutation carriers, triggering factors such as hormonal changes and commonly used drugs induce an overproduction and accumulation of toxic haem precursors in the liver. Clinically, this presents as acute attacks characterised by severe abdominal pain and a wide array of neurological and psychiatric symptoms, and, in the long-term setting, the development of primary liver cancer, hypertension and kidney failure. Treatment options are few, and therapies preventing the development of symptomatic disease and long-term complications are non-existent. Here, we provide an overview of the disorder and treatments already in use in clinical practice, in addition to other therapies under development or in the pipeline. We also introduce the pathomechanistic effects of HMBS mutations, and present and discuss emerging therapeutic options based on HMBS stabilisation and the regulation of proteostasis. These are novel mechanistic therapeutic approaches with the potential of prophylactic correction of the disease by totally or partially recovering the enzyme functionality. The present scenario appears promising for upcoming patient-tailored interventions in AIP.

Crystallographic and molecular dynamics simulation analysis of NAD synthetase from methicillin resistant Staphylococcus aureus (MRSA)

NAD synthetase (NadE) catalyzes the last step in NAD biosynthesis, transforming deamido-NAD⁺ into NAD⁺ by a two-step reaction with co-substrates ATP and amide donor ammonia. In this study, we report the crystal structure of Staphylococcus aureus NAD synthetase enzyme (saNadE) at 2.3 Å resolution. We used this structure to perform molecular dynamics simulations of apo-enzyme, enzyme-substrate (NadE with ATP and NaAD) and enzyme-intermediate complexes (NadE with NaAD-AMP) to investigate key binding interactions and explore the conformational transitions and flexibility of the binding pocket. Our results show large shift of N-terminal region in substrate bound form which is important for ATP binding. Substrates drive the correlated movement of loop regions surrounding it as well as some regions distal to the active site and stabilize them at complex state. Principal component analysis of atomic projections distinguish feasible trajectories to delineate distinct motions in enzyme-substrate to enzyme-intermediate states. Our results suggest mixed binding involving dominant induced fit and conformational selection. MD simulation extracted ensembles of NadE could potentially be utilized for in silico screening and structure based design of more effective Methicillin Resistant Staphylococcus aureus (MRSA) inhibitors.

Elucidating the potential neurotoxicity of chiral phenthoate: Molecular insight from experimental and computational studies

Chiral organophosphorus pollutants are existed ubiquitously in the ecological environment, but the enantioselective toxicities of these nerve agents to humans and their molecular bases have not been fully elucidated. Using experimental and computational approaches, this story was to explore the neurotoxic response process of the target acetylcholinesterase (AChE) to chiral phenthoate and further decipher the microscopic mechanism of such toxicological effect at the enantiomeric level. The results showed that the toxic reaction of AChE with chiral phenthoate exhibited significant enantioselectivity, and (R)-phenthoate (K＝1.486 × 10⁵ M^-1) has a bioaffinity for the nerve enzyme nearly three times that of (S)-phenthoate (K＝4.503 × 10⁴ M^-1). Dynamic research outcomes interpreted the wet experiments, and the inherent conformational flexibility of the target enzyme has a great influence on the enantioselective neurotoxicological action processes, especially reflected in the conformational changes of the three key loop regions (i.e. residues His-447, Gly-448, and Tyr-449; residues Gly-122, Phe-123, and Tyr-124; and residues Thr-75, Leu-76, and Tyr-77) around the reaction patch. This was supported by the quantitative results of conformational studies derived from circular dichroism spectroscopy (α-helix: 34.7%→30.2%/31.6%; β-sheet: 23.6%→19.5%/20.7%; turn: 19.2%→22.4%/21.9%; and random coil: 22.5%→27.9%/25.8%). Meanwhile, via analyzing the modes of toxic action and free energies, we can find that (R)-phenthoate has a strong inhibitory effect on the enzymatic activity of AChE, as compared with (S)-phenthoate, and electrostatic energy (-23.79/－17.77 kJ mol^-1) played a critical role in toxicological reactions. These points were the underlying causes of chiral phenthoate displaying different degrees of enantioselective neurotoxicity.

Structure-based discovery of novel inhibitors of Mycobacterium tuberculosis CYP121 from Indonesian natural products

Tuberculosis (TB) continues to be a serious global health threat with the emergence of multidrug-resistant tuberculosis (MDR-TB) and extremely drug-resistant tuberculosis (XDR-TB). There is an urgent need to discover new drugs to deal with the advent of drug-resistant TB variants. This study aims to find new M. tuberculosis CYP121 inhibitors by the screening of Indonesian natural products using the principle of structure-based drug design and discovery. In this work, eight natural compounds isolated from Rhoeo spathacea and Pluchea indica were selected based on their antimycobacterial activity. Derivatives compound were virtually designed from these natural molecules to improve the interaction of ligands with CYP121. Virtual screening of ligands was carried out using AutoDock Vina followed by 50 ns molecular dynamics simulation using YASARA to study the inhibition mechanism of the ligands. Two ligands, i.e., kaempferol (KAE) and its benzyl derivative (KAE3), are identified as the best CYP121 inhibitors based on their binding affinities and adherence to the Lipinski's rule. Results of molecular dynamics simulation indicate that KAE and KAE3 possess a unique inhibitory mechanism against CYP121 that is different from GGJ (control ligand). The control ligand alters the overall dynamics of the receptor, which is indicated by changes in residue flexibility away from CYP121 binding site. Meanwhile, the dynamic changes caused by the binding of KAE and KAE3 are isolated around the binding site of CYP121. These ligands can be developed for further potential biological activities.

Shared Signature Dynamics Tempered by Local Fluctuations Enables Fold Adaptability and Specificity

Recent studies have drawn attention to the evolution of protein dynamics, in addition to sequence and structure, based on the premise structure-encodes-dynamics-encodes-function. Of interest is to understand how functional differentiation is accomplished while maintaining the fold, or how intrinsic dynamics plays out in the evolution of structural variations and functional specificity. We performed a systematic computational analysis of 26,899 proteins belonging to 116 CATH superfamilies. Characterizing cooperative mechanisms and convergent/divergent features that underlie the shared/differentiated dynamics of family members required a methodology that lends itself to efficient analyses of large ensembles of proteins. We therefore introduced, SignDy, an integrated pipeline for evaluating the signature dynamics of families based on elastic network models. Our analysis confirmed that family members share conserved, highly cooperative (global) modes of motion. Importantly, our analysis discloses a subset of motions that sharply distinguishes subfamilies, which lie in a low-to-intermediate frequency regime of the mode spectrum. This regime has maximal impact on functional differentiation of families into subfamilies, while being evolutionarily conserved among subfamily members. Notably, the high-frequency end of the spectrum also reveals evolutionary conserved features across and within subfamilies; but in sharp contrast to global motions, high-frequency modes are minimally collective. Modulation of robust/conserved global dynamics by low-to-intermediate frequency fluctuations thus emerges as a versatile mechanism ensuring the adaptability of selected folds and the specificity of their subfamilies. SignDy further allows for dynamics-based categorization as a new layer of information relevant to distinctive mechanisms of action of subfamilies, beyond sequence or structural classifications.

Exploring the Ligand Binding/Unbinding Pathway by Selectively Enhanced Sampling of Ligand in a Protein-Ligand Complex

Understanding the protein-ligand binding is of fundamental biological interest and is essential for structure-based drug design. The difficulty in capturing the dynamic process, however, poses a great challenge for current experimental and theoretical simulation techniques. A selective integrated-tempering-sampling molecular dynamics (SITSMD) method offering an option for selectively enhanced sampling of the ligand in a protein-ligand complex was utilized to quantitatively illuminate the binding of benzamidine to the wild-type trypsin protease and its two mutants (S214E and S214K). The SITSMD simulations could produce consistent results as the extensive conventional MD simulation and gave additional insights into the binding pathway for the test protein-ligand complex system using significantly saved computational resource and time, indicating the potential of such a method in investigating protein-ligand binding. Additionally, the simulations identified the different roles of trypsin-benzamidine van der Waals (vdW) and electrostatic interactions in the binding: the former interaction works as the driving force for dragging the benzamidine close to the native binding pocket, and the latter interaction mainly contributes to stabilizing the benzamidine inside the pocket. The S214E mutation introduces more favorable electrostatic interactions, and as a result, both vdW and electrostatic interactions drive the benzamidine binding, lowering the binding and unbinding free energy barrier. In contrast, the S214K mutation prohibits the binding of the benzamidine to the native ligand binding pocket by introducing disliked charge-charge interactions. In summary, these findings suggest that the change in specific residues could modify the protein druggability, including the binding kinetics and thermodynamics.

Methyl-Based NMR Spectroscopy Methods for Uncovering Structural Dynamics in Large Proteins and Protein Complexes

NMR spectroscopy is particularly adept at site-specifically monitoring dynamic processes in proteins, such as protein folding, domain movements, ligand binding, and side-chain rotations. By coupling the favorable spectroscopic properties of highly dynamic side-chain methyl groups with transverse-relaxation-optimized spectroscopy (TROSY), it is now possible to routinely study such dynamic processes in high-molecular-weight proteins and complexes approaching 1 MDa. In this Perspective, we describe many elegant methyl-based NMR experiments that probe slow (second) to fast (picosecond) dynamics in large systems. To demonstrate the power of these methods, we also provide interesting examples of studies that utilized each methyl-based NMR technique to uncover functionally important dynamics. In many cases, the NMR experiments are paired with site-directed mutagenesis and/or other biochemical assays to put the dynamics and function into context. Our vision of the future of structural biology involves pairing methyl-based NMR spectroscopy with biochemical studies to advance our knowledge of the motions large proteins and macromolecular complexes use to choreograph complex functions. Such studies will be essential in elucidating the critical structural dynamics that underlie function and characterizing alterations in these processes that can lead to human disease.

Normal mode analysis as a method to derive protein dynamics information from the Protein Data Bank

Normal mode analysis (NMA) can facilitate quick and systematic investigation of protein dynamics using data from the Protein Data Bank (PDB). We developed an elastic network model-based NMA program using dihedral angles as independent variables. Compared to the NMA programs that use Cartesian coordinates as independent variables, key attributes of the proposed program are as follows: (1) chain connectivity related to the folding pattern of a polypeptide chain is naturally embedded in the model; (2) the full-atom system is acceptable, and owing to a considerably smaller number of independent variables, the PDB data can be used without further manipulation; (3) the number of variables can be easily reduced by some of the rotatable dihedral angles; (4) the PDB data for any molecule besides proteins can be considered without coarse-graining; and (5) individual motions of constituent subunits and ligand molecules can be easily decomposed into external and internal motions to examine their mutual and intrinsic motions. Its performance is illustrated with an example of a DNA-binding allosteric protein, a catabolite activator protein. In particular, the focus is on the conformational change upon cAMP and DNA binding, and on the communication between their binding sites remotely located from each other. In this illustration, NMA creates a vivid picture of the protein dynamics at various levels of the structures, i.e., atoms, residues, secondary structures, domains, subunits, and the complete system, including DNA and cAMP. Comparative studies of the specific protein in different states, e.g., apo- and holo-conformations, and free and complexed configurations, provide useful information for studying structurally and functionally important aspects of the protein.

Intrinsic dynamics study identifies two amino acids of TIMP-1 critical for its LRP-1-mediated endocytosis in neurons

The tissue inhibitor of metalloproteinases-1 (TIMP-1) exerts inhibitory activity against matrix metalloproteinases and cytokine-like effects. We previously showed that TIMP-1 reduces neurite outgrowth in mouse cortical neurons and that this cytokine-like effect depends on TIMP-1 endocytosis mediated by the low-density lipoprotein receptor-related protein-1 (LRP-1). To gain insight into the interaction between TIMP-1 and LRP-1, we considered conformational changes that occur when a ligand binds to its receptor. TIMP-1 conformational changes have been studied using biomolecular simulations, and our results provide evidence for a hinge region that is critical for the protein movement between the N- and C-terminal TIMP-1 domains. In silico mutants have been proposed on residues F12 and K47, which are located in the hinge region. Biological analyses of these mutants show that F12A or K47A mutation does not alter MMP inhibitory activity but impairs the effect of TIMP-1 on neurite outgrowth. Interestingly, these mutants bind to LRP-1 but are not endocytosed. We conclude that the intrinsic dynamics of TIMP-1 are not involved in its binding to LRP-1 but rather in the initiation of endocytosis and associated biological effects.

Dynamics, flexibility, and allostery in molecular chaperonins

The chaperonins are a family of molecular chaperones present in all three kingdoms of life. They are classified into Group I and Group II. Group I consists of the bacterial variants (GroEL) and the eukaryotic ones from mitochondria and chloroplasts (Hsp60), while Group II consists of the archaeal (thermosomes) and eukaryotic cytosolic variants (CCT or TRiC). Both groups assemble into a dual ring structure, with each ring providing a protective folding chamber for nascent and denatured proteins. Their functional cycle is powered by ATP binding and hydrolysis, which drives a series of structural rearrangements that enable encapsulation and subsequent release of the substrate protein. Chaperonins have elaborate allosteric mechanisms to regulate their functional cycle. Long-range negative cooperativity between the two rings ensures alternation of the folding chambers. Positive intra-ring cooperativity, which facilitates concerted conformational transitions within the protein subunits of one ring, has only been demonstrated for Group I chaperonins. In this review, we describe our present understanding of the underlying mechanisms and the structure-function relationships in these complex protein systems with a particular focus on the structural dynamics, allostery, and associated conformational rearrangements.

Increasing the sampling efficiency of protein conformational transition using velocity-scaling optimized hybrid explicit/implicit solvent REMD simulation

The application of temperature replica exchange molecular dynamics (REMD) simulation on protein motion is limited by its huge requirement of computational resource, particularly when explicit solvent model is implemented. In the previous study, we developed a velocity-scaling optimized hybrid explicit/implicit solvent REMD method with the hope to reduce the temperature (replica) number on the premise of maintaining high sampling efficiency. In this study, we utilized this method to characterize and energetically identify the conformational transition pathway of a protein model, the N-terminal domain of calmodulin. In comparison to the standard explicit solvent REMD simulation, the hybrid REMD is much less computationally expensive but, meanwhile, gives accurate evaluation of the structural and thermodynamic properties of the conformational transition which are in well agreement with the standard REMD simulation. Therefore, the hybrid REMD could highly increase the computational efficiency and thus expand the application of REMD simulation to larger-size protein systems.

Cooperative Conformational Regulation in N-Heterocyclic Fluorohydrins

Seven-membered N-heterocycles are flexible ring structures, and their conformational control is important to their bioactivity. Our prior work shows that stereoselective monofluorination, if installed diastereoselectively, can bias a seven-membered, substituted azepane ring to one major conformation. However, multiple fluorination may not provide as much conformational bias due to conflicting effects. Here we show in our model azepane system that fluorohydrins can confer strong conformational bias if the relative configuration of the fluorine and hydroxy substitutent is appropriate to enable cooperative conformational control.

Toward in silico biomolecular manipulation through static modes: atomic scale characterization of HIV-1 protease flexibility

Probing biomolecular flexibility with atomic-scale resolution is a challenging task in current computational biology for fundamental understanding and prediction of biomolecular interactions and associated functions. This paper makes use of the static mode method to study HIV-1 protease considered as a model system to investigate the full biomolecular flexibility at the atomic scale, the screening of active site biomechanical properties, the blind prediction of allosteric sites, and the design of multisite strategies to target deformations of interest. Relying on this single calculation run of static modes, we demonstrate that in silico predictive design of an infinite set of complex excitation fields is reachable, thanks to the storage of the static modes in a data bank that can be used to mimic single or multiatom contact and efficiently anticipate conformational changes arising from external stimuli. All along this article, we compare our results to data previously published and propose a guideline for efficient, predictive, and custom in silico experiments.

Exploring transition pathway and free-energy profile of large-scale protein conformational change by combining normal mode analysis and umbrella sampling molecular dynamics

Large-scale conformational changes of proteins are usually associated with the binding of ligands. Because the conformational changes are often related to the biological functions of proteins, understanding the molecular mechanisms of these motions and the effects of ligand binding becomes very necessary. In the present study, we use the combination of normal-mode analysis and umbrella sampling molecular dynamics simulation to delineate the atomically detailed conformational transition pathways and the associated free-energy landscapes for three well-known protein systems, viz., adenylate kinase (AdK), calmodulin (CaM), and p38α kinase in the absence and presence of respective ligands. For each protein under study, the transient conformations along the conformational transition pathway and thermodynamic observables are in agreement with experimentally and computationally determined ones. The calculated free-energy profiles reveal that AdK and CaM are intrinsically flexible in structures without obvious energy barrier, and their ligand binding shifts the equilibrium from the ligand-free to ligand-bound conformation (population shift mechanism). In contrast, the ligand binding to p38α leads to a large change in free-energy barrier (ΔΔG ≈ 7 kcal/mol), promoting the transition from DFG-in to DFG-out conformation (induced fit mechanism). Moreover, the effect of the protonation of D168 on the conformational change of p38α is also studied, which reduces the free-energy difference between the two functional states of p38α and thus further facilitates the conformational interconversion. Therefore, the present study suggests that the detailed mechanism of ligand binding and the associated conformational transition is not uniform for all kinds of proteins but correlated to their respective biological functions.

Screening and evaluation of small organic molecules as ClpB inhibitors and potential antimicrobials

Inhibition of ClpB, the bacterial representative of the heat-shock protein 100 family that is associated with virulence of several pathogens, could be an effective strategy to develop new antimicrobial agents. Using a high-throughput screening method, we have identified several compounds that bind to different conformations of ClpB and analyzed their effect on the ATPase and chaperone activities of the protein. Two of them inhibit these functional properties as well as the growth of Gram negative bacteria (E. coli), displaying antimicrobial activity under thermal or oxidative stress conditions. This activity is abolished upon deletion of ClpB, indicating that the action of these compounds is related to the stress cellular response in which ClpB is involved. Moreover, their moderate toxicity in human cell lines suggests that they might provide promising leads against bacterial growth.

Normal mode analysis based on an elastic network model for biomolecules in the Protein Data Bank, which uses dihedral angles as independent variables

We have developed a computer program, named PDBETA, that performs normal mode analysis (NMA) based on an elastic network model that uses dihedral angles as independent variables. Taking advantage of the relatively small number of degrees of freedom required to describe a molecular structure in dihedral angle space and a simple potential-energy function independent of atom types, we aimed to develop a program applicable to a full-atom system of any molecule in the Protein Data Bank (PDB). The algorithm for NMA used in PDBETA is the same as the computer program FEDER/2, developed previously. Therefore, the main challenge in developing PDBETA was to find a method that can automatically convert PDB data into molecular structure information in dihedral angle space. Here, we illustrate the performance of PDBETA with a protein-DNA complex, a protein-tRNA complex, and some non-protein small molecules, and show that the atomic fluctuations calculated by PDBETA reproduce the temperature factor data of these molecules in the PDB. A comparison was also made with elastic-network-model based NMA in a Cartesian-coordinate system.

Investigating the structural origin of trpzip2 temperature dependent unfolding fluorescence line shape based on a Markov state model simulation

Vibrationally resolved fluorescence spectra of the β-hairpin trpzip2 peptide at two temperatures as well as during a T-jump unfolding process are simulated on the basis of a combination of Markov state models and quantum chemistry schemes. The broad asymmetric spectral line shape feature is reproduced by considering the exciton-phonon couplings. The temperature dependent red shift observed in the experiment has been attributed to the state population changes of specific chromophores. Through further theoretical study, it is found that both the environment's electric field and the chromophores' geometry distortions are responsible for tryptophan fluorescence shift.

A dynamic model of long-range conformational adaptations triggered by nucleotide binding in GroEL-GroES

The molecular chaperone, GroEL, essential for correct protein folding in E. coli, is composed of 14 identical subunits organized in two interacting rings, each providing a folding chamber for non-native substrate proteins. The oligomeric assembly shows positive cooperativity within each ring and negative cooperativity between the rings. Although it is well known that ATP and long-range allosteric interactions drive the functional cycle of GroEL, an atomic resolution view of how ligand binding modulates conformational adaptations over long distances remains a major challenge. Moreover, little is known on the relation between equilibrium dynamics at physiological temperatures and the allosteric transitions in GroEL. Here we present multiple all-atom molecular dynamics simulations of the GroEL-GroES assemblies at different stages of the functional cycle. Combined with an extensive analysis of the complete set of experimentally available structures, principal component analysis and conformer plots, we provide an explicit evaluation of the accessible conformational space of unliganded GroEL. Our results suggest the presence of pre-existing conformers at the equatorial domain level, and a shift of the conformational ensemble upon ATP-binding. At the inter-ring interface the simulations capture a remarkable offset motion of helix D triggered by ATP-binding to the folding active ring. The reorientation of helix D, previously only observed upon GroES association, correlates with a change of the internal dynamics in the opposite ring. This work contributes to the understanding of the molecular mechanisms in GroEL and highlights the ability of all-atom MD simulations to model long-range structural changes and allosteric events in large systems.