Generative artificial intelligence for de novo protein design.

Engineering new molecules with desirable functions and properties has the potential to extend our ability to engineer proteins beyond what nature has so far evolved. Advances in the so-called 'de novo' design problem have recently been brought forward by developments in artificial intelligence. Generative architectures, such as language models and diffusion processes, seem adept at generating novel, yet realistic proteins that display desirable properties and perform specified functions. State-of-the-art design protocols now achieve experimental success rates nearing 20%, thus widening the access to de novo designed proteins. Despite extensive progress, there are clear field-wide challenges, for example, in determining the best in silico metrics to prioritise designs for experimental testing, and in designing proteins that can undergo large conformational changes or be regulated by post-translational modifications. With an increase in the number of models being developed, this review provides a framework to understand how these tools fit into the overall process of de novo protein design. Throughout, we highlight the power of incorporating biochemical knowledge to improve performance and interpretability.

Discovery and pharmacophoric characterization of chemokine network inhibitors using phage-display, saturation mutagenesis and computational modelling.

CC and CXC-chemokines are the primary drivers of chemotaxis in inflammation, but chemokine network redundancy thwarts pharmacological intervention. Tick evasins promiscuously bind CC and CXC-chemokines, overcoming redundancy. Here we show that short peptides that promiscuously bind both chemokine classes can be identified from evasins by phage-display screening performed with multiple chemokines in parallel. We identify two conserved motifs within these peptides and show using saturation-mutagenesis phage-display and chemotaxis studies of an exemplar peptide that an anionic patch in the first motif and hydrophobic, aromatic and cysteine residues in the second are functionally necessary. AlphaFold2-Multimer modelling suggests that the peptide occludes distinct receptor-binding regions in CC and in CXC-chemokines, with the first and second motifs contributing ionic and hydrophobic interactions respectively. Our results indicate that peptides with broad-spectrum anti-chemokine activity and therapeutic potential may be identified from evasins, and the pharmacophore characterised by phage display, saturation mutagenesis and computational modelling.

Membranome 3.0: Database of single-pass membrane proteins with AlphaFold models.

The Membranome database provides comprehensive structural information on single-pass (i.e., bitopic) membrane proteins from six evolutionarily distant organisms, including protein-protein interactions, complexes, mutations, experimental structures, and models of transmembrane α-helical dimers. We present a new version of this database, Membranome 3.0, which was significantly updated by revising the set of 5,758 bitopic proteins and incorporating models generated by AlphaFold 2 in the database. The AlphaFold models were parsed into structural domains located at the different membrane sides, modified to exclude low-confidence unstructured terminal regions and signal sequences, validated through comparison with available experimental structures, and positioned with respect to membrane boundaries. Membranome 3.0 was re-developed to facilitate visualization and comparative analysis of multiple 3D structures of proteins that belong to a specified family, complex, biological pathway, or membrane type. New tools for advanced search and analysis of proteins, their interactions, complexes, and mutations were included. The database is freely accessible at https://membranome.org.

Current structure predictors are not learning the physics of protein folding.

SUMMARY: Motivation. Predicting the native state of a protein has long been considered a gateway problem for understanding protein folding. Recent advances in structural modeling driven by deep learning have achieved unprecedented success at predicting a protein's crystal structure, but it is not clear if these models are learning the physics of how proteins dynamically fold into their equilibrium structure or are just accurate knowledge-based predictors of the final state. Results. In this work, we compare the pathways generated by state-of-the-art protein structure prediction methods to experimental data about protein folding pathways. The methods considered were AlphaFold 2, RoseTTAFold, trRosetta, RaptorX, DMPfold, EVfold, SAINT2 and Rosetta. We find evidence that their simulated dynamics capture some information about the folding pathway, but their predictive ability is worse than a trivial classifier using sequence-agnostic features like chain length. The folding trajectories produced are also uncorrelated with experimental observables such as intermediate structures and the folding rate constant. These results suggest that recent advances in structure prediction do not yet provide an enhanced understanding of protein folding. Availability. The data underlying this article are available in GitHub at https://github.com/oxpig/structure-vs-folding/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Deconvolution of conformational exchange from Raman spectra of aqueous RNA nucleosides.

Ribonucleic acids (RNAs) are key to the central dogma of molecular biology. While Raman spectroscopy holds great potential for studying RNA conformational dynamics, current computational Raman prediction and assignment methods are limited in terms of system size and inclusion of conformational exchange. Here, a framework is presented that predicts Raman spectra using mixtures of sub-spectra corresponding to major conformers calculated using classical and ab initio molecular dynamics. Experimental optimization allowed purines and pyrimidines to be characterized as predominantly syn and anti, respectively, and ribose into exchange between equivalent south and north populations. These measurements are in excellent agreement with Raman spectroscopy of ribonucleosides, and previous experimental and computational results. This framework provides a measure of ribonucleoside solution populations and conformational exchange in RNA subunits. It complements other experimental techniques and could be extended to other molecules, such as proteins and carbohydrates, enabling biological insights and providing a new analytical tool.

Revitalizing the concept of bond order through delocalization measures in real space.

Ab initio quantum chemistry is an independent source of information supplying an ever widening group of experimental chemists. However, bridging the gap between these ab initio data and chemical insight remains a challenge. In particular, there is a need for a bond order index that characterizes novel bonding patterns in a reliable manner, while recovering the familiar effects occurring in well-known bonds. In this article, through a large body of calculations, we show how the delocalization index derived from Quantum Chemical Topology (QCT) serves as such a bond order. This index is defined in a parameter-free, intuitive and consistent manner, and with little qualitative dependency on the level of theory used. The delocalization index is also able to detect the subtler bonding effects that underpin most practical organic and inorganic chemistry. We explore and connect the properties of this index and open the door for its extensive usage in the understanding and discovery of novel chemistry.