Nuclear magnetic resonance signal chemical shifts and molecular simulations: a multidisciplinary approach to modeling copper protein structures

Can multiscale simulations unravel the function of metallo-enzymes to improve knowledge-based drug discovery?

Metallo-enzymes are a large class of biomolecules promoting specialized chemical reactions. Quantum-classical quantum mechanics/molecular mechanics molecular dynamics, describing the metal site at quantum mechanics level, while accounting for the rest of system at molecular mechanics level, has an accessible time-scale limited by its computational cost. Hence, it must be integrated with classical molecular dynamics and enhanced sampling simulations to disentangle the functions of metallo-enzymes. In this review, we provide an overview of these computational methods and their capabilities. In particular, we will focus on some systems such as CYP19A1 a Fe-dependent enzyme involved in estrogen biosynthesis, and on Mg2+-dependent DNA/RNA processing enzymes/ribozymes and the spliceosome, a protein-directed ribozyme. This information may guide the discovery of drug-like molecules and genetic manipulation tools.

Solution NMR refinement of a metal ion bound protein using metal ion inclusive restrained molecular dynamics methods

Correctly calculating the structure of metal coordination sites in a protein during the process of nuclear magnetic resonance (NMR) structure determination and refinement continues to be a challenging task. In this study, we present an accurate and convenient means by which to include metal ions in the NMR structure determination process using molecular dynamics (MD) simulations constrained by NMR-derived data to obtain a realistic and physically viable description of the metal binding site(s). This method provides the framework to accurately portray the metal ions and its binding residues in a pseudo-bond or dummy-cation like approach, and is validated by quantum mechanical/molecular mechanical (QM/MM) MD calculations constrained by NMR-derived data. To illustrate this approach, we refine the zinc coordination complex structure of the zinc sensing transcriptional repressor protein Staphylococcus aureus CzrA, generating over 130 ns of MD and QM/MM MD NMR-data compliant sampling. In addition to refining the first coordination shell structure of the Zn(II) ion, this protocol benefits from being performed in a periodically replicated solvation environment including long-range electrostatics. We determine that unrestrained (not based on NMR data) MD simulations correlated to the NMR data in a time-averaged ensemble. The accurate solution structure ensemble of the metal-bound protein accurately describes the role of conformational sampling in allosteric regulation of DNA binding by zinc and serves to validate our previous unrestrained MD simulations of CzrA. This methodology has potentially broad applicability in the structure determination of metal ion bound proteins, protein folding and metal template protein-design studies.

A Candidate Ion-Retaining State in the Inward-Facing Conformation of Sodium/Galactose Symporter: Clues from Atomistic Simulations

The recent Vibrio parahaemolyticus sodium/galactose (vSGLT) symporter crystal structure captures the protein in an inward-facing substrate-bound conformation, with the sodium ion placed, by structural alignment, in a site equivalent to the Na2 site of the leucine transporter (LeuT). A recent study, based on molecular dynamics simulations, showed that the sodium ion spontaneously leaves its initial position diffusing outside vSGLT, toward the intracellular space. This suggested that the crystal structure corresponds to an ion-releasing state of the transporter. Here, using metadynamics, we identified a more stable Na(+) binding site corresponding to a putative ion-retaining state of the transporter. In addition, our simulations, consistently with mutagenesis studies, highlight the importance of D189 that, without being one of the Na(+)-coordinating residues, regulates its binding/release.