Enhanced Sampling Molecular Dynamics Simulations Research Articles

Modeling the Hydrolysis of Iron-Sulfur Clusters.

Iron-sulfur (FeS) clusters are essential metal cofactors involved in a wide variety of biological functions. Their catalytic efficiency, biosynthesis, and regulation depend on FeS stability in aqueous solution. Here, molecular modeling is used to investigate the hydrolysis of an oxidized (ferric) mononuclear FeS cluster by bare dissociation and water substitution mechanisms in neutral and acidic solution. First, approximate electronic structure descriptions of FeS reactions by density functional theory are validated against high-level wave function CCSD(T) calculations. Solvation contributions are included by an all-atom model with hybrid quantum chemical/molecular mechanical (QM/MM) potentials and enhanced sampling molecular dynamics simulations. The free energy profile obtained for FeS cluster hydrolysis indicates that the hybrid functional M06 together with an implicit solvent correction capture the most important aspects of FeS cluster reactivity in aqueous solution. Then, 20 reaction channels leading to two consecutive Fe-S bond ruptures were explored with this calibrated model. For all protonation states, nucleophilic substitution with concerted bond breaking and forming to iron is the preferred mechanism, both kinetic and thermodynamically. In neutral solution, proton transfer from water to the sulfur leaving group is also concerted. Dissociative reactions show higher barriers and will not be relevant for FeS reactivity when exposed to solvent. These hydrolysis mechanisms may help to explain the stability and catalytic mechanisms of FeS clusters of multiple sizes and proteins.

Modulating Ligand Dissociation through Methyl Isomerism in Accessory Sites: Binding of Retinol to Cellular Carriers.

Due to the poor aqueous solubility of retinoids, evolution has tuned their binding to cellular proteins to address specialized physiological roles by modulating uptake, storage, and delivery to specific targets. With the aim to disentangle the structure-function relationships in these proteins and disclose clues for engineering selective carriers, the binding mechanism of the two most abundant retinol-binding isoforms was explored by using enhanced sampling molecular dynamics simulations and surface plasmon resonance. The distinctive dynamics of the entry portal site in the holo species was crucial to modulate retinol dissociation. Remarkably, this process is controlled to a large extent by the replacement of Ile by Leu in the two isoforms, thus suggesting that fine control of ligand release can be achieved through a rigorous selection of conservative mutations in accessory sites.

Exploring Conformational Transitions and Free-Energy Profiles of Proton-Coupled Oligopeptide Transporters.

Proteins involved in peptide uptake and transport belong to the proton-coupled oligopeptide transporter (POT) family. Crystal structures of POT family members reveal a common fold consisting of two domains of six transmembrane α helices that come together to form a "V" shaped transporter with a central substrate binding site. Proton-coupled oligopeptide transporters operate through an alternate access mechanism, where the membrane transporter undergoes global conformational changes, alternating between inward-facing (IF), outward-facing (OF), and occluded (OC) states. Conformational transitions are promoted by proton and ligand binding; however, due to the absence of crystallographic models of the outward-open state, the role of H+ and ligands is still not fully understood. To provide a comprehensive picture of the POT conformational equilibrium, conventional and enhanced sampling molecular dynamics simulations of PepTst in the presence or absence of ligand and protonation were performed. Free-energy profiles of the conformational variability of PepTst were obtained from microseconds of adaptive biasing force (ABF) simulations. Our results reveal that both proton and ligand significantly change the conformational free-energy landscape. In the absence of ligand and protonation, only transitions involving IF and OC states are allowed. After protonation of the residue Glu300, the wider free-energy well for Glu300 protonated PepTst indicates a greater conformational variability relative to the apo system, and OF conformations became accessible. For the Glu300 protonated Holo-PepTst, the presence of a second free-energy minimum suggests that OF conformations are not only accessible, but also stable. The differences in the free-energy profiles demonstrate that transitions toward outward-facing conformation occur only after protonation, which is likely the first step in the mechanism of peptide transport. Our extensive ABF simulations provide a fully atomic description of all states of the transport process, offering a model for the alternating access mechanism and how protonation and ligand control the conformational changes.

Ring Puckering Landscapes of Glycosaminoglycan-Related Monosaccharides from Molecular Dynamics Simulations.

The conformational flexibility of the glycosaminoglycans (GAGs) is known to be key in their binding and biological function, for example in regulating coagulation and cell growth. In this work, we employ enhanced sampling molecular dynamics simulations to probe the ring conformations of GAG-related monosaccharides, including a range of acetylated and sulfated GAG residues. We first perform unbiased MD simulations of glucose anomers and the epimers glucuronate and iduronate. These calculations indicate that in some cases, an excess of 15 μs is required for adequate sampling of ring pucker due to the high energy barriers between states. However, by applying our recently developed msesMD simulation method (multidimensional swarm-enhanced sampling molecular dynamics), we were able to quantitatively and rapidly reproduce these ring pucker landscapes. From msesMD simulations, the puckering free energy profiles were then compared for 15 further monosaccharides related to GAGs; this includes to our knowledge the first simulation study of sulfation effects on β-GalNAc ring puckering. For the force field employed, we find that in general the calculated pucker free energy profiles for sulfated sugars were similar to the corresponding unsulfated profiles. This accords with recent experimental studies suggesting that variation in ring pucker of sulfated GAG residues is primarily dictated by interactions with surrounding residues rather than by intrinsic conformational preference. As an exception to this, however, we predict that 4-O-sulfation of β-GalNAc leads to reduced ring rigidity, with a significant lowering in energy of the 1C4 ring conformation; this observation may have implications for understanding the structural basis of the biological function of β-GalNAc-containing glycosaminoglycans such as dermatan sulfate.

Accelerated Molecular Dynamics Applied to the Peptaibol Folding Problem.

The use of enhanced sampling molecular dynamics simulations to facilitate the folding of proteins is a relatively new approach which has quickly gained momentum in recent years. Accelerated molecular dynamics (aMD) can elucidate the dynamic path from the unfolded state to the near-native state, “flattened” by introducing a non-negative boost to the potential. Alamethicin F30/3 (Alm F30/3), chosen in this study, belongs to the class of peptaibols that are 7–20 residue long, non-ribosomally synthesized, amphipathic molecules that show interesting membrane perturbing activity. The recent studies undertaken on the Alm molecules and their transmembrane channels have been reviewed. Three consecutive simulations of ~900 ns each were carried out where N-terminal folding could be observed within the first 100 ns, while C-terminal folding could only be achieved almost after 800 ns. It took ~1 μs to attain the near-native conformation with stronger potential boost which may take several μs worth of classical MD to produce the same results. The Alm F30/3 hexamer channel was also simulated in an E. coli mimicking membrane under an external electric field that correlates with previous experiments. It can be concluded that aMD simulation techniques are suited to elucidate peptaibol structures and to understand their folding dynamics.

The Highly Recurrent PP2A Aα-Subunit Mutation P179R Alters Protein Structure and Impairs PP2A Enzyme Function to Promote Endometrial Tumorigenesis.

Somatic mutation of the protein phosphatase 2A (PP2A) Aα-subunit gene PPP2R1A is highly prevalent in high-grade endometrial carcinoma. The structural, molecular, and biological basis by which the most recurrent endometrial carcinoma-specific mutation site P179 facilitates features of endometrial carcinoma malignancy has yet to be fully determined. Here, we used a series of structural, biochemical, and biological approaches to investigate the impact of the P179R missense mutation on PP2A function. Enhanced sampling molecular dynamics simulations showed that arginine-to-proline substitution at the P179 residue changes the protein's stable conformation profile. A crystal structure of the tumor-derived PP2A mutant revealed marked changes in A-subunit conformation. Binding to the PP2A catalytic subunit was significantly impaired, disrupting holoenzyme formation and enzymatic activity. Cancer cells were dependent on PP2A disruption for sustained tumorigenic potential, and restoration of wild-type Aα in a patient-derived P179R-mutant cell line restored enzyme function and significantly attenuated tumorigenesis and metastasis in vivo. Furthermore, small molecule-mediated therapeutic reactivation of PP2A significantly inhibited tumorigenicity in vivo. These outcomes implicate PP2A functional inactivation as a critical component of high-grade endometrial carcinoma disease pathogenesis. Moreover, they highlight PP2A reactivation as a potential therapeutic strategy for patients who harbor P179R PPP2R1A mutations. SIGNIFICANCE: This study characterizes a highly recurrent, disease-specific PP2A PPP2R1A mutation as a driver of endometrial carcinoma and a target for novel therapeutic development.See related commentary by Haines and Huang, p. 4009.

Parallel G-triplexes and G-hairpins as potential transitory ensembles in the folding of parallel-stranded DNA G-Quadruplexes.

Guanine quadruplexes (G4s) are non-canonical nucleic acids structures common in important genomic regions. Parallel-stranded G4 folds are the most abundant, but their folding mechanism is not fully understood. Recent research highlighted that G4 DNA molecules fold via kinetic partitioning mechanism dominated by competition amongst diverse long-living G4 folds. The role of other intermediate species such as parallel G-triplexes and G-hairpins in the folding process has been a matter of debate. Here, we use standard and enhanced-sampling molecular dynamics simulations (total length of ∼0.9 ms) to study these potential folding intermediates. We suggest that parallel G-triplex per se is rather an unstable species that is in local equilibrium with a broad ensemble of triplex-like structures. The equilibrium is shifted to well-structured G-triplex by stacked aromatic ligand and to a lesser extent by flanking duplexes or nucleotides. Next, we study propeller loop formation in GGGAGGGAGGG, GGGAGGG and GGGTTAGGG sequences. We identify multiple folding pathways from different unfolded and misfolded structures leading towards an ensemble of intermediates called cross-like structures (cross-hairpins), thus providing atomistic level of description of the single-molecule folding events. In summary, the parallel G-triplex is a possible, but not mandatory short-living (transitory) intermediate in the folding of parallel-stranded G4.

Unbinding of Translocator Protein 18 kDa (TSPO) Ligands: From in Vitro Residence Time to in Vivo Efficacy via in Silico Simulations.

Translocator protein 18 kDa (TSPO) is a validated pharmacological target for the development of new treatments for neurological disorders. N,N-Dialkyl-2-phenylindol-3-ylglyoxylamides (PIGAs) are effective TSPO modulators and potentially useful therapeutics for the treatment of anxiety, central nervous system pathologies featuring astrocyte loss, and inflammatory-based neuropathologies. For this class of compounds, no correlation exists between the TSPO binding affinity and the corresponding functional efficacy. Rather, their biological effectiveness correlates with the kinetics of the unbinding events and more specifically with the residence time (RT). So far, the structural reasons for the different recorded RT of congeneric PIGAs remain elusive. Here, to understand the different kinetics of PIGAs, their unbinding paths were studied by employing enhanced-sampling molecular dynamics simulations. Results of these studies revealed how subtle structural differences between PIGAs have a substantial effect on the unbinding energetics. In particular, during the egress from the TSPO binding site, slow-dissociating PIGAs find tight interactions with the protein LP1 region thereby determining a long RT. Further support to these findings was achieved by in vivo studies, which demonstrated how the anxiolytic effect observed for the inspected PIGAs correlated with their RT to TSPO.

Advances in enhanced sampling molecular dynamics simulations for biomolecules

Molecular dynamics simulation has emerged as a powerful computational tool for studying biomolecules as it can provide atomic insights into the conformational transitions involved in biological functions. However, when applied to complex biological macromolecules, the conformational sampling ability of conventional molecular dynamics is limited by the rugged free energy landscapes, leading to inherent timescale gaps between molecular dynamics simulations and real biological processes. To address this issue, several advanced enhanced sampling methods have been proposed to improve the sampling efficiency in molecular dynamics. In this review, the theoretical basis, practical applications, and recent improvements of both constraint and unconstrained enhanced sampling methods are summarized. Further-more, the combined utilizations of different enhanced sampling methods that take advantage of both approaches are also briefly discussed.

Impact of branching on the conformational heterogeneity of the lipopolysaccharide from Klebsiella pneumoniae: Implications for vaccine design

Resistance of Klebsiella pneumoniae (KP) to antibiotics has motivated the development of an efficacious KP human vaccine that would not be subject to antibiotic resistance. Klebsiella lipopolysaccharide (LPS) associated O polysaccharide (OPS) types have provoked broad interest as a vaccine antigen as there are only 4 that predominate worldwide (O1, O2a, O3, O5). Klebsiella O1 and O2 OPS are polygalactans that share a common D-Gal-I structure, for which a variant D-Gal-III was recently discovered. To understand the potential impact of this variability on antigenicity, a detailed molecular picture of the conformational differences associated with the addition of the D-Gal-III (1 → 4)-α-Galp branch is presented using enhanced-sampling molecular dynamics simulations. In D-Gal-I two major conformational states are observed while the presence of the 1 → 4 branch in D-Gal-III resulted in only a single dominant extended state. Stabilization of the more folded states in D-Gal-I is due to a O4-H⋯O2 hydrogen bond in the linear backbone that cannot occur in D-Gal-III as the O4 is in the Galp(1 → 4)Galp glycosidic linkage. The impact of branching in D-Gal-III also significantly decreases the accessibility of the monosaccharides in the linear backbone region of D-Gal-I, while the accessibility of the terminal D-Gal-II region of the OPS is not substantially altered. The present results suggest that a vaccine that targets both the D-Gal-I and D-Gal-III LPS can be developed by using D-Gal-III as the antigen combined with cross-reactivity experiments using the Gal-II polysaccharide to assure that this region of the LPS is the primary epitope of the antigen.

Conformational dynamics and interfacial interactions of peptide-appended pillar[5]arene water channels in biomimetic membranes.

Peptide appended pillar[5]arene (PAP) is an artificial water channel resembling biological water channel proteins, which has shown a significant potential for designing bioinspired water purification systems. Given that PAP channels need to be incorporated at a high density in membrane matrices, it is critical to examine the role of channel-channel and channel-membrane interactions in governing the structural and functional characteristics of channels. To resolve the atomic-scale details of these interactions, we have carried out atomistic molecular dynamics (MD) simulations of multiple PAP channels inserted in a lipid or a block-copolymer (BCP) membrane matrix. Classical MD simulations on a sub-microsecond timescale showed clustering of channels only in the lipid membrane, but enhanced sampling MD simulations showed thermodynamically-favorable dimerized states of channels in both lipid and BCP membranes. The dimerized configurations of channels, with an extensive buried surface area, were stabilized via interactions between the aromatic groups in the peptide arms of neighboring channels. The conformational metrics characterizing the orientational and structural changes in channels revealed a higher flexibility in the lipid membrane as opposed to the BCP membrane although hydrogen bonds between the channel and the membrane molecules were not a major contributor to the stability of channels in the BCP membrane. We also found that the channels undergo wetting/dewetting transitions in both lipid and BCP membranes with a marginally higher probability of undergoing a dewetting transition in the BCP membrane. Collectively, these results highlight the role of channel dynamics in governing channel-channel and channel-membrane interfacial interactions, and provide atomic-scale insights needed to design stable and functional biomimetic membranes for efficient separations.

Deciphering the Evolution of Cephalosporin Resistance to Ceftolozane-Tazobactam in Pseudomonas aeruginosa.

Pseudomonas aeruginosa produces a class C β-lactamase (e.g., PDC-3) that robustly hydrolyzes early generation cephalosporins often at the diffusion limit; therefore, bacteria possessing these β-lactamases are resistant to many β-lactam antibiotics. In response to this significant clinical threat, ceftolozane, a 3' aminopyrazolium cephalosporin, was developed. Combined with tazobactam, ceftolozane promised to be effective against multidrug-resistant P. aeruginosa Alarmingly, Ω-loop variants of the PDC β-lactamase (V213A, G216R, E221K, E221G, and Y223H) were identified in ceftolozane/tazobactam-resistant P. aeruginosa clinical isolates. Herein, we demonstrate that the Escherichia coli strain expressing the E221K variant of PDC-3 had the highest minimum inhibitory concentrations (MICs) against a panel of β-lactam antibiotics, including ceftolozane and ceftazidime, a cephalosporin that differs in structure largely in the R2 side chain. The kcat values of the E221K variant for both substrates were equivalent, whereas the Km for ceftolozane (341 ± 64 µM) was higher than that for ceftazidime (174 ± 20 µM). Timed mass spectrometry, thermal stability, and equilibrium unfolding studies revealed key mechanistic insights. Enhanced sampling molecular dynamics simulations identified conformational changes in the E221K variant Ω-loop, where a hidden pocket adjacent to the catalytic site opens and stabilizes ceftolozane for efficient hydrolysis. Encouragingly, the diazabicyclooctane β-lactamase inhibitor avibactam restored susceptibility to ceftolozane and ceftazidime in cells producing the E221K variant. In addition, a boronic acid transition state inhibitor, LP-06, lowered the ceftolozane and ceftazidime MICs by 8-fold for the E221K-expressing strain. Understanding these structural changes in evolutionarily selected variants is critical toward designing effective β-lactam/β-lactamase inhibitor therapies for P. aeruginosa infections.IMPORTANCE The presence of β-lactamases (e.g., PDC-3) that have naturally evolved and acquired the ability to break down β-lactam antibiotics (e.g., ceftazidime and ceftolozane) leads to highly resistant and potentially lethal Pseudomonas aeruginosa infections. We show that wild-type PDC-3 β-lactamase forms an acyl enzyme complex with ceftazidime, but it cannot accommodate the structurally similar ceftolozane that has a longer R2 side chain with increased basicity. A single amino acid substitution from a glutamate to a lysine at position 221 in PDC-3 (E221K) causes the tyrosine residue at 223 to adopt a new position poised for efficient hydrolysis of both cephalosporins. The importance of the mechanism of action of the E221K variant, in particular, is underscored by its evolutionary recurrences in multiple bacterial species. Understanding the biochemical and molecular basis for resistance is key to designing effective therapies and developing new β-lactam/β-lactamase inhibitor combinations.

Effect of Ca2+ on the promiscuous target-protein binding of calmodulin.

Calmodulin (CaM) is a calcium sensing protein that regulates the function of a large number of proteins, thus playing a crucial part in many cell signaling pathways. CaM has the ability to bind more than 300 different target peptides in a Ca2+-dependent manner, mainly through the exposure of hydrophobic residues. How CaM can bind a large number of targets while retaining some selectivity is a fascinating open question. Here, we explore the mechanism of CaM selective promiscuity for selected target proteins. Analyzing enhanced sampling molecular dynamics simulations of Ca2+-bound and Ca2+-free CaM via spectral clustering has allowed us to identify distinct conformational states, characterized by interhelical angles, secondary structure determinants and the solvent exposure of specific residues. We searched for indicators of conformational selection by mapping solvent exposure of residues in these conformational states to contacts in structures of CaM/target peptide complexes. We thereby identified CaM states involved in various binding classes arranged along a depth binding gradient. Binding Ca2+ modifies the accessible hydrophobic surface of the two lobes and allows for deeper binding. Apo CaM indeed shows shallow binding involving predominantly polar and charged residues. Furthermore, binding to the C-terminal lobe of CaM appears selective and involves specific conformational states that can facilitate deep binding to target proteins, while binding to the N-terminal lobe appears to happen through a more flexible mechanism. Thus the long-ranged electrostatic interactions of the charged residues of the N-terminal lobe of CaM may initiate binding, while the short-ranged interactions of hydrophobic residues in the C-terminal lobe of CaM may account for selectivity. This work furthers our understanding of the mechanism of CaM binding and selectivity to different target proteins and paves the way towards a comprehensive model of CaM selectivity.

Mechanism of Mg2+-Accompanied Product Release in Sugar Nucleotidyltransferases

The nucleotidyl transfer reaction, catalyzed by sugar nucleotidyltransferases (SNTs), is assisted by two active site Mg2+ ions. While studying this reaction using X-ray crystallography, we captured snapshots of the pyrophosphate (product) as it exits along a pocket. Surprisingly, one of the active site Mg2+ ions remains coordinated to the exiting pyrophosphate. This hints at the participation of Mg2+ in the process of product release, besides its role in catalyzing nucleotidyl transfer. These observations are further supported by enhanced sampling molecular dynamics simulations. Free energy computations suggest that the product release is likely to be rate limiting in SNTs, and the origin of the high free energy barrier for product release could be traced back to the "slow" conformational change of an Arg residue at the exit end of the pocket. These results establish a dual role for Mg2+, and propose a general mechanism of product release during the nucleotidyl transfer by SNTs.

Conformational Heterogeneity of the HIV Envelope Glycan Shield

To better understand the conformational properties of the glycan shield covering the surface of the HIV gp120/gp41 envelope (Env) trimer, and how the glycan shield impacts the accessibility of the underlying protein surface, we performed enhanced sampling molecular dynamics (MD) simulations of a model glycosylated HIV Env protein and related systems. Our simulation studies revealed a conformationally heterogeneous glycan shield with a network of glycan-glycan interactions more extensive than those observed to date. We found that partial preorganization of the glycans potentially favors binding by established broadly neutralizing antibodies; omission of several specific glycans could increase the accessibility of other glycans or regions of the protein surface to antibody or CD4 receptor binding; the number of glycans that can potentially interact with known antibodies is larger than that observed in experimental studies; and specific glycan conformations can maximize or minimize interactions with individual antibodies. More broadly, the enhanced sampling MD simulations described here provide a valuable tool to guide the engineering of specific Env glycoforms for HIV vaccine design.

Bacterial Outer Membrane Porins as Electrostatic Nanosieves: Exploring Transport Rules of Small Polar Molecules.

Transport of molecules through biological membranes is a fundamental process in biology, facilitated by selective channels and general pores. The architecture of some outer membrane pores in Gram-negative bacteria, common to other eukaryotic pores, suggests them as prototypes of electrostatically regulated nanosieve devices. In this study, we sensed the internal electrostatics of the two most abundant outer membrane channels of Escherichia coli, using norfloxacin as a dipolar probe in single molecule electrophysiology. The voltage dependence of the association rate constant of norfloxacin interacting with these nanochannels follows an exponential trend, unexpected for neutral molecules. We combined electrophysiology, channel mutagenesis, and enhanced sampling molecular dynamics simulations to explain this molecular mechanism. Voltage and temperature dependent ion current measurements allowed us to quantify the transversal electric field inside the channel as well as the distance where the applied potential drops. Finally, we proposed a general model for transport of polar molecules through these electrostatic nanosieves. Our model helps to further understand the basis for permeability in Gram-negative pathogens, contributing to fill in the innovation gap that has limited the discovery of effective antibiotics in the last 20 years.

Energetics and structural characterization of the "DFG-flip" conformational transition of B-RAF kinase: a SITS molecular dynamics study.

B-RAF protein kinase is a promising target to treat malignant melanoma. The kinase activity of B-RAF is regulated by a "DFG-flip" conformational transition between functional DFG-in and DFG-out states. The difficulty in resolving the activation loop in crystal structures and the even greater difficulty in experimentally capturing high-energy-level transient structures render elusive the molecular mechanism of the B-RAF functional conformational transition. Here, a homology modeling technique and an enhanced sampling molecular dynamics simulation were used to identify and energetically characterize the conformational transition pathway of B-RAF on a multi-dimensional free-energy landscape. The results reveal that the conformational transition is a two-state transition, with the evaluated free-energy barrier comparable to those of other kinds of kinases as reported in the previous literature. Hydrophobic interactions between activation loop and neighboring segments are suggested to dominate the conformational transition and determine the free-energy barrier. The detailed analysis of hydrophobic interactions involved in the conformational transition may show a suitable pathway for the development of the B-RAF inhibitor.

Stereoselective binding of agonists to the β2-adrenergic receptor: insights into molecular details and thermodynamics from molecular dynamics simulations.

The β2-adrenergic receptor (β2-AR) is one of the most studied G-protein-coupled receptors. When interacting with ligand molecules, it exhibits a binding characteristic that is strongly dependent on ligand stereoconfiguration. In particular, many experimental and theoretical studies confirmed that stereoisomers of an important β2-AR agonist, fenoterol, are associated with diverse mechanisms of binding and activation of β2-AR. The objective of the present study was to explore the stereoselective binding of fenoterol to β2-AR through the application of an advanced computational methodology based on enhanced-sampling molecular dynamics simulations and potentials of interactions tailored to investigate the stereorecognition effects. The results remain in very good, quantitative agreement with the experimental data (measured in the context of ligand-receptor affinities and their dependence on the temperature), which provides an additional validation for the applied computational protocols. Additionally, our results contribute to the understanding of stereoselective agonist binding by β2-AR. Although the significant role of the N2936.55 residue is confirmed, we additionally show that stereorecognition does not depend solely on the N293-ligand interactions; the stereoselective effects rely on the co-operation of several residues located on both the 6th and 7th transmembrane domains and on extracellular loops. The magnitude and character of the contributions of these residues may be very diverse and result in either enhancing or reducing the stereoselective effects. The same is true when considering the enthalpic and entropic contributions to the binding free energies, which also are dependent on the ligand stereoconfiguration.

Mapping transiently formed and sparsely populated conformations on a complex energy landscape.

Determining the structures, kinetics, thermodynamics and mechanisms that underlie conformational exchange processes in proteins remains extremely difficult. Only in favourable cases is it possible to provide atomic-level descriptions of sparsely populated and transiently formed alternative conformations. Here we benchmark the ability of enhanced-sampling molecular dynamics simulations to determine the free energy landscape of the L99A cavity mutant of T4 lysozyme. We find that the simulations capture key properties previously measured by NMR relaxation dispersion methods including the structure of a minor conformation, the kinetics and thermodynamics of conformational exchange, and the effect of mutations. We discover a new tunnel that involves the transient exposure towards the solvent of an internal cavity, and show it to be relevant for ligand escape. Together, our results provide a comprehensive view of the structural landscape of a protein, and point forward to studies of conformational exchange in systems that are less characterized experimentally.

Changes in the folding landscape of the WW domain provide a molecular mechanism for an inherited genetic syndrome.

WW domains are small domains present in many human proteins with a wide array of functions and acting through the recognition of proline-rich sequences. The WW domain belonging to polyglutamine tract-binding protein 1 (PQBP1) is of particular interest due to its direct involvement in several X chromosome-linked intellectual disabilities, including Golabi-Ito-Hall (GIH) syndrome, where a single point mutation (Y65C) correlates with the development of the disease. The mutant cannot bind to its natural ligand WBP11, which regulates mRNA processing. In this work we use high-field high-resolution NMR and enhanced sampling molecular dynamics simulations to gain insight into the molecular causes the disease. We find that the wild type protein is partially unfolded exchanging among multiple beta-strand-like conformations in solution. The Y65C mutation further destabilizes the residual fold and primes the protein for the formation of a disulphide bridge, which could be at the origin of the loss of function.