All-atom Replica Exchange Molecular Dynamics Research Articles

Computational investigation of the inhibitory interaction of IRF3 and SARS-CoV-2 accessory protein ORF3b

ORF3b is one of the SARS-CoV-2 accessory proteins. Previous experimental study suggested that ORF3b prevents IRF3 translocating to nucleus. However, the biophysical mechanism of ORF3b-IRF3 interaction is elusive. Here, we explored the conformation ensemble of ORF3b using all-atom replica exchange molecular dynamics simulation. Disordered ORF3b has mixed α-helix, β-turn and loop conformers. The potential ORF3b-IRF3 binding modes were searched by docking representative ORF3b conformers with IRF3, and 50 ORF3b-IRF3 complex poses were screened using molecular dynamics simulations ranging from 500 to 1000 ns. We found that ORF3b binds IRF3 predominantly on its CBP binding and phosphorylated pLxIS motifs, with CBP binding site has the highest binding affinity. The ORF3b-IRF3 binding residues are highly conserved in SARS-CoV-2. Our results provided biophysics insights into ORF3b-IRF3 interaction and explained its interferon antagonism mechanism.

In silico study on graphene quantum dots modified with various functional groups inhibiting α‑synuclein dimerization

Hypothesis: Graphene quantum dots (GQDs) with various functional groups are hypothesized to inhibit the α-synuclein (αS) dimerization, a crucial step in Parkinson’s disease pathogenesis. The potential of differently functionalized GQDs is systematically explored.Experiments: All-atom replica-exchange molecular dynamics simulations (accumulating to 75.6 μs) in explicit water were performed to study the dimerization of the αS non-amyloid component region and the influence of GQDs modified with various functional groups. Conformation ensemble, binding behavior, and free energy analysis were conducted.Findings: All studied GQDs inhibit β-sheet and backbone hydrogen bond formation in αS dimers, leading to looser oligomeric conformations. Charged GQDs severely impede the growth of extended β-sheets by providing extra contact surface. GQD binding primarily disrupts αS inter-peptide interactions through π-π stacking, CH-π interactions, and for charged GQDs, additionally through salt-bridge and hydrogen bonding interactions. GQD-COO– showed the most optimal inhibitory effect, binding mode, and intensity, which holds promise for the development of nanomedicines targeting amyloid aggregation in neurodegenerative diseases.

Binding and dimerization of PGLa peptides in anionic lipid bilayer studied by replica exchange molecular dynamics

The 21-residue PGLa peptide is well known for antimicrobial activity attributed to its ability to compromize bacterial membranes. Using all-atom explicit solvent replica exchange molecular dynamics with solute tempering, we studied PGLa binding to a model anionic DMPC/DMPG bilayer at the high peptide:lipid ratio that promotes PGLa dimerization (a two peptides per leaflet system). As a reference we used our previous simulations at the low peptide:lipid ratio (a one peptide per leaflet system). We found that the increase in the peptide:lipid ratio suppresses PGLa helical propensity, tilts the bound peptide toward the bilayer hydrophobic core, and forces it deeper into the bilayer. Surprisingly, at the high peptide:lipid ratio PGLa binding induces weaker bilayer thinning, but deeper water permeation. We explain these effects by the cross-correlations between lipid shells surrounding PGLa that leads to a much diminished efflux of DMPC lipids from the peptide proximity at the high peptide:lipid ratio. Consistent with the experimental data the propensity for PGLa dimerization was found to be weak resulting in coexistence of monomers and dimers with distinctive properties. PGLa dimers assemble via apolar criss-cross interface and become partially expelled from the bilayer residing at the bilayer-water boundary. We rationalize their properties by the dimer tendency to preserve favorable electrostatic interactions between lysine and phosphate lipid groups as well as to avoid electrostatic repulsion between lysines in the low dielectric environment of the bilayer core. PGLa homedimer interface is predicted to be distinct from that involved in PGLa-magainin heterodimers.

Influence of Protonation on the Norepinephrine Inhibiting α-Synuclein 71-82 Oligomerization.

The pathogenesis of Parkinson's disease (PD) is closely linked to the massive presence of Lewy vesicles and Lewy axons in the cytoplasm of neurons, mainly consisting of α-synuclein (αS). Norepinephrine (NE), whose secretion can be increased by exercise, has been demonstrated to prevent the fibrillation of αS and to break down the mature αS fibrils. In this work, we focus on the influence of protonation on the inhibitory ability of NE by using amyloid core fragment αS71-82 as a template. All-atom replica-exchange molecular dynamics simulations (accumulating to 33.6 μs) in explicit water were performed to explore the inhibitory effect of protonated and nonprotonated NE on αS oligomerization. Our results show that NE/NE+ can lead to a significant decrease in β-sheet content with increasing temperature, while isolated αS maintains relatively higher β-sheet conformations until 363 K, implying that both NE and NE+ can lower the critical temperature required for αS fibril decomposition. NE and NE+ also lead to the formation of less compact αS oligomers by preventing the backbone hydrogen bonds and the side-chain packing. The protonation would affect the binding affinity, interaction modes, and binding intensity of NE with αS. Interesting, NE and NE+ have a distinct binding free energy in the electrostatic and solvation terms, which mostly counter each other and produce a weak binding intensity with αS. Our work contributes to a better understanding of the inhibitory mechanism of NE and NE+ on αS oligomerization relevant to PD pathogenesis, which may provide clues for the design of antiamyloid medicine.

Binding of viral nuclear localization signal peptides to importin-α nuclear transport protein

Using all-atom replica-exchange molecular dynamics simulations, we mapped the mechanisms of binding of the nuclear localization signal (NLS) sequence from Venezuelan equine encephalitis virus (VEEV) capsid protein to importin-α (impα) transport protein. Our objective was to identify the VEEV NLS sequence fragment that confers native, experimentally resolved binding to impα as well as to study associated binding energetics and conformational ensembles. The two selected VEEV NLS peptide fragments, KKPK and KKPKKE, show strikingly different binding mechanisms. The minNLS peptide KKPK binds non-natively and nonspecifically by adopting five diverse conformational clusters with low similarity to the x-ray structure 3VE6 of NLS-impα complex. Despite the prevalence of non-native interactions, the minNLS peptide still largely binds to the impα major NLS binding site. In contrast, the coreNLS peptide KKPKKE binds specifically and natively, adopting a largely homogeneous binding ensemble with a dominant, highly native-like conformational cluster. The coreNLS peptide retains most of native binding interactions, including π-cation contacts and a tryptophan cage. While KKPK binding is governed by a complex multistate free energy landscape featuring transitions between multiple binding poses, the coreNLS peptide free energy map is simple, exhibiting a single dominant native-like bound basin. We argue that the origin of the coreNLS peptide binding specificity is several electrostatic interactions formed by the two C-terminal amino acids, Lys10 and Glu11, with impα. The coreNLS sequence is then sufficient for native binding, but none of the amino acids flanking minNLS, including Lys10 and Glu11, are strictly necessary for the native pose. Our analyses indicate that the VEEV coreNLS sequence is virtually unique among human and viral proteins interacting with impα making it a potential target for VEEV-specific inhibitors.

Dissecting the Molecular Mechanisms of the Co-Aggregation of Aβ40 and Aβ42 Peptides: A REMD Simulation Study.

The aggregation of amyloid-β protein (Aβ) into oligomers and amyloid fibrils is closely related to Alzheimer's disease (AD). Aβ40 and Aβ42, as two most prominent isoforms of Aβ peptides, can cross-interact with each other and form co-aggregates, which affect the progression of the disease. However, the molecular determinants underlying Aβ40 and Aβ42 cross-interaction and the structural details of their co-oligomers remain elusive. Herein, we performed all-atom explicit-solvent replica exchange molecular dynamics simulations on Aβ40-Aβ42 heterogeneous and Aβ40/Aβ42 homogeneous dimer systems to dissect the co-aggregation mechanisms of the two isoforms. Our results show that the interpeptide main-chain interaction of Aβ40-Aβ42 is stronger than that of Aβ40-Aβ40 and Aβ42-Aβ42. The positions of hotspot residues in heterodimers and homodimers display high similarity, implying similar molecular recognition sites for both cross-interaction and self-interaction. Contact maps of Aβ40-Aβ42 heterodimers reveal that residue pairs crucial for cross-interaction are mostly located in the C-terminal hydrophobic regions of Aβ40 and Aβ42 peptides. Conformational analysis shows that Aβ40 and Aβ42 monomers can co-assemble into β-sheet-rich heterodimers with shorter β-sheets than those in homodimers, which is decremental to monomer addition. Similar molecular recognition sites and β-sheet distribution of Aβ40 and Aβ42 peptides are observed in heterodimers and homodimers, which may provide the molecular basis for the two isoforms' co-aggregation and cross-seeding. Our work dissects the co-aggregation mechanisms of Aβ40 and Aβ42 peptides at the atomic level, which will help for in-depth understanding of the cross-talk between the two Aβ isoforms and the pathogenesis of AD.

Binding of Venezuelan Equine Encephalitis Virus Inhibitors to Importin-α Receptors Explored with All-Atom Replica Exchange Molecular Dynamics.

Although Venezuelan equine encephalitis virus (VEEV) is a life-threatening pathogen with a capacity for epidemic outbreaks, there are no FDA-approved VEEV antivirals for humans. VEEV cytotoxicity is partially attributed to the formation of a tetrameric complex between the VEEV capsid protein, the nuclear import proteins importin-α and importin-β, and the nuclear export protein CRM1, which together block trafficking through the nuclear pore complex. Experimental studies have identified small molecules from the CL6662 scaffold as potential inhibitors of the viral nuclear localization signal (NLS) sequence binding to importin-α. However, little is known about the molecular mechanism of CL6662 inhibition. To address this issue, we employed all-atom replica exchange molecular dynamics simulations to probe, in atomistic detail, the binding mechanism of CL6662 ligands to importin-α. Three ligands, including G281-1485 and two congeners with varying hydrophobicities, were considered. We investigated the distribution of ligand binding poses, their locations, and ligand specificities measured by the strength of binding interactions. We found that G281-1485 binds nonspecifically without forming well-defined binding poses throughout the NLS binding site. Binding of the less hydrophobic congener becomes strongly on-target with respect to the NLS binding site but remains nonspecific. However, a more hydrophobic congener is a strongly specific binder and the only ligand out of three to form a well-defined binding pose, while partially overlapping with the NLS binding site. On the basis of free energy estimates, we argue that all three ligands weakly compete with the viral NLS sequence for binding to importin-α in an apparent compromise to preserve host NLS binding. We further show that all-atom replica exchange binding simulations are a viable tool for studying ligands binding nonspecifically without forming well-defined binding poses.

Lysine Acetylation Changes the Mechanism of Aβ25-35 Peptide Binding and Dimerization in the DMPC Bilayer.

The impact of Lys28 acetylation on Alzheimer's Aβ peptide binding to the lipid bilayer has not been previously studied, either experimentally or computationally. To probe this common post-translational modification, we performed all-atom replica exchange molecular dynamics simulations targeting binding and aggregation of acetylated acAβ25-35 peptide within the DMPC bilayer. Using the unmodified Aβ25-35 studied previously as a reference, our results can be summarized as follows. First, Lys28 acetylation strengthens the Aβ25-35 hydrophobic moment and consequently promotes the helical structure across the peptide extending it into the N-terminus. Second, because Lys28 acetylation disrupts electrostatic contact between Lys28 and lipid phosphate groups, it reduces the binding affinity of acAβ25-35 peptides to the DMPC bilayer. Accordingly, although acetylation preserves the bimodal binding featuring a preferred inserted state and a less probable surface bound state, it decreases the stability of the former. Third, acetylation promotes acAβ25-35 aggregation and eliminates monomers as thermodynamically viable species. More importantly, acAβ25-35 retains as the most thermodynamically stable the inserted dimer with unique head-to-tail helical aggregation interface. However, due to enhanced helix structure, this dimer state becomes less stable and is less likely to propagate into higher order aggregates. Thus, acetylation is predicted to facilitate the formation of low-molecular-weight oligomers. Other post-translational modifications, including phosphorylation and oxidation, reduce helical propensity and have divergent impact on aggregation. Consequently, acetylation, when considered in its totality, has distinct consequences on Aβ25-35 binding and aggregation in the lipid bilayer.

Mechanisms of Binding of Antimicrobial Peptide PGLa to DMPC/DMPG Membrane.

PGLa belongs to a class of antimicrobial peptides showing strong affinity to anionic bacterial membranes. Using all-atom explicit solvent replica exchange molecular dynamics with solute tempering, we studied binding of PGLa to a model anionic dimyristoylphosphatidylcholine/dimyristoylphosphatidylglycerol (DMPC/DMPG) bilayer. Due to a strong hydrophobic moment, PGLa upon binding adopts a helical structure and two distinct bound states separated by a significant free energy barrier. In these states, the C-terminus helix is either surface bound or inserted into the bilayer, whereas the N-terminus remains anchored in the bilayer. Analysis of the free energy landscape indicates that the transition between the two states involves a C-terminus helix rotation permitting the peptide to preserve the interactions between cationic Lys amino acids and anionic lipid phosphorus groups. We calculated the free energy of PGLa binding and showed that it is mostly governed by the balance between desolvation of PGLa positive charges and formation of electrostatic PGLa-lipid interactions. PGLa binding induces minor bilayer thinning but causes pronounced lipid redistribution resulting from an influx of DMPG lipids into the binding footprint and efflux of DMPC lipids. Our in silico results rationalize the S-state detected in NMR experiments.

Met35 Oxidation Hinders Aβ25-35 Peptide Aggregation within the Dimyristoylphosphatidylcholine Bilayer.

Using all-atom explicit solvent replica exchange molecular dynamics simulations, we studied the aggregation of oxidized (ox) Aβ25-35 peptides into dimers mediated by the zwitterionic dimyristoylphosphatidylcholine (DMPC) lipid bilayer. By comparing oxAβ25-35 aggregation with that observed for reduced and phosphorylated Aβ25-35 peptides, we elucidated plausible impact of post-translational modifications on cytotoxicity of Aβ peptides involved in Alzheimer's disease. We found that Met35 oxidation reduces helical propensity in oxAβ25-35 peptides bound to the lipid bilayer and enhances backbone fluctuations. These factors destabilize the wild-type head-to-tail dimer interface and lower the aggregation propensity. Met35 oxidation diversifies aggregation pathways by adding monomeric species to the bound conformational ensemble. The oxAβ25-35 dimer becomes partially expelled from the DMPC bilayer and as a result inflicts limited disruption to the bilayer structure compared to wild-type Aβ25-35. Interestingly, the effect of Ser26 phosphorylation is largely opposite, as it preserves the wild-type head-to-tail aggregation interface and strengthens, not weakens, aggregation propensity. The differing effects can be attributed to the sequence locations of these post-translational modifications, since in contrast to Ser26 phosphorylation, Met35 oxidation directly affects the wild-type C-terminal aggregation interface. A comparison with experimental data is provided.

Phytochemicals from Ayurvedic plants as potential medicaments for ovarian cancer: an in silico analysis.

Ovarian cancer is one of the highly prominent gynecological malignancies after breast cancer. Although myriad literature is available, there is no specific biomarker available for the personalized treatment strategy. The unavailability of effective drug therapy for ovarian cancer calls for an urgent push in its development from the multidisciplinary scientific community. Indian Ayurvedic medicine pharmacology is widely appreciated and accepted for its immense healthcare benefits. Bioinformatics and cheminformatics approaches can be effectively used to screen phytochemicals present in the Indian Ayurvedic plants against ovarian cancer target receptors. Recent studies discern that POTE, a cancer-testis antigen (CTA) family, plays a crucial role in the proliferation and progression of cancers including ovarian cancer. Specifically, POTEE paralog has been observed to be hypermethylated in ovarian cancer. This study undertakes an in silico analysis of Indian Ayurvedic plants for their anticancer efficacy against ovarian cancer proliferation target receptor POTEE. Structures of 100 phytochemicals from 11 Ayurvedic plants were screened with ADME criteria, and qualified phytochemicals were subjected to molecular docking and interaction analysis. Only 6 phytochemicals having a high affinity to the target receptor (POTEE) were then subjected to an all-atom replica exchange molecular dynamics simulation for 50ns. Binding affinities of 6 phytochemicals cedeodarin, deodarin, hematoxylin, matairesinol, quercetin, and taxifolin with POTEE were -8.1, -7.7, -7.7, -7.9, -8.0, and - 7.7kcal/mol, respectively, and their RMSD were recorded as zero. This study concludes that phytochemicals present in Indian Ayurvedic plants namely Cedrus deodara and Asparagus racemosus possess inhibitory effects against ovarian cancer proliferation receptor POTEE.

Conformational Ensemble of TteAdoCbl Riboswitch Provides Stable Structural Elements for Conformation Selection and Population Shift in Cobalamin Recognition.

Cobalamin riboswitch is a cis-regulatory element widely found in the 5'-UTRs of the vitamin B12-associated genes in bacteria, resulting in modulation and production of a particular protein. Thermoanaerobacter tengcongensis (Tte) AdoCbl riboswitches are the largest of the known riboswitches with 210 nucleotides, partially due to its long peripheral P6-extension, which enable high affinity of AdoCbl. Two structural elements, T-loop/T-looplike motif and kissing loop are key to the global folding of the RNA. While the structure of the TteAdoCbl riboswitch complex is known, we still do not understand the structure and conformation before AdoCbl ligand recognition. In order to delineate the conformational changes and the stabilities of long-range interactions, we have performed extensive all-atom replica-exchange molecular dynamics simulations of the TteAdoCbl riboswitch with a total simulation time of 2296 ns. We found that both the T-loop/T-looplike motif and kissing loop are very stable with ligand binding. The gating conformation changes of P6-extension allow the ligand to bind to the preorganized kissing loop binding pocket. The T-loop/T-looplike motif has much more hydrogen bonds than observed in TteAdoCbl riboswitch complex crystal structure, indicating an allosteric response of the T-loop/T-looplike motif. Our study demonstrated that the conformational ensemble of TteAdoCbl riboswitch provides stable structural elements for conformation selection and population shift in cobalamin recognition.

Partitioning of Aβ Peptide Fragments into Blood-Brain Barrier Mimetic Bilayer.

We used all-atom replica-exchange umbrella sampling molecular dynamics simulations to investigate the partitioning of the charged tetrapeptide KLVF and its neutral apolar counterpart VVIA into the blood-brain barrier (BBB)-mimetic bilayer. Our findings allowed us to reconstruct the partitioning mechanism for these two Aβ peptide fragments. Despite dissimilar sequences, their permeation shares significant common features. Computations of free energies and permeabilities show that partitioning of both peptides is highly unfavorable, ruling out passive transport. The peptides experience multiple rotational transitions within the bilayer and typically cause considerable lipid disorder and bilayer thinning. Near the bilayer midplane, they lose almost entirely their solvation shell and the interactions with the lipid headgroups. The peptides cause complex reorganization within the proximal bilayer region. Upon insertion, they induce striking cholesterol influx reversed by its depletion and the influx of DMPC when the peptides reach the midplane. The differences in partitioning mechanisms are due to the much higher polarity of KLVF peptide, the permeation of which is more unfavorable and which exclusively assumes vertical orientations within the bilayer. In contrast, VVIA positions itself flat between the leaflets, causing minor disorder and even thickening of the BBB-mimetic bilayer. Due to the high density of the cholesterol-rich BBB bilayer, the unfavorable work associated with the peptide insertion provides a significant, but not dominant, contribution to the partition free energy, which is still governed by dehydration and loss of peptide-headgroup interactions. Comparison with experiments indicates that KLVF and VVIA permeation is similar to that of proline tetrapeptide, mannitol, or cimetidine, all of which exhibit no passive transport.

Effects of Ions and Small Compounds on the Structure of Aβ42 Monomers

Aggregation of amyloid-β (Aβ) proteins in the brain is a hallmark of Alzheimer's disease. This phenomenon can be promoted or inhibited by adding small molecules to the solution where Aβ is embedded. These molecules affect the ensemble of conformations sampled by Aβ monomers even before aggregation starts. Here, we perform extensive all-atom replica exchange molecular dynamics (REMD) simulations to provide a comparative study of the ensemble of conformations sampled by Aβ42 monomers in solutions that promote (i.e., aqueous solution containing NaCl) and inhibit (i.e., aqueous solutions containing scyllo-inositol or 4-aminophenol) aggregation. Simulations performed in pure water are used as our reference. We find that secondary-structure content is only affected in an antagonistic manner by promoters and inhibitors at the C-terminus and the central hydrophilic core. Moreover, the end of the C-terminus binds more favorably to the central hydrophobic core region of Aβ42 in NaCl adopting a type of strand-loop-strand structure that is disfavored by inhibitors. Nonpolar residues that form the dry core of larger aggregates of Aβ42 (e.g., PDB ID 2BEG) are found at close proximity in these strand-loop-strand structures, suggesting that their formation could play an important role in initiating nucleation. In the presence of inhibitors, the C-terminus binds the central hydrophilic core with a higher probability than in our reference simulation. This sensitivity of the C-terminus, which is affected in an antagonistic manner by inhibitors and promoters, provides evidence for its critical role in accounting for aggregation.

Methionine Oxidation Changes the Mechanism of A\u03b2 Peptide Binding to the DMPC Bilayer

Using all-atom explicit solvent replica exchange molecular dynamics simulations with solute tempering, we study the effect of methionine oxidation on Aβ10–40 peptide binding to the zwitterionic DMPC bilayer. By comparing oxidized and reduced peptides, we identified changes in the binding mechanism caused by this modification. First, Met35 oxidation unravels C-terminal helix in the bound peptides. Second, oxidation destabilizes intrapeptide interactions and expands bound peptides. We explain these outcomes by the loss of amphiphilic character of the C-terminal helix due to oxidation. Third, oxidation “polarizes” Aβ binding to the DMPC bilayer by strengthening the interactions of the C-terminus with lipids while largely releasing the rest of the peptide from bilayer. Fourth, in contrast to the wild-type peptide, oxidized Aβ induces significantly smaller bilayer thinning and drop in lipid density within the binding footprint. These observations are the consequence of mixing oxidized peptide amino acids with lipids promoted by enhanced Aβ conformational fluctuations. Fifth, methionine oxidation reduces the affinity of Aβ binding to the DMPC bilayer by disrupting favorable intrapeptide interactions upon binding, which offset the gains from better hydration. Reduced binding affinity of the oxidized Aβ may represent the molecular basis for its reduced cytotoxicity.

Investigating the Influence of Arginine Dimethylation on Nucleosome Dynamics Using All-Atom Simulations and Kinetic Analysis.

The dimethylation of Arg at the 42nd position (R42me2) of H3 histone, a post-translational modification (PTM) in nucleosomes close to the DNA entry/exit region, showed controversial gene regulations. To address this discrepancy, we performed comprehensive all-atom replica-exchange molecular dynamics simulations with and without a single PTM, either symmetric (R42me2s) or asymmetric (R42me2a) dimethylation. Together with a kinetics analysis, our simulations showed that DNA at the entry/exit region in the R42me2a nucleosome adopts a relatively more open conformation than that in the unmodified nucleosome, whereas R42me2s exhibits significantly weaker or even negligible effects on DNA dynamics and structures, which may provide clues of the discrepancy of gene regulation by R42me2. Our approach will be useful to study the mechanism of nucleosome dynamical change induced by a subtle modification.

How Reactive are Druggable Cysteines in Protein Kinases?

Targeted covalent inhibitors (TCIs) have been successfully developed as high-affinity and selective inhibitors of enzymes of the protein kinase family. These drugs typically act by undergoing an electrophilic addition with an active-site cysteine residue, so design of a TCI begins with the identification of a "druggable" cysteine. These electrophilic additions generally require deprotonation of the thiol to form a reactive anionic thiolate, so the acidity of the residue is a critical factor. Few experimental measurements of the p Ka's of druggable cysteines have been reported, so computational prediction could prove to be very important in selecting reactive cysteine targets. Here we report the computed p Ka's of druggable cysteines in selected protein kinases that are of clinical relevance for targeted therapies. The p Ka's of the cysteines were calculated using advanced computational methods based on all-atom replica-exchange thermodynamic integration molecular dynamics simulations in explicit solvent. We found that the acidities of druggable cysteines within protein kinases are diverse and elevated, indicating enormous differences in their reactivity. Constant-pH molecular dynamics simulations were also performed on selected protein kinases, and the results confirmed this varied range in the acidities of druggable cysteines. Many of these active-site cysteines have low exposure to solvent molecules, elevating their p Ka values. Electrostatic interactions with nearby anionic residues also elevate the p Ka's of cysteine residues in the active site. The results suggest that some cysteine residues within kinase binding sites will be slow to react with a TCI because of their low acidity. Several oncogenic kinase mutations were also modeled and found to have p Ka's similar to that of the wild-type kinase.

A spin-1 representation for dual-funnel energy landscapes.

The interconversion between the left- and right-handed helical folds of a polypeptide defines a dual-funneled free energy landscape. In this context, the funnel minima are connected through a continuum of unfolded conformations, evocative of the classical helix-coil transition. Physical intuition and recent conjectures suggest that this landscape can be mapped by assigning a left- or right-handed helical state to each residue. We explore this possibility using all-atom replica exchange molecular dynamics and an Ising-like model, demonstrating that the energy landscape architecture is at odds with a two-state picture. A three-state model-left, right, and unstructured-can account for most key intermediates during chiral interconversion. Competing folds and excited conformational states still impose limitations on the scope of this approach. However, the improvement is stark: Moving from a two-state to a three-state model decreases the fit error from 1.6 kBT to 0.3 kBT along the left-to-right interconversion pathway.

Coupled Mutations-Enabled Glycerol Transportation in an Aquaporin Z Mutant.

Aquaporins are transmembrane channel proteins with key function being transportation of water or other small substrates. Escherichia coli Aqp Z transports water molecules only, whereas Glp F is permeable to glycerol. It is intriguing to explore the possibility to induce glycerol permeability in Aqp Z by targeted mutations. The Aqp Z mutants with mutated selectivity filter (SF) residues exhibit poor permeability for both glycerol and water. For addressing the complexity of protein systems, pair correlation information in protein sequence analyses is instructive to identify residues that are coupled by coevolution and motion. In this study, we analyze the correlation between residues and unravel the clustering patterns of coupled residues, beyond SF residues, in aquaglyceroporins (AQGPs). The identified coupled motifs are proposed to be sequenced into aquaporin (Aqp Z) to introduce glycerol permeability. These residues are located in the vicinity of SF region, C-loop, and M6–M7 linkage domain. Significant enlargement of SF pore size of the proposed Aqp Z mutant is observed by an all-atom replica exchange molecular dynamics simulation, which is critical to facilitate considerable glycerol passage as characterized in calculated free-energy landscapes. Clearly, the hidden connections among residues play crucial roles in water/glycerol selectivity. In contrast, single-site mutation-based scheme may even lead to undesirable effects in AQGPs, such as the blocking of water transportation by aromatic π-stacked gate. As demonstrated in this work, the pair correlation analysis guided rational mutagenesis provides a feasible strategy to modulate proteins’ functions.

Cholesterol Changes the Mechanisms of Aβ Peptide Binding to the DMPC Bilayer.

Using isobaric-isothermal all-atom replica-exchange molecular dynamics (REMD) simulations, we investigated the equilibrium binding of Aβ10-40 monomers to the zwitterionic dimyristoylphosphatidylcholine (DMPC) bilayer containing cholesterol. Our previous REMD simulations, which studied binding of the same peptide to the cholesterol-free DMPC bilayer, served as a control, against which we measured the impact of cholesterol. Our findings are as follows. First, addition of cholesterol to the DMPC bilayer partially expels the Aβ peptide from the hydrophobic core and promotes its binding to bilayer polar headgroups. Using thermodynamic and energetics analyses, we argued that Aβ partial expulsion is not related to cholesterol-induced changes in lateral pressure within the bilayer but is caused by binding energetics, which favors Aβ binding to the surface of the densely packed cholesterol-rich bilayer. Second, cholesterol has a protective effect on the DMPC bilayer structure against perturbations caused by Aβ binding. More specifically, cholesterol reduces bilayer thinning and overall depletion of bilayer density beneath the Aβ binding footprint. Third, we found that the Aβ peptide contains a single cholesterol binding site, which involves hydrophobic C-terminal amino acids (Ile31-Val36), Phe19, and Phe20 from the central hydrophobic cluster, and cationic Lys28 from the turn region. This binding site accounts for about 76% of all Aβ-cholesterol interactions. Because cholesterol binding site in the Aβ10-40 peptide does not contain the GXXXG motif featured in cholesterol interactions with the transmembrane domain C99 of the β-amyloid precursor protein, we argued that the binding mechanisms for Aβ and C99 are distinct reflecting their different conformations and positions in the lipid bilayer. Fourth, cholesterol sharply reduces the helical propensity in the bound Aβ peptide. As a result, cholesterol largely eliminates the emergence of helical structure observed upon Aβ transition from a water environment to the cholesterol-free DMPC bilayer. We explain this effect by the formation of hydrogen bonds between cholesterol and the Aβ backbone, which prevent helix formation. Taken together, we expect that our simulations will advance understanding of a molecular-level mechanism behind the role of cholesterol in Alzheimer's disease.