Applied Mathematics
Vol.09 No.10(2018), Article ID:87974,40 pages
10.4236/am.2018.910076

Global Geometrical Constraints on the Shape of Proteins and Their Influence on Allosteric Regulation

Naoto Morikawa

Genocript, Zama, Japan

Copyright © 2018 by author and Scientific Research Publishing Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).

http://creativecommons.org/licenses/by/4.0/

Received: September 14, 2018; Accepted: October 21, 2018; Published: October 24, 2018

ABSTRACT

Proteins are the workhorse molecules of the cell, which are obtained by folding long chains of amino acids. Since not all shapes are obtained as a folded chain of amino acids, there should be global geometrical constraints on the shape. Moreover, since the function of a protein is largely determined by its shape, constraints on the shape should have some influence on its interaction with other proteins. In this paper, we consider global geometrical constraints on the shape of proteins. Using a mathematical toy model, in which proteins are represented as closed chains of tetrahedrons, we have identified not only global geometrical constraints on the shape of proteins, but also their influence on protein interactions. As an example, we show that a garlic-bulb like structure appears as a result of the constraints. Regarding the influence of global geometrical constraints on interactions, we consider their influence on the structural coupling of two distal sites in allosteric regulation. We then show the inseparable relationship between global geometrical constraints and protein interactions; i.e. they are different sides of the same coin. This finding could be important for the understanding of the basic mechanisms of allosteric regulation of protein functions.

Keywords:

Differential Geometry, Protein Structure, Discrete Mathematics, Allosteric Regulation, Simultaneous Equations for Shape

1. Introduction

In this paper, we consider global geometrical constraints on the shape of proteins, using the mathematical toy model of proteins proposed in [1] . Proteins are the workhorse molecules of the cell, which are obtained as a complex of folded chains of amino acids. Since the function of proteins depends primarily on their shape, structural studies are essential for understanding proteins. In our approach, protein molecules are represented as a complex of closed trajectories of tetrahedrons. Then, the surface of proteins is obtained as the intersection of a pair of four-dimensional cones [2] . Interactions between proteins are defined (or mimicked) as “fusion and fission” of closed trajectories.

Previously, two types of geometrical constraints are known in the study of protein structures. One is a set of constraints on the backbone conformation due to collisions between atoms [3] . The backbone conformation is determined by torsion angle pairs ( ϕ , ψ ) along the backbone, and their allowed values are shown in the Ramachandran map [4] . The other is a set of constraints on relative distances between certain pairs of atoms, which are obtained from either physical experiments or theoretical estimates. The determination of protein structures which satisfy a set of constraints on inter-atomic distances, known as the distance geometry problem, is an important problem in structural biology [5] .

In virology, another type of geometrical constraints, the symmetry of the virus structure, is also considered. Viruses are metastable macromolecular assemblies composed of the viral genome enclosed within the protein shells, called viral capsids [6] . Virus capsids are highly specific assemblies that are formed from a large number of often identical subunits. Formulated in [7] is a set of structural constraints on the subunit arrangements, using an extension of the underlying symmetry group. On the other hand, [8] finds that some viruses allow their representation as two-dimensional monohedral tilings of a bound surface, where each tile represents a subunit. Note that viral molecules consist of separeated parts. Protein molecules are obtained by folding a chain of linked parts and it is impossible to describe the shape of proteins by symmetry alone nor to describe their surface by tiling of basic subunits.

What we will consider below are global constraints on the shape of a complex of folded chains of basic blocks, such as triangles and tetrahedrons. One of the advantages of our model is the correspondence between “the shape of molecules” and “interaction between molecules”. Since a protein’s function is largely determined by its shape, constraints on the shape of a protein should have some influence on its interaction with other proteins. In our model, the geometrical constraints on the shape of a molecule correspond to the constraints on the interaction between three molecules, such as allosteric regulations. In the section before the conclusion, we will explain the correspondence between geometrical constraints on the shape and allosteric regulations using an example. An introduction to allosteric regulation is also given there.

Finally, Genocript (http://www.genocript.com) is the one-man bio-venture started by Naoto Morikawa in 2000 which is developing software tools for protein structure analysis.

2. Discrete Differential Geometry of Triangles

Now, let us consider the case of closed trajectories of triangles to explain the basic ideas behind our approach. For detailed description, see [1] and [2] .

In the following, the coordinates of points in the N-dimensional Eucledean space E N ( N = 3 or 4) are represented by a monomial in N indeterminates x 0 , x 1 , , x N 1 for space saving purposes. For example, point ( l , m , n ) 3 E 3 is represented by x 0 l x 1 m x 2 n , where denotes the set of all integers. ( 0,0,0 ) , ( 0, m , n ) , ( l ,0, n ) , and ( l , m ,0 ) are denoted by 1, x 1 m x 2 n , x 0 1 x 2 n , and x 0 1 x 1 m , respectively. Let p = x 0 l x 1 m x 2 n . Then, points ( l + k , m , n ) , ( l , m + k , n ) and ( l , m , n + k ) are represented by monomials p x 0 k , p x 1 k and p x 2 k , respectively. Note that x i x j = x j x i for all pairs of i and j.

2.1. Flows of Triangles

Flows of triangles are defined using unit cubes in E3. As shown in Figure 1(a), unit cubes are piled-up in the direction of ( 1, 1, 1 ) in E3, where each of the three upper faces is divided into two triangles by the vertical diagonal (thick line). Then, a flow of triangles is obtained along the diagonals (Figure 1(b)). That is, the piled-up cubes form a mountain range-like structure and the vertical diagonals on its surface determine a flow of “slant” triangles on the slope.

As an example, let us consider the unit cube with the eight corner points 0, x 0 , x 0 x 1 , x 1 , x 2 , x 0 x 2 , x 0 x 1 x 2 , and x 1 x 2 (Figure 1(c)). Let P 0 = 1 , P 1 = x 0 , P 2 = x 0 x 1 , and P 3 = x 1 . Then, the upper face P 0 P 1 P 2 P 3 is divided into two “slant” triangles P 0 P 1 P 2 and P 0 P 3 P 2 . The triangle flow goes down (or up) along the edge P 0 P 2 at P 0 P 1 P 2 and P 0 P 3 P 2 .

In the following, we give the mathematical definition of the mountain range-like structure and the associated flow of triangles.

Definition 1 (Standard Lattice) The three-dimensional standard lattice L3 is the three-dimensional lattice generated by three vectors ( 1,0,0 ) , ( 0,1,0 ) , and ( 0,0,1 ) . Using the monomial representation, L3 is defined by

L 3 : = { x 0 l x 1 m x 2 n | l , m , n } E 3 .

Let p 0 , p 1 , and p 2 L 3 . We denote the convex hull of p 0 , p 1 , and p 2 by [ p 0 , p 1 , p 2 ] , i.e.,

[ p 0 , p 1 , p 2 ] : = { p 0 a p 1 b p 2 c | a , b , c , a , b , c 0 , a + b + c = 1 } ,

where is the set of all real numbers. The group of all permutations of the three-element set { 0,1,2 } is denoted by Sym3. For example, ρ ( 0 ) = 2 , ρ ( 1 ) = 1 , and ρ ( 2 ) = 0 for ρ = ( 02 ) S y m 3 , where ρ is written in cyclic notation.

Definition 2 (Slant Triangles) Let a L 3 and ρ S y m 3 . A slant triangle a [ x ρ ( 0 ) x ρ ( 1 ) ] is the triangle defined by three points a, a x ρ ( 0 ) , and a x ρ ( 0 ) x ρ ( 1 ) , i.e.,

a [ x ρ ( 0 ) x ρ ( 1 ) ] : = [ a , a x ρ ( 0 ) , a x ρ ( 0 ) x ρ ( 1 ) ] .

Figure 1. Flow of triangles: (a) A mountain range-like structure obtained by piling up unit cubes in the direction of ( 1, 1, 1 ) , whose peaks are P a = ( 2 , 0 , 0 ) , P b = ( 1 , 0 , 1 ) , P c = ( 1 , 2 , 1 ) , P d = ( 1 , 3 , 0 ) , and P e = ( 3 , 2 , 1 ) . The diagonal edges of “slant” triangles are drawn with thick lines. Shown above is the top view of the structure; (b) The flow of triangles determined by the mountain range-like structure of (a); (c) A unit cube and its top view (above), where P 0 = ( 0 , 0 , 0 ) , P 1 = ( 1 , 0 , 0 ) , P 2 = ( 1 , 1 , 0 ) , and P 3 = ( 0 , 1 , 0 ) .

The line segment joining vertex a and vertex a x ρ ( 0 ) x ρ ( 1 ) is called the diagonal edge of the slant triangle. The set S 2 of all slant triangles is defined by:

S 2 : = { a [ x ρ ( 0 ) x ρ ( 1 ) ] : a L 3 , ρ S y m 3 } .

Example 1 In the case of Figure 1(c),

P 0 P 1 P 2 = [ x 0 x 1 ] ( a = 1 and ρ = ( 012 ) ) , P 0 P 3 P 2 = [ x 1 x 0 ] ( a = 1 and ρ = ( 102 ) ) .

Their diagonal edges are the line segment P 0 P 2 .

By abuse of notation, we denote the vectors ( 0,1,1 ) , ( 1,0,1 ) , and ( 1,1,0 ) by the monomial x 1 x 2 , x 0 x 2 , and x 0 x 1 respectively in the following definition.

Definition 3 (Gradient of Slant Triangles) Let s = a [ x ρ ( 0 ) x ρ ( 1 ) ] S 3 , the gradient Ds of s is defined by

D s : = x ρ ( 0 ) x ρ ( 1 ) .

Example 2 In the case of Example 1, the slope of P 0 P 1 P 2 is given by

D ( P 0 P 1 P 2 ) = D [ x 0 x 1 ] = x 0 x 1 .

Flows of slant triangles along the diagonal edges are defined as follows.

Definition 4 (Local Trajectories of Slant Triangles) Let s S 2 . The local trajectory of slant triangles at s is a set of three consecutive slant triangles, consisting of s and two adjacent slant triangles which do not include the diagonal edge of s. By patching “consistent” local trajectories together, we will obtain a flow of slant triangles as shown in Figure 1(a).

Let s = a [ x ρ ( 0 ) x ρ ( 1 ) ] S 2 . The local trajectory of slant triangles at s is either

{ s D D , s , s U U } or { s D D , s , s U D } or { s D U , s , s U U } or { s D U , s , s U D } ,

where

{ s D D : = a x ρ ( 0 ) [ x ρ ( 1 ) x ρ ( 0 ) ] S 2 , s D U : = a x ρ ( 0 ) x ρ ( 2 ) 1 [ x ρ ( 2 ) x ρ ( 1 ) ] S 2 , s U D : = a [ x ρ ( 0 ) x ρ ( 2 ) ] S 2 , s U U : = a x ρ ( 1 ) 1 [ x ρ ( 1 ) x ρ ( 0 ) ] S 2

(Figure 2(a)).

A flow of “flat” triangles is defined on the hyperplane H 2 D using the projection π p t of E3 onto H 2 D , where

{ H 2 D : = { x 0 l x 1 m x 2 n | l , m , n , l + m + n = 0 } E 3 , π p t : E 3 H 2 D , π p t ( x 0 l x 1 m x 2 n ) : = x 0 ( 2 l m n ) / 3 x 1 ( l + 2 m n ) / 3 x 2 ( l m + 2 n ) / 3 .

Definition 5 (Flat Triangles) Let s = a [ x ρ ( 0 ) x ρ ( 1 ) ] S 2 . The projection π H of s on H 2 D is defined by

π H ( s ) : = [ π p t ( a ) , π p t ( a x ρ ( 0 ) ) , π p t ( a x ρ ( 0 ) x ρ ( 1 ) ) ] H 2 D .

π H ( s ) is called a flat triangle. The line segment joining π p t ( a ) and π p t ( a x ρ ( 0 ) x ρ ( 1 ) ) is called the diagonal edge of the flat triangle (Figure 1(c)). The set B 2 of all flat triangles on H 2 D is defined by:

B 2 : = { π H ( s ) | s S 2 } .

Remark π H ( s ) is a projection of triangles and π p t ( p ) is a projection of points.

By projecting slant triangles onto H 2 D , we obtain a two-dimensional flow of flat triangles on H 2 D .

Definition 6 (Local Trajectories of Flat Triangles) Let t B 2 . The local trajectory of flat triangles at t is a projection image of a local trajectory of slant triangles at s S 2 by π H , where π H ( s ) = t . That is, there exists some s S 2 such that π H ( s ) = t and the local trajectory at s is given by either

{ π H ( s D D ) , π H ( s ) , π H ( s U U ) } or { π H ( s D D ) , π H ( s ) , π H ( s U D ) } or { π H ( s D U ) , π H ( s ) , π H ( s U U ) } or { π H ( s D U ) , π H ( s ) , π H ( s U D ) } ,

where the definition of s D D , s D U , s U D , and s U U are given above (immediately after Definition 4). Note that π H ( s D D ) = π H ( s D U ) and π H ( s U D ) = π H ( s U U ) (Figure 2(a)). By patching “consistent” local trajectories together, we will obtain a flow of flat triangles as shown in Figure 1(a).

2.2. Vector Fields of Triangles

As shown in Figure 1(a), a mountain range-like structure induces a flow of triangles on B 2 . We can define a “tangent space” structure on the space B 2 of flat triangles, where each flat triangle assume one of the three gradient vectors x 1 x 2 , x 0 x 2 , and x 0 x 1 .

Figure 2. Slant triangles: (a) The four local trajectories of slant triangles at s B 2 . Shown above is the top view of the trajectory; (b) Positional relationship of the slant triangles of π H 1 ( [ x 1 x 2 ] ) and the 2-face { ( x 1 x 2 ) ( x 0 x 2 ) m 2 ( x 0 x 1 ) n | 0 m , n } of the cotangent cone C o n e * { ( x 1 x 2 ) ( x 0 x 2 ) 2 } , where slant triangles shown are A: x 0 2 x 1 2 x 2 2 [ x 1 x 2 ] , B: x 0 x 1 2 x 2 2 [ x 0 x 1 ] , C: x 0 x 1 2 x 2 [ x 2 x 0 ] , D: x 0 x 1 x 2 [ x 1 x 2 ] , E: x 1 x 2 [ x 0 x 1 ] , F: x 1 [ x 2 x 0 ] , and G: [ x 1 x 2 ] . Shown above is a schematic diagram of the relationship between the three vertices of a slant triangle and a 2-face of C o n e * { ( x 1 x 2 ) ( x 0 x 2 ) 2 } . In the diagram, the diagonal edges of slant triangles are drawn with thick line, where the diagonal edges on the 2-face are colored black and the others are colored grey. Note that all slant triangles are projected onto the same flat triangle π H ( [ x 1 x 2 ] ) by π H . Triangles A, B, and C are included in C o n e * { ( x 1 x 2 ) ( x 0 x 2 ) 2 } . Triangle D intersects the 2-face of the cotangent cone. Triangles E, F,and G are located outside the cotangent cone.

Definition 7 (Tangent Space) The tangent space T B 2 on B 2 is defined by

{ T B 2 : = B 2 × { x 1 x 2 , x 0 x 2 , x 0 x 1 } , π : T B 2 B 2 , π ( ( t , x i x j ) ) : = t .

Let t B 2 . The tangent space at t is denoted by T B 2 [ t ] . Note that there exists a one-to-one correspondence

T B 2 [ t ] ~ { x 1 x 2 , x 0 x 2 , x 0 x 1 } .

An inverse function of the projection π , i.e., a vector field of T B 2 on B 2 , is induced by a mountain range-like structure consisting of piled-up unit cubes.

Definition 8 (Tangent Cones) Let A L 3 , the three-dimensional tangent cone C o n e A is defined by

C o n e A : = { a x 0 l x 1 m x 2 n | a A ,0 l , m , n } L 3 .

Roughly speaking, C o n e A is the triangular cone whose top vertices are given by A L 3 . The set of all the top vertices of a tangent cone c is denoted by t o p ( c ) . In general, t o p ( C o n e A ) A .

Definition 9 ( t t o p ( c ) ) Let c be a three-dimensional tangent cone. The peaks on the boundary of c is defined by

t t o p ( c ) : = { p t o p ( c ) | i { 0,1,2 } s .t . p ( x i ) N C o n e ( t o p ( c ) \ { p } ) for N } ,

where denotes the set of all natural numbers.

Example 3 Let c be the tangent cone corresponds to the mountain range-like structure of Figure 1(a). Then,

{ c = C o n e { x 0 2 , x 0 x 2 , x 0 1 x 1 2 x 2 , x 0 x 1 3 , x 0 3 x 1 2 x 2 1 } , t o p ( c ) = t t o p ( c ) = { x 0 2 , x 0 x 2 , x 0 1 x 1 2 x 2 , x 0 x 1 3 , x 0 3 x 1 2 x 2 1 } .

Definition 10 ( p t ( c ) ) Let c be a three-dimensional tangent cone. The surface lattice points p t ( c ) of c is the set of all the L3 lattice points on the surface of c, i.e.,

p t ( c ) : = { a x i l x j m | a t o p ( c ) , { i , j } { 0,1,2 } , 0 l , m s .t . a x i l x j m C o n e { b x 0 x 1 x 2 } for b t o p ( c ) } L 3 .

Example 4 The surface lattice points of the three-dimensional tangent cone C o n e { 1 } are given by

p t ( C o n e { 1 } ) = U l = 0 U m = 0 U n = 0 L 3 ,

where

{ U l = 0 : = { x 1 m x 2 n | 0 m , n } , U m = 0 : = { x 0 l x 2 n | 0 l , n } , U n = 0 : = { x 0 l x 1 m | 0 l , m } .

That is, p t ( C o n e { 1 } ) is the union of the three 2-faces of C o n e { 1 } .

Definition 11 ( d S c ) Let c be a three-dimensional tangent cone. Then, d S c is the set of all the slant triangles on the surface of c, i.e.

d S c : = { s S 2 | all the vertices of s are included in p t ( c ) } .

Definition 12 (Vector Fields) Let c be a three-dimensional tangent cone. The vector field V c induced by c on B 2 is defined by

V c ( t ) : = D s ( t B 2 ) ,

where s d S c such that t = π H ( s ) (Note that s is uniquely determined for each t B 2 ). The flow of triangles determined by V c is called the flow of triangles induced by c.

Example 5 In the case of Figure 1(c),

V c ( π ( [ x 0 x 1 ] ) ) = D [ x 0 x 1 ] = x 0 x 1 ,

where c = C o n e { 1 } (Recall that P 0 P 1 P 2 = [ x 0 x 1 ] ).

2.3. Contour of Closed Trajectories of Triangles

The ridge lines of tangent cones are given by three vectors ( 1,0,0 ) , ( 0,1,0 ) , and ( 0,0,1 ) . To compute the contour of the region of H 2 D swept by a set of closed trajectories of flat triangles, we will consider another type of triangular cones whose ridge lines are given by the slopes of slant triangles, i.e., ( 0,1,1 ) , ( 1,0,1 ) , and ( 1,1,0 ) .

Definition 13 (Conjugate Lattice) The conjugate lattice L 3 is the three-dimensional lattice generated by three vectors ( 0,1,1 ) , ( 1,0,1 ) , and ( 1,1,0 ) . Using the monomial representation, L 3 is given by

L 3 : = { x 0 m + n x 1 l + n x 2 l + m | l , m , n } L 3 .

Definition 14 (Cotangent Cones) Let A L 3 . The three-dimensional cotangent cone C o n e * A is defined by

C o n e * A : = { a ( x 1 x 2 ) l ( x 0 x 2 ) m ( x 0 x 1 ) n | a A , 0 l , m , n } L * 3 .

We denote the top vertices of a cotangent cone c by t o p ( c ) . Note that we define cotangent cones only for A L 3 .

For a given cotangent cone, we can put a “roof” on the cone.

Definition 15 (Cotangent Roofs) Let A L 3 . The three-dimensional cotangent roof R o o f * A is defined by

R o o f * A : = { p L 3 | N s .t . N > 0 and p ( x 1 x 2 ) N , p ( x 0 x 2 ) N , p ( x 0 x 1 ) N C o n e * A } .

Roughly speaking, R o o f * A is obtained by putting as many unit cubes of L 3 as possible on C o n e * A .

Example 6 In the case of Figure 3(a),

R o o f * t t o p ( c 0 ) = R o o f * { x 0 2 , x 0 x 2 , x 0 1 x 1 2 x 2 , x 0 x 1 3 , x 0 3 x 1 2 x 2 1 } = C o n e * { x 0 x 2 1 , 1 , x 0 1 x 1 } .

Definition 16 ( p t ( w ) ) . Let w be a three-dimensional cotangent cone. The surface lattice points p t ( w ) of w is the set of all the L 3 lattice points on the surface of w, i.e.,

p t ( w ) : = { a ( e 3 / x i ) l ( e 3 / x j ) m | a t o p ( w ) , { i , j } { 0,1,2 } ,0 l , m , a ( e 3 / x i ) l ( e 3 / x j ) m C o n e * { b e 3 2 } for b t o p ( w ) } L 3 ,

where e 3 : = x 0 x 1 x 2 (For example, e 3 / x 1 = x 0 x 2 and e 3 2 = x 0 2 x 1 2 x 2 2 ).

Example 7 The surface lattice points of the three-dimensional cotangent cone C o n e * { 1 } are given by

p t ( C o n e * { 1 } ) = V l = 0 V m = 0 V n = 0 L 3 ,

where

{ V l = 0 : = { ( x 0 x 2 ) m ( x 0 x 1 ) n | 0 m , n } , V m = 0 : = { ( x 1 x 2 ) l ( x 0 x 1 ) n | 0 l , n } , V n = 0 : = { ( x 1 x 2 ) l ( x 0 x 2 ) m | 0 l , m } .

Figure 3. The closed trajectory of Figure 1 (a) and the associated contour pair: (a) Tangent cone c 0 = C o n e { P a , P b , P c , P d , P e } and the associated cotangent roof w c ( c 0 ) = C o n e * { K a , K b , K c } ; (b) Inverted cotangent roof I R o o f * ϕ t ( c 0 ) = I C o n e * { K d , K e } ; (c) The region R c ( w , i v ) for the contour pair ( w , i v ) = ( C o n e * { K a , K b , K c } , I C o n e * { K d , K e } ) associated with c 0 . In the figure, P a = x 0 2 , P b = x 0 x 2 , P c = x 0 1 x 1 2 x 2 , P d = x 0 x 1 3 , P e = x 0 3 x 1 2 x 2 1 , P f = x 0 4 x 1 2 , P g = x 0 2 x 1 x 2 , P h = x 0 x 1 2 x 2 3 , P i = x 0 x 1 4 x 2 , P j = x 0 2 x 1 3 x 2 , P k = x 0 2 x 1 2 , P l = x 0 3 x 1 3 , K a = x 0 x 2 1 , K b = 1 , K c = x 0 1 x 1 , K d = x 0 4 x 1 3 x 2 , and K e = x 0 3 x 1 4 x 2 3 .

Let c be a tangent cone and w be a cotangent cone. Then, we can divide all the slant triangles of the flow induced by c into three groups: 1) inside w, 2) outside w, and 3) on the surface of w (Figure 2(b)). In particular, we can compute the contour of closed trajectories induced by a tangent cone using a cotangent cone as shown below.

Definition 17 ( w c ( c 0 ) ) . Let c 0 be a three-dimensional tangent cone. Suppose that t t o p ( c 0 ) L 3 . The three-dimensional cotangent roof w c ( c 0 ) associated with c 0 is defined by

w c ( c 0 ) : = R o o f * t t o p ( c 0 ) .

Definition 18 ( ϕ t ( c 0 ) ) . Let c 0 be a three-dimensional tangent cone. Suppose that t t o p ( c 0 ) L 3 . The contour vertices ϕ t ( c 0 ) with respect to c 0 is defined by

ϕ t ( c 0 ) : = p t ( c 0 ) p t ( w c ( c 0 ) ) L 3 .

That is, ϕ t ( c 0 ) is the set of L 3 lattice points on the intersection of the surface of c 0 and the surface of w c ( c 0 ) .

Definition 19 ( Φ t ( c 0 ) ) . Let c 0 be a three-dimensional tangent cone. Suppose that t t o p ( c 0 ) L 3 . Let Φ t ( c 0 ) E 3 be the polygonal line obtained by joining the adjacent L 3 lattice points of ϕ t ( c 0 ) . Since all the points of ϕ t ( c 0 ) are on the surface of w c ( c 0 ) , the points of ϕ t ( c 0 ) are connected along the surface of the associated cotangent roof. Note that Φ t ( c 0 ) forms a closed polygonal line if t o p ( c 0 ) is finite (Figure 3(a)).

Definition 20 ( m t ( c 0 ) ) . Let c 0 be a three-dimensional tangent cone. The one-dimensional surface mesh m t ( c 0 ) with respect to c 0 is defined by

m t ( c 0 ) : = π p t ( Φ t ( c 0 ) ) H 2 D .

We also call m t ( c 0 ) the (one-dimensional) contour with respect to c 0 .

Definition 21 Let c 0 be a three-dimensional tangent cone. Set

{ I N t ( c 0 ) : = { s d S c 0 | s is contained inside w c ( c 0 ) } , O U T t ( c 0 ) : = { s d S c 0 | s is contained outside w c ( c 0 ) } , B D t ( c 0 ) : = { s d S c 0 | s intersects with the surface of w c ( c 0 ) } .

Theorem 1 Let c be a three-dimensional tangent cone. Suppose that t o p ( c ) is finite and t t o p ( c ) L 3 . Then, w c ( c ) divides all the slant triangles of d S c (i.e., slant triangles on the surface of c) into two groups: inside the roof and outside the roof. That is, B D t ( c ) = and π H ( I N t ( c ) ) exactly corresponds to the region swept by all the closed trajectories of V c .

Proof. Φ t ( c ) forms a closed polygonal line because ϕ t ( c ) contains only a finite number of points. Note that Φ t ( c ) consists of the diagonal edges of slant triangles on the surface of c. Since flows of slant triangles go along the diagonal edge at each slant triangle, there is no slant triangle crossing Φ t ( c ) . That is, Φ t ( c ) divides the flow of the slant triangles on the surface of c into two parts: the inside Φ t ( c ) and the outside Φ t ( c ) . □

Definition 22 ( R t ( c ) ) . Let c be a three-dimensional tangent cone. Let V c be the vector field induced by c. We define the region R t ( c ) of H 2 D by

R t ( c ) : = s I N t ( c ) π H ( s ) H 2 D .

By Theorem 1, R t ( c ) corresponds to the region of H 2 D swept by all the closed trajectories of V c if t o p ( c ) is finite and t t o p ( c ) L 3 .

We can compute the contour of R t ( c ) instantly.

Corollary 1 (“Contour” of R t ( c ) ). Let c a three-dimensional tangent cone. Suppose that t o p ( c ) is finite and t t o p ( c ) L 3 . Then, m t ( c ) gives the “contour” of R t ( c ) .

Proof. It follows immediately from the theorem. □

Example 8 In the case of Figure 3(a), the closed polygonal line m t ( c 0 ) consists of 12 vertices and 12 line segments, where

{ c 0 = C o n e { P a , P b , P c , P d , P e } , w c ( c 0 ) = C o n e * { K a , K b , K c } , ϕ t ( c 0 ) = { P a , P g , P b , P b h , P h , P h c , P c , P c i , P i , P d , P j , P k , P l , P c , P f , P f a } ,

and P b h = x 0 x 1 x 2 2 , P h c = x 1 2 x 2 , P c i = x 1 3 x 2 , and P f a = x 0 3 x 1 .

2.4..Constraints on the Contour of Closed Trajectories

Let c be a three-dimensional tangent cone. We have computed the contour of R t ( c ) using c and the associated cotangent roof w c ( c ) (the upper row of Figure 4). Now, we will compute regions of H 2 D without using tangent cones (the lower row of Figure 4).

Definition 23 (Inverted Cotangent Cones) Let A L 3 . The three-dimensional inverted cotangent cone I C o n e * A is defined by

I C o n e * A : = { a ( x 1 x 2 ) l ( x 0 x 2 ) m ( x 0 x 1 ) n | a A , 0 l , m , n } L * 3 .

We denote the top vertices of an inverted cotangent cone i v by i t o p ( i v ) .

Definition 24 (Inverted Cotangent Roofs) Let A L 3 . The three-dimensional inverted cotangent roof I R o o f * A is defined by

I R o o f * A : = { p L 3 | N s .t . N < 0 and p ( x 1 x 2 ) N , p ( x 0 x 2 ) N , p ( x 0 x 1 ) N I C o n e * A } .

Example 9 In the case of Figure 3(b),

I R o o f * ϕ t ( c ) = I R o o f * { x 0 4 x 1 2 , x 0 x 1 2 x 2 3 , x 0 x 1 4 x 2 , x 0 3 x 1 3 } = I C o n e * { x 0 4 x 1 3 x 2 , x 0 3 x 1 4 x 2 3 } .

Definition 25 ( p t ( i v ) ) . Let iv be a three-dimensional inverted cotangent cone. The surface lattice points p t ( i v ) of iv is the set of all the L 3 lattice points on the surface of iv, i.e.,

p t ( i v ) : = { a ( e 3 / x i ) l ( e 3 / x j ) m | a i t o p ( i v ) , { i , j } { 0,1,2 } ,0 l , m , a ( e 3 / x i ) l ( e 3 / x j ) m I C o n e * { b e 3 2 } for b i t o p ( i v ) } L * 3 .

Recall that e 3 = x 0 x 1 x 2 . For example, e 3 / x 1 = x 0 x 2 and e 3 2 = x 0 2 x 1 2 x 2 2 .

Figure 4. Schematic diagram showing the procedure for computing the contour of a region of H 2 D . The upper row shows the procedure for a region R t ( c ) specified by a tangent cone c. The lower row shows the procedure for a region R c ( w , i v ) specified by a contour pair ( w , i v ) .

Definition 26 (Contour Pairs) Let w be a three-dimensional cotangent cone. Let iv be a three-dimensional inverted cotangent cone. A pair ( w , i v ) of w and iv is called a three-dimensional contour pair.

Definition 27 ( ϕ c ( w , i v ) ) . Let ( w , i v ) be a three-dimensional contour pair. The contour vertices ϕ c ( w , i v ) with respect to ( w , i v ) is defined by

ϕ c ( w , i v ) : = p t ( w ) p t ( i v ) L 3 .

That is, ϕ c ( w , i v ) is the set of L 3 lattice points on the intersection of the surface of w and the surface of iv.

Definition 28 ( Φ c ( w , i v ) ) . Let ( w , i v ) be a three-dimensional contour pair. Let Φ c ( w , i v ) E 3 be the polygonal line obtained by joining the adjacent L 3 lattice points of ϕ c ( w , i v ) (Figure 3(c)). Since all the points of ϕ c ( w , i v ) are on the surface of w (or iv), the points of ϕ c ( w , i v ) are connected along the surface of w (or iv). Note that ϕ c ( w , i v ) forms a closed polygonal line if t o p ( w ) or i t o p ( i v ) is finite (Figure 3(c)).

Definition 29 ( m c ( w , i v ) ) . Let ( w , i v ) be a three-dimensional contour pair. The one-dimensional surface mesh m c ( w , i v ) with respect to ( w , i v ) is defined by

m c ( w , i v ) : = π p t ( Φ c ( w , i v ) ) H 2 D .

We also call m c ( w , i v ) the (one-dimensional) contour with respect to ( w , i v ) .

Remark Note that we have two types of one-dimensional surface meshes, i.e., the contour m t ( c ) with respect to a tangent cone c and the contour m c ( w , i v ) with respect to a contour pair ( w , i v ) .

Definition 30 ( R c ( w , i v ) ) . Let ( w , i v ) be a three-dimensional contour pair. R c ( w , i v ) H 2 D is the region enclosed by m c ( w , i v ) (Figure 3(c)).

Definition 31 ( ( w c ( c 0 ) , i v c ( c 0 ) ) ) . Let c 0 be a three-dimensional tangent cone. Suppose that t o p ( c 0 ) is finite and t t o p ( c 0 ) L 3 . The three-dimensional contour pair ( w c ( c 0 ) , i v c ( c 0 ) ) associated with c 0 is a pair of the associated cotangent roof w c ( c 0 ) and the inverted cotangent roof i v c ( c 0 ) which is defined by

i v c ( c 0 ) : = I R o o f * ϕ t ( c 0 ) .

Theorem 2 Let c be a three-dimensional tangent cone. Suppose that t o p ( c ) is finite and t t o p ( c ) L 3 . Then,

R t ( c ) = R c ( w c ( c ) , i v c ( c ) ) .

Proof. Recall that

ϕ t ( c ) = p t ( c ) p t ( w c ( c ) ) , ϕ c ( w c ( c ) , i v c ( c ) ) = p t ( w c ( c ) ) p t ( i v c ( c ) ) .

Note that i v c ( c ) ϕ t ( c ) , which implies p t ( i v c ( c ) ) ϕ t ( c ) because the slope inclination of cotangent cones is steeper than that of tangent cones. On the other hand, p t ( w c ( c ) ) ϕ t ( c ) by definition. Therefore, ϕ c ( w c ( c ) , i v c ( c ) ) ϕ t ( c ) . Suppose that ϕ c ( w c ( c ) , i v c ( c ) ) ϕ t ( c ) . Let p be an L 3 lattice point of ϕ c ( w c ( c ) , i v c ( c ) ) not included in ϕ t ( c ) . Then, p p t ( i v c ( c ) ) and p p t ( c ) . That is, p resides on the part of the surface of i v c ( c ) which is expanded by the “roof” operation. But the expanded part of i v c ( c ) is strictly contained in w c ( c ) , i.e., p p t ( w c ( c ) ) , which is a contradiction. Therefore, ϕ t ( c ) = ϕ c ( w c ( c ) , i v c ( c ) ) . In particular, R t ( c ) = R c ( w c ( c ) , i v c ( c ) ) . □

Example 10 In the case of Figure 3,

{ c 0 = C o n e { P a , P b , P c , P d , P e } , w c ( c 0 ) = C o n e * { K a , K b , K c } , i v c ( c 0 ) = I R o o f * { P a , P b , P c , P d , P e , P f , P g , P h , P i , P j , P k , P l } = I R o o f * { P f , P h , P i , P k } = I C o n e * { K d , K e } .

In Theorem 2, we have computed the “contour” of R t ( c ) for a give tangent cone c (Figure 5(a)). Now, we will compute the “contour” of R c ( w , i v ) for a given contour pair ( w , i v ) (Figure 5(b)).

Definition 32 Let ( w , i v ) be a three-dimensional contour pair. The tangent cone c t ( w , i v ) associated with ( w , i v ) is defined by

c t ( w , i v ) : = C o n e ϕ c ( w , i v ) .

Note that t o p ( c t ( w , i v ) ) L 3 .

Theorem 3 Let ( w , i v ) be a three-dimensional contour pair. Then,

R c ( w , i v ) R t ( c t ( w , i v ) ) .

Proof. Suppose that q H 2 D such that q R t ( c t ( w , i v ) ) and q R c ( w , i v ) . Then, there exists an L 3 lattice point p of ϕ t ( c t ( w , i v ) ) (i.e., points on the contour) such that π p t ( p ) R c ( w , i v ) ( π p t ( p ) is on the contour of R t ( c t ( w , i v ) ) .) It follows that p is outside iv because p p t ( w ) . However, the part of the surface of w not included in iv is strictly contained in c. Therefore, p p t ( c t ( w , i v ) ) , which is a contradiction. □

Remark In general, c t ( w , i v ) ϕ c ( w , i v ) does not imply p t ( c t ( w , i v ) ) ϕ c ( w , i v ) because the slope inclination of tangent cones is gentler than that of cotangent cones. Therefore, there exists a contour pair ( w , i v ) s.t.

ϕ t ( c t ( w , i v ) ) = p t ( c t ( w , i v ) ) p t ( w c ( c t ( w , i v ) ) ) = p t ( c t ( w , i v ) ) p t ( w ) ϕ c ( w , i v ) .

In particular, ϕ t ( c t ( w , i v ) ) ϕ c ( w , i v ) , i.e., R t ( c t ( w , i v ) ) R c ( w , i v ) .

Starting with some definitions, we will consider the correspondence between the two types of regions of H 2 D , i.e., R t ( c ) s and R c ( w , i v ) s.

Definition 33 Sets of there-dimensional cones are defined by

Figure 5. Maps between R t , n and R c , n and the relevant maps between sets of n + 1 -dimensional cones ( n = 2 or 3): (a) and the relevant maps; (b) τ n : R c , n R t , n and the relevant maps. The SECT denotes the set of all self-eclipsed closed trajectory complexes.

{ T C 2 : = { C o n e A | A L 3 and A is finite } , C C 2 : = { C o n e * A | A L 3 and A is finite } , I C C 2 : = { I C o n e * A | A L 3 and A is finite } .

Note that t o p ( c ) L 3 for c T C 2 .

Definition 34 ( R t ,2 and R c ,2 ) Sets of regions of H 2 D are defined by

{ R t ,2 : = { R t ( c ) | c T C 2 } , R c ,2 : = { R c ( w , i v ) | ( w , i v ) C C 2 × I C C 2 } .

Definition 35 ( ι 2 and τ 2 ) Maps between regions of H 2 D are defined by

It is not difficult to show that the maps are well-defined. By Theorem 2, we have

ι 2 ( R t ( c ) ) = R t ( c ) .

See Figure 5 for the correlation between the relevant maps.

Theorem 4 ι 2 is not surjective.

Proof. Let us consider the case of Figure 6(a). We have

{ w 0 = C o n e * { K a } , i v 0 = I C o n e * { K b } , ϕ c ( w 0 , i v 0 ) = { P a , P a b , P b , P c , P c d , P d , P d e , P e , P f , P f a } ,

where K a = x 0 1 x 2 1 , K b = x 0 2 x 1 4 x 2 2 , P a = 1 , P a b = x 1 x 2 , P b = x 1 2 x 2 2 , P c = x 0 1 x 1 2 x 2 , P c d = x 1 3 x 2 , P d = x 0 x 1 4 x 2 , P d e = x 0 x 1 3 , P e = x 0 x 1 2 x 2 1 , P f = x 0 2 x 1 2 , and P f a = x 0 x 1 . Then, R c ( w 0 , i v 0 ) forms a flattened hexagon as shown in the figure.

Figure 6. Self-eclipsed close trajectories: (a) Contour pair ( w 0 , i v 0 ) and R c ( w 0 , i v 0 ) ; (b) The tangent cone c 0 : = c c ( w 0 , i v 0 ) and the self-eclipsed region R t ( c 0 ) ; (c) The contour pair ( w c ( c 0 ) , i v c ( c 0 ) ) associated with c 0 .

On the other hand, the tangent cone c t ( w 0 , i v 0 ) and its associated cotangent roof is given by

{ c t ( w 0 , i v 0 ) = C o n e { P a , P b , P c , P d , P e , P f } = C o n e { P a , P c , P e } , w c ( c t ( w 0 , i v 0 ) ) = R o o f * { P a , P c , P e } = C o n e { K a } , ϕ t ( c t ( w 0 , i v 0 ) ) = { P a , P a b , P b , P c , P g , P h , P i , P e , P f , P f a }

(Figure 6(b)). Then, the region R t ( c t ( w 0 , i v 0 ) ) gets dented on the bottom, where P g = x 1 3 x 2 , P h = x 1 2 , and P i = x 0 x 1 3 . That is,

τ 2 ( R c ( w 0 , i v 0 ) ) R c ( w 0 , i v 0 ) .

It follows immediately that R c ( w 0 , i v 0 ) ι 2 ( R t ,2 ) .

Corollary 2 (Self-eclipse of Rc(w, iv)) There exists a contour pair ( w , i v ) such that τ 2 ( R c ( w , i v ) ) R c ( w , i v ) .

That is, not all closed polygonal lines defined by contour pairs correspond to a closed trajectory of triangles induced by a tangent cone. In other words, there exist global geometrical constraints on the contour of closed trajectories of triangles.

Definition 36 (Self-eclipsed closed trajectory complexes) A contour pair ( w , i v ) is called a self-eclipsed contour pair if τ 2 ( R c ( w , i v ) ) R c ( w , i v ) . The complex of closed trajectories of triangles induced by a three-dimensional tangent cone c is called a self-eclipsed closed trajectory complex (abbreviated as SECT) if there exists a self-eclipsed contour pair ( w , i v ) such that c = c t ( w , i v ) .

Example 11 In the case of Figure 6(b),

τ 2 ( R c ( C o n e * { K a } , I C o n e * { K b } ) ) = R t ( C o n e { P a , P c , P e } ) = R c ( C o n e * { K a } , I C o n e * { K c , K d } ) R c ( C o n e * { K a } , I C o n e * { K b } ) ,

where K c = x 0 2 x 1 3 x 2 and K d = x 0 x 1 3 x 2 2 (Figure 6(c)). That is, the closed trajectory of Figure 6(b) is a self-eclipsed closed trajectory complex.

In the next section, we will consider geometrical constraints on the shape of closed trajectories of tetrahedrons. As an example, it will be shown that a garlic bulb-like structure appears as a result of the constraints, where a flattened dodecahedron gets dented on the bottom and has vertical linear grooves on the side (Figure 11).

3. Mathematical Toy Model of Protein Molecules

3.1. Flows of Tetrahedrons

Now let us consider the case of tetrahedrons. To define a flow of tetrahedrons, we use unit cubes in the four-dimensional Euclidean space E4. By piling up unit cubes in the direction of ( 1, 1, 1, 1 ) , we will obtain a flow of “slant” tetrahedrons as in the case of flows of triangles.

Definition 37 (Standard Lattice) The four-dimensional standard lattice L4 is defined by

L 4 : = { x 0 l x 1 m x 2 n x 3 k | l , m , n , k } E 4 .

Let p 0 , p 1 , p 2 , and p 3 L 4 . We denote the convex hull of four points p 0 , p 1 , p 2 , and p 3 by [ p 0 , p 1 , p 2 , p 3 ] , i.e.,

[ p 0 , p 1 , p 2 , p 3 ] : = { p 0 a p 1 b p 2 c p 3 d | 0 a , b , c , d , a + b + c + d = 1 } ,

We denote the group of all permutations of the four-element set { 0,1,2,3 } by S y m 4 .

Definition 38 (Slant Tetrahedrons) Let a L 4 and ρ S y m 4 . A slant tetrahedron a [ x ρ ( 0 ) x ρ ( 1 ) x ρ ( 2 ) ] is the tetrahedron defined by four points a, a x ρ ( 0 ) , a x ρ ( 0 ) x ρ ( 1 ) , and a x ρ ( 0 ) x ρ ( 1 ) x ρ ( 2 ) , i.e.,

a [ x ρ ( 0 ) x ρ ( 1 ) x ρ ( 2 ) ] : = [ a , a x ρ ( 0 ) , a x ρ ( 0 ) x ρ ( 1 ) , a x ρ ( 0 ) x ρ ( 1 ) x ρ ( 2 ) ] .

The line segment joining vertex a and vertex a x ρ ( 0 ) x ρ ( 1 ) x ρ ( 2 ) , i.e., the cube-diagonal, is called the diagonal edge of the slant tetrahedron. Then, the four upper faces of each unit cube are divided into six tetrahedrons along the diagonal edge as shown in Example 12. The set S 3 of all slant tetrahedrons is defined by:

S 3 : = { a [ x ρ ( 0 ) x ρ ( 1 ) x ρ ( 2 ) ] | a L 4 , ρ S y m 4 } .

Example 12 Shown in Figure 7(a) is a four-dimensional unit cube at the origin P 1 . The upper face P 1 P y P y z P z P x P x y P x y z P x x of the cube is divided into six tetrahedrons along the cube-diagonal P 1 P x y z :

{ P 1 P x P x y P x y z = [ x 0 x 1 x 2 ] = [ 1 , x 0 , x 0 x 1 , x 0 x 1 x 2 ] , P 1 P y P x y P x y z = [ x 1 x 0 x 2 ] = [ 1 , x 1 , x 1 x 0 , x 1 x 0 x 2 ] , P 1 P y P y z P x y z = [ x 1 x 2 x 0 ] = [ 1 , x 1 , x 1 x 2 , x 1 x 2 x 0 ] , P 1 P z P y z P x y z = [ x 2 x 1 x 0 ] = [ 1 , x 2 , x 2 x 1 , x 2 x 1 x 0 ] , P 1 P z P x z P x y z = [ x 2 x 0 x 1 ] = [ 1 , x 2 , x 2 x 0 , x 2 x 0 x 1 ] , P 1 P x P x z P x y z = [ x 0 x 2 x 1 ] = [ 1 , x 0 , x 0 x 2 , x 0 x 2 x 1 ] .

Figure 7. Slant tetrahedrons: (a) A four-dimensional unit cube shown in the Schlegel diagram (below) and its projection image on hyperplane H 3 D (above), where P x l y m z n w k = x 0 l x 1 m x 2 n x 3 k L 4 and Q x l y m z n w k = π p t ( x 0 l x 1 m x 2 n x 3 k ) H 3 D ; (b) All the local trajectories at s = P 1 P x P x y P x y z B 3 . The diagonal edges of tetrahedrons are drawn with thick lines.

Definition 39 (Gradient of Slant Tetrahedrons) Let s = a [ x ρ ( 0 ) x ρ ( 1 ) x ρ ( 2 ) ] S 4 . The gradient Ds of s is defined by

D s : = x ρ ( 0 ) x ρ ( 1 ) x ρ ( 2 ) .

By abuse of notation, we denote vector ( l , m , n , k ) by monomial x 0 l x 1 m x 2 n x 3 k in the definition.

Example 13 In the case of Example 12,

D ( P 1 P x P x y P x y z ) = D [ x 0 x 1 x 2 ] = x 0 x 1 x 2 .

Flows of slant tetrahedrons along the diagonal edges are defined as follows.

Definition 40 (Local Trajectories of Slant Tetrahedrons) Let s S 3 . A local trajectory of slant tetrahedrons at s is a set of three consecutive slant tetrahedrons, consisting of s and two adjacent slant tetrahedrons which do not include the diagonal edge of s. By patching “consistent” local trajectories together, we obtain a flow of slant tetrahedrons.

Let s = a [ x ρ ( 0 ) x ρ ( 1 ) x ρ ( 2 ) ] S 3 . Then, the local trajectory of slant tetrahedrons at s is either

{ s D D , s , s U U } or { s D D , s , s U D } or { s D U , s , s U U } or { s D U , s , s U D } ,

where

{ s D D : = a x ρ ( 0 ) [ x ρ ( 1 ) x ρ ( 2 ) x ρ ( 0 ) ] S 3 , s D U : = a x ρ ( 0 ) x ρ ( 3 ) 1 [ x ρ ( 3 ) x ρ ( 1 ) x ρ ( 2 ) ] S 3 , s U D : = a [ x ρ ( 0 ) x ρ ( 1 ) x ρ ( 3 ) ] S 3 , s U U : = a x ρ ( 2 ) 1 [ x ρ ( 2 ) x ρ ( 0 ) x ρ ( 1 ) ] S 3 .

Example 14 In the case of Example 13, the four local trajectories at P 1 P x P x y P x y z ( = [ x 0 x 1 x 2 ] ) are shown in Figure 7(b).

A flow of “flat” tetrahedrons is defined on the hyperplane H 3 D using the projection π p t of E 4 onto H 3 D , where

{ H 3 D : = { x 0 l x 1 m x 2 n x 3 k | l , m , n , k , l + m + n + k = 0 } E 4 , π p t : E 4 H 3 D , π p t ( x 0 l x 1 m x 2 n x 3 k ) : = x 0 ( 3 x y z w ) / 4 x 1 ( x + 3 y z w ) / 4 x 2 ( x y + 3 z w ) / 4 x 3 ( x y z + 3 w ) / 4 .

Definition 41 (Flat Tetrahedrons) Let s = a [ x ρ ( 0 ) x ρ ( 1 ) x ρ ( 2 ) ] S 3 . The projection π H of s on H 3 D is defined by

π H ( s ) : = [ π p t ( a ) , π p t ( a x ρ ( 0 ) ) , π p t ( a x ρ ( 0 ) x ρ ( 1 ) ) , π p t ( a x ρ ( 0 ) x ρ ( 1 ) x ρ ( 2 ) ) ] H 3 D .

π H ( s ) is called a flat tetrahedron. The line segment joining π p t ( a ) and π p t ( a x ρ ( 0 ) x ρ ( 1 ) x ρ ( 2 ) ) is called the diagonal edge of the flat tetrahedron. The set B 3 of all flat tetrahedrons on H 3 D is defined by:

B 3 : = { π H ( s ) | s S 3 } .

Example 15 In the case of Example 12 (Figure 7(a)), the projection image Q 1 Q x Q x y Q y Q z Q x z Q x y z Q y z of the upper face P 1 P x P x y P y P z P x z P x y z P y z by π L is divided into six flat tetrahedrons:

{ Q 1 Q x Q x y Q x y z = π H ( [ x 0 x 1 x 2 ] ) , Q 1 Q y Q x y Q x y z = π H ( [ x 1 x 0 x 2 ] ) , Q 1 Q y Q y z Q x y z = π H ( [ x 1 x 2 x 0 ] ) , Q 1 Q z Q y z Q x y z = π H ( [ x 2 x 1 x 0 ] ) , Q 1 Q z Q x z Q x y z = π H ( [ x 2 x 0 x 1 ] ) , Q 1 Q x Q x z Q x y z = π H ( [ x 0 x 2 x 1 ] ) .

Remark. Note that each tetrahedron has two long edges and four short edges, where the diagonal edge correspond to a short edge. Flows of tetrahedrons go along the diagonal edge at each tetrahedron.

By projecting slant tetrahedrons onto H 3 D , we obtain a three-dimensional flow of flat tetrahedrons on H 3 D .

Definition 42 (Local Trajectories of Flat Tetrahedrons) The local trajectory of flat tetrahedrons at t B 3 is a projection image of a local trajectory of slant tetrahedrons at s S 3 by π H , where π H ( s ) = t . The local trajectory at s is given by either

{ π H ( s D D ) , π H ( s ) , π H ( s U U ) } or { π H ( s D D ) , π H ( s ) , π H ( s U D ) } or { π H ( s D U ) , π H ( s ) , π H ( s U U ) } or { π H ( s D U ) , π H ( s ) , π H ( s U D ) } .

Note that π H ( s D D ) = π H ( s D U ) and π H ( s U D ) = π H ( s U U ) . By patching “consistent” local trajectories together, we obtain a flow of flat tetrahedrons as shown in Figure 8.

Figure 8. Flows of tetrahedrons: (a) Closed trajectories of flat tetrahedrons induced by c 0 = C o n e { Q x y z , Q x z w , Q x y w } , where Q x y z = x 0 x 1 x 2 , Q x z w = x 0 x 2 x 3 , and Q x y w = x 0 x 1 x 3 . Thick polygonal lines (black and grey) indicate the diagonal edges of flat tetrahedrons, where the black line indicate the polygonal line passing through Q x y z , Q x z w , and Q x y w . Note that there exist infinitely many closed trajectories of length six and length twelve; (b) Closed trajectories of flat tetrahedrons induced by c 1 = C o n e { Q x y z , Q x z w , Q x y w , Q y z w } , where Q y z w = x 1 x 2 x 3 . Thick polygonal lines (black and grey) indicate the diagonal edges of flat tetrahedrons. Black lines indicate the polygonal lines passing through Q x y z , Q x z w , Q x y w , or Q y z w . Grey polygonal lines correspond to closed trajectories of length six and length twelve.

3.2. Vector Fields of Tetrahedrons

The tangent space on the space B 3 of flat tetrahedrons is defined in the same way as the tangent space T B 2 on B 2 .

Definition 43 (Tangent Space) The tangent space T B 3 on B 3 is defined by

{ T B 3 : = B 3 × { x 1 x 2 x 3 , x 0 x 2 x 3 , x 0 x 1 x 3 , x 0 x 1 x 2 } , π : T B 3 B 3 , π ( ( t , x i x j x k ) ) : = t .

Let t B 3 . The tangent space at t is denoted by T B 3 [ t ] . Note that there exists a one-to-one correspondence

T B 3 [ t ] ~ { x 1 x 2 x 3 , x 0 x 2 x 3 , x 0 x 1 x 3 , x 0 x 1 x 2 } .

Tangent cones are also defined similarly for A L 4 .

Definition 44 (Tangent Cones) Let A L 4 . The four-dimensional tangent cone C o n e A is defined by

C o n e A : = { a x 0 l x 1 m x 2 n x 3 k | a A , 0 l , m , n , k } L 4 .

The set of all the top vertices of a cone c is denoted by t o p ( c ) .

Definition 45 t t o p ( c ) . Let c be a four-dimensional tangent cone. The peaks on the boundary of c is defined by

t t o p ( c ) : = { p t o p ( c ) | i { 0,1,2,3 } s .t . p ( x i ) N C o n e { t o p ( c ) \ { p } } for N } .

Example 16 In the case of Figure 8,

t t o p ( c 0 ) = { Q x y z , Q x z w , Q x y w } , t t o p ( c 1 ) = { Q x y z , Q x z w , Q x y w , Q y z w } .

Definition 46 ( p t ( c ) ) . Let c be a four-dimensional tangent cone. The surface lattice points p t ( c ) of c is the set of all the L4 lattice points included in the 3-faces of c, i.e.,

p t ( c ) : = { a x i l x j m x k n | a t o p ( c ) , { i , j , k } { 0,1,2,3 } ,0 l , m , k , a x i l x j m x k n C o n e { b x 0 x 1 x 2 x 3 } for b t o p ( c ) } L 4 .

Example 17 The surface lattice points of the four-dimensional tangent cone C o n e { 1 } are given by

p t ( C o n e { 1 } ) = U l = 0 U m = 0 U n = 0 U k = 0 L 4 ,

where

{ U l = 0 : = { x 1 m x 2 n x 3 k | 0 m , n , k } , U m = 0 : = { x 0 l x 2 n x 3 k | 0 l , n , k } , U n = 0 : = { x 0 l x 1 m x 3 k | 0 l , m , k } , U k = 0 : = { x 0 l x 1 m x 2 n | 0 l , m , n } .

That is, p t ( C o n e { 1 } ) is the union of the four 3-faces of C o n e { 1 } .

Definition 47 (dSc) Let c be a four-dimensional tangent cone. Then, d S c is the set of all the slant tetrahedrons included in the 3-faces of c, i.e.

d S c : = { a [ x ρ ( 0 ) x ρ ( 1 ) x ρ ( 2 ) ] S 3 | alltheverticesareincludedin p t ( c ) } .

Definition 48 (Vector Fields) Let c be a four-dimensional tangent cone. The vector field V c induced by c on B 3 is defined by

V c ( t ) : = D s ( t B 3 ) ,

where s d S c π H 1 ( t ) (Note that s is uniquely determined for each t B 3 ). The flow of tetrahedrons determined by V c is called the flow of tetrahedrons induced by c.

Example 18 In the case of Figure 7(a),

V c ( Q 1 Q x Q x y Q x y z ) = D ( P 1 P x P x y P x y z ) = x 0 x 1 x 2 ,

where c = C o n e { 1 } , Q 1 Q x Q x y Q x y z = π H ( P 1 P x P x y P x y z ) , and P 1 P x P x y P x y z = [ x 0 x 1 x 2 ] .

Unlike the case of flows of triangles, infinitely many closed trajectories are induced by a tangent cone.

Example 19 Shown in Figure 8(a) is the closed trajectories of the flow induced by C o n e { x 0 x 1 x 2 , x 0 x 2 x 3 , x 0 x 1 x 3 } . Two types of closed trajectories, one is length 6 and the other is length 12, are alternately stacked infinitely.

Example 20 By putting another top vertex x 1 x 2 x 3 on the tangent cone of Figure 8(a), we obtain a decomposition of a rhombic dodecahedron into four closed trajectories of tetrahedrons (Figure 8(b)). Then, each triplet of the four top vertices of C o n e { x 0 x 1 x 2 , x 0 x 2 x 3 , x 0 x 1 x 3 , x 1 x 2 x 3 } induces infinitely many closed trajectories outside the rhombic dodecahedron (grey polygonal lines).

3.3. Shape of Closed Trajectories of Tetrahedrons

To compute the surface (2-faces) of the region of H 3 D swept by a set of closed trajectories of flat tetrahedrons, we will consider another type of cones whose ridge lines are given by the “slopes” of slant tetrahedrons, i.e., ( 0,1,1,1 ) , ( 1,0,1,1 ) , ( 1,1,0,1 ) , and ( 1,1,1,0 ) .

Definition 49 (Conjugate Lattice) The four-dimensional conjugate lattice L 4 is defined by

L 4 : = { x 0 m + n + k x 1 l + n + k x 2 l + m + k x 3 l + m + n | l , m , n , k } L 4 .

Definition 50 (Cotangent Cones) Let A L 4 . The four-dimensional cotangent cone C o n e * A is defined by

C o n e * A : = { a ( x 1 x 2 x 3 ) l ( x 0 x 2 x 3 ) m ( x 0 x 1 x 3 ) n ( x 0 x 1 x 2 ) k | a A , 0 l , m , n , k } L 4 .

We denote the top vertices of a cotangent cone c by t o p ( c ) .

Definition 51 (Cotangent Roofs) Let A L 4 . The four-dimensional cotangent roof R o o f * A is defined by

R o o f * A : = { p L 4 | N s .t . N > 0 and p ( x 1 x 2 x 3 ) N , p ( x 0 x 2 x 3 ) N , p ( x 0 x 1 x 3 ) N , p ( x 0 x 1 x 2 ) N C o n e * A } .

Example 21 In the case of Figure 8,

R o o f * t t o p ( c 0 ) = R o o f * { x 0 x 1 x 2 , x 0 x 2 x 3 , x 0 x 1 x 3 } = C o n e * { x 0 x 1 x 2 , x 0 x 2 x 3 , x 0 x 1 x 3 } ,

R o o f * t t o p ( c 1 ) = R o o f * { x 0 x 1 x 2 , x 0 x 2 x 3 , x 0 x 1 x 3 , x 1 x 2 x 3 } = C o n e * { 1 } .

Definition 52 ( p t ( w ) ) . Let w be a four-dimensional cotangent cone. The surface lattice points p t ( w ) of w is the set of all the L 4 lattice points included in the 3-faces of w, i.e.,

p t ( w ) : = { a ( e 4 / x i ) l ( e 4 / x j ) m ( e 4 / x k ) n | a t o p ( w ) , { i , j , k } { 0,1,2,3 } , 0 l , m , n , a ( e 4 / x i ) l ( e 4 / x j ) m ( e 4 / x k ) n C o n e * { b e 4 3 } for b t o p ( w ) } L * 4 ,

where e 4 : = x 0 x 1 x 2 x 3 (For example, e 4 / x 1 = x 0 x 2 x 3 and e 4 3 = x 0 3 x 1 3 x 2 3 x 3 3 ).

Example 22 The surface lattice points of the four-dimensional cotangent cone C o n e * { 1 } are given by

p t ( C o n e * { 1 } ) = V l = 0 V m = 0 V n = 0 V k = 0 L 4 ,

where

{ V l = 0 : = { ( x 0 x 2 x 3 ) m ( x 0 x 1 x 3 ) n ( x 0 x 1 x 2 ) k | 0 m , n , k } , V m = 0 : = { ( x 1 x 2 x 3 ) l ( x 0 x 1 x 3 ) n ( x 0 x 1 x 2 ) k | 0 l , n , k } , V n = 0 : = { ( x 1 x 2 x 3 ) l ( x 0 x 2 x 3 ) m ( x 0 x 1 x 2 ) k | 0 l , m , k } , V k = 0 : = { ( x 1 x 2 x 3 ) l ( x 0 x 2 x 3 ) m ( x 0 x 1 x 3 ) n | 0 l , m , n } .

As in the case of slant triangles, all the slant tetrahedrons of the flow induced by a tangent cone are divided into three groups by a cotangent cone: 1) inside the cotangent cone, 2) outside the cotangent cone, and3) on the surface (3-faces) of the cotangent cone (Figure 9). Unlike the case of triangles, multiple types of slant tetrahedrons are on the surface (3-faces) of the cotangent cone (tetrahedrons B, C, D, E, F in Figure 9). Among them, only two types of tetrahedrons (B and F) flow through the surface (2-faces).

Definition 53 ( w c ( c 0 ) ) . Let c 0 be a four-dimensional tangent cone. Suppose that t t o p ( c 0 ) L 4 . The four-dimensional cotangent roof w c ( c 0 ) associated with c 0 is defined by

w c ( c 0 ) : = R o o f * t t o p ( c 0 ) .

Definition 54 ( ϕ t ( c 0 ) ) . Let c 0 be a four-dimensional tangent cone. Suppose that t t o p ( c 0 ) L 4 . The surface vertices ϕ t ( c 0 ) with respect to c 0 is defined by

ϕ t ( c 0 ) : = p t ( c 0 ) p t ( w c ( c 0 ) ) L 4 .

That is, ϕ t ( c 0 ) is the set of L 4 lattice points on the intersection of the 3-faces of c 0 and the 3-faces of w c ( c 0 ) .

Definition 55 ( Φ t ( c 0 ) ) . Let c 0 be a four-dimensional tangent cone. Suppose that t t o p ( c 0 ) L 4 . Let Φ t ( c 0 ) E 4 be the set of the polygonal lines obtained by joining the adjacent L 4 lattice points of surface vertices ϕ t ( c 0 ) . Since all the points of ϕ t ( c 0 ) are included in the 3-faces of w c ( c 0 ) , the points of ϕ t ( c 0 ) are connected along the 3-faces of the associated cotangent roof.

Definition 56 ( m t ( c 0 ) ) . Let c 0 be a four-dimensional tangent cone. The two-dimensional surface mesh m t ( c 0 ) with respect to c 0 is defined by

m t ( c 0 ) : = π p t ( Φ t ( c 0 ) ) H 3 D .

We also call m t ( c 0 ) (two-dimensional) surface vein with respect to c 0 .

Definition 57 Let c 0 be a four-dimensional tangent cone. Set

{ I N t ( c 0 ) : = { s d S c 0 | s iscontainedinside w c ( c 0 ) } , O U T t ( c 0 ) : = { s d S c 0 | s iscontainedoutside w c ( c 0 ) } , B D t ( c 0 ) : = { s d S c 0 | s intersectswiththe3-facesof w c ( c 0 ) } .

Definition 58 ( R t ( c 0 ) ) . Let c 0 be a four-dimensional tangent cone. We define the region R t ( c 0 ) of H 3 D by

Figure 9. Positional relationship of slant triangles of π H 1 ( [ x 0 x 1 x 2 ] ) and the 3-face { ( x 0 x 2 x 3 ) m ( x 0 x 1 x 3 ) n ( x 0 x 1 x 2 ) k | 0 m , n , k } of the cotangent cone C o n e * { 1 } . Slant triangles shown are A: x 0 2 x 1 2 x 2 x 3 [ x 2 x 3 x 0 ] , B: x 0 2 x 1 x 2 x 3 [ x 1 x 2 x 3 ] , C: x 0 x 1 x 2 x 3 [ x 0 x 1 x 2 ] , D: x 0 x 1 x 2 [ x 3 x 0 x 1 ] , E: x 0 x 1 [ x 2 x 3 x 0 ] , F: x 0 [ x 1 x 2 x 3 ] , and G: [ x 0 x 1 x 2 ] . Shown above is a schematic diagram of the relationship between the four vertices of a slant tetrahedron and the 3-face of C o n e * { 1 } . In the figure, the diagonal edges of slant tetrahedrons are drawn with thick line, where the diagonal edges included in the 3-face are colored black and the others are colored grey. Tetrahedron A is included in C o n e * { 1 } . Tetrahedrons B, C, D, E, and F intersect the 3-face of C o n e * { 1 } . Tetrahedron G is located outside C o n e * { 1 } . Note that the diagonal edges of tetrahedrons B and F cross the 3-face of C o n e * { 1 } .

R t ( c 0 ) : = s I N t ( c 0 ) π H ( s ) H 3 D .

In the case of flows of triangles, R t ( c ) corresponds to all the closed trajectories of V c if t o p ( c ) is finite and t t o p ( c ) L 3 (Theorem 1). In the case of flows of tetrahedrons, we have the following result.

Theorem 5 There exist a four-dimensional tangent cone c such that R t ( c ) does not contain all the closed trajectories of V c .

Proof. For example, set c 0 = C o n e * { x 0 x 1 x 2 , x 0 x 2 x 3 , x 0 x 1 x 3 } (Figure 8(a)). Then, V c 0 induces infinitely many closed trajectories of tetrahedrons. However, we can not construct a cotangent cone which covers all the closed trajectories because w c ( c 0 ) = c 0 (Note that more than three vertices are required to construct a “roof” on a tangent cone.) In particular, R t ( c 0 ) does not contain all the closed trajectories of V c 0 . □

As in the case of flows of triangles, we can compute R t ( c ) as the intersection of w c ( c ) and c if t t o p ( c ) L 4 .

Example 23 In the case of Figure 8(b) (or Figure10(a)),

{ c 1 = C o n e { x 0 x 1 x 2 , x 0 x 2 x 3 , x 0 x 1 x 3 , x 1 x 2 x 3 } , w c ( c 1 ) = R o o f * { x 0 x 1 x 2 , x 0 x 2 x 3 , x 0 x 1 x 3 , x 1 x 2 x 3 } = C o n e * { 1 } , ϕ t ( c 1 ) = { P y z w , P x z w , P x y w , P x y z , P x 2 y 2 z w , P x 2 y z 2 w , P x 2 y z w 2 , P x y 2 z 2 w , P x y 2 z w 2 , P x y z 2 w 2 } ,

where P x l y m z n w k = x 0 l x 1 m x 2 n x 3 k L 4 . Then,

Figure 10. Rhombic dodecahedrons: (a) Four closed trajectories of the flow induced by c 1 = C o n e { x 0 x 1 x 2 , x 0 x 2 x 3 , x 0 x 1 x 3 , x 1 x 2 x 3 } (See also Figure 8 (b)); (b) The closed trajectory of length 24 induced by c 2 = C o n e { x 0 x 2 , x 1 x 2 , x 1 x 3 } ; (c) The closed trajectory of length 24 induced by c 3 = C o n e { x 0 x 2 , x 1 x 2 , x 2 x 3 , x 0 x 1 x 3 } . Thick polygonal lines (black and grey) indicate the diagonal edges of flat tetrahedrons, where the black lines indicate the polygonal lines included in w c ( c i ) ( i = 1 , 2 , 3 ).

R t ( c 1 ) is the rhombic dodecahedron surrounded by the black thick lines, which consists of four closed trajectories of length six, i.e., consists of 24 tetrahedrons.

All the 24 tetrahedrons of R t ( c 1 ) are type D of Figure 9:

1) Q 1 is inside w c ( c 1 ) ,

2) Q x 2 y z w , Q x y 2 z w , Q x y z 2 w and Q x y z w 2 are outside w c ( c 1 ) ,

3) all the diagonal edges of the tetrahedrons p t ( w c ( c 1 ) ) ,

where Q x l y m z n w k = π p t ( x 0 l x 1 m x 2 n x 3 k ) H 3 D . Other closed trajectories of the vector field are outside w c ( c 1 ) .

In the case of flows of tetrahedrons, we should also consider the case of t t o p ( c ) L 4 .

Example 24 Putting more unit cubes on the tangent cone c 1 of Figure10(a), we obtain rhombic dodecaherons consisting of a closed trajectory of length 24 (Figure 10(b) and Figure 10(c)). However, we can not compute the shape of the rhombic dodecaherons using cotangent roofs. Cotangent roofs are not defined because

t t o p ( c i ) L 4 ( i = 2 , 3 ) .

Remark In the three-dimensional case, t t o p ( c ) L 3 implies the existence of loopholes on the contour, i.e., the existence of the triangles of type D (Figure 2(b)). On the other hand, in the four-dimensional case, t t o p ( c ) L 4 implies not only the existence of loopholes but also dents and bulges on the surface. That is, the tetrahedrons of type B and F correspond to loopholes, the tetrahedrons of type C to dents, and the tetrahedrons of type E to bulges (Figure 9).

To define a four-dimensional cotangent roof R o o f * A for any A L 4 , we consider a set of the “closest L 4 lattice points to a” for each a A .

Definition 59 (STAND) Let a L 4 . The stand of a is defined by

s t a n d ( a ) : = { { a } if d e g ( a ) 0 m o d 3, { a x i | 0 i 3 } if d e g ( a ) 2 m o d 3, { a x i x j | 0 i , j 3 } if d e g ( a ) 1 m o d 3 ,

where d e g ( x 0 l x 1 m x 2 n x 3 k ) : = l + m + n + k . Note that s t a n d ( a ) L 4 .

Let A L 4 . The stand of A is defined by

S T A N D ( A ) : = a A s t a n d ( a ) L 4 .

Definition 60 (Extended Cotangent Roofs) Let A L 4 such that A L 4 . The four-dimensional extended cotangent roof E r o o f * A is defined by

E r o o f * A : = R o o f * S T A N D ( A ) .

Definition 61 (Extended w c ( c 0 ) ) Let c 0 be a four-dimensional tangent cone. The four-dimensional (extended) cotangent roof w c ( c 0 ) associated with c 0 is defined by

w c ( c 0 ) : = E r o o f * t t o p ( c 0 ) .

Note that

E r o o f * t t o p ( c 0 ) = R o o f * t t o p (c0)

if t t o p ( c 0 ) L 4 .

Definition 62 (Extended ϕ t ( c 0 ) ) Let c 0 be a four-dimensional tangent cone. The (extended) surface vertices ϕ t ( c 0 ) with respect to c 0 is defined by

ϕ t ( c 0 ) : = p t ( c 0 ) p t ( w c ( c 0 ) ) L 4 .

Example 25 In the case of Figure 10(b),

{ c 2 = C o n e { x 0 x 2 , x 1 x 2 , x 1 x 3 } , w c ( c 2 ) = E r o o f * { x 0 x 2 , x 1 x 2 , x 1 x 3 } = C o n e * { 1 } , ϕ t ( c 2 ) = { P x z w , P x 2 y z w 2 , P x y w , P x 2 y 2 z w , P x y z , P x y z 2 w 2 , P y z w } .

Then, R t ( c 2 ) is the rhombic dodecahedron surrounded by the black and grey thick lines, which consists of a closed trajectories of length 24.

Note that x 0 x 2 , x 1 x 2 , and x 1 x 3 are outside C o n e * { 1 } . The rhombic dodecahedron consists of not only type D but also type E tetrahedrons of Figure 9. In the figure, the tetrahedron with the grey diagonal edge are type E. The tetrahedrons with the black diagonal edge are type D.

Example 26 In the case of Figure 10(c),

{ c 3 = C o n e { x 0 x 2 , x 1 x 2 , x 2 x 3 , x 0 x 1 x 3 } , w c ( c 3 ) = E r o o f * { x 0 x 2 , x 1 x 2 , x 2 x 3 , x 0 x 1 x 3 } = C o n e * { 1 } , ϕ t ( c 3 ) = { P x z w , P x 2 y z w 2 , P x y w , P x 2 y 2 z w , P x y z , P x y 2 z w 2 , P y z w } .

Then, R t ( c 3 ) is the rhombic dodecahedron surrounded by the black and grey thick lines, which consists of a closed trajectories of length 24. In the figure, the tetrahedrons with the grey diagonal edge are type E. The tetrahedrons with the black diagonal edge are type D.

3.4. Constraints on the Shape of Protein Molecules

In this paper, we consider the shape of complexes of closed trajectories of tetrahedrons as a simplified geometrical model of protein molecules. As in the case of flows of triangles, we will specify the shape of regions of H 3 D using a pair of cotangent cones.

Definition 63 (Inverted Cotangent Cones) Let A L 4 . The four-dimensional inverted cotangent cone I C o n e * A is defined by

I C o n e * A : = { a ( x 1 x 2 x 3 ) l ( x 0 x 2 x 3 ) m ( x 0 x 1 x 3 ) n ( x 0 x 1 x 2 ) k | a A , 0 l , m , n , k } L * 4 .

We denote the top vertices of an inverted cotangent cone iv by i t o p ( i v ) .

Definition 64 (Inverted Cotangent Roofs) Let A L 4 . The four-dimensional inverted cotangent roof I R o o f * A is defined by

I R o o f * A : = { p L 4 | N s .t . N < 0 and p ( x 1 x 2 x 3 ) N , p ( x 0 x 2 x 3 ) N , p ( x 0 x 1 x 3 ) N , p ( x 0 x 1 x 2 ) N I C o n e * A } .

Let A L 4 . The four-dimensional extended inverted cotangent roof E I r o o f * A is defined by

E I r o o f * A : = I R o o f * S T A N D ( A ) L 4 .

Example 27 In the case of Figure 10(a),

I R o o f * t t o p ( c 1 ) = I R o o f * { x 0 x 1 x 2 , x 0 x 2 x 3 , x 0 x 1 x 3 , x 1 x 2 x 3 } = I C o n e * { x 0 3 x 1 3 x 2 3 x 3 3 } .

Definition 65 ( p t ( i v ) ) . Let iv be a four-dimensional inverted cotangent cone. The surface lattice points p t ( i v ) of iv is the set of all the L 4 lattice points included in the 3-faces of iv, i.e.,

p t ( i v ) : = { a ( e 4 / x i ) l ( e 4 / x j ) m ( e 4 / x k ) n | a i t o p ( i v ) , { i , j , k } { 0,1,2,3 } , 0 l , m , n , a ( e 4 / x i ) l ( e 4 / x j ) m ( e 4 / x k ) n I C o n e * { b e 4 3 } for b i t o p ( i v ) } L 4 .

Recall that e 4 = x 0 x 1 x 2 x 3 . For example, e 4 / x 2 = x 0 x 1 x 3 and e 4 3 = x 0 3 x 1 3 x 2 3 x 3 3 .

Definition 66 (Surface Pairs) Let w be a three-dimensional cotangent cone. Let iv be a three-dimensional inverted cotangent cone. A pair ( w , i v ) of w and iv is called a four-dimensional surface pair.

Definition 67 ( ϕ c ( w , i v ) ) . Let ( w , i v ) be a four-dimensional surface pair. The surface vertices ϕ c ( w , i v ) with respect to ( w , i v ) is defined by

ϕ c ( w , i v ) : = p t ( w ) p t ( i v ) L 4 .

That is, ϕ c ( w , i v ) is the set of L 4 lattice points on the intersection of the 3-faces of w and the 3-faces of iv.

Definition 68 ( Φ c ( w , i v ) ) . Let ( w , i v ) be a four-dimensional surface pair. Let Φ c ( w , i v ) E 3 be the set of the polygonal lines obtained by joining the adjacent L 4 lattice points of ϕ c ( w , i v ) . That is, Φ ( w , i v ) is the set of all the diagonal edges (of slat tetrahedrons) whose end points are included in ϕ c ( w , i v ) . Since all the points of ϕ c ( w , i v ) are on the 3-faces of w (or iv), the points of ϕ c ( w , i v ) are connected along the 3-face of w (or iv).

Definition 69 ( m c ( w , i v ) ) . Let ( w , i v ) be a four-dimensional surface pair. The two-dimensional surface mesh m c ( w , i v ) with respect to ( w , i v ) is defined by

m c ( w , i v ) : = π p t ( Φ c ( w , i v ) ) H 3 D .

We also call m c ( w , i v ) (two-dimensional) surface vein with respect to ( w , i v ) .

Definition 70 ( ( w c ( c 0 ) , i v c ( c 0 ) ) ) . Let c 0 be a four-dimensional tangent cone. The four-dimensional surface pair ( w c ( c 0 ) , i v c ( c 0 ) ) associated with c 0 is a pair of the (extended) cotangent roof w c ( c 0 ) associated with c 0 and the (extended) inverted cotangent roof i v c ( c 0 ) associated with c 0 which is defined by

i v c ( c 0 ) : = E I r o o f * t t o p ( c 0 ) .

Definition 71 ( R c ( w , i v ) ) . Let ( w , i v ) be a four-dimensional surface pair. The region R c ( w , i v ) of H 3 D is the region covered by the union of all the closed trajectories of flat tetrahedrons surrounded by m c ( w , i v ) .

Example 28 In the case of Figure 10(a),

{ w c ( c 1 ) = R o o f * { x 0 x 1 x 2 , x 0 x 2 x 3 , x 0 x 1 x 3 , x 1 x 2 x 3 } = C o n e * { 1 } , i v c ( c 1 ) = I R o o f * { x 0 x 1 x 2 , x 0 x 2 x 3 , x 0 x 1 x 3 , x 1 x 2 x 3 } = I C o n e * { x 0 3 x 1 3 x 2 3 x 3 3 } , ϕ c ( w c ( c 1 ) , i v c ( c 1 ) ) = { P y z w , P x z w , P x y w , P x y z , P x 2 y 2 z w , P x 2 y z 2 w , P x 2 y z w 2 , P x y 2 z 2 w , P x y 2 z w 2 , P x y z 2 w 2 } .

Then, the surface vein m c ( w c ( c 1 ) , i v c ( c 1 ) ) corresponds to the set of the diagonal edges colored black. R c ( w c ( c 1 ) , i v c ( c 1 ) ) is the rhombic dodecahedron consisting of four closed trajectories of length four.

Example 29 In the case of Figure 10(b),

{ w c ( c 2 ) = E r o o f * { x 0 x 2 , x 1 x 2 , x 1 x 3 } = C o n e * { 1 } , i v c ( c 2 ) = E I r o o f * { x 0 x 2 , x 1 x 2 , x 1 x 3 } = I C o n e * { x 0 3 x 1 3 x 2 3 x 3 3 } , ϕ c ( w c ( c 2 ) , i v c ( c 2 ) ) = { P x z w , P x 2 y z w 2 , P x y w , P x 2 y 2 z w , P x y z , P x y z 2 w 2 , P y z w } .

Then, the surface vein m c ( w c ( c 2 ) , i v c ( c 2 ) ) corresponds to the set of the diagonal edges colored black. Unlike the case of flows of triangles, some diagonal edges (colored grey) are not included in m c ( w c ( c 2 ) , i v c ( c 2 ) ) . R c ( w c ( c 2 ) , i v c ( c 2 ) ) is the rhombic dodecahedron consisting of a closed trajectories of length 24. In particular,

R c ( w c ( c 2 ) , i v c ( c 2 ) ) = R c ( w c ( c 1 ) , i v c ( c 1 ) ) .

Example 30 In the case of Figure 10(c),

{ w c ( c 3 ) = E r o o f * { x 0 x 2 , x 1 x 2 , x 2 x 3 , x 0 x 1 x 3 } = C o n e * { 1 } , i v c ( c 3 ) = E I r o o f * { x 0 x 2 , x 1 x 2 , x 2 x 3 , x 0 x 1 x 3 } = I C o n e * { x 0 3 x 1 3 x 2 3 x 3 3 } , ϕ c ( w c ( c 3 ) , i v c ( c 3 ) ) = { P x z w , P x 2 y z w 2 , P x y w , P x 2 y 2 z w , P x y z , P x y 2 z w 2 , P y z w } .

Then, the surface vein m c ( w c ( c 3 ) , i v c ( c 3 ) ) corresponds to the set of the diagonal edges colored black. The diagonal edges colored grey are not included in m c ( w c ( c 3 ) , i v c ( c 3 ) ) . R c ( w c ( c 3 ) , i v c ( c 3 ) ) is also the rhombic dodecahedron consisting of another closed trajectories of length 24. In particular,

R c ( w c ( c 3 ) , i v c ( c 3 ) ) = R c ( w c ( c 1 ) , i v c ( c 1 ) ) .

Example 31 In the case of Figure 8(a),

{ w c ( c 0 ) = E r o o f * { x 0 x 1 x 2 , x 0 x 2 x 3 , x 0 x 1 x 3 } = C o n e * { x 0 x 1 x 2 , x 0 x 2 x 3 , x 0 x 1 x 3 } , i v c ( c 0 ) = E I r o o f * { x 0 x 1 x 2 , x 0 x 2 x 3 , x 0 x 1 x 3 } = I C o n e * { x 0 x 1 x 2 , x 0 x 2 x 3 , x 0 x 1 x 3 } , ϕ c ( w c ( c 0 ) , i v c ( c 0 ) ) = { P x y z , P x 2 y z 2 w , P x z w , P x 2 y z w 2 , P x y w , P x 2 y 2 z w } .

Then, the surface vein m c ( w c ( c 0 ) , i v c ( c 0 ) ) is empty.

In the above examples, a tangent cone c is given first. Then, we compute the surface vein m c ( w c ( c ) , i v c ( c ) ) for ( w c ( c ) , i v c ( c ) ) (Figure 5(a)). Now, let us consider the case where a surface pair ( w , i v ) is given first (Figure 5(b)).

Definition 72 Let ( w , i v ) be a four-dimensional surface pair. The tangent cone c t ( w , i v ) associated with ( w , i v ) is defined by

c t ( w , i v ) : = C o n e ϕ c ( w , i v ) .

Starting with some definitions, we will consider the vector field induced by the tangent cone c t ( w , i v ) associated with the given surface pair ( w , i v ) .

Definition 73 Sets of four-dimensional cones are defined by

T C 3 : = { C o n e A | A L 4 and A isfinite } , C C 3 : = { C o n e * A | A L 4 and A isfinite } , I C C 3 : = { I C o n e * A | A L 4 and A isfinite } .

Definition 74 ( R t ,3 and R c ,3 ) Sets of regions on H 2 D are defined by

{ R t ,3 : = { R t ( c ) | c T C 3 } , R c ,3 : = { R c ( w , i v ) | ( w , i v ) C C 3 × I C C 3 } .

Definition 75 ( ι 3 and τ 3 ) Maps between regions on H 3 D are defined by

Theorem 6 ι 3 is not surjective.

Proof. Let us consider the case of Figure 11(a), where

{ ( w s , i v s ) = ( C o n e * { P e } , I C o n e * { P f } ) R c ( w s , i v s ) = aflatteneddodecahedron .

On the other hand, shown in Figure 11(b) are the closed trajectories of τ 3 ( R c ( w s , i v s ) ) = R t ( c t ( w s , i v s ) ) , where

{ c t ( w s , i v s ) = C o n e { P a , P b , P c , P d } , w c ( c t ( w s , i v s ) ) = C o n e * { P e } , R t ( c t ( w s , i v s ) ) = agarlicbulb-likestructure .

Note that Q p dose not appear because P p is buried beneath the surface of c t ( w s , i v s ) due to the difference of slope inclination between tangent cones and cotangent cones. The flattened dodecahedron of (a) gets dented on the bottom and has vertical linear grooves on the side. That is, we obtain a garlic bulb-like structure as a result of the constraints.

Therefore,

τ 3 ( R c ( w s , i v s ) ) R c ( w s , i v s ) .

It follows immediately that R c ( w s , i v s ) ι 3 ( R t ,3 ) . □

Corollary 3 (Self-eclipse of R c ( w , i v ) ) There exists a four-dimensional surface pair ( w , i v ) such that τ 3 ( R c ( w , i v ) ) R c ( w , i v ) .

That is, there exist global geometrical constraints on the shape of complexes of closed trajectories of tetrahedrons.

Definition 76 (Self-eclipsed protein molecules (Toy model)) A surface pair ( w , i v ) is called a self-eclipsed surface pair if τ 3 ( R c ( w , i v ) ) R c ( w , i v ) . The complex of closed trajectories induced by a tangent cone c is called a self-eclipsed protein molecule if there exists a self-eclipsed surface pair ( w , i v ) such that c = c t ( w , i v ) .

Example 32 In the case of the garlic bulb-like structure obtained in the proof of Theorem 6 (Figure 11(a) and Figure 11(b)),

τ 3 ( R c ( C o n e * { P e } , I C o n e * { P f } ) ) = R t ( C o n e { P a , P b , P c , P d } ) = R c ( C o n e * { P e } , I C o n e * { P g , P h , P i } ) R c ( C o n e * { P e } , I C o n e * { P f } ) .

Example 33 Shown in Figure 11(c) is a complex of closed trajectories of R t ( c 1 ) , where

{ c 1 = C o n e { P b , P c , P d , P j , P k , P l , P m , P n , P o } , w c ( c 1 ) = C o n e * { P e } , R t ( c 1 ) = agarlicbulb-likestructure .

Figure 11. Self-eclipsed protein molecules: (a) The surface vein m c ( w s , i v s ) with respect to ( w s , i v s ) . m c ( w s , i v s ) is indicated by the black thick polygonal lines, where ( w s , i v s ) = ( C o n e * { P e } , I C o n e * { P f } ) ; (b) Closed trajectories of V c t ( w s , i v s ) . Shown above is the top view. Shown below is the bottom view. m t ( c t ( w s , i v s ) ) is indicated by the black thick polygonal lines. The grey thick polygonal lines are the diagonal edges of type E tetrahedrons of Figure 9; (c) Closed trajectories of V c 1 , where c 1 = C o n e { P b , P c , P d , P j , P k , P l , P m , P n , P o } . m t ( c 1 ) is indicated by the black thick polygonal lines. Note that all the tetrahedrons are type D. In the figure, Q x = π p t ( P x ) ( x = a , b , c , , p ), where P a = x 0 2 x 1 0 x 2 2 x 3 2 , P b = x 0 3 x 1 2 x 2 1 x 3 3 , P c = x 0 3 x 1 2 x 2 3 x 3 1 , P d = x 0 1 x 1 2 x 2 3 x 3 3 , P e = x 0 1 x 1 0 x 2 1 x 3 1 , P f = x 0 6 x 1 6 x 2 6 x 3 6 , P g = x 0 5 x 1 4 x 2 5 x 3 4 , P h = x 0 4 x 1 4 x 2 5 x 3 5 , P i = x 0 5 x 1 4 x 2 4 x 3 5 , P j = x 0 3 x 1 1 x 2 2 x 3 3 , P k = x 0 2 x 1 1 x 2 3 x 3 3 , P l = x 0 3 x 1 1 x 2 3 x 3 2 , P m = x 0 3 x 1 2 x 2 2 x 3 2 , P n = x 0 2 x 1 2 x 2 2 x 3 3 , P o = x 0 2 x 1 2 x 2 3 x 3 2 , and P p = x 0 5 x 1 6 x 2 5 x 3 5 ( Q p is the diametrically opposite point to Q a ).

R t ( c 1 ) consists of 13 closed trajectories of length six, which sweep the same region as R t ( c t ( w s , i v s ) ) considered in the proof of Theorem 6.

On the other hand,

{ ( w c ( c 1 ) , i v c ( c 1 ) ) = ( C o n e * { P e } , I C o n e * { P q , P r , P s } ) , R c ( w c ( c 1 ) , i v c ( c 1 ) ) = R t ( c 1 ) ,

where P q = x 0 5 x 1 4 x 2 5 x 3 4 , P r = x 0 4 x 1 4 x 2 5 x 3 5 , and P s = x 0 5 x 1 4 x 2 4 x 3 5 .

Then, c t ( w c ( c 1 ) , i v c ( c 1 ) ) = c 1 and

τ 3 ( R c ( w c ( c 1 ) , i v c ( c 1 ) ) ) = R t ( c t ( w c ( c 1 ) , i v c ( c 1 ) ) ) = R t ( c 1 ) = R c ( w c ( c 1 ) , i v c ( c 1 ) ) .

That is, The complex of closed trajectories included in R c ( w c ( c 1 ) , i v c ( c 1 ) ) is not a self-eclipsed protein molecule.

4. Systems of Simultaneous Equations for Shape

Protein molecules interact each other and form an intermediate complex to perform their function. In high-throughput proteomics, proteins are characterized using a interaction network between proteins and intermediate protein complexes. Since the function of a protein is primarily determined by the three-dimensional shape, it is the shape of proteins that is characterized by the interaction network.

In this section, we will consider “interaction” between closed trajectories of triangles as a simplified geometrical description of protein interactions. Despite its simplicity, the closed trajectory model of protein interaction gives a novel geometrical interpretation of the difference between direct interactions of two proteins and cooperative interactions of three proteins (such as allosteric regulation).

4.1. Fusion and Fission of Closed Trajectories of Triangles

We have seen in the previous sections that vector fields of triangles are associated with three-dimensional tangent cones. Here we will define “fusion and fission” of closed trajectories of triangles using the tangent cone structure. For the sake of simplicity, we only consider the case of flows of triangles.

Let V c a and V c b be two vector fields of triangles induced by three-dimensional tangent cones c a and c b respectively. Then, the vector field V c a can be obtained from the other V c b by “putting unit cubes on” and/or “taking unit cubes from” the tangent cone c b . Suppose that the closed trajectories of V c a and V c b sweep the same region, i.e., R t ( c a ) = R t ( c b ) (See Definition 22). Then, the two vector fields V c a and V c b give two different decompositions of the same region into a set of closed trajectories of triangles.

Definition 77 (Addition of closed trajectories) Addition is defined between sets of all the closed trajectories of vector fields. Given two three-dimensional tangent cones c a and c b . Let { m i | 0 i N a } and { n j | 0 j N b } be the set of all the closed trajectories of the vector fields V c a and V c b respectively, where N a and N b is the numbers of the closed trajectories. Then, addition of closed trajectories of triangles is defined by

0 i N a m i : = 0 j N a n j if R t ( c a ) = R t ( c b ) .

Example 34 Shown in Figure 12(a) is a flow of triangles consisting of 36 closed trajectories of length six and infinitely many open trajectories of various lengths, where each closed trajectory sweeps a hexagonal region. By putting unit cubes on the associated tangent cone, we obtain another decomposition of the same region into a set of closed trajectories as shown in Figure 12(b). Then, we have

0 i 35 x i = 0 j 5 m j + 0 i 3 x i + x 4 + x 8 + x 9 + x 14 + x 15 + x 20 + x 21 + x 22 + x 26 + x 27 + x 28 + x 31 + 32 i 35 x i .

By removing the common terms from both sides, we obtain

x 5 + x 6 + x 7 + x 10 + x 11 + x 12 + x 13 + x 16 + x 17 + x 18 + x 19 + x 23 + x 24 + x 25 + x 29 + x 30 = m 0 + m 1 + m 2 + m 3 + m 4 + m 5 .

Figure 12. Addition of closed trajectories: (a) A flow of triangles consisting of 36 closed trajectories of length six. Shown below is the corresponding tangent cone (top view); (b) Closed trajectories of triangles obtained by putting unit cubes on the tangent cone of (a); (c) Closed trajectories of triangles obtained by putting unit cubes on the tangent cone of (b).

Note that m0 is the closed trajectory given in Figure 1. That is, we obtained m0 as a result of “fusion and fission” of 16 hexagons.

Example 35 By putting more unit cubes on the tangent cone of Figure 12(b), we obtain another equation

m 0 + m 1 + m 2 + m 3 + m 4 + m 5 + x 15 + x 22 + x 28 = m 6 .

Then, m0 appears as a “factor” of a longer closed trajectory m6.

In the above examples, closed trajectories xis (of length six) are given first. Then, m0 is obtained as a result of interactions of the xis. The challenge we propose is to give a set of equations of m0 on variables xis first, and solve the system of simultaneous equations.

Open Problem 1 (Simultaneous equations for shape) Let { x 0 , x 1 , , x N } be a finite set of closed trajectories of length six, i.e., hexagons. Let { m 0 , m 1 , , m M } be a finite set of closed trajectories of length longer than six. Suppose that we are given a finite set of addition equations with respect to xis and mjs:

f k ( x 0 , x 1 , , x N , m 0 , m 1 , , m M ) = g k ( x 0 , x 1 , , x N , m 0 , m 1 , , m M ) , ( 0 k K ) ,

where fks and gks are finite sets of terms separated by addition sign (i.e., addition expressions with coefficients one).

Find three-dimensional tangent cones for the variables xis that make the addition equations true, where xis are assigned the closed trajectory induced by the corresponding tangent cone. Then, mis are obtained as intermediate products of the interactions between xis.

Since the interaction of closed trajectories is primarily determined by their contours, it is their contours which are characterized by a set of simultaneous equations. In the case of Open Problem 1, a closed trajectory m0 is characterized using interactions between closed trajectories of length six and other closed trajectories. Therefore, the set of equations is nothing but a specification of the shape of m0 if m0 is uniquely determined.

4.2. Allosteric Regulation of Interactions

Now let us consider the difference between direct interactions of two proteins and cooperative interactions of three proteins (such as allosteric regulation). In our closed trajectory model, allosteric regulation corresponds to the complex of self-eclipsed closed trajectories of triangles (Definition 36). We will start with a brief introduction to allosteric regulation.

4.2.1. Introduction to Allosteric Regulation

In biological systems, all proteins bind to other molecules to carry out their functions. For example, enzymes bind to one or more reactant molecules to catalyze chemical reactions in our body. The region on the surface to which other molecules bind is called the active site.

The binding of a molecule at an active site is often controlled by the binding of another molecule at a distant site other than the active site. This type of regulation of protein function is called allosteric regulation. The distant site is called an allosteric site.

Allosteric regulation, which is known as “the second secret of life”’, second only to the genetic code [9] [10] , is ubiquitous in biological processes. But we still lack general understanding of the mechanisms underlying the coupling between allosteric and active sites [11] [12] .

Allosteric regulation is typically triggered by the binding of a small molecule, but also triggered by the binding of another protein. When proteins bind to other molecules or proteins, changes in conformation and/or dynamics occur within the protein. Classically, allosteric regulation was considered to be induced through a change in conformation of the protein. Today, it is believed that allostery can take place through a change in the dynamic fluctuations (i.e., internal motions and vibrations) of the protein even without obvious conformational changes [13] [14] [15] .

Currently almost all the drugs modify the actions of proteins by directly binding to their active sites. On the other hand, gaining increasing attention recently in drug discovery is another type of drugs, called allosteric drugs, which bind to the allosteric sites on their target proteins [16] [17] [18] . This is because allosteric drugs have several advantages over traditional drugs, such as higher specificity, fewer side effects, and lower toxicity.

However, allosteric drug discovery is more challenging than traditional drug discovery due to difficulties in identification of allosteric sites, prediction of drug modulatory effects, and others. For example, allosteric sites may have features we are not yet aware of because of our insufficient understanding of how coupling between the active site and the allosteric site occurs.

In the past ten years, various computational approaches have been developed for identification and characterization of allosteric sites as the first step in allosteric drug discovery [19] [20] [21] . In a static geometry-based model, protein structures are transformed into residue interaction graphs (RIGs), where amino acid residues are graph nodes and their interactions are the graph edges [22] . Then, allosteric regulation and communication are characterized using a local centrality measure (local closeness) and other newly developed quantifiable measures (binding leverage and leverage coupling). In normal mode analysis (NMA)-based models, proteins are represented by a set of Ca atoms interconnected by a network of elastic springs. Then, the structural fluctuations of a protein are decomposed into harmonic orthogonal modes and the long-range nature of allosteric communication is often well-described by low-frequency modes. However, Molecular dynamics (MD) remains the standard computational tool for structural analysis when structures are available.

4.2.2. The SECT Model of Allosteric Regulation

In this paper, we propose a novel geometrical interpretation of the long-distance regulation of protein interactions (with no conformational change). In particular, we consider how the coupling between active and allosteric sites occurs using the simplified structural description. Analysis of this model has allowed us to characterize a novel geometrical aspect of the structural coupling between active and allosteric sites.

Definition 78 (The SECT model of allosteric regulation) The Self-Eclipsed Closed Trajectory (SECT) model is a simplified geometrical interpretation of protein interactions, where

Ÿ Protein molecules correspond to a complex of closed trajectories of triangles,

Ÿ Protein interactions are represented as additions of the corresponding complexes of closed trajectories,

Ÿ Protein molecules with allosteric sites correspond to a complex of self-eclipsed closed trajectories,

Ÿ An allosteric site is the region of the contour eclipsed by the active site.

Note that the SECT model is a purely theoretical model. The SECT model was devised in the process of searching the definition formula of the shape of proteins. In the model, the constraint on the contour and the coupling between two sites are two sides of the same coin as shown in the example below. Roughly speaking, active and allosteric sites are closely tied to each other as an entity and its shadow.

Remark. Recall that (slant) triangles are flowing on the slopes of a “tangent” cone. On the other hand, the contour of a closed trajectory is drawn on the slopes of a “cotangent” cone. Since the “slope inclination” of tangent cones is gentler than that of cotangent cones, the existence of a local trajectory of triangles at one region (active site) of the contour affect the existence of another local trajectory at another region (allosteric site) of the contour.

Example 36 Shown in Figure 13 is the interactions between the three closed trajectories m0, m2, and x22 of Figure 12(b). The interaction between m0 and x22

Figure 13. Allosteric regulation and complexes of self-eclipsed closed trajectories: (a) An active site and an allosteric site of a self-eclipsed closed trajectory m0. Shown below is the corresponding three-dimensional tangent cone (top view); (b) Interaction of m0 and x22, i.e., m 0 + x 22 = m 6 + m 7 ; (c) Interaction of m0, m2, and x22, i.e., m 0 + x 22 + m 2 = m 8 .

is obtained by putting two unit cubes on the tangent cone of (a) (Figure 13(b)). We also obtain the interaction of the three closed trajectories m0, m2, and x22 by putting one more unit cube on the tangent cone of (b) (Figure 13(c)). However, m0 and m2 do not interact without the binding of x22 because of the overlap of the slopes of the tangent cone, i.e., self-eclipse. That is,

{ m 0 + x 22 = m 6 + m 7 , m 0 + m 2 = m 0 + m 2 , ( nointeractionwithout x 22 ) m 0 + x 22 + m 2 = m 8 .

In this case, x22 activates the interaction between m0 and m2.

5. Conclusions

We have proposed a novel simplified geometrical description of the shape of protein molecules and their interactions. Using the model, we have identified not only global geometrical constraints on the shape of proteins, but also their influence on protein interactions. As an example of the global constraints, a “garlic-bulb like structure” was shown. As an example of their influence on interactions, the structural coupling between active and allosteric sites was considered. In particular, our model gives a novel geometrical interpretation of the long-distance regulation of protein interactions, which could be important for the understanding of the basic mechanisms of allosteric regulation of protein functions.

As for future research questions, we have already proposed an open problem in the text, i.e., the problem of simultaneous equations for shape. Just as the function (i.e., shape) of a protein is determined by its interaction with other molecules, the shape of a complex of closed trajectories may be determined uniquely by its interaction (i.e., fusion and fission) with other closed trajectories.

Another open problem is about the relationship between the shape of a complex of closed trajectories of tetrahedrons and the triangle flow induced on the surface. How far can we learn about the three-dimensional shape from the two-dimensional surface triangle flow? To answer the question, we should consider a patchwork of locally overlapping three-dimensional cones that cover the surface of the complex because the surface triangle flow has “singular points”.

From the viewpoint of computer engineering, development of a “protein description language” will be required if proteins are to be designed automatically. The model of the closed trajectories of tetrahedrons (or n-simplices) is expected to give the semantics of the language as in the case of mathematical semantics and the lambda calculus [23] .

Conflicts of Interest

The author declares no conflicts of interest regarding the publication of this paper.

Cite this paper

Morikawa, N. (2018) Global Geometrical Constraints on the Shape of Proteins and Their Influence on Allosteric Regulation. Applied Mathematics, 9, 1116-1155. https://doi.org/10.4236/am.2018.910076

References

  1. 1. Morikawa, N. (2014) Discrete Differential Geometry of n-Simplices and Protein Structure Analysis. Applied Mathematics, 5, 2458-2463. https://doi.org/10.4236/am.2014.516237

  2. 2. Morikawa, N. (2017) Discrete Differential Geometry and the Structural Study of Protein Complexes. Open Journal of Discrete Mathematics, 7, 148-164. https://doi.org/10.4236/ojdm.2017.73014

  3. 3. Saitoh, S., Nakai, T. and Nishikawa, K. (1993) A Geometrical Constraint Approach for Reproducing the Native Backbone Conformation of a Protein. Proteins, 15, 191-204. https://doi.org/10.1002/prot.340150209

  4. 4. Ramachandran, G.N., Ramakrishnan, C. and Sasisekharan, V. (1963) Stereochemistry of Polypeptide Chain Configurations. Journal of Molecular Biology, 7, 95-99. https://doi.org/10.1016/S0022-2836(63)80023-6

  5. 5. Lund, O., Hansen, J., Brunak, S. and Bohr, J. (1996) Relationship between Protein Structure and Geometrical Constraints. Protein Science, 5, 2217-2225. https://doi.org/10.1002/pro.5560051108

  6. 6. Prasad, B.V. and Schmid, M.F. (2012) Principles of Virus Structural Organization. Advances in Experimental Medicine and Biology, 726, 17-47. https://doi.org/10.1007/978-1-4614-0980-9_3

  7. 7. Keef, T., Wardman, J.P., Ranson, N.A., Stockley, P.G. and Twarock, R. (2013) Structural Constraints on the Three-Dimensional Geometry of Simple Viruses: Case Studies of a New Predictive Tool. Acta Crystallographica, Section A, A69, 140-150. https://doi.org/10.1107/S0108767312047150

  8. 8. Mannige, R.V. and Brooks 3rd, C.L. (2008) Tilable Nature of Virus Capsids and the Role of Topological Constraints in Natural Capsid Design. Physical Review E, Statistical, Nonlinear, and Soft Matter Physics, 77, Article ID: 051902. https://doi.org/10.1103/PhysRevE.77.051902

  9. 9. Monod, J. (1977) Chance and Necessity: Essay on the Natural Philosophy of Modern Biology. Penguin Books Ltd., London.

  10. 10. Edelstein, S.J. (2013) Allosteric Interactions after 50 Years. Journal of Molecular Biology, 425, 1391-1395. https://doi.org/10.1016/j.jmb.2013.03.026

  11. 11. Laskowski, R.A., Gerick, F. and Thornton, J.M. (2009) The Structural Basis of Allosteric Regulation in Proteins. FEBS Letters, 583, 1692-1698. https://doi.org/10.1016/j.febslet.2009.03.019

  12. 12. Liu, J. and Nussinov, R. (2016) Allostery: An Overview of Its History, Concepts, Methods, and Applications. PLoS Computational Biology, 12, e1004966. https://doi.org/10.1371/journal.pcbi.1004966

  13. 13. Cooper, A. and Dryden, D.T. (1984) Allostery without Conformational Change. European Biophysics Journal, 11, 103-109. https://doi.org/10.1007/BF00276625

  14. 14. Kalodimos, C.G. (2011) NMR Reveals Novel Mechanisms of Protein Activity Regulation. Protein Science, 20, 773-782. https://doi.org/10.1002/pro.614

  15. 15. Motlagh, H.N., Wrabl, J.O., Li, J. and Hilser, V.J. (2014) The Ensemble Nature of Allostery. Nature, 508, 331-339. https://doi.org/10.1038/nature13001

  16. 16. Nussinov, R. and Tsai, C.J. (2013) Allostery in Disease and in Drug Discovery. Cell, 153, 293-305. https://doi.org/10.1016/j.cell.2013.03.034

  17. 17. Guarnera, E. and Berezovsky, I.N. (2016) Allosteric Sites: Remote Control in Regulation of Protein Activity. Current Opinion in Structural Biology, 37, 1-8. https://doi.org/10.1016/j.sbi.2015.10.004

  18. 18. Huang, W., Nussinov, R. and Zhang, J. (2017) Tools for Allosteric Drug Discovery: Site Identification and Focus Library Design. Methods in Molecular Biology, 1529, 439-446. https://doi.org/10.1007/978-1-4939-6637-0_23

  19. 19. Lu, S., Huang, W. and Zhang, J. (2014) Recent Computational Advances in the Identification of Allosteric Sites in Proteins. Drug Discovery Today, 19, 1595-1600. https://doi.org/10.1016/j.drudis.2014.07.012

  20. 20. Wagner, J.R., Lee, C.T., Durrant, J.D., Malmstrom, R.D., Feher, V.A. and Amaro, R.E. (2016) Emerging Computational Methods for the Rational Discovery of Allosteric Drugs. Chemical Reviews, 116, 6370-6390. https://doi.org/10.1021/acs.chemrev.5b00631

  21. 21. Greener, G.G. and Sternberg, M.J. (2018) Structure-Based Prediction of Protein Allostery. Current Opinion in Structural Biology, 50, 1-8. https://doi.org/10.1016/j.sbi.2017.10.002

  22. 22. Mitternacht, S. and Berezovsky, I.N. (2011) A Geometry-Based Generic Predictor for Catalytic and Allosteric Sites. Protein Engineering Design and Selection, 24, 405-409. https://doi.org/10.1093/protein/gzq115

  23. 23. Scott, D.S. and Strachey, C. (1971) Towards a Mathematical Semantics for Computer Languages, Technical Monograph PRG-6. Programming Research Group, Oxford University Computing Laboratory, Oxford.

List of Symbols

a [ x ρ ( 0 ) x ρ ( 1 ) ]

slant triangle defined by three points a, a x ρ ( 0 ) and a x ρ ( 0 ) x ρ ( 1 ) (Def.2)

a [ x ρ ( 0 ) x ρ ( 1 ) x ρ ( 2 ) ]

slant tetrahedron defined by four points a, a x ρ ( 0 ) , a x ρ ( 0 ) x ρ ( 1 ) and a x ρ ( 0 ) x ρ ( 1 ) x ρ ( 2 ) (Def.38)

B N

set of all flat triangles/tetrahedrons (Def.5, Def.41)

B D t ( C o n e A )

set of slant triangles/tetrahedrons of d S ( C o n e A ) (Def.21, Def.57)

C C N

set of N + 1 -dimensional tangent cones (Def.33, Def.73)

C o n e A

three/four-dimensional tangent cone (Def.8, Def.44)

C o n e * A

three/four-dimensional cotangent cone (Def.14, Def.50)

( C o n e * A , I C o n e * B )

three/four-dimensional contour/surface pair (Def.26, Def.66)

c t ( C o n e * A , I C o n e * B )

three/four-dimensional tangent cone associated with ( C o n e * A , I C o n e * B )

(Def.32, Def.72)

p t ( C o n e A )

surface lattice points of C o n e A (Def.10, Def.46)

p t ( C o n e * A )

surface lattice points of C o n e * A (Def.16, Def.52)

p t ( I C o n e * A )

surface lattice points of I C o n e * A (Def.25, Def.65)

t t o p ( C o n e A )

peaks on the boundary of C o n e A (Def.9, Def.45)

D s

gradient of slant triangle/tetrahedron s (Def.3, Def.39)

d S ( C o n e A )

set of all the slant triangles/tetrahedrons on the surface of C o n e A

(Def.11, Def.47)

E I r o o f * A

four-dimensional extended inverted cotangent cone (Def.64)

E N

N-dimensional Eucledean space

E r o o f * A

four-dimensional extended cotangent roof (Def.60)

H N D

hyperplane in E N + 1 (above Def.5, above Def.41)

I N t ( C o n e A )

set of slant triangles/tetrahedrons of d S ( C o n e A ) (Def.21, Def.57)

I C C N

set of N + 1 -dimensional tangent cones (Def.33, Def.73)

I C o n e * A

three/four-dimensional inverted cotangent cone (Def.23, Def.63)

ι N

map from R t , N to R c , N (Def.35, Def.75)

I R o o f * A

three/four-dimensional inverted cotangent roof (Def.24, Def.64)

i t o p ( I C o n e * A )

set of all the top vertices of I C o n e * A (Def.23, Def.63)

i w c ( C o n e A )

three/four-dimensional inverted cotangent roof associated with C o n e A

(Def.31, Def.70)

L N

N-dimensional standard lattice (Def.1, Def.37)

L N

N-dimensional conjugate lattice (Def.13, Def.49)

m c ( C o n e * A , I C o n e * B )

contour/surface vein with respect to ( C o n e * A , I C o n e * B ) (Def.29, Def.69)

m t ( C o n e A )

contour/surface vein with respect to C o n e A (Def.20, Def.56)

set of all natural numbers

O U T t ( C o n e A )

set of slant triangles/tetrahedrons of d S ( C o n e A ) (Def.21, Def.57)

[ p 0 , p 1 , p 2 ]

convex hull of points p0, p1 and p2 (above Def.2)

[ p 0 , p 1 , p 2 , p 3 ]

convex hull of points p0, p1, p2 and p3 (above Def.38)

ϕ c ( C o n e * A , I C o n e * B )

contour/surface vertices with respect to ( C o n e * A , I C o n e * B )

(Def.27, Def.67)

φ t ( C o n e A )

contour/surface vertices with respect to C o n e A (Def.18, Def.54)

φ t ( C o n e A )

extended surface vertices with respect to C o n e A (Def.62)

Φ c ( C o n e * A , I C o n e * B )

polygonal line (or set of the polygonal lines) obtained by joining the adjacent L N lattice points of ϕ c ( C o n e * A , I C o n e * B ) (Def.28, Def.68)

Φ t ( C o n e A )

polygonal line (or set of the polygonal lines) obtained by joining the

adjacent L N lattice points of ϕ t ( C o n e A ) (Def.19, Def.55)

π

projection of T B N onto B N (Def.7, Def.43)

π H

projection of S N onto H N D (Def.5, Def.41)

π p t

projection of E N + 1 onto H N D (above Def.5, above Def.41)

+

addition of closed trajectories (Def.77)

set of all real numbers

R c ( C o n e * A , I C o n e * B )

region of H N D (Def.30, Def.71)

R t ( C o n e A )

region of H N D (Def.22, Def.58)

R o o f * A

three/four-dimensional cotangent roof (Def.15, Def.51)

R c , N

region of H N D (Def.34, Def.74)

R t , N

region of H N D (Def.34, Def.74)

S E C T

self-eclipsed closed trajectory complexes (Def.36)

S E C T

model of allosteric regulation (Def.78)

Self-eclipsed contour/surface pair

(Def.36, Def.76)

Self-eclipsed protein molecule

(Def.76)

S N

set of all slant triangles/tetrahedrons (Def.2, Def.38)

s D D

slant triangle/tetrahedron adjacent to s (Def.4, Def.40)

s D U

slant triangle/tetrahedron adjacent to s (Def.4, Def.40)

s U D

slant triangle/tetrahedron adjacent to s (Def.4, Def.40)

s U U

slant triangle/tetrahedron adjacent to s (Def.4, Def.40)

s t a n d ( a )

stand of a (Def.59)

S T A N D ( A )

stand of A (Def.59)

S y m N

group of all permutations of the N-elements set { 0,1, , N 1 }

τ N

map from R c , N to R t , N (Def.35, Def.75)

T B N

tangent space on B N (Def.7, Def.43)

T C N

set of N + 1 -dimensional tangent cones (Def.33, Def.73)

t o p ( C o n e A )

set of all the top vertices of C o n e A (Def.8, Def.44)

V C o n e A

vector field induce by C o n e A on B N (Def.12, Def.48)

w c ( C o n e A )

three/four-dimensional cotangent roof associated with C o n e A

(Def.17, Def.53)

w c ( C o n e A )

four-dimensional extended cotangent roof associated with C o n e A

(Def.61)

x 0 l x 1 m x 2 n

point or vector ( l , m , n ) in E 3

x 0 l x 1 m x 2 n x 3 k

point or vector ( l , m , n , k ) in E 4

set of all integers