Fast evaluation of protein dynamics from deficient 15N relaxation data

AlphaFold and what is next: bridging functional, systems and structural biology

INTRODUCTION

The DeepMind's AlphaFold (AF) has revolutionized biomedical research by providing both experts and non-experts with an invaluable tool for predicting protein structures. However, while AF is highly effective for predicting structures of rigid and globular proteins, it is not able to fully capture the dynamics, conformational variability, and interactions of proteins with ligands and other biomacromolecules.

AREAS COVERED

In this review, we present a comprehensive overview of the latest advancements in 3D model predictions for biomacromolecules using AF. We also provide a detailed analysis its of strengths and limitations, and explore more recent iterations, modifications, and practical applications of this strategy. Moreover, we map the path forward for expanding the landscape of AF toward predicting structures of every protein and peptide in the proteome in the most physiologically relevant form. This discussion is based on an extensive literature search performed using PubMed and Google Scholar.

EXPERT OPINION

While significant progress has been made to enhance AF's modeling capabilities, we argue that a combined approach integrating both various in silico and in vitro methods will be most beneficial for the future of structural biology, bridging the gaps between static and dynamic features of proteins and their functions.

Structural insights into the DNA-binding mechanism of BCL11A: The integral role of ZnF6

The transcription factor BCL11A is a critical regulator of the switch from fetal hemoglobin (HbF: α2γ2) to adult hemoglobin (HbA: α2β2) during development. BCL11A binds at a cognate recognition site (TGACCA) in the γ-globin gene promoter and represses its expression. DNA-binding is mediated by a triple zinc finger domain, designated ZnF456. Here, we report comprehensive investigation of ZnF456, leveraging X-ray crystallography and NMR to determine the structures in both the presence and absence of DNA. We delve into the dynamics and mode of interaction with DNA. Moreover, we discovered that the last zinc finger of BCL11A (ZnF6) plays a different role compared to ZnF4 and 5, providing a positive entropic contribution to DNA binding and γ-globin gene repression. Comprehending the DNA binding mechanism of BCL11A opens avenues for the strategic, structure-based design of novel therapeutics targeting sickle cell disease and β-thalassemia.

Pitfalls in measurements of R1 relaxation rates of protein backbone 15N nuclei

The dynamics of the backbone and side-chains of protein are routinely studied by interpreting experimentally determined 15N spin relaxation rates. R1(15N), the longitudinal relaxation rate, reports on fast motions and encodes, together with the transverse relaxation R2, structural information about the shape of the molecule and the orientation of the amide bond vectors in the internal diffusion frame. Determining error-free 15N longitudinal relaxation rates remains a challenge for small, disordered, and medium-sized proteins. Here, we show that mono-exponential fitting is sufficient, with no statistical preference for bi-exponential fitting up to 800 MHz. A detailed comparison of the TROSY and HSQC techniques at medium and high fields showed no statistically significant differences. The least error-prone DD/CSA interference removal technique is the selective inversion of amide signals while avoiding water resonance. The exchange of amide with solvent deuterons appears to affect the rate R1 of solvent-exposed amides in all fields tested and in each DD/CSA interference removal technique in a statistically significant manner. In summary, the most accurate R1(15N) rates in proteins are achieved by selective amide inversion, without the addition of D2O. Importantly, at high magnetic fields stronger than 800 MHz, when non-mono-exponential decay is involved, it is advisable to consider elimination of the shortest delays (typically up to 0.32 s) or bi-exponential fitting.

Structural Insights into the DNA-Binding Mechanism of BCL11A: The Integral Role of ZnF6

The transcription factor BCL11A is a critical regulator of the switch from fetal hemoglobin (HbF: α 2 γ 2 ) to adult hemoglobin (HbA: α 2 β 2 ) during development. BCL11A binds at a cognate recognition site (TGACCA) in the γ-globin gene promoter and represses its expression. DNA-binding is mediated by a triple zinc finger domain, designated ZnF456. Here, we report comprehensive investigation of ZnF456, leveraging X-ray crystallography and NMR to determine the structures in both the presence and absence of DNA. We delve into the dynamics and mode of interaction with DNA. Moreover, we discovered that the last zinc finger of BCL11A (ZnF6) plays a special role in DNA binding and γ-globin gene repression. Our findings help account for some rare γ-globin gene promoter mutations that perturb BCL11A binding and lead to increased HbF in adults (hereditary persistence of fetal hemoglobin). Comprehending the DNA binding mechanism of BCL11A opens avenues for the strategic, structure-based design of novel therapeutics targeting sickle cell disease and β-thalassemia.

Increased slow dynamics defines ligandability of BTB domains

Efficient determination of protein ligandability, or the propensity to bind small-molecules, would greatly facilitate drug development for novel targets. Ligandability is currently assessed using computational methods that typically consider the static structural properties of putative binding sites or by experimental fragment screening. Here, we evaluate ligandability of conserved BTB domains from the cancer-relevant proteins LRF, KAISO, and MIZ1. Using fragment screening, we discover that MIZ1 binds multiple ligands. However, no ligands are uncovered for the structurally related KAISO or LRF. To understand the principles governing ligand-binding by BTB domains, we perform comprehensive NMR-based dynamics studies and find that only the MIZ1 BTB domain exhibits backbone µs-ms time scale motions. Interestingly, residues with elevated dynamics correspond to the binding site of fragment hits and recently defined HUWE1 interaction site. Our data argue that examining protein dynamics using NMR can contribute to identification of cryptic binding sites, and may support prediction of the ligandability of novel challenging targets.

NMR as a "Gold Standard" Method in Drug Design and Discovery

Studying disease models at the molecular level is vital for drug development in order to improve treatment and prevent a wide range of human pathologies. Microbial infections are still a major challenge because pathogens rapidly and continually evolve developing drug resistance. Cancer cells also change genetically, and current therapeutic techniques may be (or may become) ineffective in many cases. The pathology of many neurological diseases remains an enigma, and the exact etiology and underlying mechanisms are still largely unknown. Viral infections spread and develop much more quickly than does the corresponding research needed to prevent and combat these infections; the present and most relevant outbreak of SARS-CoV-2, which originated in Wuhan, China, illustrates the critical and immediate need to improve drug design and development techniques. Modern day drug discovery is a time-consuming, expensive process. Each new drug takes in excess of 10 years to develop and costs on average more than a billion US dollars. This demonstrates the need of a complete redesign or novel strategies. Nuclear Magnetic Resonance (NMR) has played a critical role in drug discovery ever since its introduction several decades ago. In just three decades, NMR has become a "gold standard" platform technology in medical and pharmacology studies. In this review, we present the major applications of NMR spectroscopy in medical drug discovery and development. The basic concepts, theories, and applications of the most commonly used NMR techniques are presented. We also summarize the advantages and limitations of the primary NMR methods in drug development.

A Conformational Switch in the Zinc Finger Protein Kaiso Mediates Differential Readout of Specific and Methylated DNA Sequences

Recognition of the epigenetic mark 5-methylcytosine (mC) at CpG sites in DNA has emerged as a novel function of many eukaryotic transcription factors (TFs). It remains unclear why the sequence specificity of these TFs differs for CpG-methylated motifs and consensus motifs. Here, we dissect the structural and dynamic basis for this differential DNA binding specificity in the human zinc finger TF Kaiso, which exhibits high affinity for two consecutive mCpG sites in variable contexts and also for a longer, sequence-specific Kaiso binding site (KBS). By integrating structural analysis and DNA binding studies with targeted protein mutagenesis and nucleotide substitutions, we identify distinct mechanisms for readout of methylated and KBS motifs by Kaiso. We show that a key glutamate residue (E535), critical for mCpG site recognition, adopts different conformations in complexes with specific and methylated DNA. These conformational differences, together with intrinsic variations in DNA flexibility and/or solvation at TpG versus mCpG sites, contribute to the different DNA affinity and sequence specificity. With methylated DNA, multiple direct contacts between E535 and the 5' mCpG site dominate the binding affinity, allowing for tolerance of different flanking DNA sequences. With KBS, Kaiso employs E535 as part of an indirect screen of the 5' flanking sequence, relying on key tyrosine-DNA interactions to stabilize an optimal DNA conformation and select against noncognate sites. These findings demonstrate how TFs use conformational adaptation and exploit variations in DNA flexibility to achieve distinct DNA readout outcomes and target a greater variety of regulatory and epigenetic sites than previously appreciated.

Characterization of Internal Protein Dynamics and Conformational Entropy by NMR Relaxation

Recent studies suggest that the fast timescale motion of methyl-bearing side chains may play an important role in mediating protein activity. These motions have been shown to encapsulate the residual conformational entropy of the folded state that can potentially contribute to the energetics of protein function. Here, we provide an overview of how to characterize these motions using nuclear magnetic resonance (NMR) spin relaxation methods. The strengths and limitations of several techniques are highlighted in order to assist with experimental design. Particular emphasis is placed on the practical aspects of sample preparation, data collection, data fitting, and statistical analysis. Additionally, discussion of the recently refined "entropy meter" is presented and its use in converting NMR observables to conformational entropy is illustrated. Taken together, these methods should yield new insights into the complex interplay between structure and dynamics in protein function.

Simple and Robust Study of Backbone Dynamics of Crystalline Proteins Employing 1H-15N Dipolar Coupling Dispersion

We report a new solid-state multidimensional NMR approach based on the cross-polarization with variable-contact pulse sequence [ Paluch , P. ; Pawlak , T. ; Amoureux , J.-P. ; Potrzebowski , M. J. J. Magn. Reson. 233 , 2013 , 56 ], with ¹H inverse detection and very fast magic angle spinning (ν_R = 60 kHz), dedicated to the measurement of local molecular motions of ¹H-¹⁵N vectors. The introduced three-dimensional experiments, ¹H-¹⁵N-¹H and hCA(N)H, are particularly useful for the study of molecular dynamics of proteins and other complex structures. The applicability and power of this methodology have been revealed by employing as a model sample the GB-1 small protein doped with Na₂CuEDTA. The results clearly prove that the dispersion of ¹H-¹⁵N dipolar coupling constants well correlates with higher order structure of the protein. Our approach complements the conventional studies and offers a fast and reasonably simple method.