Multiscale modeling of nucleosome dynamics

Conformational Dynamics of the Nucleosomal Histone H2B Tails Revealed by Molecular Dynamics Simulations

Epigenetic modifications of histone N-terminal tails play a critical role in regulating the chromatin structure and biological processes such as transcription and DNA repair. One of the key post-translational modifications (PTMs) is the acetylation of lysine residues on histone tails. Epigenetic modifications are ubiquitous in the development of diseases, such as cancer and neurological disorders. Histone H2B tails are critical regulators of nucleosome dynamics, biological processes, and certain diseases. Here, we report all-atomistic molecular dynamics (MD) simulations of the nucleosome to demonstrate that acetylation of the histone tails changes their conformational space and interaction with DNA. We perform simulations of H2B tails, critical regulators of gene regulation, in both the lysine-acetylated (ACK) and unacetylated wild type (WT) states. To explore the effects of salt concentration, we use two different NaCl concentrations to perform simulations at microsecond time scales. Salt can modulate the effects of electrostatic interactions between the DNA phosphate backbone and histone tails. Upon acetylation, H2B tails shift their secondary structure helical propensity. The number of contacts between the DNA and the H2B tail decreases. We characterize the conformational dynamics of the H2B tails by principal component analysis (PCA). The ACK tails become more compact at increased salt concentrations, but conformations from the WT tails display the most contacts with DNA at both salt concentrations. Mainly, H2B acetylation may increase the DNA accessibility for regulatory proteins to bind, which can aid in gene regulation and NCP stability.

Structural studies of catalytic peptides using molecular dynamics simulations

Many self-assembling peptides can form amyloid like structures with different sizes and morphologies. Driven by non-covalent interactions, their aggregation can occur through distinct pathways. Additionally, they can bind metal ions to create enzyme like active sites that allow them to catalyze diverse reactions. Due to the non-crystalline nature of amyloids, it is quite challenging to elucidate their structures using experimental spectroscopic techniques. In this aspect, molecular dynamics (MD) simulations provide a useful tool to derive structures of these macromolecules in solution. They can be further validated by comparing with experimentally measured structural parameters. However, these simulations require a multi-step process starting from the selection of the initial structure to the analysis of MD trajectories. There are multiple force fields, parametrization protocols, equilibration processes, software and analysis tools available for this process. Therefore, it is complicated for non-experts to select the most relevant tools and perform these simulations effectively. In this chapter, a systematic methodology that covers all major aspects of modeling of catalytic peptides is provided in a user-friendly manner. It will be helpful for researchers in this critical area of research.

Computational insights into the cross-talk between medin and Aβ: implications for age-related vascular risk factors in Alzheimer's disease

The aggregation of medin forming aortic medial amyloid is linked to arterial wall degeneration and cerebrovascular dysfunction. Elevated levels of arteriolar medin are correlated with an increased presence of vascular amyloid-β (Aβ) aggregates, a hallmark of Alzheimer's disease (AD) and vascular dementia. The cross-interaction between medin and Aβ results in the formation of heterologous fibrils through co-aggregation and cross-seeding processes both in vitro and in vivo. However, a comprehensive molecular understanding of the cross-interaction between medin and Aβ-two intrinsically disordered proteins-is critically lacking. Here, we employed atomistic discrete molecular dynamics simulations to systematically investigate the self-association, co-aggregation and also the phenomenon of cross-seeding between these two proteins. Our results demonstrated that both Aβ and medin were aggregation prone and their mixture tended to form β-sheet-rich hetero-aggregates. The formation of Aβ-medin hetero-aggregates did not hinder Aβ and medin from recruiting additional Aβ and medin peptides to grow into larger β-sheet-rich aggregates. The β-barrel oligomer intermediates observed in the self-aggregations of Aβ and medin were also present during their co-aggregation. In cross-seeding simulations, preformed Aβ fibrils could recruit isolated medin monomers to form elongated β-sheets. Overall, our comprehensive simulations suggested that the cross-interaction between Aβ and medin may contribute to their pathological aggregation, given the inherent amyloidogenic tendencies of both medin and Aβ. Targeting medin, therefore, could offer a novel therapeutic approach to preserving brain function during aging and AD by improving vascular health.

In Silico Analysis of the Effect of Hydrastis canadensis on Controlling Breast Cancer

Background and Objectives: Breast cancer is a significant type of cancer among women worldwide. Studies have reported the anti-carcinogenic activity of Hydrastis Canadensis (Goldenseal) in cancer cell lines. Hydrastis Canadensis could help eliminate toxic substances due to its anti-cancer, anti-inflammatory, and other properties. The design phase includes the identification of potential and effective molecules through modern computational techniques. Objective: This work aims to study Hydrastis Canadensis's effect in controlling hormone-independent breast cancer through in-silico analysis. Materials and Methods: The preliminary screening of reported phytochemicals includes biomolecular networking. Identifying functionally relevant phytochemicals and the respective target mutations/genes leads to selecting 3D proteins of the desired mutations being considered the target. Interaction studies have been conducted using docking. The kinetic and thermodynamic stability of complexes was studied through molecular dynamic simulation and MM-PBSA/GBSA analysis. Pharmacodynamic and pharmacokinetic features have been predicted. The mechanism-wise screening, functional enrichment, and interactional studies suggest that canadaline and Riboflavin effectively interact with the target proteins. Results: Hydrastis Canadensis has been identified as the effective formulation containing all these constituents. The phytoconstituents; Riboflavin and Canadensis showed good interaction with the targets of hormone-independent breast cancer. The complexes were found to be kinetically and thermodynamically stable. Conclusions: Hydrastis Canadensis has been identified as effective in controlling 'hormone-independent or basal-like breast cancer' followed by 'hormone-dependent breast cancer: Luminal A' and Luminal B.

Combining molecular dynamics simulations and scoring method to computationally model ubiquitylated linker histones in chromatosomes

The chromatin in eukaryotic cells plays a fundamental role in all processes during a cell's life cycle. This nucleoprotein is normally tightly packed but needs to be unpacked for expression and division. The linker histones are critical for such packaging processes and while most experimental and simulation works recognize their crucial importance, the focus is nearly always set on the nucleosome as the basic chromatin building block. Linker histones can undergo several modifications, but only few studies on their ubiquitylation have been conducted. Mono-ubiquitylated linker histones (HUb), while poorly understood, are expected to influence DNA compaction. The size of ubiquitin and the globular domain of the linker histone are comparable and one would expect an increased disorder upon ubiquitylation of the linker histone. However, the formation of higher order chromatin is not hindered and ubiquitylation of the linker histone may even promote gene expression. Structural data on chromatosomes is rare and HUb has never been modeled in a chromatosome so far. Descriptions of the chromatin complex with HUb would greatly benefit from computational structural data. In this study we generate molecular dynamics simulation data for six differently linked HUb variants with the help of a sampling scheme tailored to drive the exploration of phase space. We identify conformational sub-states of the six HUb variants using the sketch-map algorithm for dimensionality reduction and iterative HDBSCAN for clustering on the excessively sampled, shallow free energy landscapes. We present a highly efficient geometric scoring method to identify sub-states of HUb that fit into the nucleosome. We predict HUb conformations inside a nucleosome using on-dyad and off-dyad chromatosome structures as reference and show that unbiased simulations of HUb produce significantly more fitting than non-fitting HUb conformations. A tetranucleosome array is used to show that ubiquitylation can even occur in chromatin without too much steric clashes.

An understanding of coronavirus and exploring the molecular dynamics simulations to find promising candidates against the Mpro of nCoV to combat the COVID-19: A systematic review

The first infection case of new coronavirus was reported at the end of 2019 and after then, the cases are reported in all nations across the world in a very short period. Further, the regular news of mutations in the virus has made life restricted with appropriate behavior. To date, a new strain (Omicron and its new subvariant Omicron XE) has brought fear amongst us due to a higher trajectory of increase in the number of cases. The researchers thus started giving attention to this viral infection and discovering drug-like candidates to cure the infections. Finding a drug for any viral infection is not an easy task and takes plenty of time. Therefore, computational chemistry/bioinformatics is followed to get promising molecules against viral infection. Molecular dynamics (MD) simulations are being explored to get drug candidates in a short period. The molecules are screened via molecular docking, which provides preliminary information which can be further verified by molecular dynamics (MD) simulations. To understand the change in structure, MD simulations generated several trajectories such as root mean square deviation (RMSD), root mean square fluctuation (RMSF), hydrogen bonding, and radius of gyration for the main protease (Mpro) of the new coronavirus (nCoV) in the presence of small molecules. Additionally, change in free energy for the formation of complex of Mpro of nCoV with the small molecule can be determined by applying molecular mechanics with generalized born and surface area solvation (MM-GBSA). Thus, the promising molecules can be further explored for clinical trials to combat coronavirus disease-19 (COVID-19).

Nucleosome destabilization by polyamines

The roles and molecular interactions of polyamines (PAs) in the nucleus are not fully understood. Here their effect on nucleosome stability, a key regulatory factor in eukaryotic gene control, is reported, as measured in agarose embedded nuclei of H2B-GFP expressor HeLa cells. Nucleosome stability was assessed by quantitative microscopy [1,2] in situ, in close to native state of chromatin, preserving the nucleosome constrained topology of the genomic DNA. A robust destabilizing effect was observed in the millimolar concentration range in the case of spermine, spermidine as well as putrescine, which was strongly pH and salt concentration-dependent, and remained significant also at neutral pH. The integrity of genomic DNA was not affected by PA treatment, excluding DNA break-elicited topological relaxation as a factor in destabilization. The binding of PAs to DNA was demonstrated by the displacement of ethidium bromide, both from deproteinized nuclear halos and from plasmid DNA. The possibility that DNA methylation patterns may be influenced by PA levels is contemplated in the context of gene expression and DNA methylation correlations identified in the NCI-60 panel-based CellMiner database: methylated loci in subsets of high-ODC1 cell lines and the dependence of PER3 DNA methylation on PA metabolism.

A potential peptide inhibitor of SARS-CoV-2 S and human ACE2 complex

The disease COVID-19 has caused heavy socio-economic burden and there is immediate need to control it. The disease is caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) virus. The viral entry into human cell depends on the attachment of spike (S) protein via its receptor binding domain (RBD) to human cell receptor angiotensin-converting enzyme 2 (hACE2). Thus, blocking the virus attachment to hACE2 could serve as potential therapeutics for viral infection. We have designed a peptide inhibitor (ΔABP-α2) targeting the RBD of S protein using in-silico approach. Docking studies and computed affinities suggested that peptide inhibitor binds at the RBD with ∼95-fold higher affinity than hACE2. Molecular dynamics (MD) simulation confirms the stable binding of inhibitor to hACE2. Immunoinformatics studies suggest non-immunogenic and non-toxic nature of peptide. Thus, the proposed peptide could serve as potential blocker for viral attachment.

Communicated by Ramaswamy H. Sarma

Identification of zebrafish fumarate hydratase active site by molecular docking and simulation studies

Fumarate hydratase (FH), one of the members of TCA cycle, acts as a catalyte for the synthesis of malate from fumarate. FH has been proposed to play as a tumour suppressor leading to the pathogenicity of leiomyomas, renal cell carcinoma and paraganglioma. Mutations in the active site of FH lead to alteration in the protein structure. Similarly, binding of several chemical inhibitors to the active site also leads to the disruption of protein structural integrity thereby leading to protein dysfunction. Therefore, in order to address this mechanism leading to cancer, the binding efficiency of potential human FH inhibitor citrate to zebrafish fh has been extensively analysed in this study by molecular docking and simulation experiments followed by quantification of fumarate hydratase enzyme activity to validate and confirm the findings. Molecular docking revealed stronger interaction of zebrafish fh protein with inhibitor citrate when compared to natural substrate fumarate. Study on the dynamics of docked structures further confirmed that citrate was found to possess more binding affinity than fumarate. In vitro biochemical analysis also revealed concentration dependent potential inhibitory effect of citrate on zebrafish fh, thus confirming the findings of the in-silico experiments.Communicated by Ramaswamy H. Sarma.

Data-driven computational analysis of allosteric proteins by exploring protein dynamics, residue coevolution and residue interaction networks

BACKGROUND

Computational studies of allosteric interactions have witnessed a recent renaissance fueled by the growing interest in modeling of the complex molecular assemblies and biological networks. Allosteric interactions in protein structures allow for molecular communication in signal transduction networks.

METHODS

In this work, we performed a large scale comprehensive and multi-faceted analysis of >300 diverse allosteric proteins and complexes with allosteric modulators. By modeling and exploring coarse-grained dynamics, residue coevolution, and residue interaction networks for allosteric proteins, we have determined unifying molecular signatures shared by allosteric systems.

RESULTS

The results of this study have suggested that allosteric inhibitors and allosteric activators may differentially affect global dynamics and network organization of protein systems, leading to diverse allosteric mechanisms. By using structural and functional data on protein kinases, we present a detailed case study that that included atomic-level analysis of coevolutionary networks in kinases bound with allosteric inhibitors and activators.

CONCLUSIONS

We have found that coevolutionary networks can form direct communication pathways connecting functional regions and can recapitulate key regulatory sites and interactions responsible for allosteric signaling in the studied protein systems. The results of this computational investigation are compared with the experimental studies and reveal molecular signatures of known regulatory hotspots in protein kinases.

GENERAL SIGNIFICANCE

This study has shown that allosteric inhibitors and allosteric activators can have a different effect on residue interaction networks and can exploit distinct regulatory mechanisms, which could open up opportunities for probing allostery and new drug combinations with broad range of activities.

Protobiotic Systems Chemistry Analyzed by Molecular Dynamics

Systems chemistry has been a key component of origin of life research, invoking models of life's inception based on evolving molecular networks. One such model is the graded autocatalysis replication domain (GARD) formalism embodied in a lipid world scenario, which offers rigorous computer simulation based on defined chemical kinetics equations. GARD suggests that the first pre-RNA life-like entities could have been homeostatically-growing assemblies of amphiphiles, undergoing compositional replication and mutations, as well as rudimentary selection and evolution. Recent progress in molecular dynamics has provided an experimental tool to study complex biological phenomena such as protein folding, ligand-receptor interactions, and micellar formation, growth, and fission. The detailed molecular definition of GARD and its inter-molecular catalytic interactions make it highly compatible with molecular dynamics analyses. We present a roadmap for simulating GARD's kinetic and thermodynamic behavior using various molecular dynamics methodologies. We review different approaches for testing the validity of the GARD model by following micellar accretion and fission events and examining compositional changes over time. Near-future computational advances could provide empirical delineation for further system complexification, from simple compositional non-covalent assemblies towards more life-like protocellular entities with covalent chemistry that underlies metabolism and genetic encoding.

Application of molecular dynamics simulation to improve the theoretical prediction for collisional cross section of aromatic compounds with long alkyl chains in crude oils

RATIONALE

Molecular dynamics (MD) simulations with finite temperature were performed to improve the theoretical prediction of collisional cross section (CCS) values, especially for aromatic compounds containing long alkyl chains.

METHODS

In this study, the CCS values of 11 aromatic compounds with long alkyl chains were calculated by MD simulations while considering internal energy at 300, 500, and 700 K, and the results were compared with experimentally determined values.

RESULTS

The CCS values calculated at higher energies showed better agreement with the experimental values. Polycyclic aromatic hydrocarbons (PAHs) such as pentacene and benz[b]anthracene were also investigated, and better agreement between the theoretical and experimental results was observed when higher temperature (or higher internal energy) was considered.

CONCLUSIONS

The data presented in this study show that the internal degrees of freedom of ions must be considered to accurately predict the CCS values of aromatic compounds with a flexible structure measured by ion mobility mass spectrometry.

Molecular Mechanism of Spontaneous Nucleosome Unraveling

Meters of DNA wrap around histone proteins to form nucleosomes and fit inside the micron-diameter nucleus. For the genetic information encoded in the DNA to become available for transcription, replication, and repair, the DNA-histone assembly must be disrupted. Experiment has indicated that the outer stretches of nucleosomal DNA "breathe" by spontaneously detaching from and reattaching to the histone core. Here, we report direct observation of spontaneous DNA breathing in atomistic molecular dynamics simulations, detailing a microscopic mechanism of the DNA breathing process. According to our simulations, the outer stretches of nucleosomal DNA detach in discrete steps involving 5 or 10 base pairs, with the detachment process being orchestrated by the motion of several conserved histone residues. The inner stretches of nucleosomal DNA are found to be more stably associated with the histone core by more abundant nonspecific DNA-protein contacts, providing a microscopic interpretation of nucleosome unraveling experiments. The CG content of nucleosomal DNA is found to anticorrelate with the extent of unwrapping, supporting the possibility that AT-rich segments may signal the start of transcription by forming less stable nucleosomes.

Biophysical simulations and structure-based modeling of residue interaction networks in the tumor suppressor proteins reveal functional role of cancer mutation hotspots in molecular communication

In the current study, we have combined molecular simulations and energetic analysis with dynamics-based network modeling and perturbation response scanning to determine molecular signatures of mutational hotspot residues in the p53, PTEN, and SMAD4 tumor suppressor proteins. By examining structure, energetics and dynamics of these proteins, we have shown that inactivating mutations preferentially target a group of structurally stable residues that play a fundamental role in global propagation of dynamic fluctuations and mediating allosteric interaction networks. Through integration of long-range perturbation dynamics and network-based approaches, we have quantified allosteric potential of residues in the studied proteins. The results have revealed that mutational hotspot sites often correspond to high centrality mediating centers of the residue interaction networks that are responsible for coordination of global dynamic changes and allosteric signaling. Our findings have also suggested that structurally stable mutational hotpots can act as major effectors of allosteric interactions and mutations in these positions are typically associated with severe phenotype. Modeling of shortest inter-residue pathways has shown that mutational hotspot sites can also serve as key mediating bridges of allosteric communication in the p53 and PTEN protein structures. Multiple regression models have indicated that functional significance of mutational hotspots can be strongly associated with the network signatures serving as robust predictors of critical regulatory positions responsible for loss-of-function phenotype. The results of this computational investigation are compared with the experimental studies and reveal molecular signatures of mutational hotspots, providing a plausible rationale for explaining and localizing disease-causing mutations in tumor suppressor genes.

Singular value decomposition for the correlation of atomic fluctuations with arbitrary angle

Many proteins exhibit a critical property called allostery, which enables intra-molecular transmission of information between distal sites. Microscopically, allosteric response is closely related to correlated atomic fluctuations. Conventional correlation analysis correlates the atomic fluctuations at two sites by taking the dot product (DP) between the fluctuations, which accounts only for the parallel and antiparallel components. Here, we present a singular value decomposition (SVD) method that analyzes the correlation coefficient of fluctuation dynamics with an arbitrary angle between the correlated directions. In a model allosteric system, the second PDZ domain (PDZ2) in the human PTP1E protein, approximately one third of the strong correlations have near-perpendicular directions, which are underestimated in the conventional method. The discrimination becomes more prominent for residue pairs with larger separation. The results of the proposed SVD method are more consistent with the experimentally determined PDZ2 dynamics than those of conventional method. In addition, the SVD method improved the prediction accuracy of the allosteric sites in a dataset of 23 known allosteric monomer proteins. The proposed method may inspire extended investigation not only into allostery, but also into protein dynamics and drug design.

Force-extension behavior of DNA in the presence of DNA-bending nucleoid associated proteins

Interactions between nucleoid associated proteins (NAPs) and DNA affect DNA polymer conformation, leading to phenomena such as concentration dependent force-extension behavior. These effects, in turn, also impact the local binding behavior of the protein, such as high forces causing proteins to unbind, or proteins binding favorably to locally bent DNA. We develop a coarse-grained NAP-DNA simulation model that incorporates both force- and concentration-dependent behaviors, in order to study the interplay between NAP binding and DNA conformation. This model system includes multi-state protein binding and unbinding, motivated by prior work, but is now dependent on the local structure of the DNA, which is related to external forces acting on the DNA strand. We observe the expected qualitative binding behavior, where more proteins are bound at lower forces than at higher forces. Our model also includes NAP-induced DNA bending, which affects DNA elasticity. We see semi-quantitative matching of our simulated force-extension behavior to the reported experimental data. By using a coarse-grained simulation, we are also able to look at non-equilibrium behaviors, such as dynamic extension of a DNA strand. We stretch a DNA strand at different rates and at different NAP concentrations to observe how the time scales of the system (such as pulling time and unbinding time) work in concert. When these time scales are similar, we observe measurable rate-dependent changes in the system, which include the number of proteins bound and the force required to extend the DNA molecule. This suggests that the relative time scales of different dynamic processes play an important role in the behavior of NAP-DNA systems.

Modulation of nucleosomal DNA accessibility via charge-altering post-translational modifications in histone core

BACKGROUND

Controlled modulation of nucleosomal DNA accessibility via post-translational modifications (PTM) is a critical component to many cellular functions. Charge-altering PTMs in the globular histone core-including acetylation, phosphorylation, crotonylation, propionylation, butyrylation, formylation, and citrullination-can alter the strong electrostatic interactions between the oppositely charged nucleosomal DNA and the histone proteins and thus modulate accessibility of the nucleosomal DNA, affecting processes that depend on access to the genetic information, such as transcription. However, direct experimental investigation of the effects of these PTMs is very difficult. Theoretical models can rationalize existing observations, suggest working hypotheses for future experiments, and provide a unifying framework for connecting PTMs with the observed effects.

RESULTS

A physics-based framework is proposed that predicts the effect of charge-altering PTMs in the histone core, quantitatively for several types of lysine charge-neutralizing PTMs including acetylation, and qualitatively for all phosphorylations, on the nucleosome stability and subsequent changes in DNA accessibility, making a connection to resulting biological phenotypes. The framework takes into account multiple partially assembled states of the nucleosome at the atomic resolution. The framework is validated against experimentally known nucleosome stability changes due to the acetylation of specific lysines, and their effect on transcription. The predicted effect of charge-altering PTMs on DNA accessibility can vary dramatically, from virtually none to a strong, region-dependent increase in accessibility of the nucleosomal DNA; in some cases, e.g., H4K44, H2AK75, and H2BK57, the effect is significantly stronger than that of the extensively studied acetylation sites such H3K56, H3K115 or H3K122. Proximity to the DNA is suggestive of the strength of the PTM effect, but there are many exceptions. For the vast majority of charge-altering PTMs, the predicted increase in the DNA accessibility should be large enough to result in a measurable modulation of transcription. However, a few possible PTMs, such as acetylation of H4K77, counterintuitively decrease the DNA accessibility, suggestive of the repressed chromatin. A structural explanation for the phenomenon is provided. For the majority of charge-altering PTMs, the effect on DNA accessibility is simply additive (noncooperative), but there are exceptions, e.g., simultaneous acetylation of H4K79 and H3K122, where the combined effect is amplified. The amplification is a direct consequence of the nucleosome-DNA complex having more than two structural states. The effect of individual PTMs is classified based on changes in the accessibility of various regions throughout the nucleosomal DNA. The PTM's resulting imprint on the DNA accessibility, "PTMprint," is used to predict effects of many yet unexplored PTMs. For example, acetylation of H4K44 yields a PTMprint similar to the PTMprint of H3K56, and thus acetylation of H4K44 is predicted to lead to a wide range of strong biological effects.

CONCLUSION

Charge-altering post-translational modifications in the relatively unexplored globular histone core may provide a precision mechanism for controlling accessibility to the nucleosomal DNA.

Mapping allosteric linkage to channel gating by extracellular domains in the human epithelial sodium channel

The epithelial sodium channel (ENaC) mediates sodium absorption in lung, kidney, and colon epithelia. Channels in the ENaC/degenerin family possess an extracellular region that senses physicochemical changes in the extracellular milieu and allosterically regulates the channel opening. Proteolytic cleavage activates the ENaC opening, by the removal of specific segments in the finger domains of the α- and γ ENaC-subunits. Cleavage causes perturbations in the extracellular region that propagate to the channel gate. However, it is not known how the channel structure mediates the propagation of activation signals through the extracellular sensing domains. Here, to identify the structure-function determinants that mediate allosteric ENaC activation, we performed MD simulations, thiol modification of residues substituted by cysteine, and voltage-clamp electrophysiology recordings. Our simulations of an ENaC heterotetramer, α₁βα₂γ, in the proteolytically cleaved and uncleaved states revealed structural pathways in the α-subunit that are responsible for ENaC proteolytic activation. To validate these findings, we performed site-directed mutagenesis to introduce cysteine substitutions in the extracellular domains of the α-, β-, and γ ENaC-subunits. Insertion of a cysteine at the α-subunit Glu⁵⁵⁷ site, predicted to stabilize a closed state of ENaC, inhibited ENaC basal activity and retarded the kinetics of proteolytic activation by 2-fold. Our results suggest that the lower palm domain of αENaC is essential for ENaC activation. In conclusion, our integrated computational and experimental approach suggests key structure-function determinants for ENaC proteolytic activation and points toward a mechanistic model for the allosteric communication in the extracellular domains of the ENaC/degenerin family channels.

Engineering extrinsic disorder to control protein activity in living cells

Optogenetic and chemogenetic control of proteins has revealed otherwise inaccessible facets of signaling dynamics. Here, we use light- or ligand-sensitive domains to modulate the structural disorder of diverse proteins, thereby generating robust allosteric switches. Sensory domains were inserted into nonconserved, surface-exposed loops that were tight and identified computationally as allosterically coupled to active sites. Allosteric switches introduced into motility signaling proteins (kinases, guanosine triphosphatases, and guanine exchange factors) controlled conversion between conformations closely resembling natural active and inactive states, as well as modulated the morphodynamics of living cells. Our results illustrate a broadly applicable approach to design physiological protein switches.

Nucleosome Histone Tail Conformation and Dynamics: Impacts of Lysine Acetylation and a Nearby Minor Groove Benzo[a]pyrene-Derived Lesion

Histone tails in nucleosomes play critical roles in regulation of many biological processes, including chromatin compaction, transcription, and DNA repair. Moreover, post-translational modifications, notably lysine acetylation, are crucial to these functions. While the tails have been intensively studied, how the structures and dynamics of tails are impacted by the presence of a nearby bulky DNA lesion is a frontier research area, and how these properties are impacted by tail lysine acetylation remains unexplored. To obtain molecular insight, we have utilized all atom 3 μs molecular dynamics simulations of nucleosome core particles (NCPs) to determine the impact of a nearby DNA lesion, 10S (+)-trans-anti-B[a]P-N²-dG-the major adduct derived from the procarcinogen benzo[a]pyrene-on H2B tail behavior in unacetylated and acetylated states. We similarly studied lesion-free NCPs to investigate the normal properties of the H2B tail in both states. In the lesion-free NCPs, charge neutralization upon lysine acetylation causes release of the tail from the DNA. When the lesion is present, it stably engulfs part of the nearby tail, impairing the interactions between DNA and tail. With the tail in an acetylated state, the lesion still interacts with part of it, although unstably. The lesion's partial entrapment of the tail should hinder the tail from interacting with other nucleosomes, and other proteins such as acetylases, deacetylases, and acetyl-lysine binding proteins, and thus disrupt critical tail-governed processes. Hence, the lesion would impede tail functions modulated by acetylation or deacetylation, causing aberrant chromatin structures and impaired biological transactions such as transcription and DNA repair.

The role of histone tails in the nucleosome: a computational study

Histone tails play an important role in gene transcription and expression. We present here a systematic computational study of the role of histone tails in the nucleosome, using replica exchange molecular dynamics simulations with an implicit solvent model and different well-established force fields. We performed simulations for all four histone tails, H4, H3, H2A, and H2B, isolated and with inclusion of the nucleosome. The results confirm predictions of previous theoretical studies for the secondary structure of the isolated tails but show a strong dependence on the force field used. In the presence of the entire nucleosome for all force fields, the secondary structure of the histone tails is destabilized. Specific contacts are found between charged lysine and arginine residues and DNA phosphate groups and other binding sites in the minor and major DNA grooves. Using cluster analysis, we found a single dominant configuration of binding to DNA for the H4 and H2A histone tails, whereas H3 and H2B show multiple binding configurations with an equal probability. The leading stabilizing contribution for those binding configurations is the attractive interaction between the positively charged lysine and arginine residues and the negatively charged phosphate groups, and thus the resulting charge neutralization. Finally, we present results of molecular dynamics simulations in explicit solvent to confirm our conclusions. Results from both implicit and explicit solvent models show that large portions of the histone tails are not bound to DNA, supporting the complex role of these tails in gene transcription and expression and making them possible candidates for binding sites of transcription factors, enzymes, and other proteins.

Controlling Allosteric Networks in Proteins

Allosteric transition, defined as conformational changes induced by ligand binding, is one of the fundamental properties of proteins. Allostery has been observed and characterized in many proteins, and has been recently utilized to control protein function via regulation of protein activity. Here, we review the physical and evolutionary origin of protein allostery, as well as its importance to protein regulation, drug discovery, and biological processes in living systems. We describe recently developed approaches to identify allosteric pathways, connected sets of pairwise interactions that are responsible for propagation of conformational change from the ligand-binding site to a distal functional site. We then present experimental and computational protein engineering approaches for control of protein function by modulation of allosteric sites. As an example of application of these approaches, we describe a synergistic computational and experimental approach to rescue the cystic-fibrosis-associated protein cystic fibrosis transmembrane conductance regulator, which upon deletion of a single residue misfolds and causes disease. This example demonstrates the power of allosteric manipulation in proteins to both elucidate mechanisms of molecular function and to develop therapeutic strategies that rescue those functions. Allosteric control of proteins provides a tool to shine a light on the complex cascades of cellular processes and facilitate unprecedented interrogation of biological systems.

Coarse-Grained Protein Model Coupled with a Coarse-Grained Water Model: Molecular Dynamics Study of Polyalanine-Based Peptides

The coupling of a coarse-grained (CG) protein model with the CG water model developed by Marrink et al. (J. Phys. Chem. B 2004, 108, 750) is presented. The model was used in the molecular dynamics studies of Ac-(Ala)6-Xaa-(Ala)7-NHMe, Xaa = Ala, Leu, Val, and Gly. A Gly mutation in the middle of polyalanine is found to destabilize the helix and stabilize the hairpin by favoring a type-II' turn and probably to speed up hairpin folding. The simulations allow us to derive thermodynamic parameters of, in particular, the helical propensities (s) of amino acids in these polyalanine-based peptides. The calculated s values are 1.18 (Ala), 0.84 (Leu), 0.30 (Val), and <0.02 (Gly) at 291 K, in excellent agreement with experimental values (R(2)=0.970). Analyses using a structural approach method show that the helical propensity difference of these amino acids mainly comes from solvation effect. Leu and Val have lower helical propensities than Ala mainly because the larger side chains shield the solvation of helical structures, while Gly has a much poorer helical propensity mainly due to the much better solvation for the coil structures than for the helical structures. Overall, the model is at least about 10(2) times faster than current all-atom MD methods with explicit solvent.

Toward a Coarse-Grained Protein Model Coupled with a Coarse-Grained Solvent Model: Solvation Free Energies of Amino Acid Side Chains

Recently, we reported that molecular dynamics (MD) simulations using a coarse-grained (CG) peptide model coupled with a CG water model are able to reproduce many of the structural and thermodynamic features of short peptides with nonpolar side chains at 10(3) times the normal speed (JCTC, 2007, 3, 2146-2161). To further develop a CG protein model for MD simulations, we systematically parametrized the side chains of all 20 naturally occurring amino acids. We developed the parameters by fitting the dihedral potentials of 13 small molecules, the densities and self-solvation free energies of liquids of eight organic molecules, and the hydration free energies of 35 small organic molecules. In a set of 11 classes of compounds (105 in total) including alkanes, alcohols, ethers, ketones/aldehydes, amines, amides, aromatics, carboxylic acids, sulfides/thiols, alkyl ammoniums, and carboxylate ions, the average error in the calculated hydration free energies compared with experimental results is about 1.4 kJ/mol. The average error in the calculated transfer free energies of the 19 side-chain analogues of amino acids from cyclohexane to water is about 2.2 kJ/mol. These results are comparable to the results of all-atom models.

Molecular dynamics simulations: advances and applications

Molecular dynamics simulations have evolved into a mature technique that can be used effectively to understand macromolecular structure-to-function relationships. Present simulation times are close to biologically relevant ones. Information gathered about the dynamic properties of macromolecules is rich enough to shift the usual paradigm of structural bioinformatics from studying single structures to analyze conformational ensembles. Here, we describe the foundations of molecular dynamics and the improvements made in the direction of getting such ensemble. Specific application of the technique to three main issues (allosteric regulation, docking, and structure refinement) is discussed.

REMD Simulations Reveal the Dynamic Profile and Mechanism of Action of Deleterious, Rescuing, and Stabilizing Perturbations to NBD1 from CFTR

Cystic Fibrosis (CF) is a lethal, genetic disease caused by mutations to the CFTR chloride channel. The most common CF causing mutation is the deletion of F508 from the first Nucleotide Binding Domain (F508del-NBD1). This mutation leads to a thermally unstable domain and a misfolded, nonfunctioning CFTR. Replica Exchange MD simulations were used to simulate seven NBD1 constructs including wt and F508del-NBD1 both alone and in the presence of known rescuing mutations as well as F508del-NBD1 in complex with a known small (ligand) stabilizer. Analyzing the resulting trajectories suggests that differences in the biochemical properties of the constructs result from local and coupled differences in their dynamic profiles. A comparative analysis of these profiles as well as of the resulting trajectories reveals how the different perturbations exert their deleterious, rescuing, and stabilizing effects on NBD1. These simulations may therefore be useful for the design and mechanism-of-action analysis of new NBD1 stabilizers.

Programming of Essential Hypertension: What Pediatric Cardiologists Need to Know

Hypertension is recognized as one of the major contributing factors to cardiovascular disease, but its etiology remains incompletely understood. Known genetic and environmental influences can only explain a small part of the variability in cardiovascular disease risk. The missing heritability is currently one of the most important challenges in blood pressure and hypertension genetics. Recently, some promising approaches have emerged that move beyond the DNA sequence and focus on identification of blood pressure genes regulated by epigenetic mechanisms such as DNA methylation, histone modification and microRNAs. This review summarizes information on gene-environmental interactions that lead toward the developmental programming of hypertension with specific reference to epigenetics and provides pediatricians and pediatric cardiologists with a more complete understanding of its pathogenesis.

Partial Unwrapping and Histone Tail Dynamics in Nucleosome Revealed by Coarse-Grained Molecular Simulations

Nucleosomes, basic units of chromatin, are known to show spontaneous DNA unwrapping dynamics that are crucial for transcriptional activation, but its structural details are yet to be elucidated. Here, employing a coarse-grained molecular model that captures residue-level structural details up to histone tails, we simulated equilibrium fluctuations and forced unwrapping of single nucleosomes at various conditions. The equilibrium simulations showed spontaneous unwrapping from outer DNA and subsequent rewrapping dynamics, which are in good agreement with experiments. We found several distinct partially unwrapped states of nucleosomes, as well as reversible transitions among these states. At a low salt concentration, histone tails tend to sit in the concave cleft between the histone octamer and DNA, tightening the nucleosome. At a higher salt concentration, the tails tend to bound to the outer side of DNA or be expanded outwards, which led to higher degree of unwrapping. Of the four types of histone tails, H3 and H2B tail dynamics are markedly correlated with partial unwrapping of DNA, and, moreover, their contributions were distinct. Acetylation in histone tails was simply mimicked by changing their charges, which enhanced the unwrapping, especially markedly for H3 and H2B tails.

Nucleosomes, folding units of chromatin, wrap DNA about 1.75 turns and provide bottlenecks for transcription. Recent experiments showed that nucleosomes are not rigid but dynamic, showing spontaneous and partial unwrapping which is thus important for transcriptional activation. Experimentally, however, one cannot directly watch DNA unwrapping at high resolution. On the other hand, molecular dynamics simulations have high spatio-temporal resolution and thus can be powerful and complementary to experiments. Here, we put forward coarse-grained modeling of protein-DNA interactions at residue-level resolution, which is rather generic and thus can be applied to any protein-DNA complexes. By this method, we could reveal spontaneous and salt-concentration dependent partial unwrapping of DNA from nucleosomes. In addition to consistency with single molecule experiments, the simulation showed multiple and distinct intermediate states of unwrapping. Interestingly, partial unwrapping of DNA is correlated with certain parts of histone tail dynamics. Deleting positive charges in histone tails that mimics histone acetylation facilitated partial unwrapping, most significantly for H3 and H2B.

Rational Coupled Dynamics Network Manipulation Rescues Disease-Relevant Mutant Cystic Fibrosis Transmembrane Conductance Regulator

Many cellular functions necessary for life are tightly regulated by protein allosteric conformational change, and correlated dynamics between protein regions has been found to contribute to the function of proteins not previously considered allosteric. The ability to map and control such dynamic coupling would thus create opportunities for the extension of current therapeutic design strategy. Here, we present an approach to determine the networks of residues involved in the transfer of correlated motion across a protein, and apply our approach to rescue disease-causative mutant cystic fibrosis transmembrane regulator (CFTR) ion channels, ΔF508 and ΔI507, which together constitute over 90% of cystic fibrosis cases. We show that these mutations perturb dynamic coupling within the first nucleotide-binding domain (NBD1), and uncover a critical residue that mediates trans-domain coupled dynamics. By rationally designing a mutation to this residue, we improve aberrant dynamics of mutant CFTR as well as enhance surface expression and function of both mutants, demonstrating the rescue of a disease mutation by rational correction of aberrant protein dynamics.

The histone H3 N-terminal tail: a computational analysis of the free energy landscape and kinetics

Extensive molecular dynamics simulations and Markov State models are used to characterize the free energy landscape and kinetics of the histone H3 N-terminal tail, which plays a critical role in regulating chromatin dynamics and gene activity.

Close encounters with DNA

Over the past ten years, the all-atom molecular dynamics method has grown in the scale of both systems and processes amenable to it and in its ability to make quantitative predictions about the behavior of experimental systems. The field of computational DNA research is no exception, witnessing a dramatic increase in the size of systems simulated with atomic resolution, the duration of individual simulations and the realism of the simulation outcomes. In this topical review, we describe the hallmark physical properties of DNA from the perspective of all-atom simulations. We demonstrate the amazing ability of such simulations to reveal the microscopic physical origins of experimentally observed phenomena. We also discuss the frustrating limitations associated with imperfections of present atomic force fields and inadequate sampling. The review is focused on the following four physical properties of DNA: effective electric charge, response to an external mechanical force, interaction with other DNA molecules and behavior in an external electric field.

Rational design of a ligand-controlled protein conformational switch

Design of a regulatable multistate protein is a challenge for protein engineering. Here we design a protein with a unique topology, called uniRapR, whose conformation is controlled by the binding of a small molecule. We confirm switching and control ability of uniRapR in silico, in vitro, and in vivo. As a proof of concept, uniRapR is used as an artificial regulatory domain to control activity of kinases. By activating Src kinase using uniRapR in single cells and whole organism, we observe two unique phenotypes consistent with its role in metastasis. Activation of Src kinase leads to rapid induction of protrusion with polarized spreading in HeLa cells, and morphological changes with loss of cell-cell contacts in the epidermal tissue of zebrafish. The rational creation of uniRapR exemplifies the strength of computational protein design, and offers a powerful means for targeted activation of many pathways to study signaling in living organisms.

Structural role of RKS motifs in chromatin interactions: a molecular dynamics study of HP1 bound to a variably modified histone tail

The current understanding of epigenetic signaling assigns a central role to post-translational modifications that occur in the histone tails. In this context, it has been proposed that methylation of K9 and phosphorylation of S10 in the tail of histone H3 represent a binary switch that controls its reversible association to heterochromatin protein 1 (HP1). To test this hypothesis, we performed a comprehensive molecular dynamics study in which we analyzed a crystallographically defined complex that involves the HP1 chromodomain and an H3 tail peptide. Microsecond-long simulations show that the binding of the trimethylated K9 H3 peptide in the aromatic cage of HP1 is only slightly affected by S10 phosphorylation, because the modified K9 and S10 do not interact directly with one another. Instead, the phosphate group of S10 seems to form a persistent intramolecular salt bridge with R8, an interaction that can provoke a major structural change and alter the hydrogen-bonding regime in the H3-HP1 complex. These observations suggest that interactions between adjacent methyl-lysine and phosphoserine side chains do not by themselves provide a binary switch in the H3-HP1 system, but arginine-phosphoserine interactions, which occur in both histones and nonhistone proteins in the context of a conserved RKS motif, are likely to serve a key regulatory function.

Conformational changes relevant to channel activity and folding within the first nucleotide binding domain of the cystic fibrosis transmembrane conductance regulator

Deletion of Phe-508 (F508del) in the first nucleotide binding domain (NBD1) of the cystic fibrosis transmembrane conductance regulator (CFTR) leads to defects in folding and channel gating. NMR data on human F508del NBD1 indicate that an H620Q mutant, shown to increase channel open probability, and the dual corrector/potentiator CFFT-001 similarly disrupt interactions between β-strands S3, S9, and S10 and the C-terminal helices H8 and H9, shifting a preexisting conformational equilibrium from helix to coil. CFFT-001 appears to interact with β-strands S3/S9/S10, consistent with docking simulations. Decreases in T(m) from differential scanning calorimetry with H620Q or CFFT-001 suggest direct compound binding to a less thermostable state of NBD1. We hypothesize that, in full-length CFTR, shifting the conformational equilibrium to reduce H8/H9 interactions with the uniquely conserved strands S9/S10 facilitates release of the regulatory region from the NBD dimerization interface to promote dimerization and thereby increase channel open probability. These studies enabled by our NMR assignments for F508del NBD1 provide a window into the conformational fluctuations within CFTR that may regulate function and contribute to folding energetics.

Allosteric modulation balances thermodynamic stability and restores function of ΔF508 CFTR

Most cystic fibrosis is caused by a deletion of a single residue (F508) in CFTR (cystic fibrosis transmembrane conductance regulator) that disrupts the folding and biosynthetic maturation of the ion channel protein. Progress towards understanding the underlying mechanisms and overcoming the defect remains incomplete. Here, we show that the thermal instability of human ΔF508 CFTR channel activity evident in both cell-attached membrane patches and planar phospholipid bilayers is not observed in corresponding mutant CFTRs of several non-mammalian species. These more stable orthologs are distinguished from their mammalian counterparts by the substitution of proline residues at several key dynamic locations in first N-terminal nucleotide-binding domain (NBD1), including the structurally diverse region, the γ-phosphate switch loop, and the regulatory insertion. Molecular dynamics analyses revealed that addition of the prolines could reduce flexibility at these locations and increase the temperatures of unfolding transitions of ΔF508 NBD1 to that of the wild type. Introduction of these prolines experimentally into full-length human ΔF508 CFTR together with the already recognized I539T suppressor mutation, also in the structurally diverse region, restored channel function and thermodynamic stability as well as its trafficking to and lifetime at the cell surface. Thus, while cellular manipulations that circumvent its culling by quality control systems leave ΔF508 CFTR dysfunctional at physiological temperature, restoration of the delicate balance between the dynamic protein's inherent stability and channel activity returns a near-normal state.

Coarse-grain Protein Models

Coarse-graining is a powerful approach for modeling biomolecules that, over the last few decades, has been extensively applied to proteins. Coarse-grain models offer access to large systems and to slow processes without becoming computationally unmanageable. In addition, they are very versatile, enabling both the protein representation and the energy function to be adapted to the biological problem in hand. This review concentrates on modeling soluble proteins and their assemblies. It presents an overview of the coarse-grain representations, of the associated interaction potentials, and of the optimization procedures used to define them. It then shows how coarse-grain models have been used to understand processes involving proteins, from their initial folding to their functional properties, their binary interactions, and the assembly of large complexes.

Energetic and structural basis for activation of the epithelial sodium channel by matriptase

Limited proteolysis, accomplished by endopeptidases, is a ubiquitous phenomenon underlying the regulation and activation of many enzymes, receptors, and other proteins synthesized as inactive precursors. Serine proteases make up one of the largest and most conserved families of endopeptidases involved in diverse cellular activities, including wound healing, blood coagulation, and immune responses. Heteromeric α,β,γ-epithelial sodium channels (ENaC) associated with diseases like cystic fibrosis and Liddle's syndrome are irreversibly stimulated by membrane-anchored proteases (MAPs) and furin-like convertases. Matriptase/channel activating protease-3 (CAP3) is one of the several MAPs that potently activate ENaC. Despite identification of protease cleavage sites, the basis for the enhanced susceptibility of α- and γ-ENaC to proteases remains elusive. Here, we elucidate the energetic and structural bases for activation of ENaC by CAP3. We find a region near the γ-ENaC furin site that has previously not been identified as a critical cleavage site for CAP3-mediated stimulation. We also report that CAP3 mediates cleavage of ENaC at basic residues downstream of the furin site. Our results indicate that surface proteases alone are sufficient to fully activate uncleaved ENaC and explain how ENaC in epithelia expressing surface-active proteases can appear refractory to soluble proteases. Our results support a model in which proteases prime ENaC for activation by cleaving at the furin site, and cleavage at downstream sites is accomplished by membrane surface proteases or extracellular soluble proteases. On the basis of our results, we propose a dynamics-driven "anglerfish" mechanism that explains less stringent sequence requirements for substrate recognition and cleavage by matriptase than by furin.

Unwrapping of nucleosomal DNA ends: a multiscale molecular dynamics study

To permit access to DNA-binding proteins involved in the control and expression of the genome, the nucleosome undergoes structural remodeling including unwrapping of nucleosomal DNA segments from the nucleosome core. Here we examine the mechanism of DNA dissociation from the nucleosome using microsecond timescale coarse-grained molecular dynamics simulations. The simulations exhibit short-lived, reversible DNA detachments from the nucleosome and long-lived DNA detachments not reversible on the timescale of the simulation. During the short-lived DNA detachments, 9 bp dissociate at one extremity of the nucleosome core and the H3 tail occupies the space freed by the detached DNA. The long-lived DNA detachments are characterized by structural rearrangements of the H3 tail including the formation of a turn-like structure at the base of the tail that sterically impedes the rewrapping of DNA on the nucleosome surface. Removal of the H3 tails causes the long-lived detachments to disappear. The physical consistency of the CG long-lived open state was verified by mapping a CG structure representative of this state back to atomic resolution and performing molecular dynamics as well as by comparing conformation-dependent free energies. Our results suggest that the H3 tail may stabilize the nucleosome in the open state during the initial stages of the nucleosome remodeling process.

Modelling chromatin structure and dynamics: status and prospects

The packaging of genomic DNA into chromatin in the eukaryotic cell nucleus demands extensive compaction. This requires attractive nucleosome-nucleosome interactions to overcome repulsion between the negatively charged DNA segments as well as other constraints. At the same time, DNA must be dynamically accessible to the cellular machinery that operates on it. Recent progress in the experimental characterisation of the higher order structure and dynamics of well-defined chromatin fibres has stimulated the attempts at theoretical description of chromatin and the nucleosome. Here we review the present status of chromatin modelling, with particular emphasis on coarse-grained computer simulation models, the role of electrostatic interactions, and discuss future perspectives in the field.

Dissecting DNA-histone interactions in the nucleosome by molecular dynamics simulations of DNA unwrapping

The nucleosome complex of DNA wrapped around a histone protein octamer organizes the genome of eukaryotes and regulates the access of protein factors to the DNA. We performed molecular dynamics simulations of the nucleosome in explicit water to study the dynamics of its histone-DNA interactions. A high-resolution histone-DNA interaction map was derived that revealed a five-nucleotide periodicity, in which the two DNA strands of the double helix made alternating contacts. On the 100-ns timescale, the histone tails mostly maintained their initial positions relative to the DNA, and the spontaneous unwrapping of DNA was limited to 1-2 basepairs. In steered molecular dynamics simulations, external forces were applied to the linker DNA to investigate the unwrapping pathway of the nucleosomal DNA. In comparison with a nucleosome without the unstructured N-terminal histone tails, the following findings were obtained: 1), Two main barriers during unwrapping were identified at DNA position ±70 and ±45 basepairs relative to the central DNA basepair at the dyad axis. 2), DNA interactions of the histone H3 N-terminus and the histone H2A C-terminus opposed the initiation of unwrapping. 3), The N-terminal tails of H2A, H2B, and H4 counteracted the unwrapping process at later stages and were essential determinants of nucleosome dynamics. Our detailed analysis of DNA-histone interactions revealed molecular mechanisms for modulating access to nucleosomal DNA via conformational rearrangements of its structure.

Nuclear remodeling as a mechanism for genomic instability in cancer

This chapter focuses on the three-dimensional organization of the nucleus in normal, early genomically unstable, and tumor cells. A cause-consequence relationship is discussed between nuclear alterations and the resulting genomic rearrangements. Examples are presented from studies on conditional Myc deregulation, experimental tumorigenesis in mouse plasmacytoma, nuclear remodeling in Hodgkin's lymphoma, and in adult glioblastoma. A model of nuclear remodeling is proposed for cancer progression in multiple myeloma. Current models of nuclear remodeling are described, including our model of altered nuclear architecture and the onset of genomic instability.

Discrete molecular dynamics distinguishes nativelike binding poses from decoys in difficult targets

Virtual screening is one of the major tools used in computer-aided drug discovery. In structure-based virtual screening, the scoring function is critical to identifying the correct docking pose and accurately predicting the binding affinities of compounds. However, the performance of existing scoring functions has been shown to be uneven for different targets, and some important drug targets have proven especially challenging. In these targets, scoring functions cannot accurately identify the native or near-native binding pose of the ligand from among decoy poses, which affects both the accuracy of the binding affinity prediction and the ability of virtual screening to identify true binders in chemical libraries. Here, we present an approach to discriminating native poses from decoys in difficult targets for which several scoring functions failed to correctly identify the native pose. Our approach employs Discrete Molecular Dynamics simulations to incorporate protein-ligand dynamics and the entropic effects of binding. We analyze a collection of poses generated by docking and find that the residence time of the ligand in the native and nativelike binding poses is distinctly longer than that in decoy poses. This finding suggests that molecular simulations offer a unique approach to distinguishing the native (or nativelike) binding pose from decoy poses that cannot be distinguished using scoring functions that evaluate static structures. The success of our method emphasizes the importance of protein-ligand dynamics in the accurate determination of the binding pose, an aspect that is not addressed in typical docking and scoring protocols.

Crystal structure of the plant epigenetic protein arginine methyltransferase 10

Protein arginine methyltransferase 10 (PRMT10) is a type I arginine methyltransferase that is essential for regulating flowering time in Arabidopsis thaliana. We present a 2.6 Å resolution crystal structure of A. thaliana PRMT 10 (AtPRMT10) in complex with a reaction product, S-adenosylhomocysteine. The structure reveals a dimerization arm that is 12-20 residues longer than PRMT structures elucidated previously; as a result, the essential AtPRMT10 dimer exhibits a large central cavity and a distinctly accessible active site. We employ molecular dynamics to examine how dimerization facilitates AtPRMT10 motions necessary for activity, and we show that these motions are conserved in other PRMT enzymes. Finally, functional data reveal that the 10 N-terminal residues of AtPRMT10 influence substrate specificity, and that enzyme activity is dependent on substrate protein sequences distal from the methylation site. Taken together, these data provide insights into the molecular mechanism of AtPRMT10, as well as other members of the PRMT family of enzymes. They highlight differences between AtPRMT10 and other PRMTs but also indicate that motions are a conserved element of PRMT function.