How to Determine Accurate Conformational Ensembles by Metadynamics Metainference: A Chignolin Study Case.
Ontology highlight
ABSTRACT:
SUBMITTER: Paissoni C
PROVIDER: S-EPMC8187852 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:A major challenge in the development of antibody biotherapeutics is their tendency to aggregate. One root cause for aggregation is exposure of hydrophobic surface regions to the solvent. Many current techniques predict the relative aggregation propensity of antibodies via precalculated scales for the hydrophobicity or aggregation propensity of single amino acids. However, those scales cannot describe the nonadditive effects of a residue's surrounding on its hydrophobicity. Therefore, they are inherently limited in their ability to describe the impact of subtle differences in molecular structure on the overall hydrophobicity. Here, we introduce a physics-based approach to describe hydrophobicity in terms of the hydration free energy using grid inhomogeneous solvation theory (GIST). We apply this method to assess the effects of starting structures, conformational sampling, and protonation states on the hydrophobicity of antibodies. Our results reveal that high-quality starting structures, i.e., crystal structures, are crucial for the prediction of hydrophobicity and that conformational sampling can compensate errors introduced by the starting structure. On the other hand, sampling of protonation states only leads to good results when combined with high-quality structures, whereas it can even be detrimental otherwise. We conclude by pointing out that a single static homology model may not be adequate for predicting hydrophobicity.
| S-EPMC7820740 | biostudies-literature
Project description:Earliest events in the aggregation process, such as single molecule reconfiguration, are extremely important and the most difficult to characterize in experiments. To this end, we have used well-tempered bias exchange metadynamics simulations to determine the equilibrium ensembles of an insulin molecule under amyloidogenic conditions of low pH and high temperature. A bin-based clustering method that uses statistics accumulated in bias exchange metadynamics trajectories was employed to construct a detailed thermodynamic and kinetic model of insulin folding. The highest lifetime, lowest free-energy ensemble identified consisted of native conformations adopted by a folded insulin monomer in solution, namely, the R-, the Rf-, and the T-states of insulin. The lowest free-energy structure had a root mean square deviation of only 0.15 nm from native x-ray structure. The second longest-lived metastable state was an unfolded, compact monomer with little similarity to the native structure. We have identified three additional long-lived, metastable states from the bin-based model. We then carried out an exhaustive structural characterization of metastable states on the basis of tertiary contact maps and per-residue accessible surface areas. We have also determined the lowest free-energy path between two longest-lived metastable states and confirm earlier findings of non-two-state folding for insulin through a folding intermediate. The ensemble containing the monomeric intermediate retained 58% of native hydrophobic contacts, however, accompanied by a complete loss of native secondary structure. We have discussed the relative importance of nativelike versus nonnative tertiary contacts for the folding transition. We also provide a simple measure to determine the importance of an individual residue for folding transition. Finally, we have compared and contrasted this intermediate with experimental data obtained in spectroscopic, crystallographic, and calorimetric measurements during early stages of insulin aggregation. We have also determined stability of monomeric insulin by incubation at a very low concentration to isolate protein-protein interaction effects.
| S-EPMC5406382 | biostudies-literature
Project description:A detailed description of the events ruling ligand/protein interaction and an accurate estimation of the drug affinity to its target is of great help in speeding drug discovery strategies. We have developed a metadynamics-based approach, named funnel metadynamics, that allows the ligand to enhance the sampling of the target binding sites and its solvated states. This method leads to an efficient characterization of the binding free-energy surface and an accurate calculation of the absolute protein-ligand binding free energy. We illustrate our protocol in two systems, benzamidine/trypsin and SC-558/cyclooxygenase 2. In both cases, the X-ray conformation has been found as the lowest free-energy pose, and the computed protein-ligand binding free energy in good agreement with experiments. Furthermore, funnel metadynamics unveils important information about the binding process, such as the presence of alternative binding modes and the role of waters. The results achieved at an affordable computational cost make funnel metadynamics a valuable method for drug discovery and for dealing with a variety of problems in chemistry, physics, and material science.
| S-EPMC3631651 | biostudies-literature
Project description:Molecular docking algorithms suggest possible structures for molecular complexes. They are used to model biological function and to discover potential ligands. A present challenge for docking algorithms is the treatment of molecular flexibility. Here, the rigid body program, DOCK, is modified to allow it to rapidly fit multiple conformations of ligands. Conformations of a given molecule are pre-calculated in the same frame of reference, so that each conformer shares a common rigid fragment with all other conformations. The ligand conformers are then docked together, as an ensemble, into a receptor binding site. This takes advantage of the redundancy present in differing conformers of the same molecule. The algorithm was tested using three organic ligand protein systems and two protein-protein systems. Both the bound and unbound conformations of the receptors were used. The ligand ensemble method found conformations that resembled those determined in X-ray crystal structures (RMS values typically less than 1.5 A). To test the method's usefulness for inhibitor discovery, multi-compound and multi-conformer databases were screened for compounds known to bind to dihydrofolate reductase and compounds known to bind to thymidylate synthase. In both cases, known inhibitors and substrates were identified in conformations resembling those observed experimentally. The ligand ensemble method was 100-fold faster than docking a single conformation at a time and was able to screen a database of over 34 million conformations from 117,000 molecules in one to four CPU days on a workstation.
| S-EPMC2143983 | biostudies-other
Project description:We propose a multiscale exploration method to characterize the conformational space populated by a protein at equilibrium. The method efficiently obtains a large set of equilibrium conformations in two stages: first exploring the entire space at a coarse-grained level of detail, then narrowing a refined exploration to selected low-energy regions. The coarse-grained exploration periodically adds all-atom detail to selected conformations to ensure that the search leads to regions which maintain low energies in all-atom detail. The second stage reconstructs selected low-energy coarse-grained conformations in all-atom detail. A low-dimensional energy landscape associated with all-atom conformations allows focusing the exploration to energy minima and their conformational ensembles. The lowest energy ensembles are enriched with additional all-atom conformations through further multiscale exploration. The lowest energy ensembles obtained from the application of the method to three different proteins correctly capture the known functional states of the considered systems.
| S-EPMC3164158 | biostudies-literature
Project description:DNA methyltransferases play crucial roles in establishing and maintenance of DNA methylation, which is an important epigenetic mark. Flipping the target cytosine out of the DNA helical stack and into the active site of protein provides DNA methyltransferases with an opportunity to access and modify the genetic information hidden in DNA. To investigate the conversion process of base flipping in the HhaI methyltransferase (M.HhaI), we performed different molecular simulation approaches on M.HhaI-DNA-S-adenosylhomocysteine ternary complex. The results demonstrate that the nonspecific binding of DNA to M.HhaI is initially induced by electrostatic interactions. Differences in chemical environment between the major and minor grooves determine the orientation of DNA. Gln237 at the target recognition loop recognizes the GCGC base pair from the major groove side by hydrogen bonds. In addition, catalytic loop motion is a key factor during this process. Our study indicates that base flipping is likely to be an "induced-fit" process. This study provides a solid foundation for future studies on the discovery and development of mechanism-based DNA methyltransferases regulators.
| S-EPMC4090504 | biostudies-literature
Project description:Multidimensional heteronuclear NMR approaches can provide nearly complete sequential signal assignments of isotopically enriched biomolecules. The availability of assignments together with measurements of spin relaxation rates, residual spin interactions, J-couplings and chemical shifts provides information at atomic resolution about internal dynamics on timescales ranging from ps to ms, both in solution and in the solid state. However, due to the complexity of biomolecules, it is not possible to extract a unique atomic-resolution description of biomolecular motions even from extensive NMR data when many conformations are sampled on multiple timescales. For this reason, powerful computational approaches are increasingly applied to large NMR data sets to elucidate conformational ensembles sampled by biomolecules. In the past decade, considerable attention has been directed at an important class of biomolecules that function by binding to a wide variety of target molecules. Questions of current interest are: "Does the free biomolecule sample a conformational ensemble that encompasses the conformations found when it binds to various targets; and if so, on what time scale is the ensemble sampled?" This article reviews recent efforts to answer these questions, with a focus on comparing ensembles obtained for the same biomolecules by different investigators. A detailed comparison of results obtained is provided for three biomolecules: ubiquitin, calmodulin and the HIV-1 trans-activation response RNA.
| S-EPMC4325279 | biostudies-literature
Project description:Intrinsically disordered proteins (IDPs) are important for signaling and regulatory pathways. In contrast to folded proteins, they sample a diverse conformational space. IDPs have residue ranges within a sequence that have been referred to as molecular recognition features (MoRFs). A MoRF can be viewed as contiguous residues exhibiting a conformational disorder that become ordered upon binding to another protein or ligand. In this work, we introduce a structural characterization of MoRFs based on entropy and mutual information (MI). In this view, a MoRF is a set of contiguous residues that exhibit a large entropy (from rotameric residue sampling) and large MI, the latter indicating a dependence among the residues' rotameric sampling comprising the MoRF. The methodology is first applied to a number of ubiquitin ensembles that were obtained based on nuclear magnetic resonance experiments. One is a denatured Ub ensemble that has a large entropy for various unitSizes (number of contiguous residues) but essentially zero MI, indicting no dependence among the residue rotamer sampling. Another ensemble does exhibit extensive regions along the sequence where there are MoRFs centered on nonsecondary structure regions. The MoRFs are present for unitSizes 2-10. That a substantial number of MoRFs are present in Ub strongly suggests a conformational selection mechanism for this protein. Two additional ensembles for the cyclin-dependent kinase inhibitor Sic1 and for the amyloid protein α-synuclein, which have been shown to be IDPs, are also analyzed. Both exhibit MoRF-like character.
| S-EPMC6644992 | biostudies-literature
Project description:Casein kinase CK2 is an essential enzyme in higher organisms, catalyzing the transfer of the ? phosphate from ATP to serine and threonine residues on protein substrates. In a number of animal tumors, CK2 activity has been shown to escape normal cellular control, making it a potential target for cancer therapy. Several crystal structures of human CK2 have been published with different conformations for the CK2? catalytic subunit. This variability reflects a high flexibility for two regions of CK2?: the interdomain hinge region, and the glycine-rich loop (p-loop). Here, we present a computational study simulating the equilibrium between three conformations involving these regions. Simulations were performed using well-tempered metadynamics combined with a path collective variables approach. This provides a reference pathway describing the conformational changes being studied, based on analysis of free energy surfaces. The free energies of the three conformations were found to be close and the paths proposed had low activation barriers. Our results indicate that these conformations can exist in water. This information should be useful when designing inhibitors specific to one conformation.
| S-EPMC4026775 | biostudies-literature
Project description:Conformational changes in nucleic acids play a key role in the way genetic information is stored, transferred, and processed in living cells. Here, we describe new approaches that employ a broad range of experimental data, including NMR-derived chemical shifts and residual dipolar couplings, small-angle X-ray scattering, and computational approaches such as molecular dynamics simulations to determine ensembles of DNA and RNA at atomic resolution. We review the complementary information that can be obtained from diverse sets of data and the various methods that have been developed to combine these data with computational methods to construct ensembles and assess their uncertainty. We conclude by surveying RNA and DNA ensembles determined using these methods, highlighting the unique physical and functional insights obtained so far.
| S-EPMC4748180 | biostudies-other