Exploring the helix-coil transition via all-atom equilibrium ensemble simulations

Interactions of HP1 Bound to H3K9me3 Dinucleosome by Molecular Simulations and Biochemical Assays

Heterochromatin protein 1 (HP1), associated with heterochromatin formation, recognizes an epigenetically repressive marker, trimethylated lysine 9 in histone H3 (H3K9me3), and generally contributes to long-term silencing. How HP1 induces heterochromatin is not fully understood. Recent experiments suggested that not one, but two nucleosomes provide a platform for this recognition. Integrating previous and new biochemical assays with computational modeling, we provide near-atomic structural models for HP1 binding to the dinucleosomes. We found that the dimeric HP1α tends to bind two H3K9me3s that are in adjacent nucleosomes, thus bridging two nucleosomes. We identified, to our knowledge, a novel DNA binding motif in the hinge region that is specific to HP1α and is essential for recognizing the H3K9me3 sites of two nucleosomes. An HP1 isoform, HP1γ, does not easily bridge two nucleosomes in extended conformations because of the absence of the above binding motif and its shorter hinge region. We propose a molecular mechanism for chromatin structural changes caused by HP1.

Unfolding of α-helical 20-residue poly-glutamic acid analyzed by multiple runs of canonical molecular dynamics simulations

Elucidating the molecular mechanism of helix-coil transitions of short peptides is a long-standing conundrum in physical chemistry. Although the helix-coil transitions of poly-glutamic acid (PGA) have been extensively studied, the molecular details of its unfolding process still remain unclear. We performed all-atom canonical molecular dynamics simulations for a 20-residue PGA, over a total of 19 μs, in order to investigate its helix-unfolding processes in atomic resolution. Among the 28 simulations, starting with the α-helical conformation, all showed an unfolding process triggered by the unwinding of terminal residues, rather than by kinking and unwinding of the middle region of the chain. The helix-coil-helix conformation which is speculated by the previous experiments was not observed. Upon comparison between the N- and C-termini, the latter tended to be unstable and easily unfolded. While the probabilities of helix elongation were almost the same among the N-terminal, middle, and C-terminal regions of the chain, unwinding of the helix was enriched at the C-terminal region. The turn and 3₁₀-helix conformations were kinetic intermediates in the formation and deformation of α-helix, consistent with the previous computational studies for Ala-based peptides.

Molecular mechanics of coiled coils loaded in the shear geometry

Coiled coils are important nanomechanical building blocks in biological and biomimetic materials. A mechanistic molecular understanding of their structural response to mechanical load is essential for elucidating their role in tissues and for utilizing and tuning these building blocks in materials applications. Using a combination of single-molecule force spectroscopy (SMFS) and steered molecular dynamics (SMD) simulations, we have investigated the mechanics of synthetic heterodimeric coiled coils of different length (3-4 heptads) when loaded in shear geometry. Upon shearing, we observe an initial rise in the force, which is followed by a constant force plateau and ultimately strand separation. The force required for strand separation depends on the coiled coil length and the applied loading rate, suggesting that coiled coil shearing occurs out of equilibrium. This out-of-equilibrium behaviour is determined by a complex structural response which involves helix uncoiling, uncoiling-assisted sliding of the helices relative to each other in the direction of the applied force as well as uncoiling-assisted dissociation perpendicular to the force axis. These processes follow a hierarchy of timescales with helix uncoiling being faster than sliding and sliding being faster than dissociation. In SMFS experiments, strand separation is dominated by uncoiling-assisted dissociation and occurs at forces between 25-45 pN for the shortest 3-heptad coiled coil and between 35-50 pN for the longest 4-heptad coiled coil. These values are highly similar to the forces required for shearing apart short double-stranded DNA oligonucleotides, reinforcing the potential role of coiled coils as nanomechanical building blocks in applications where protein-based structures are desired.

Controlling anticancer drug mediated G-quadruplex formation and stabilization by a molecular container

Controlling of ligand mediated G-quadruplex DNA (GQ-DNA) formation and stabilization is an important and challenging aspect due to its active involvement in many biologically important processes such as DNA replication, transcription, etc. Here, we have demonstrated that topotecan (TPT), a potential anticancer drug, can instigate the formation and stabilization of GQ-DNA (H24 → GQ-DNA) in the absence of Na⁺/K⁺ ions via circular dichroism, fluorescence, NMR, UV melting and molecular dynamics (MD) simulation studies. The primary binding mode of TPT to GQ was found to be stacking at the terminal rather than binding to the groove. We have also reverted this conformational transition (GQ-DNA → H24) using a molecular container, cucurbit[7]uril (CB7), by means of the translocation of the drug (TPT) from GQ-DNA to its nanocavity. Importantly, we have carried out the detection of these conformational transitions using the fluorescence color switch of the drug, which is more direct and simple than some of the other methods that involve sophisticated and complex detection techniques.

Dual binding in cohesin-dockerin complexes: the energy landscape and the role of short, terminal segments of the dockerin module

The assembly of the polysaccharide degradating cellulosome machinery is mediated by tight binding between cohesin and dockerin domains. We have used an empirical model known as FoldX as well as molecular mechanics methods to determine the free energy of binding between a cohesin and a dockerin from Clostridium thermocellum in two possible modes that differ by an approximately 180° rotation. Our studies suggest that the full-length wild-type complex exhibits dual binding at room temperature, i.e., the two modes of binding have comparable probabilities at equilibrium. The ability to bind in the two modes persists at elevated temperatures. However, single-point mutations or truncations of terminal segments in the dockerin result in shifting the equilibrium towards one of the binding modes. Our molecular dynamics simulations of mechanical stretching of the full-length wild-type cohesin-dockerin complex indicate that each mode of binding leads to two kinds of stretching pathways, which may be mistakenly taken as evidence of dual binding.

On-surface synthesis: a promising strategy toward the encapsulation of air unstable ultra-thin 2D materials

2D black phosphorus (BP) and transition metal chalcogenides (TMCs) have beneficial electronic, optical, and physical properties at the few-layer limit. However, irreversible degradation of exfoliated or chemical vapor deposition-grown ultrathin BP and TMCs like GaSe via oxidation under ambient conditions limits their applications. Herein, the on-surface growth of an oxidation-resistant 2D thin film of a metal coordination polymer is demonstrated by multiscale simulations. We show that the preparation of such heterostructures can be conducted in solution, in which pristine BP and GaSe present better stability than in an air environment. Our calculations reveal that the interaction between the polymer layer and 2D materials is dominated by van der Waals forces; thus, the electronic properties of pristine BP and GaSe are well preserved. Meanwhile, the isolation from oxygen and water can be achieved by monolayer polymers, due to the nature of their close-packed layers. Our facile strategy for enhancing the environmental stability of ultrathin materials is expected to accelerate efforts to implement 2D materials in electronic and optoelectronic applications.

Charged N-terminus of Influenza Fusion Peptide Facilitates Membrane Fusion

Cleavage of hemagglutinin precursor (HA0) by cellular proteases results in the formation of two subunits, HA1 and HA2. The N-terminal fragment of HA2, named a fusion peptide (HAfp), possess a charged, amine N-terminus. It has been shown that the N-terminus of HAfp stabilizes the structure of a helical hairpin observed for a 23-amino acid long peptide (HAfp1-23), whose larger activity than HAfp1-20 has been demonstrated recently. In this paper, we analyze the effect of N-terminal charge on peptide-mediated fusion efficiency and conformation changes at the membrane interface by comparison with the corresponding N-acetylated peptides of 20- and 23-amino acid lengths. We found that higher fusogenic activities of peptides with unmodified amino termini correlates with their ability to form helical hairpin structures oriented perpendicularly to the membrane plane. Molecular dynamics simulations showed that acetylated peptides adopt open and surface-bound conformation more often, which induced less disorder of the phospholipid chains, as compared to species with unmodified amino termini.

Incorporation of transuranium elements: coordination of Cm(iii) to human serum transferrin

Structure determination of Cm(iii)-transferrin by a combined spectroscopic and theoretical approach gives insight into the biochemical behaviour of incorporated actinides.

A molecular dynamics model for glycosylphosphatidyl-inositol anchors: “flop down” or “lollipop”?

Computational model for GPI anchors tested in DMPC and POPC bilayers. The free anchor rarely occurs as an erected “lollipop-like” conformation, it rather “flops down” onto the bilayer surface. Yet an attached protein (here green fluorescent protein) exhibits extensive orientational flexibility due to the phospho-ethanolamine linker.

Exploitation of the Large-Area Basal Plane of MoS2 and Preparation of Bifunctional Catalysts through On-Surface Self-Assembly

The development of nonprecious electrochemical catalysts for water splitting is a key step to achieve a sustainable energy supply for the future. Molybdenum disulfide (MoS2) has been extensively studied as a promising low-cost catalyst for hydrogen evolution reaction (HER), whereas HER is only catalyzed at the edge for pristine MoS2, leaving a large area of basal plane useless. Herein, on-surface self-assembly is demonstrated to be an effective, facile, and damage-free method to take full advantage of the large ratio surface of MoS2 for HER by using multiscale simulations. It is found that as supplement of edge sites of MoS2, on-MoS2 M(abt)2 (M = Ni, Co; abt = 2-aminobenzenethiolate) owns high HER activity, and the self-assembled M(abt)2 monolayers on MoS2 can be obtained through a simple liquid-deposition method. More importantly, on-surface self-assembly provides potential application for overall water splitting once the self-assembled systems prove to be of both HER and oxygen evolution reaction activities, for example, on-MoS2 Co(abt)2. This work opens up a new and promising avenue (on-surface self-assembly) toward the full exploitation of the basal plane of MoS2 for HER and the preparation of bifunctional catalysts for overall water splitting.

Molecular Dynamics Simulations of the [2Fe-2S] Cluster-Binding Domain of NEET Proteins Reveal Key Molecular Determinants That Induce Their Cluster Transfer/Release

The NEET proteins are a novel family of iron-sulfur proteins characterized by an unusual three cysteine and one histidine coordinated [2Fe-2S] cluster. Aberrant cluster release, facilitated by the breakage of the Fe-N bond, is implicated in a variety of human diseases, including cancer. Here, the molecular dynamics in the multi-microsecond timescale, along with quantum chemical calculations, on two representative members of the family (the human NAF-1 and mitoNEET proteins), show that the loss of the cluster is associated with a dramatic decrease in secondary and tertiary structure. In addition, the calculations provide a mechanism for cluster release and clarify, for the first time, crucial differences existing between the two proteins, which are reflected in the experimentally observed difference in the pH-dependent cluster reactivity. The reliability of our conclusions is established by an extensive comparison with the NMR data of the solution proteins, in part measured in this work.

Three-dimensional Modelling of the Voltage-gated Sodium Ion Channel from Anopheles gambiae Reveals Spatial Clustering of Evolutionarily Conserved Acidic Residues at the Extracellular Sites

BACKGROUND

The eukaryotic voltage-gated sodium channel(e-Nav) is a large asymmetric transmembrane protein with important functions concerning neurological function. No structure has been resolved at high resolution for this protein.

METHODS

A homology model of the transmembrane and extracellular regions of an Anopheles gambiae para-like channel with emphasis on the pore entrance has been constructed, based upon the templates provided by a prokaryotic sodium channel and a potassium two-pore channel. The latter provides a template for the extracellular regions, which are located above the entrance to the pore, which is likely to open at a side of a dome formed by these loops.

RESULTS

A model created with this arrangement shows a structure similar to low-resolution cryoelectron microscope images of a related structure. The pore entrance also shows favorable electrostatic interface.

CONCLUSION

Residues responsible for the negative charge around the pore have been traced in phylogeny to highlight their importance. This model is intended for the study of pore-blocking toxins.

Uncovering GPCR and G Protein Function by Protein Structure Network Analysis

Protein structure network (PSN) analysis is one of the graph theory-based approaches currently used for investigating structural communication in biomolecular systems. Information on the system's dynamics can be provided by atomistic molecular dynamics (MD) simulations or coarse grained elastic network models paired with normal mode analysis (ENM-NMA). This chapter reports on selected applications of PSN analysis to uncover the structural communication in G protein coupled receptors (GPCRs) and G proteins. Strategies to highlight changes in structural communication caused by mutations, ligand and protein binding are described. Conserved amino acids, sites of misfolding mutations, or ligands acting as functional switches tend to behave as hubs in the native structure networks. Densely linked regions in the protein structure graphs could be identified as playing central roles in protein stability and function. Changes in the communication pathway fingerprints depending on the bound ligand or following amino acid mutation could be highlighted as well. A bridge between misfolding and misrouting could be established in rhodopsin mutants linked to inherited blindness. The analysis of native network perturbations by misfolding mutations served to infer key structural elements of protein responsiveness to small chaperones with implications for drug discovery.

HtrA1 activation is driven by an allosteric mechanism of inter-monomer communication

The human protease family HtrA is responsible for preventing protein misfolding and mislocalization, and a key player in several cellular processes. Among these, HtrA1 is implicated in several cancers, cerebrovascular disease and age-related macular degeneration. Currently, HtrA1 activation is not fully characterized and relevant for drug-targeting this protease. Our work provides a mechanistic step-by-step description of HtrA1 activation and regulation. We report that the HtrA1 trimer is regulated by an allosteric mechanism by which monomers relay the activation signal to each other, in a PDZ-domain independent fashion. Notably, we show that inhibitor binding is precluded if HtrA1 monomers cannot communicate with each other. Our study establishes how HtrA1 trimerization plays a fundamental role in proteolytic activity. Moreover, it offers a structural explanation for HtrA1-defective pathologies as well as mechanistic insights into the degradation of complex extracellular fibrils such as tubulin, amyloid beta and tau that belong to the repertoire of HtrA1.

Ensemble simulations: folding, unfolding and misfolding of a high-efficiency frameshifting RNA pseudoknot

Massive all-atom molecular dynamics simulations were conducted across a distributed computing network to study the folding, unfolding, misfolding and conformational plasticity of the high-efficiency frameshifting double mutant of the 26 nt potato leaf roll virus RNA pseudoknot. Our robust sampling, which included over 40 starting structures spanning the spectrum from the extended unfolded state to the native fold, yielded nearly 120 μs of cumulative sampling time. Conformational microstate transitions on the 1.0 ns to 10.0 μs timescales were observed, with post-equilibration sampling providing detailed representations of the conformational free energy landscape and the complex folding mechanism inherent to the pseudoknot motif. Herein, we identify and characterize two alternative native structures, three intermediate states, and numerous misfolded states, the latter of which have not previously been characterized via atomistic simulation techniques. While in line with previous thermodynamics-based models of a general RNA folding mechanism, our observations indicate that stem-strand-sequence-separation may serve as an alternative predictor of the order of stem formation during pseudoknot folding. Our results contradict a model of frameshifting based on structural rigidity and resistance to mechanical unfolding, and instead strongly support more recent studies in which conformational plasticity is identified as a determining factor in frameshifting efficiency.

Simple MD-based model for oxidative folding of peptides and proteins

Significant strides have been recently made to fold peptides and small proteins in silico using MD simulations. However, facilities are currently lacking to include disulfide bonding in the MD models of protein folding. To address this problem, we have developed a simple empirical protocol to model formation of disulfides, which is perturbation-free, retains the same speed as conventional MD simulations and allows one to control the reaction rate. The new protocol has been tested on 15-aminoacid peptide guanylin containing four cysteine residues; the net simulation time using Amber ff14SB force field was 61 μs. The resulting isomer distribution is in qualitative agreement with experiment, suggesting that oxidative folding of guanylin in vitro occurs under kinetic control. The highly stable conformation of the so-called isomer 2(B) has been obtained for full-length guanylin, which is significantly different from the poorly ordered structure of the truncated peptide PDB ID 1GNB. In addition, we have simulated oxidative folding of guanylin within the 94-aminoacid prohormone proguanylin. The obtained structure is in good agreement with the NMR coordinates 1O8R. The proposed modeling strategy can help to explore certain fundamental aspects of protein folding and is potentially relevant for manufacturing of synthetic peptides and recombinant proteins.

Modulation of the interaction between human P450 3A4 and B. megaterium reductase via engineered loops

Chimerogenesis involving cytochromes P450 is a successful approach to generate catalytically self-sufficient enzymes. However, the connection between the different functional modules should allow a certain degree of flexibility in order to obtain functional and catalytically efficient proteins. We previously applied the molecular Lego approach to develop a chimeric P450 3A4 enzyme linked to the reductase domain of P450 BM3 (BMR). Three constructs were designed with the connecting loop containing no glycine, 3 glycine or 5 glycine residues and showed a different catalytic activity and coupling efficiency. Here we investigate how the linker affects the ability of P450 3A4 to bind substrates and inhibitors. We measure the electron transfer rates and the catalytic properties of the enzyme also in the presence of ketoconazole as inhibitor. The data show that the construct 3A4-5GLY-BMR with the longest loop better retains the binding ability and cooperativity for testosterone, compared to P450 3A4. In both 3A4-3GLY-BMR and 3A4-5GLY-BMR, the substrate induces an increase in the first electron transfer rate and a shorter lag phase related to a domain rearrangements, when compared to the construct without Gly. These data are consistent with docking results and secondary structure predictions showing a propensity to form helical structures in the loop of the 3A4-BMR and 3A4-3GLY-BMR. All three chimeras retain the ability to bind the inhibitor ketoconazole and show an IC₅₀ comparable with those reported for the wild type protein. This article is part of a Special Issue entitled: Cytochrome P450 biodiversity and biotechnology, edited by Erika Plettner, Gianfranco Gilardi, Luet Wong, Vlada Urlacher, Jared Goldstone.

Acetylation within the First 17 Residues of Huntingtin Exon 1 Alters Aggregation and Lipid Binding

Huntington's disease (HD) is a genetic neurodegenerative disorder caused by an expanded polyglutamine (polyQ) domain near the N-terminus of the huntingtin (htt) protein. Expanded polyQ leads to htt aggregation. The first 17 amino acids (Nt(17)) in htt comprise a lipid-binding domain that undergoes a number of posttranslational modifications that can modulate htt toxicity and subcellular localization. As there are three lysines within Nt(17), we evaluated the impact of lysine acetylation on htt aggregation in solution and on model lipid bilayers. Acetylation of htt-exon1(51Q) and synthetic truncated htt-exon 1 mimicking peptides (Nt(17)-Q35-P10-KK) was achieved using a selective covalent label, sulfo-N-hydroxysuccinimide (NHSA). With this treatment, all three lysine residues (K6, K9, and K15) in Nt(17) were significantly acetylated. N-terminal htt acetylation retarded fibril formation in solution and promoted the formation of larger globular aggregates. Acetylated htt also bound lipid membranes and disrupted the lipid bilayer morphology less aggressively compared with the wild-type. Computational studies provided mechanistic insights into how acetylation alters the interaction of Nt(17) with lipid membranes. Our results highlight that N-terminal acetylation influences the aggregation of htt and its interaction with lipid bilayers.

Mechanistic insights on the reduction of glutathione disulfide by protein disulfide isomerase

We explore the enzymatic mechanism of the reduction of glutathione disulfide (GSSG) by the reduced a domain of human protein disulfide isomerase (hPDI) with atomistic resolution. We use classical molecular dynamics and hybrid quantum mechanics/molecular mechanics calculations at the mPW1N/6-311+G(2d,2p):FF99SB//mPW1N/6-31G(d):FF99SB level. The reaction proceeds in two stages: (i) a thiol-disulfide exchange through nucleophilic attack of the Cys53-thiolate to the GSSG-disulfide followed by the deprotonation of Cys56-thiol by Glu47-carboxylate and (ii) a second thiol-disulfide exchange between the Cys56-thiolate and the mixed disulfide intermediate formed in the first step. The Gibbs activation energy for the first stage was 18.7 kcal·mol^-1, and for the second stage, it was 7.2 kcal·mol^-1, in excellent agreement with the experimental barrier (17.6 kcal·mol^-1). Our results also suggest that the catalysis by protein disulfide isomerase (PDI) and thiol-disulfide exchange is mostly enthalpy-driven (entropy changes below 2 kcal·mol^-1 at all stages of the reaction). Hydrogen bonds formed between the backbone of His55 and Cys56 and the Cys56-thiol result in an increase in the Gibbs energy barrier of the first thiol-disulfide exchange. The solvent plays a key role in stabilizing the leaving glutathione thiolate formed. This role is not exclusively electrostatic, because an explicit inclusion of several water molecules at the density-functional theory level is a requisite to form the mixed disulfide intermediate. In the intramolecular oxidation of PDI, a transition state is only observed if hydrogen bond donors are nearby the mixed disulfide intermediate, which emphasizes that the thermochemistry of thiol-disulfide exchange in PDI is influenced by the presence of hydrogen bond donors.

The structural basis of a high affinity ATP binding ε subunit from a bacterial ATP synthase

The ε subunit from bacterial ATP synthases functions as an ATP sensor, preventing ATPase activity when the ATP concentration in bacterial cells crosses a certain threshold. The R103A/R115A double mutant of the ε subunit from thermophilic Bacillus PS3 has been shown to bind ATP two orders of magnitude stronger than the wild type protein. We use molecular dynamics simulations and free energy calculations to derive the structural basis of the high affinity ATP binding to the R103A/R115A double mutant. Our results suggest that the double mutant is stabilized by an enhanced hydrogen-bond network and fewer repulsive contacts in the ligand binding site. The inferred structural basis of the high affinity mutant may help to design novel nucleotide sensors based on the ε subunit from bacterial ATP synthases.

Molecular Dynamics Insights into Water-Parylene C Interface: Relevance of Oxygen Plasma Treatment for Biocompatibility

Solid-water interfaces play a vital role in biomaterials science because they provide a natural playground for most biochemical reactions and physiological processes. In the study, fully atomistic molecular dynamics simulations were performed to investigate interactions between water molecules and several surfaces modeling for unmodified and modified parylene C surfaces. The introduction of -OH, -CHO, and -COOH to the surface and alterations in their coverage significantly influence the energetics of interactions between water molecules and the polymer surface. The theoretical studies were complemented with experimental measurements of contact angle, surface free energy, and imaging of osteoblast cells adhesion. Both MD simulations and experiments demonstrate that the optimal interface, in terms of biocompatibility, is obtained when 60% of native -Cl groups of parylene C surface is exchanged for -OH groups. By exploring idealized models of bare and functionalized parylene C, we obtained a unique insight into molecular interactions at the water-polymer interface. The calculated values of interaction energy components (electrostatic and dispersive) correspond well with the experimentally determined values of surface free energy components (polar and dispersive), revealing their optimal ratio for cells adhesion. The results are discussed in the context of controllable tuning and functionalization of implant polymeric coating toward improved biocompatibility.

Nuclear quantum effects in a HIV/cancer inhibitor: The case of ellipticine

Ellipticine is a natural product that is currently being actively investigated for its inhibitory cancer and HIV properties. Here we use path-integral molecular dynamics coupled with excited state calculations to characterize the role of nuclear quantum effects on the structural and electronic properties of ellipticine in water, a common biological solvent. Quantum effects collectively enhance the fluctuations of both light and heavy nuclei of the covalent and hydrogen bonds in ellipticine. In particular, for the ellipticine-water system, where the proton donor and acceptor have different proton affinities, we find that nuclear quantum effects (NQEs) strengthen both the strong and the weak H bonds. This is in contrast to what is observed for the cases where the proton affinity of the donors and acceptors is same. These structural fluctuations cause a significant red-shift in the absorption spectra and an increase in the broadening, bringing it into closer agreement with the experiments. Our work shows that nuclear quantum effects alter both qualitatively and quantitatively the optical properties of this biologically relevant system and highlights the importance of the inclusion of these effects in the microscopic understanding of their optical properties. We propose that isotopic substitution will produce a blue shift and a reduction in the broadening of the absorption peak.

Microscopic nucleation and propagation rates of an alanine-based α-helix

An infrared temperature-jump (T-jump) study by Huang et al. (Proc. Natl. Acad. Sci. U. S. A., 2002, 99, 2788-2793) showed that the conformational relaxation kinetics of an alanine-based α-helical peptide depend not only on the final temperature (Tf) but also on the initial temperature (Ti) when Tf is fixed. Their finding indicates that the folding free energy landscape of this peptide is non-two-state like, allowing for the population of conformational ensembles with different helical lengths and relaxation times in the temperature range of the experiment. Because α-helix folding involves two fundamental events, nucleation and propagation, the results of Huang et al. thus present a unique opportunity to determine their rate constants - a long-sought goal in the study of the helix-coil transition dynamics. Herein, we capitalize on this notion and develop a coarse-grained kinetic model to globally fit the thermal unfolding curve and T-jump kinetic traces of this peptide. Using this strategy, we are able to explicitly determine the microscopic rate constants of the kinetic steps encountered in the nucleation and propagation processes. Our results reveal that the time taken to form one α-helical turn (i.e., an α-helical segment with one helical hydrogen bond) is about 315 ns, whereas the time taken to elongate this nucleus by one residue (or backbone unit) is 5.9 ns, depending on the position of the residue.

Passivation of Black Phosphorus via Self-Assembled Organic Monolayers by van der Waals Epitaxy

An effective passivation approach to protect black phosphorus (BP) from degradation based on multi-scale simulations is proposed. The self-assembly of perylene-3,4,9,10-tetracarboxylic dianhydride monolayers via van der Waals epitaxy on BP does not break the original electronic properties of BP. The passivation layer thickness is only 2 nm. This study opens up a new pathway toward fine passivation of BP.

Structural Determinants of Constitutive Activation of Gα Proteins: Transducin as a Paradigm

Heterotrimeric guanine nucleotide-binding proteins (Gα proteins) are intracellular nanomachines deputed to signal transduction. The switch-on process requires the release of bound GDP from a site at the interface between GTPase and helical domains. Nucleotide release is catalyzed by G protein Coupled Receptors (GPCRs). Here we investigate the functional dynamics of wild type (WT) and six constitutively active mutants (CAMs) of the Gα protein transducin (Gt) by combining atomistic molecular dynamics (MD) simulations with Maxwell-Demod discrete MD (MDdMD) simulations of the receptor-catalyzed transition between GDP-bound and nucleotide-free states. Compared to the WT, Gt CAMs increase the overall fluctuations of nucleotide and its binding site. This is accompanied by weakening of native links involving GDP, α1, the G boxes, β1-β3, and α5. Collectively, constitutive activation by the considered mutants seems to associate with weakening of the interfaces between α5 and the surrounding portions and the interface between GTPase (G) and helical (H) domains. These mutational effects associate with increases in the overall fluctuations of the G and H domains, which reflect on the collective motions of the protein. Gt CAMs, with prominence to G56P, T325A, and F332A, prioritize collective motions of the H domain overlapping with the collective motions associated with receptor-catalyzed nucleotide release. In spite of different local perturbations, the mechanisms of nucleotide exchange catalyzed by activating mutations and by receptor are expected to employ similar molecular switches in the nucleotide binding site and to share the detachment of the H domain from the G domain.

Reaction mechanism and reaction coordinates from the viewpoint of energy flow

Reaction coordinates are of central importance for correct understanding of reaction dynamics in complex systems, but their counter-intuitive nature made it a daunting challenge to identify them. Starting from an energetic view of a reaction process as stochastic energy flows biased towards preferred channels, which we deemed the reaction coordinates, we developed a rigorous scheme for decomposing energy changes of a system, both potential and kinetic, into pairwise components. The pairwise energy flows between different coordinates provide a concrete statistical mechanical language for depicting reaction mechanisms. Application of this scheme to the C7eq → C7ax transition of the alanine dipeptide in vacuum revealed novel and intriguing mechanisms that eluded previous investigations of this well studied prototype system for biomolecular conformational dynamics. Using a cost function developed from the energy decomposition components by proper averaging over the transition path ensemble, we were able to identify signatures of the reaction coordinates of this system without requiring any input from human intuition.

Ensemble Molecular Dynamics of a Protein-Ligand Complex: Residual Inhibitor Entropy Enhances Drug Potency in Butyrylcholinesterase

Butyrylcholinesterase is a key enzyme that catalyzes the hydrolysis of the neurotransmitter acetylcholine and shows an increased activity in patients suffering from Alzheimer's disease (AD), making this enzyme a primary target in treating AD. Central to this problem, and to similar scenarios involving biomolecular recognition, is our understanding of the nature of the protein-ligand complex. The butyrylcholinesterase enzyme was studied via all-atom, explicit solvent, ensemble molecular dynamics simulations sans inhibitor and in the presence of three dialkyl phenyl phosphate inhibitors of known potency to a cumulative sampling of over 40 μs. Following the relaxation of these ensembles to conformational equilibria, binding modes for each inhibitor were identified. While classical models, which assume significant reduction in protein and ligand conformational entropies, continue to be favored in contemporary studies, our observations contradict those assumptions: bound ligands occupy many conformational states, thereby stabilizing the complex, while also promoting protein flexibility.

Dual Active Site in the Endolytic Transglycosylase gp144 of Bacteriophage phiKZ

Lytic transglycosylases are abundant peptidoglycan lysing enzymes that degrade the heteropolymers of bacterial cell walls in metabolic processes or in the course of a bacteriophage infection. The conventional catalytic mechanism of transglycosylases involves only the Glu or Asp residue. Endolysin gp144 of Pseudomonas aeruginosa bacteriophage phiKZ belongs to the family of Gram-negative transglycosylases with a modular composition and C-terminal location of the catalytic domain. Glu115 of gp144 performs the predicted role of a catalytic residue. However, replacement of this residue does not completely eliminate the activity of the mutant protein. Site-directed mutagenesis has revealed the participation of Tyr197 in the catalytic mechanism, as well as the presence of a second active site involving Glu178 and Tyr147. The existence of the dual active site was supported by computer modeling and monitoring of the molecular dynamics of the changes in the conformation and surface charge distribution as a consequence of point mutations.

Accurate Three States Model for Amino Acids with Two Chemically Coupled Titrating Sites in Explicit Solvent Atomistic Constant pH Simulations and pK(a) Calculations

Correct protonation of titratable groups in biomolecules is crucial for their accurate description by molecular dynamics simulations. In the context of constant pH simulations, an additional protonation degree of freedom is introduced for each titratable site, allowing the protonation state to change dynamically with changing structure or electrostatics. Here, we extend previous approaches for an accurate description of chemically coupled titrating sites. A second reaction coordinate is used to switch between two tautomeric states of an amino acid with chemically coupled titratable sites, such as aspartate (Asp), glutamate (Glu), and histidine (His). To this aim, we test a scheme involving three protonation states. To facilitate charge neutrality as required for periodic boundary conditions and Particle Mesh Ewald (PME) electrostatics, titration of each respective amino acid is coupled to a "water" molecule that is charged in the opposite direction. Additionally, a force field modification for Amber99sb is introduced and tested for the description of carboxyl group protonation. Our three states model is tested by titration simulations of Asp, Glu, and His, yielding a good agreement, reproducing the correct geometry of the groups in their different protonation forms. We further show that the ion concentration change due to the neutralizing "water" molecules does not significantly affect the protonation free energies of the titratable groups, suggesting that the three states model provides a good description of biomolecular dynamics at constant pH.

A benchmark for reaction coordinates in the transition path ensemble

The molecular mechanism of a reaction is embedded in its transition path ensemble, the complete collection of reactive trajectories. Utilizing the information in the transition path ensemble alone, we developed a novel metric, which we termed the emergent potential energy, for distinguishing reaction coordinates from the bath modes. The emergent potential energy can be understood as the average energy cost for making a displacement of a coordinate in the transition path ensemble. Where displacing a bath mode invokes essentially no cost, it costs significantly to move the reaction coordinate. Based on some general assumptions of the behaviors of reaction and bath coordinates in the transition path ensemble, we proved theoretically with statistical mechanics that the emergent potential energy could serve as a benchmark of reaction coordinates and demonstrated its effectiveness by applying it to a prototypical system of biomolecular dynamics. Using the emergent potential energy as guidance, we developed a committor-free and intuition-independent method for identifying reaction coordinates in complex systems. We expect this method to be applicable to a wide range of reaction processes in complex biomolecular systems.

The Low Barrier Hydrogen Bond in the Photoactive Yellow Protein: A Vacuum Artifact Absent in the Crystal and Solution

There has been considerable debate on the existence of a low-barrier hydrogen bond (LBHB) in the photoactive yellow protein (PYP). The debate was initially triggered by the neutron diffraction study of Yamaguchi et al. ( Proc. Natl. Acad. Sci., U. S. A. , 2009 , 106 , 440 - 444 ) who suggested a model in which a neutral Arg52 residue triggers the formation of the LBHB in PYP. Here, we present an alternative model that is consistent within the error margins of the Yamaguchi structure factors. The model explains an increased hydrogen bond length without nuclear quantum effects and for a protonated Arg52. We tested both models by calculations under crystal, solution, and vacuum conditions. Contrary to the common assumption in the field, we found that a single PYP in vacuum does not provide an accurate description of the crystal conditions but instead introduces strong artifacts, which favor a LBHB and a large ¹H NMR chemical shift. Our model of the crystal environment was found to stabilize the two Arg52 hydrogen bonds and crystal water positions for the protonated Arg52 residue in free MD simulations and predicted an Arg52 pK_a upshift with respect to PYP in solution. The crystal and solution environments resulted in almost identical ¹H chemical shifts that agree with NMR solution data. We also calculated the effect of the Arg52 protonation state on the LBHB in 3D nuclear equilibrium density calculations. Only the charged crystal structure in vacuum supports a LBHB if Arg52 is neutral in PYP at the previously reported level of theory ( J. Am. Chem. Soc. , 2014 , 136 , 3542 - 3552 ). We attribute the anomalies in the interpretation of the neutron data to a shift of the potential minimum, which does not involve nuclear quantum effects and is transferable beyond the Yamaguchi structure.

Genetic Algorithm Managed Peptide Mutant Screening: Optimizing Peptide Ligands for Targeted Receptor Binding

This study demonstrates the utility of genetic algorithms to search exceptionally large and otherwise intractable mutant libraries for sequences with optimal binding affinities for target receptors. The Genetic Algorithm Managed Peptide Mutant Screening (GAMPMS) program was used to search an α-conotoxin (α-CTx) MII mutant library of approximately 41 billion possible peptide sequences for those exhibiting the greatest binding affinity for the α₃β₂-nicotinic acetylcholine receptor (nAChR) isoform. A series of top resulting peptide ligands with high sequence homology was obtained, with each mutant having an estimated ΔG_bind approximately double that of the potent native α-CTx MII ligand. A consensus sequence from the top GAMPMS results was subjected to more rigorous binding free energy calculations by molecular dynamics and compared to α-CTx MII and other related variants for binding with α₃β₂-nAChR. In this study, the efficiency of GAMPMS to substantially reduce the sample population size through evolutionary selection criteria to produce ligands with higher predicted binding affinity is demonstrated.

Mechanism of Phosphorylation-Induced Folding of 4E-BP2 Revealed by Molecular Dynamics Simulations

Site-specific phosphorylation of an intrinsically disordered protein, eIF4E-binding protein isoform 2 (4E-BP2), can suppress its native function by folding it into a four-stranded β-sheet, but the mechanism of this phosphorylation-induced folding is unclear. In this work, we use all-atom molecular dynamics simulations to investigate both the folded and unfolded states of 4E-BP2 under different phosphorylation states of T37 and T46. The results show that the phosphorylated forms of both T37 and T46 play important roles in stabilizing the folded structure, especially for the β-turns and the sequestered binding motif. The phosphorylated residues not only guide the folding of the protein through several intermediate states but also affect the conformational distribution of the unfolded ensemble. Significantly, the phosphorylated residues can function as nucleation sites for the folding of the protein by forming certain local structures that are stabilized by hydrogen bonding involving the phosphate group. The region around phosphorylated T46 appears to fold before that around phosphorylated T37. These findings provide new insight into the intricate effects of protein phosphorylation.

Conformational Dynamics and Protein-Substrate Interaction of ABC Transporter BtuCD at the Occluded State Revealed by Molecular Dynamics Simulations

ATP-binding cassette (ABC) transporters are ubiquitous in all three kingdoms of life and are implicated in many clinically relevant physiological processes. They couple the energy released by ATP hydrolysis to facilitate substrate translocation across cell membranes. The crystal structures of type II ABC importers have revealed their unique transmembrane domain architecture consisting of 10 transmembrane helices and their structurally conserved nucleotide-binding domains among all ABC transporters. However, molecular details of the interactions between the importers and their substrate remain largely elusive. Taking vitamin B₁₂ importer BtuCD as an exemplar of type II importers, we investigated the dynamics of its occluded state and the detailed protein-substrate interactions using molecular dynamics simulation. Our trajectories show that the importer accommodates the substrate through a nonspecific binding mode as the substrate undergoes evident vertical and tilt motions inside the translocation cavity. Extensive hydrogen bond and hydrophobic interactions were observed between the substrate and the importer; however, most of these interactions are weak, with <38% occurrence. The presence of substrate leads to enlargement of the translocation cavity, especially at its cytoplasmic end, which may activate cytoplasmic regions and probably facilitate the transportation. The perturbations caused by periplasmic binding protein and nucleotides were also investigated. The study provides deeper insight into the translocation mechanism of BtuCD.

ATP Hydrolysis Induced Conformational Changes in the Vitamin B12 Transporter BtuCD Revealed by MD Simulations

ATP binding cassette (ABC) transporters utilize the energy of ATP hydrolysis to uni-directionally transport substrates across cell membrane. ATP hydrolysis occurs at the nucleotide-binding domain (NBD) dimer interface of ABC transporters, whereas substrate translocation takes place at the translocation pathway between the transmembrane domains (TMDs), which is more than 30 angstroms away from the NBD dimer interface. This raises the question of how the hydrolysis energy released at NBDs is "transmitted" to trigger the conformational changes at TMDs. Using molecular dynamics (MD) simulations, we studied the post-hydrolysis state of the vitamin B12 importer BtuCD. Totally 3-μs MD trajectories demonstrate a predominantly asymmetric arrangement of the NBD dimer interface, with the ADP-bound site disrupted and the ATP-bound site preserved in most of the trajectories. TMDs response to ATP hydrolysis by separation of the L-loops and opening of the cytoplasmic gate II, indicating that hydrolysis of one ATP could facilitate substrate translocation by opening the cytoplasmic end of translocation pathway. It was also found that motions of the L-loops and the cytoplasmic gate II are coupled with each other through a contiguous interaction network involving a conserved Asn83 on the extended stretch preceding TM3 helix plus the cytoplasmic end of TM2/6/7 helix bundle. These findings entail a TMD-NBD communication mechanism for type II ABC importers.

Adsorption of Synthetic Cationic Polymers on Model Phospholipid Membranes: Insight from Atomic-Scale Molecular Dynamics Simulations

Although synthetic cationic polymers represent a promising class of effective antibacterial agents, the molecular mechanisms behind their antimicrobial activity remain poorly understood. To this end, we employ atomic-scale molecular dynamics simulations to explore adsorption of several linear cationic polymers of different chemical structure and protonation (polyallylamine (PAA), polyethylenimine (PEI), polyvinylamine (PVA), and poly-l-lysine (PLL)) on model bacterial membranes (4:1 mixture of zwitterionic phosphatidylethanolamine (PE) and anionic phosphatidylglycerol (PG) lipids). Overall, our findings show that binding of polycations to the anionic membrane surface effectively neutralizes its charge, leading to the reorientation of water molecules close to the lipid/water interface and to the partial release of counterions to the water phase. In certain cases, one has even an overcharging of the membrane, which was shown to be a cooperative effect of polymer charges and lipid counterions. Protonated amine groups of polycations are found to interact preferably with head groups of anionic lipids, giving rise to formation of hydrogen bonds and to a noticeable lateral immobilization of the lipids. While all the above findings are mostly defined by the overall charge of a polymer, we found that the polymer architecture also matters. In particular, PVA and PEI are able to accumulate anionic PG lipids on the membrane surface, leading to lipid segregation. In turn, PLL whose charge twice exceeds charges of PVA/PEI does not induce such lipid segregation due to its considerably less compact architecture and relatively long side chains. We also show that partitioning of a polycation into the lipid/water interface is an interplay between its protonation level (the overall charge) and hydrophobicity of the backbone. Therefore, a possible strategy in creating highly efficient antimicrobial polymeric agents could be in tuning these polycation's properties through proper combination of protonated and hydrophobic blocks.

Electric Field Keeps Chromophore Planar and Produces High Yield Fluorescence in Green Fluorescent Protein

The green fluorescent protein and its designed variants fluoresce efficiently. Because the isolated chromophore is not fluorescent in a practical sense, it is apparent that the protein environment plays a crucial role in its efficiency. Because of various obstacles in studying excited state dynamics of complex systems, however, the detailed mechanism of this efficiency enhancement is not yet clearly elucidated. Here, by adopting excited state nonadiabatic molecular dynamics simulations together with an interpolated quantum chemical potential model of the chromophore, we find that the strong electric field from the protein matrix contributes dominantly to the motional restriction of the chromophore. The delay in twisting motion subsequently obstructs the nonradiative decay that competes with fluorescence, leading naturally to an enhancement in light-emitting efficiency. Surprisingly, steric constraints make only a minor contribution to these aspects. Through residue specific analyses, we identify a group of key residues that control the excited state behavior. Testing a series of mutant GFPs with different brightnesses also supports the view regarding the importance of protein electrostatics. Our findings may provide a useful guide toward designing new fluorescent chemical systems in the future.

Role of sequence evolution and conformational dynamics in the substrate specificity and oligomerization mode of thymidylate kinases

Thymidylate kinase (TMK) is a key enzyme for the synthesis of DNA, making it an important target for the development of anticancer, antibacterial, and antiparasitic drugs. TMK homologs exhibit significant variations in sequence, residue conformation, substrate specificity, and oligomerization mode. However, the influence of sequence evolution and conformational dynamics on its quaternary structure and function has not been studied before. Based on extensive sequence and structure analyses, our study detected several non-conserved residues which are linked by co-evolution and are implicated in the observed variations in flexibility, oligomeric assembly, and substrate specificity among the homologs. These lead to differences in the pattern of interactions at the active site in TMKs of different specificity. The method was further tested on TMK from Sulfolobus tokodaii (StTMK) which has substantial differences in sequence and structure compared to other TMKs. Our analyses pointed to a more flexible dTMP-binding site in StTMK compared to the other homologs. Binding assays proved that the protein can accommodate both purine and pyrimidine nucleotides at the dTMP binding site with comparable affinity. Additionally, the residues responsible for the narrow specificity of Brugia malayi TMK, whose three-dimensional structure is unavailable, were detected. Our study provides a residue-level understanding of the differences observed among TMK homologs in previous experiments. It also illustrates the correlation among sequence evolution, conformational dynamics, oligomerization mode, and substrate recognition in TMKs and detects co-evolving residues that affect binding, which should be taken into account while designing novel inhibitors.