Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userIntegral equation models for solvent in macromolecular crystals
358
0
0.0
(
0
)
Authors
Jonathon G. Gray
Ask ChatGPT about the research
No Arabic abstract
Solvent can occupy up to ~70% of macromolecular crystals and hence having models that predict solvent distributions in periodic systems could improve in the interpretation of crystallographic data. Yet there are few implicit solvent models applicable to periodic solutes while crystallographic structures are commonly solved assuming a flat solvent model. Here we present a newly-developed periodic version of the 3D-RISM integral equation method that is able to solve for efficiently and describe accurately water and ions distributions in periodic systems; the code can compute accurate gradients that can be used in minimizations or molecular dynamics simulations. The new method includes an extension of the OZ equation needed to yield charge neutrality for charged solutes which requires an additional contribution to the excess chemical potential that has not been previously identified; this is an important consideration for nucleic acids or any other charged system where most or all of the counter- and co-ions are part of the disordered solvent. We present of several calculations of protein, RNA and small molecule crystals to show that X-ray scattering intensities and solvent structure predicted by the periodic 3D-RISM solvent model are in closer agreement with experiment than are intensities computed using the default flat solvent model in the refmac5 or phenix refinement programs, with the greatest improvement in the 2 to 4 {AA} range. Prospects for incorporating integral equation models into crystallographic refinement are discussed.
rate research
Read More
We present a protocol for the fully automated construction of quantum mechanical-(QM)-classical hybrid models by extending our previously reported approach on self-parametrizing system-focused atomistic models (SFAM) J. Chem. Theory Comput. 2020, 16, 1646]. In this QM/SFAM approach, the size and composition of the QM region is evaluated in an automated manner based on first principles so that the hybrid model describes the atomic forces in the center of the QM region accurately. This entails the automated construction and evaluation of differently sized QM regions with a bearable computational overhead that needs to be paid for automated validation procedures. Applying SFAM for the classical part of the model eliminates any dependence on pre-existing parameters due to its system-focused quantum mechanically derived parametrization. Hence, QM/SFAM is capable of delivering a high fidelity and complete automation. Furthermore, since SFAM parameters are generated for the whole system, our ansatz allows for a convenient re-definition of the QM region during a molecular exploration. For this purpose, a local re-parametrization scheme is introduced, which efficiently generates additional classical parameters on the fly when new covalent bonds are formed (or broken) and moved to the classical region.
Transitions between different conformational states are ubiquitous in proteins, being involved in signaling, catalysis and other fundamental activities in cells. However, modeling those processes is extremely difficult, due to the need of efficiently exploring a vast conformational space in order to seek for the actual transition path for systems whose complexity is already high in the steady states. Here we report a strategy that simplifies this task attacking the complexity on several sides. We first apply a minimalist coarse-grained model to Calmodulin, based on an empirical force field with a partial structural bias, to explore the transition paths between the apo-closed state and the Ca-bound open state of the protein. We then select representative structures along the trajectory based on a structural clustering algorithm and build a cleaned-up trajectory with them. We finally compare this trajectory with that produced by the online tool MinActionPath, by minimizing the action integral using a harmonic network model, and with that obtained by the PROMPT morphing method, based on an optimal mass transportation-type approach including physical constraints. The comparison is performed both on the structural and energetic level, using the coarse-grained and the atomistic force fields upon reconstruction. Our analysis indicates that this method returns trajectories capable of exploring intermediate states with physical meaning, retaining a very low computational cost, which can allow systematic and extensive exploration of the multi-stable proteins transition pathways.
Recent studies of the hydration of micro- and nanoscale solutes have demonstrated a strong {it coupling} between hydrophobic, dispersion and electrostatic contributions, a fact not accounted for in current implicit solvent models. We present a theoretical formalism which accounts for coupling by minimizing the Gibbs free energy with respect to a solvent volume exclusion function. The solvent accessible surface is output of our theory. Our method is illustrated with the hydration of alkane-assembled solutes on different length scales, and captures the strong sensitivity to the particular form of the solute-solvent interactions in agreement with recent computer simulations.
In this work, a systematic protocol is proposed to automatically parametrize implicit solvent models with polar and nonpolar components. The proposed protocol utilizes the classical Poisson model or the Kohn-Sham density functional theory (KSDFT) based polarizable Poisson model for modeling polar solvation free energies. For the nonpolar component, either the standard model of surface area, molecular volume, and van der Waals interactions, or a model with atomic surface areas and molecular volume is employed. Based on the assumption that similar molecules have similar parametrizations, we develop scoring and ranking algorithms to classify solute molecules. Four sets of radius parameters are combined with four sets of charge force fields to arrive at a total of 16 different parametrizations for the Poisson model. A large database with 668 experimental data is utilized to validate the proposed protocol. The lowest leave-one-out root mean square (RMS) error for the database is 1.33k cal/mol. Additionally, five subsets of the database, i.e., SAMPL0-SAMPL4, are employed to further demonstrate that the proposed protocol offers some of the best solvation predictions. The optimal RMS errors are 0.93, 2.82, 1.90, 0.78, and 1.03 kcal/mol, respectively for SAMPL0, SAMPL1, SAMPL2, SAMPL3, and SAMPL4 test sets. These results are some of the best, to our best knowledge.
Two numerical methods are used to evaluate the relativistic spectrum of the two-centre Coulomb problem (for the $H_{2}^{+}$ and $Th_{2}^{179+}$ diatomic molecules) in the fixed nuclei approximation by solving the single particle time-independent Dirac equation. The first one is based on a min-max principle and uses a two-spinor formulation as a starting point. The second one is the Rayleigh-Ritz variational method combined with kinematically balanced basis functions. Both methods use a B-spline basis function expansion. We show that accurate results can be obtained with both methods and that no spurious states appear in the discretization process.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
