N4BP1 functions as a dimerization-dependent linear ubiquitin reader which regulates TNF signalling

Signalling through TNFR1 modulates proinflammatory gene transcription and programmed cell death, and its impairment causes autoimmune diseases and cancer. NEDD4-binding protein 1 (N4BP1) is a critical suppressor of proinflammatory cytokine production that acts as a regulator of innate immune signalling and inflammation. However, our current understanding about the molecular properties that enable N4BP1 to exert its suppressive potential remain limited. Here, we show that N4BP1 is a novel linear ubiquitin reader that negatively regulates NFκB signalling by its unique dimerization-dependent ubiquitin-binding module that we named LUBIN. Dimeric N4BP1 strategically positions two non-selective ubiquitin-binding domains to ensure preferential recognition of linear ubiquitin. Under proinflammatory conditions, N4BP1 is recruited to the nascent TNFR1 signalling complex, where it regulates duration of proinflammatory signalling in LUBIN-dependent manner. N4BP1 deficiency accelerates TNFα-induced cell death by increasing complex II assembly. Under proapoptotic conditions, caspase-8 mediates proteolytic processing of N4BP1, resulting in rapid degradation of N4BP1 by the 26 S proteasome, and acceleration of apoptosis. In summary, our findings demonstrate that N4BP1 dimerization creates a novel type of ubiquitin reader that selectively recognises linear ubiquitin which enables the timely and coordinated regulation of TNFR1-mediated inflammation and cell death.

The Legionella collagen-like protein employs a unique binding mechanism for the recognition of host glycosaminoglycans

Bacterial adhesion is a fundamental process which enables colonisation of niche environments and is key for infection. However, in Legionella pneumophila, the causative agent of Legionnaires' disease, these processes are not well understood. The Legionella collagen-like protein (Lcl) is an extracellular peripheral membrane protein that recognises sulphated glycosaminoglycans (GAGs) on the surface of eukaryotic cells, but also stimulates bacterial aggregation in response to divalent cations. Here we report the crystal structure of the Lcl C-terminal domain (Lcl-CTD) and present a model for intact Lcl. Our data reveal that Lcl-CTD forms an unusual dynamic trimer arrangement with a positively charged external surface and a negatively charged solvent exposed internal cavity. Through Molecular Dynamics (MD) simulations, we show how the GAG chondroitin-4-sulphate associates with the Lcl-CTD surface via unique binding modes. Our findings show that Lcl homologs are present across both the Pseudomonadota and Fibrobacterota-Chlorobiota-Bacteroidota phyla and suggest that Lcl may represent a versatile carbohydrate binding mechanism.

Abnormal myosin post-translational modifications and ATP turnover time associated with human congenital myopathy-related RYR1 mutations

AIM

Conditions related to mutations in the gene encoding the skeletal muscle ryanodine receptor 1 (RYR1) are genetic muscle disorders and include congenital myopathies with permanent weakness, as well as episodic phenotypes such as rhabdomyolysis/myalgia. Although RYR1 dysfunction is the primary mechanism in RYR1-related disorders, other downstream pathogenic events are less well understood and may include a secondary remodeling of major contractile proteins. Hence, in the present study, we aimed to investigate whether congenital myopathy-related RYR1 mutations alter the regulation of the most abundant contractile protein, myosin.

METHODS

We used skeletal muscle tissues from five patients with RYR1-related congenital myopathy and compared those with five controls and five patients with RYR1-related rhabdomyolysis/myalgia. We then defined post-translational modifications on myosin heavy chains (MyHCs) using LC/MS. In parallel, we determined myosin relaxed states using Mant-ATP chase experiments and performed molecular dynamics (MD) simulations.

RESULTS

LC/MS revealed two additional phosphorylations (Thr1309-P and Ser1362-P) and one acetylation (Lys1410-Ac) on the β/slow MyHC of patients with congenital myopathy. This method also identified six acetylations that were lacking on MyHC type IIa of these patients (Lys35-Ac, Lys663-Ac, Lys763-Ac, Lys1171-Ac, Lys1360-Ac, and Lys1733-Ac). MD simulations suggest that modifying myosin Ser1362 impacts the protein structure and dynamics. Finally, Mant-ATP chase experiments showed a faster ATP turnover time of myosin heads in the disordered-relaxed conformation.

CONCLUSIONS

Altogether, our results suggest that RYR1 mutations have secondary negative consequences on myosin structure and function, likely contributing to the congenital myopathic phenotype.

Binding pocket dynamics along the recovery stroke of human β-cardiac myosin

The druggability of small-molecule binding sites can be significantly affected by protein motions and conformational changes. Ligand binding, protein dynamics and protein function have been shown to be closely interconnected in myosins. The breakthrough discovery of omecamtiv mecarbil (OM) has led to an increased interest in small molecules that can target myosin and modulate its function for therapeutic purposes (myosin modulators). In this work, we use a combination of computational methods, including steered molecular dynamics, umbrella sampling and binding pocket tracking tools, to follow the evolution of the OM binding site during the recovery stroke transition of human β-cardiac myosin. We found that steering two internal coordinates of the motor domain can recapture the main features of the transition and in particular the rearrangements of the binding site, which shows significant changes in size, shape and composition. Possible intermediate conformations were also identified, in remarkable agreement with experimental findings. The differences in the binding site properties observed along the transition can be exploited for the future development of conformation-selective myosin modulators.

Comparative study of binding pocket structure and dynamics in cardiac and skeletal myosin

The development of small molecule myosin modulators has seen an increased effort in recent years due to their possible use in the treatment of cardiac and skeletal myopathies. Omecamtiv mecarbil (OM) is the first-in-class cardiac myotrope and the first to enter clinical trials. Its selectivity toward slow/beta-cardiac myosin lies at the heart of its function; however, little is known about the underlying reasons for selectivity to this isoform as opposed to other closely related ones such as fast-type skeletal myosins. In this work, we compared the structure and dynamics of the OM binding site in cardiac and in fasttype IIa skeletal myosin to identify possible reasons for OM selectivity. We found that the different shape, size, and composition of the binding pocket in skeletal myosin directly affects the binding mode and related affinity of OM, which is potentially a result of weaker interactions and less optimal molecular recognition. Moreover, we identified a side pocket adjacent to the OM binding site that shows increased accessibility in skeletal myosin compared with the cardiac isoform. These findings could pave the way to the development of skeletal-selective compounds that can target this region of the protein and potentially be used to treat congenital myopathies where muscle weakness is related to myosin loss of function.

Reconstruction of ARNT PAS-B Unfolding Pathways by Steered Molecular Dynamics and Artificial Neural Networks

Several experimental studies indicated that large conformational changes, including partial domain unfolding, have a role in the functional mechanisms of the basic helix loop helix Per/ARNT/SIM (bHLH-PAS) transcription factors. Recently, single-molecule atomic force microscopy (AFM) revealed two distinct pathways for the mechanical unfolding of the ARNT PAS-B. In this work we used steered molecular dynamics simulations to gain new insights into this process at an atomistic level. To reconstruct and classify pathways sampled in multiple simulations, we designed an original approach based on the use of self-organizing maps (SOMs). This led us to identify two types of unfolding pathways for the ARNT PAS-B, which are in good agreement with the AFM findings. Analysis of average forces mapped on the SOM revealed a stable conformation of the PAS-B along one pathway, which represents a possible structural model for the intermediate state detected by AFM. The approach here proposed will facilitate the study of other signal transmission mechanisms involving the folding/unfolding of PAS domains.

Heart Failure Drug Modifies the Intrinsic Dynamics of the Pre-Power Stroke State of Cardiac Myosin

Omecamtiv mecarbil (OM), currently investigated for the treatment of heart failure, is the first example of a new class of drugs (cardiac myotropes) that can modify muscle contractility by directly targeting sarcomeric proteins. Using atomistic molecular dynamics simulations, we show that the binding of OM to the pre-power stroke state of cardiac myosin inhibits the functional motions of the protein and potentially affects Pi release from the nucleotide binding site. We also show that the changes in myosin ATPase activity induced by a set of OM analogues can be predicted from their relative affinity to the pre-power stroke state compared to the near rigor one, indicating that conformational selectivity plays an important role in determining the activity of these compounds.

Structure and functional analysis of the Legionella pneumophila chitinase ChiA reveals a novel mechanism of metal-dependent mucin degradation

Chitinases are important enzymes that contribute to the generation of carbon and nitrogen from chitin, a long chain polymer of N-acetylglucosamine that is abundant in insects, fungi, invertebrates and fish. Although mammals do not produce chitin, chitinases have been identified in bacteria that are key virulence factors in severe respiratory, gastrointestinal and urinary diseases. However, it is unclear how these enzymes are able to carry out this dual function. Legionella pneumophila is the causative agent of Legionnaires' disease, an often-fatal pneumonia and its chitinase ChiA is essential for the survival of L. pneumophila in the lung. Here we report the first atomic resolution insight into the pathogenic mechanism of a bacterial chitinase. We derive an experimental model of intact ChiA and show how its N-terminal region targets ChiA to the bacterial surface after its secretion. We provide the first evidence that L. pneumophila can bind mucins on its surface, but this is not dependent on ChiA. This demonstrates that additional peripheral mucin binding proteins are also expressed in L. pneumophila. We also show that the ChiA C-terminal chitinase domain has novel Zn2+-dependent peptidase activity against mammalian mucin-like proteins, namely MUC5AC and the C1-esterase inhibitor, and that ChiA promotes bacterial penetration of mucin gels. Our findings suggest that ChiA can facilitate passage of L. pneumophila through the alveolar mucosa, can modulate the host complement system and that ChiA may be a promising target for vaccine development.

Binding Pockets in Proteins Induced by Mechanical Stress

We report the first observation of a pocket that opens as a result of a mechanical force applied to an Ig-like domain from the cardiac muscle. This previously unseen mechanism of pocket formation is revealed by molecular dynamics simulations under force. Preliminary investigations show that this "mechano-pocket" is potentially druggable and could be found in other domains from the same fold family, suggesting the existence of a general mechanism of pocket formation under mechanical stress.

In silico identification of rescue sites by double force scanning

Motivation

A deleterious amino acid change in a protein can be compensated by a second-site rescue mutation. These compensatory mechanisms can be mimicked by drugs. In particular, the location of rescue mutations can be used to identify protein regions that can be targeted by small molecules to reactivate a damaged mutant.

Results

We present the first general computational method to detect rescue sites. By mimicking the effect of mutations through the application of forces, the double force scanning (DFS) method identifies the second-site residues that make the protein structure most resilient to the effect of pathogenic mutations. We tested DFS predictions against two datasets containing experimentally validated and putative evolutionary-related rescue sites. A remarkably good agreement was found between predictions and experimental data. Indeed, almost half of the rescue sites in p53 was correctly predicted by DFS, with 65% of remaining sites in contact with DFS predictions. Similar results were found for other proteins in the evolutionary dataset.

Availability and implementation

The DFS code is available under GPL at https://fornililab.github.io/dfs/.

Supplementary information

Supplementary data are available at Bioinformatics online.

Allosteric modulation of cardiac myosin dynamics by omecamtiv mecarbil

New promising avenues for the pharmacological treatment of skeletal and heart muscle diseases rely on direct sarcomeric modulators, which are molecules that can directly bind to sarcomeric proteins and either inhibit or enhance their activity. A recent breakthrough has been the discovery of the myosin activator omecamtiv mecarbil (OM), which has been shown to increase the power output of the cardiac muscle and is currently in clinical trials for the treatment of heart failure. While the overall effect of OM on the mechano-chemical cycle of myosin is to increase the fraction of myosin molecules in the sarcomere that are strongly bound to actin, the molecular basis of its action is still not completely clear. We present here a Molecular Dynamics study of the motor domain of human cardiac myosin bound to OM, where the effects of the drug on the dynamical properties of the protein are investigated for the first time with atomistic resolution. We found that OM has a double effect on myosin dynamics, inducing a) an increased coupling of the motions of the converter and lever arm subdomains to the rest of the protein and b) a rewiring of the network of dynamic correlations, which produces preferential communication pathways between the OM binding site and distant functional regions. The location of the residues responsible for these effects suggests possible strategies for the future development of improved drugs and the targeting of specific cardiomyopathy-related mutations.

Using Local States To Drive the Sampling of Global Conformations in Proteins

Conformational changes associated with protein function often occur beyond the time scale currently accessible to unbiased molecular dynamics (MD) simulations, so that different approaches have been developed to accelerate their sampling. Here we investigate how the knowledge of backbone conformations preferentially adopted by protein fragments, as contained in precalculated libraries known as structural alphabets (SA), can be used to explore the landscape of protein conformations in MD simulations. We find that (a) enhancing the sampling of native local states in both metadynamics and steered MD simulations allows the recovery of global folded states in small proteins; (b) folded states can still be recovered when the amount of information on the native local states is reduced by using a low-resolution version of the SA, where states are clustered into macrostates; and (c) sequences of SA states derived from collections of structural motifs can be used to sample alternative conformations of preselected protein regions. The present findings have potential impact on several applications, ranging from protein model refinement to protein folding and design.

Towards the identification of the allosteric Phe-binding site in phenylalanine hydroxylase

The enzyme phenylalanine hydroxylase (PAH) is defective in the inherited disorder phenylketonuria. PAH, a tetrameric enzyme, is highly regulated and displays positive cooperativity for its substrate, Phe. Whether Phe binds to an allosteric site is a matter of debate, despite several studies worldwide. To address this issue, we generated a dimeric model for Phe-PAH interactions, by performing molecular docking combined with molecular dynamics simulations on human and rat wild-type sequences and also on a human G46S mutant. Our results suggest that the allosteric Phe-binding site lies at the dimeric interface between the regulatory and the catalytic domains of two adjacent subunits. The structural and dynamical features of the site were characterized in depth and described. Interestingly, our findings provide evidence for lower allosteric Phe-binding ability of the G46S mutant than the human wild-type enzyme. This also explains the disease-causing nature of this mutant.

Anatomy of protein disorder, flexibility and disease-related mutations

Integration of protein structural information with human genetic variation and pathogenic mutations is essential to understand molecular mechanisms associated with the effects of polymorphisms on protein interactions and cellular processes. We investigate occurrences of non-synonymous SNPs in ordered and disordered protein regions by systematic mapping of common variants and disease-related SNPs onto these regions. We show that common variants accumulate in disordered regions; conversely pathogenic variants are significantly depleted in disordered regions. These different occurrences of pathogenic and common SNPs can be attributed to a negative selection on random mutations in structurally highly constrained regions. New approaches in the study of quantitative effects of pathogenic-related mutations should effectively account for all the possible contexts and relative functional constraints in which the sequence variation occurs.

Protein-protein interaction networks studies and importance of 3D structure knowledge

Protein-protein interaction networks (PPINs) are a powerful tool to study biological processes in living cells. In this review, we present the progress of PPIN studies from abstract to more detailed representations. We will focus on 3D interactome networks, which offer detailed information at the atomic level. This information can be exploited in understanding not only the underlying cellular mechanisms, but also how human variants and disease-causing mutations affect protein functions and complexes' stability. Recent studies have used structural information on PPINs to also understand the molecular mechanisms of binding partner selection. We will address the challenges in generating 3D PPINs due to the restricted number of solved protein structures. Finally, some of the current use of 3D PPINs will be discussed, highlighting their contribution to the studies in genotype-phenotype relationships and in the optimization of targeted studies to design novel chemical compounds for medical treatments.

Specialized Dynamical Properties of Promiscuous Residues Revealed by Simulated Conformational Ensembles

The ability to interact with different partners is one of the most important features in proteins. Proteins that bind a large number of partners (hubs) have been often associated with intrinsic disorder. However, many examples exist of hubs with an ordered structure, and evidence of a general mechanism promoting promiscuity in ordered proteins is still elusive. An intriguing hypothesis is that promiscuous binding sites have specific dynamical properties, distinct from the rest of the interface and pre-existing in the protein isolated state. Here, we present the first comprehensive study of the intrinsic dynamics of promiscuous residues in a large protein data set. Different computational methods, from coarse-grained elastic models to geometry-based sampling methods and to full-atom Molecular Dynamics simulations, were used to generate conformational ensembles for the isolated proteins. The flexibility and dynamic correlations of interface residues with a different degree of binding promiscuity were calculated and compared considering side chain and backbone motions, the latter both on a local and on a global scale. The study revealed that (a) promiscuous residues tend to be more flexible than nonpromiscuous ones, (b) this additional flexibility has a higher degree of organization, and (c) evolutionary conservation and binding promiscuity have opposite effects on intrinsic dynamics. Findings on simulated ensembles were also validated on ensembles of experimental structures extracted from the Protein Data Bank (PDB). Additionally, the low occurrence of single nucleotide polymorphisms observed for promiscuous residues indicated a tendency to preserve binding diversity at these positions. A case study on two ubiquitin-like proteins exemplifies how binding promiscuity in evolutionary related proteins can be modulated by the fine-tuning of the interface dynamics. The interplay between promiscuity and flexibility highlighted here can inspire new directions in protein–protein interaction prediction and design methods.

GSATools: analysis of allosteric communication and functional local motions using a structural alphabet

Motivation: GSATools is a free software package to analyze conformational ensembles and to detect functional motions in proteins by means of a structural alphabet. The software integrates with the widely used GROMACS simulation package and can generate a range of graphical outputs. Three applications can be supported: (i) investigation of the conformational variability of local structures; (ii) detection of allosteric communication; and (iii) identification of local regions that are critical for global functional motions. These analyses provide insights into the dynamics of proteins and allow for targeted design of functional mutants in theoretical and experimental studies.

Availability: The C source code of the GSATools, along with a set of pre-compiled binaries, is freely available under GNU General Public License from http://mathbio.nimr.mrc.ac.uk/wiki/GSATools.

Contact:alessandro.pandini@kcl.ac.uk or jkleinj@nimr.mrc.ac.uk

Supplementary information:Supplementary data are available at Bioinformatics online.

Decrypting Prion Protein Conversion into a β-Rich Conformer by Molecular Dynamics

Prion diseases are fatal neurodegenerative diseases characterized by the formation of β-rich oligomers and the accumulation of amyloid fibrillar deposits in the central nervous system. Understanding the conversion of the cellular prion protein into its β-rich polymeric conformers is fundamental to tackling the early stages of the development of prion diseases. In this paper, we have identified unfolding and refolding steps critical to the conversion into a β-rich conformer for different constructs of the ovine prion protein by molecular dynamics simulations. By combining our results with in vitro experiments, we show that the folded C-terminus of the ovine prion protein is able to recurrently undergo a drastic conformational change by displacement of the H1 helix, uncovering of the H2H3 domain, and formation of persistent β-sheets between H2 and H3 residues. The observed β-sheets refold toward the C-terminus exposing what we call a "bending region" comprising residues 204-214. This is strikingly coincident with the region harboring mutations determining the fate of the prion oligomerization process. The β-rich intermediate is used here for the construction of a putative model for the assembly into an oligomeric aggregate. The results presented here confirm the importance of the H2H3 domain for prion oligomer formation and therefore its potential use as molecular target in the design of novel prion inhibitors.

Detection of allosteric signal transmission by information-theoretic analysis of protein dynamics

Allostery offers a highly specific way to modulate protein function. Therefore, understanding this mechanism is of increasing interest for protein science and drug discovery. However, allosteric signal transmission is difficult to detect experimentally and to model because it is often mediated by local structural changes propagating along multiple pathways. To address this, we developed a method to identify communication pathways by an information-theoretical analysis of molecular dynamics simulations. Signal propagation was described as information exchange through a network of correlated local motions, modeled as transitions between canonical states of protein fragments. The method was used to describe allostery in two-component regulatory systems. In particular, the transmission from the allosteric site to the signaling surface of the receiver domain NtrC was shown to be mediated by a layer of hub residues. The location of hubs preferentially connected to the allosteric site was found in close agreement with key residues experimentally identified as involved in the signal transmission. The comparison with the networks of the homologues CheY and FixJ highlighted similarities in their dynamics. In particular, we showed that a preorganized network of fragment connections between the allosteric and functional sites exists already in the inactive state of all three proteins.

In silico phosphorylation of the autoinhibited form of p47(phox): insights into the mechanism of activation

Activation of the multicomponent enzyme NADPH oxidase requires the interaction between the tandem SH3 domain of the cytosolic subunit p47(phox) and the cytoplasmic tail of membrane-bound p22(phox). In the resting state, p47(phox) exists in an autoinhibited conformation stabilized by intramolecular contacts between the SH3 domains and an adjacent polybasic region. Phosphorylation of three serine residues, Ser(303), Ser(304), and Ser(328) within this polybasic region has been shown to be sufficient for the disruption of the intramolecular interactions thereby inducing an active state of p47(phox). This active conformation is accessible to the cytoplasmic tail of p22(phox) and initiates the formation of the membrane-bound functional enzyme complex. Molecular dynamics simulations reveal insights in the mechanism of activation of the autoinhibited form of p47(phox) by in silico phosphorylation, of the three serine residues, Ser(303), Ser(304), and Ser(328). The simulations highlight the major collective coordinates generating the opening and the closing of the two SH3 domains and the residues that cause the unmasking of the p22(phox) binding site.

Structural alphabets derived from attractors in conformational space

Background

The hierarchical and partially redundant nature of protein structures justifies the definition of frequently occurring conformations of short fragments as 'states'. Collections of selected representatives for these states define Structural Alphabets, describing the most typical local conformations within protein structures. These alphabets form a bridge between the string-oriented methods of sequence analysis and the coordinate-oriented methods of protein structure analysis.

Results

A Structural Alphabet has been derived by clustering all four-residue fragments of a high-resolution subset of the protein data bank and extracting the high-density states as representative conformational states. Each fragment is uniquely defined by a set of three independent angles corresponding to its degrees of freedom, capturing in simple and intuitive terms the properties of the conformational space. The fragments of the Structural Alphabet are equivalent to the conformational attractors and therefore yield a most informative encoding of proteins. Proteins can be reconstructed within the experimental uncertainty in structure determination and ensembles of structures can be encoded with accuracy and robustness.

Conclusions

The density-based Structural Alphabet provides a novel tool to describe local conformations and it is specifically suitable for application in studies of protein dynamics.

Accurate description of nitrogenous base flexibility in classical molecular dynamics simulations of nucleotides bound to proteins

Proteins can induce significant distortions in planar cyclic compounds upon binding, in particular in nucleotide-enzyme complexes. An accurate representation of the ring flexibility is thus desirable when modeling these systems through classical force fields, especially when deformations are supposed to be involved in the catalytic mechanism. In this study, we use a newly developed general procedure to determine sets of dihedral parameters for planar cycles that accurately reproduce their out-of-plane normal modes as determined at the quantum mechanical (QM) level. The optimization allows the deviation from reference data to be reduced for the pyrimidine bases to values comparable to the accuracy of the QM data. Furthermore, the influence of the description of ring flexibility in protein-ligand interactions is assessed through molecular dynamics simulations of the complex between uridine and the pyrimidine-specific nucleoside hydrolase YeiK using the AMBER force field. The differences in ligand-protein interactions emerging from different parameter sets are also discussed with respect to existing biochemical data.

On the suitability of strictly localized orbitals for hybrid QM/MM calculations

In the QM/MM method we have developed (LSCF/MM), the QM and the MM parts are held together by means of strictly localized bonding orbitals (SLBOs). Generally these SLBOs are derived from localized bond orbitals (LBOs) that undergo tails deletion, resulting in a nonpredictable change of their properties. An alternative set of SLBOs is provided by the extremely localized molecular orbitals (ELMOs) approach, where the orbitals are rigorously localized on some prefixed atoms without tails on the other atoms of the molecule. A comparative study of SLBOs arising from various localization schemes and ELMOs is presented to test the reliability and the transferability of these functions within the Local Self-Consistent Field (LSCF) framework. Two types of chemical bonds were considered: C--C and C--O single bonds. The localized functions are obtained on the ethane and the methanol molecules, and are tested on beta-alanine and diethyl ether molecules. Moreover, the various protonation forms of beta-alanine have been investigated to illustrate how well the polarity variation of the chemical bond can be handled throughout a chemical process. At last, rotation energy profiles around C--C and C--O bonds are reproduced for butane and fluoromethanol. Energetic, geometric, as well as electronic factors all indicate that ELMO functions are much more transferable from one molecule to another, leading to results closer to the usual SCF reference than any other calculations involving any other localized orbitals. When the shape of the orbital is the most important factor then ELMO functions will perform as well as any other localized orbital.

Optimal virtual orbitals to relax wave functions built up with transferred extremely localized molecular orbitals

Extremely localized molecular orbitals (ELMOs), namely orbitals strictly localized on molecular fragments, are easily transferable from one molecule to another one. Hence, they provide a natural way to set up the electronic structure of large molecules using a data base of orbitals obtained from model molecules. However, this procedure obviously increases the energy with respect to a traditional MO calculation. To gain accuracy, it is important to introduce a partial electron delocalization. This can be carried out by defining proper optimal virtual orbitals that supply an efficient set for nonorthogonal configurations to be employed in VB-like expansions.